samtools过滤多重比对reads

samtools过滤多重比对reads

• samtools

samtools过滤多重比对reads

samtools过滤多重比对reads的命令如下：

samtools view -h -q 1 -F 4 -F 256 big.sort.dedup.bam |grep -v XA:Z |grep -v SA:Z |samtools view -b - > big.sort.dedup.unique.bam

其中参数：

-q 1：过滤掉比对质量小于1的reads；

-F 4：过滤掉没有比对上的reads；

-F 256：过滤掉比对上的多次的reads。