diff --git a/2019_GDB/rules/ANALYSIS_RULES b/2019_GDB/rules/ANALYSIS_RULES index 7d9fd4d491cb8cd09ae1bd0a0136173835387dd8..0a213daff69633deaa521867056c435ddcee0265 100755 --- a/2019_GDB/rules/ANALYSIS_RULES +++ b/2019_GDB/rules/ANALYSIS_RULES @@ -158,7 +158,7 @@ rule casc: os.path.join(RESULTS_DIR, "assembly/{assembly}.fa") output: cas1=os.path.join(RESULTS_DIR, "analysis/crispr/casc/{assembly}_casc_output/{assembly}.results.txt"), - cas2=os.path.join(RESULTS_DIR, "analysis/crispr/casc/{assembly}_casc_output/casc_CRISPR_output.txt") +# cas2=os.path.join(RESULTS_DIR, "analysis/crispr/casc/{assembly}_casc_output/casc_CRISPR_output.txt") # cas3=directory(os.path.join(RESULTS_DIR, "analysis/crispr/casc/{assembly}_casc_output")) log: os.path.join(RESULTS_DIR, "analysis/crispr/casc/{assembly}_casc_output/casc.log") threads: config["mmseq2"]["threads"] @@ -169,8 +169,8 @@ rule casc: export PATH={config[CASC][PATH]} &&\ export PERL5LIB={config[CASC][PERL5LIB]} &&\ casc -i {input} -o $(dirname {output.cas1}) -n 12 --conservative &&\ - echo $(basename -s ".results.txt" {output.cas1}) >> {output.cas2} &&\ - cat {output.cas1}) | awk '$8 == "true" {{sum += $2}} END {{print sum}}' >> {output.cas2} &&\ + echo $(basename -s ".results.txt" {output.cas1}) >> results/analysis/crispr/casc/casc_CRISPR_output.txt &&\ + cat {output.cas1} | awk '$8 == "true" {{sum += $2}} END {{print sum}}' >> results/analysis/crispr/casc/casc_CRISPR_output.txt &&\ date) &> >(tee {log}) """ diff --git a/rules/ANALYSIS_RULES b/rules/ANALYSIS_RULES index 7d9fd4d491cb8cd09ae1bd0a0136173835387dd8..0a213daff69633deaa521867056c435ddcee0265 100755 --- a/rules/ANALYSIS_RULES +++ b/rules/ANALYSIS_RULES @@ -158,7 +158,7 @@ rule casc: os.path.join(RESULTS_DIR, "assembly/{assembly}.fa") output: cas1=os.path.join(RESULTS_DIR, "analysis/crispr/casc/{assembly}_casc_output/{assembly}.results.txt"), - cas2=os.path.join(RESULTS_DIR, "analysis/crispr/casc/{assembly}_casc_output/casc_CRISPR_output.txt") +# cas2=os.path.join(RESULTS_DIR, "analysis/crispr/casc/{assembly}_casc_output/casc_CRISPR_output.txt") # cas3=directory(os.path.join(RESULTS_DIR, "analysis/crispr/casc/{assembly}_casc_output")) log: os.path.join(RESULTS_DIR, "analysis/crispr/casc/{assembly}_casc_output/casc.log") threads: config["mmseq2"]["threads"] @@ -169,8 +169,8 @@ rule casc: export PATH={config[CASC][PATH]} &&\ export PERL5LIB={config[CASC][PERL5LIB]} &&\ casc -i {input} -o $(dirname {output.cas1}) -n 12 --conservative &&\ - echo $(basename -s ".results.txt" {output.cas1}) >> {output.cas2} &&\ - cat {output.cas1}) | awk '$8 == "true" {{sum += $2}} END {{print sum}}' >> {output.cas2} &&\ + echo $(basename -s ".results.txt" {output.cas1}) >> results/analysis/crispr/casc/casc_CRISPR_output.txt &&\ + cat {output.cas1} | awk '$8 == "true" {{sum += $2}} END {{print sum}}' >> results/analysis/crispr/casc/casc_CRISPR_output.txt &&\ date) &> >(tee {log}) """