| 2 | 1/1 | 返回列表 |
| 查看: 1002 | 回復(fù): 1 | |||
| 本帖產(chǎn)生 1 個 BioEPI ,點擊這里進行查看 | |||
xiezhancnu木蟲 (小有名氣)
|
[求助]
vcf格式疑問,想請教大神! 已有1人參與
|
||
|
各位好: 請問一下 vcf格式中 QUAL這一列 我查了一下是這樣的: QUAL:Phred格式(Phred_scaled)的質(zhì)量值,表示在該位點存在variant的可能性;該值越高,則variant的可能性越大;計算方法:Phred值 = -10 * log (1-p) p為variant存在的概率; 通過計算公式可以看出值為10的表示錯誤概率為0.1,該位點為variant的概率為90%。 我想問一下 這里p是怎么算出來的? 我一直對這個vcf文件中QUAL是怎么來的不太明白 想請教 金幣不是問題,解釋清楚了,要多少給多少,呵呵~~ |
至尊木蟲 (著名寫手)
|
先看文獻http://www.ncbi.nlm.nih.gov/pubmed/9521921 然后我的理解和你稍微不同,phred不是測variant,是測測序序列的準(zhǔn)確性的。 像你說的phred=10,就是在那個位點得到的堿基有90%的可能性是準(zhǔn)確的 這個P的計算有點像machine learning,phred程序中保存了大量的數(shù)據(jù)。這些數(shù)據(jù)來源于測已知序列,記錄了測序波形和精度與正確堿基判斷之間的關(guān)系。然后通過對當(dāng)前波形和精度的判斷,在已知數(shù)據(jù)中找匹配值,而得到P。 |
| 2 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 286求調(diào)劑 +5 | Faune 2026-03-21 | 5/250 |
|
|---|---|---|---|---|
|
[考研] 299求調(diào)劑 +3 | 某某某某位 2026-03-21 | 3/150 |
|
|
[考研] 一志愿華中科技大學(xué)071000,求調(diào)劑 +3 | 沿岸有貝殼6 2026-03-21 | 3/150 |
|
|
[考研] 能源材料化學(xué)課題組招收碩士研究生8-10名 +5 | 脫穎而出 2026-03-16 | 15/750 |
|
|
[考研] 297求調(diào)劑 +9 | 戲精丹丹丹 2026-03-17 | 9/450 |
|
|
[考研] 材料 336 求調(diào)劑 +3 | An@. 2026-03-18 | 4/200 |
|
|
[考研] 材料專碩英一數(shù)二306 +7 | z1z2z3879 2026-03-18 | 7/350 |
|
|
[考研] 321求調(diào)劑 +9 | 何潤采123 2026-03-18 | 11/550 |
|
|
[考研] 一志愿蘇州大學(xué)材料求調(diào)劑,總分315(英一) +5 | sbdksD 2026-03-19 | 5/250 |
|
|
[考研] 一志愿 南京航空航天大學(xué)大學(xué) ,080500材料科學(xué)與工程學(xué)碩 +5 | @taotao 2026-03-20 | 5/250 |
|
|
[考研] 工科材料085601 279求調(diào)劑 +7 | 困于星晨 2026-03-17 | 9/450 |
|
|
[考博] 招收博士1-2人 +3 | QGZDSYS 2026-03-18 | 3/150 |
|
|
[考研] 求調(diào)劑 +3 | 暗涌afhb 2026-03-16 | 3/150 |
|
|
[考研] 收復(fù)試調(diào)劑生 +4 | 雨后秋荷 2026-03-18 | 4/200 |
|
|
[考研] 312求調(diào)劑 +8 | 陌宸希 2026-03-16 | 9/450 |
|
|
[考研]
|
胡辣湯放糖 2026-03-15 | 6/300 |
|
|
[考研] 考研化學(xué)學(xué)碩調(diào)劑,一志愿985 +4 | 張vvvv 2026-03-15 | 6/300 |
|
|
[考研] 290求調(diào)劑 +3 | p asserby. 2026-03-15 | 4/200 |
|
|
[考研] 318求調(diào)劑 +3 | Yanyali 2026-03-15 | 3/150 |
|
|
[考研] 288求調(diào)劑 +4 | 奇點0314 2026-03-14 | 4/200 |
|