We read every piece of feedback, and take your input very seriously.
To see all available qualifiers, see our documentation.
1 parent 80d4ec3 commit 7469280Copy full SHA for 7469280
pipeline_consensus/Snakefile
@@ -31,7 +31,7 @@ rule call_consensus:
31
"{out_dir}/consensus_sequences/{sample}.fa"
32
shell:
33
"""
34
- samtools mpileup -A -Q 0 -d 300000 {input} | ivar consensus -p {output} -n N -m 10
+ samtools mpileup -A -Q 0 -d 300000 {input} | ivar consensus -p {output} -m 10
35
36
37
rule trim_reads:
@@ -44,7 +44,7 @@ rule trim_reads:
44
tmp="{out_dir}/trimmed_bams/{sample}.trimmed.bam"
45
46
47
- ivar trim -i {input} -b {params.bed} -p {params.tmp}
+ ivar trim -e -i {input} -b {params.bed} -p {params.tmp}
48
samtools sort -T {wildcards.sample}.trim -o {output} {params.tmp}
49
rm {params.tmp}
50
0 commit comments