就是指唯一比對的reads
現在人們已經開始避免使用unique mapped reads這個概念了,而轉向使用mapq值來保留高質量的比對結果。因為mapq值反應了一組比對結果發生的可能性,MapQ = -10 log10(P), 比如結果為10,那就是1/10的概率會出現這個比對結果,如果我們認為0.05%是一個小概率的話,那個mapq值為15就可以用於篩選了, 如果認為0.01%是個小概率的話,mapq值為20就可以用於篩選了。但是人們往往從30這個值開始試起(1/1000的概率),如果它的篩選結果符合你的測序要求,就可以使用它了。如果不行,可以適當的調整這個篩選值。
samtools view -bhS -q 30 input.sam > output.bam