MAPQ 一定是255嗎?

不一定,不同開發(fā)者的定義不同。

早上爬樓看到這個問題:

各位有個基礎(chǔ)問題我不是很明白,SAM文件中的MAPQ map質(zhì)量等于255究竟是什么意思,網(wǎng)上查了好多,包括官網(wǎng)的解釋都是not available,但是像STAR這種軟件給出的卻是unique map,那么not available和unique map之間有什么關(guān)系嗎?還是每個軟件對MAPQ數(shù)值的定義本身就不同?

我們知道,在cellranger中:

MAPQ adjustment
For reads that align to a single exonic locus but also align to 1 or more non-exonic loci, the exonic locus is prioritized and the read is considered to be confidently mapped to the exonic locus with MAPQ 255.

于是,一頓搜索操作,整理成ppt:


https://support.10xgenomics.com/single-cell-gene-expression/software/pipelines/latest/algorithms/overview
http://www.acgt.me/blog/2014/12/16/understanding-mapq-scores-in-sam-files-does-37-42
https://www.biostars.org/p/178387/
https://www.biostars.org/p/372666/

最后編輯于
?著作權(quán)歸作者所有,轉(zhuǎn)載或內(nèi)容合作請聯(lián)系作者
【社區(qū)內(nèi)容提示】社區(qū)部分內(nèi)容疑似由AI輔助生成,瀏覽時請結(jié)合常識與多方信息審慎甄別。
平臺聲明:文章內(nèi)容(如有圖片或視頻亦包括在內(nèi))由作者上傳并發(fā)布,文章內(nèi)容僅代表作者本人觀點,簡書系信息發(fā)布平臺,僅提供信息存儲服務(wù)。

友情鏈接更多精彩內(nèi)容