Phred Quality score(Q score)可用來評估次世代定序的鹼基品質,Q30 所表示的鹼基定序錯誤概率為何?
詳細解析
本題觀念
次世代定序(NGS)在單一次反應中會產生數以百萬計的鹼基讀取結果,因此需要一套客觀且可量化的品質評估工具,來判斷每一個鹼基的定序可信度,Phred quality score(Q score)正是為此目的所設計的對數化品質分數。醫事檢驗師在執行次世代定序相關檢驗、判讀定序報告或撰寫品管紀錄時,都必須清楚掌握Q score與鹼基定序錯誤機率之間的換算關係,這正是本題所要考核的核心概念。
Phred quality score最早應用於Sanger定序技術的層析圖判讀,其數學定義為 Q = -10 × log10(P),其中P代表該鹼基被錯誤定序的機率;將公式反推可得到 P = 10^(-Q/10)。這種以對數尺度呈現誤差機率的設計,讓極小的機率差異能用簡單的整數分數清楚表示:Q值每增加10,代表定序錯誤機率降低為原本的十分之一,也就是準確度提升一個數量級。這套換算關係後來被延伸應用到次世代定序平台(如Illumina)的品質報告系統,成為業界通用的定序品質溝通語言。
選項分析
A:1/10。將P=1/10代入公式可得 Q = -10 × log10(1/10) = 10,對應的是Q10而非Q30,代表正確率僅約90%,屬於品質偏低的等級,不符合本題所問,不正確。
B:1/100。將P=1/100代入公式可得 Q = -10 × log10(1/100) = 20,對應的是Q20而非Q30,正確率約99%,雖優於Q10但仍低於Q30所代表的品質水準,不正確。
C:1/1000。將P=1/1000代入公式可得 Q = -10 × log10(1/1000) = 30,與題目所問的Q30完全吻合,代表正確率約99.9%,也就是每定序1000個鹼基平均出現1個錯誤,此為Illumina等業界公認的次世代定序資料品質「黃金標準」門檻,正確,為本題答案。
D:1/10000。將P=1/10000代入公式可得 Q = -10 × log10(1/10000) = 40,對應的是比Q30更高等級的Q40,代表正確率達99.99%,並非Q30所指的錯誤機率,不正確。
答案解析
本題答案為C。依Phred quality score的標準公式 Q = -10 × log10(P),將Q=30代入可反推得到 P = 10^(-30/10) = 10^(-3) = 1/1000,也就是每定序1000個鹼基,平均會出現1個錯誤鹼基,正確率達99.9%。Illumina與GATK(Broad Institute)等次世代定序技術文件皆採用相同的定義與換算方式,並將Q30訂為次世代定序資料可信賴使用的最低品質門檻。在實務判讀上,若某個鹼基位點的Q值低於30,代表其定序錯誤機率高於千分之一,後續進行變異點呼叫(variant calling)等下游分析時應特別留意該位點資料的可靠度,必要時需以其他方法(如Sanger定序)進行確認。
核心知識點
Phred quality score是評估定序品質最核心的量化工具,其對應關係應熟記:Q10對應錯誤機率1/10(正確率90%)、Q20對應1/100(正確率99%)、Q30對應1/1000(正確率99.9%)、Q40對應1/10000(正確率99.99%)。Q值每提升10,錯誤機率即降低十倍,這種對數關係讓即使是極微小的品質差異,也能以簡單整數清楚呈現。醫事檢驗師在執行次世代定序為基礎的遺傳診斷或基因體定序檢驗時,應熟記Q30作為業界公認的最低可信門檻,並能依公式在Q值與錯誤機率之間相互換算,以利判讀定序報告品質與撰寫品管紀錄。