diff --git a/Snakefile b/Snakefile index b29ddc3fedc7f97d36ca830e7222e4b9c093e7e2..c8d28aa2b166a36b6abcda7283ae17867e23441f 100644 --- a/Snakefile +++ b/Snakefile @@ -330,7 +330,7 @@ rule mapped_basenum: output: num=out_path("{sample}/bams/{sample}.mapped.basenum") conda: "envs/samtools.yml" - shell: "samtools view -F 4 {input.bam} | wc -c > {output.num}" + shell: "samtools view -F 4 {input.bam} | cut -f10 | wc -c > {output.num}" rule unique_num: