unique mapped reads


就是指唯一比對的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

參考:關於map當中的unique mapped reads問題


免責聲明!

本站轉載的文章為個人學習借鑒使用,本站對版權不負任何法律責任。如果侵犯了您的隱私權益,請聯系本站郵箱yoyou2525@163.com刪除。



 
粵ICP備18138465號   © 2018-2025 CODEPRJ.COM