@@ -10,10 +10,6 @@ def get_extension_path() -> Path:
1010 return Path (__file__ ).parent .resolve ()
1111
1212
13- def cenote_output () -> Path :
14- return VIRUS_FP / "cenote_taker" / "{sample}.fasta"
15-
16-
1713rule all_cenote_taker :
1814 input :
1915 expand (
@@ -119,7 +115,7 @@ rule align_virus_reads:
119115 output :
120116 temp (VIRUS_FP / "alignments" / "{sample}.sam" ),
121117 params :
122- index = str (cenote_output () ),
118+ index = str (VIRUS_FP / "cenote_taker" / "filtered" / "{sample}.fasta" ),
123119 threads : 6
124120 conda :
125121 "envs/sbx_cenote_taker.yml"
@@ -137,7 +133,7 @@ rule process_virus_alignment:
137133 sorted = temp (VIRUS_FP / "alignments" / "{sample}.sorted.bam" ),
138134 bai = temp (VIRUS_FP / "alignments" / "{sample}.sorted.bam.bai" ),
139135 params :
140- target = str (cenote_output () ),
136+ target = str (VIRUS_FP / "cenote_taker" / "filtered" / "{sample}.fasta" ),
141137 conda :
142138 "envs/sbx_cenote_taker.yml"
143139 container :
@@ -170,7 +166,7 @@ rule virus_mpileup:
170166 input :
171167 bam = VIRUS_FP / "alignments" / "{sample}.sorted.bam" ,
172168 idx = VIRUS_FP / "alignments" / "{sample}.sorted.bam.bai" ,
173- contigs = cenote_output () ,
169+ contigs = VIRUS_FP / "cenote_taker" / "filtered" / "{sample}.fasta" ,
174170 output :
175171 VIRUS_FP / "alignments" / "{sample}.mpileup" ,
176172 conda :
@@ -282,7 +278,7 @@ rule virus_coverage_per_gene:
282278 output :
283279 tsv = VIRUS_FP / "alignments" / "{sample}.gene_coverage.tsv" ,
284280 params :
285- contigs = cenote_output ( ),
281+ contigs = str ( VIRUS_FP / "cenote_taker" / "filtered" / "{sample}.fasta" ),
286282 conda :
287283 "envs/sbx_cenote_taker.yml"
288284 container :
0 commit comments