用samtools进行过滤,选择MapQ非0的多映射读取
具体代码示例:
samtools view -bq 1 file.bam > output.bam
其中,-b表示输出BAM格式文件,-q 1表示只选择MapQ值大于等于1的读取。
上一篇:BWA-mem和Sambamba读取组行错误
下一篇:BWA: 检测短读序列中的变异