Skip to content
Snippets Groups Projects
Commit 1735ea0d authored by Susheel Busi's avatar Susheel Busi
Browse files

fixed my (SSB) silly typos - facepalming

parent 51c5727f
No related branches found
No related tags found
1 merge request!67WIP: Checkpoint snakefile
......@@ -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})
"""
......
......@@ -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})
"""
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment