It looks like you're new here. If you want to get involved, click one of these buttons!

- 7.4K All Categories
- 206 Announcements
- 6.5K Ask the GATK team
- 235 GATK Documentation Guide
- 30 Tutorials
- 43 FAQs
- 11 Solutions to Problems
- 33 Methods and Algorithms
- 19 Dictionary
- 13 GATK 4 Alpha
- 1 Pipelining with Queue
- 30 Developer Zone
- 55 Archive
- 454 Tools related to GATK
- 135 MuTect v1
- 11 ReCapSeg
- 42 Oncotator
- 266 GenomeSTRiP

You can opt in to receive email notifications, for example when your questions get answered or when there are new announcements, by following the instructions given here.

Wrap blocks of code, error messages and BAM/VCF snippets--especially content with hashes (#)--with lines with three backticks ( ``` ) each to make a code block as demonstrated here.

Picard 2.10.2 is now available. As of 2.10.0, Picard supports NovaSeq CBCL data. Download and read release notes at **https://github.com/broadinstitute/picard/releases**.

**GATK4-BETA.2** is here. That's TWO, as in the second beta release. Be sure to read about the known issues before test driving. See Article#9881 to start and https://github.com/broadinstitute/gatk/blob/master/README.md for details.

Sheila
Broad InstituteMember, Broadie, Moderator

You may have noticed that a lot of the scores that are output by the GATK are in Phred scale. The Phred scale was originally used to represent base quality scores emitted by the Phred program in the early days of the Human Genome Project (see this Wikipedia article for more historical background). Now they are widely used to represent probabilities and confidence scores in other contexts of genome science.

In the context of sequencing, Phred-scaled quality scores are used to represent how confident we are in the assignment of each base call by the sequencer.

In the context of variant calling, Phred-scaled quality scores can be used to represent many types of probabilities. The most commonly used in GATK is the QUAL score, or variant quality score. It is used in much the same way as the base quality score: the variant quality score is a Phred-scaled estimate of how confident we are that the variant caller correctly identified that a given genome position displays variation in at least one sample.

In today’s sequencing output, by convention, most useable Phred-scaled base quality scores range from 2 to 40, with some variations in the range depending on the origin of the sequence data (see the FASTQ format documentation for details). However, Phred-scaled quality scores in general can range anywhere from 0 to infinity. A **higher score** indicates a higher probability that a particular decision is **correct**, while conversely, a **lower score** indicates a higher probability that the decision is **incorrect**.

The Phred quality score (Q) is logarithmically related to the error probability (E).

$$ Q = -10 \log E $$

So we can interpret this score as an estimate of **error**, where the error is *e.g.* the probability that the base is called **incorrectly** by the sequencer, but we can also interpret it as an estimate of **accuracy**, where the accuracy is *e.g.* the probability that the base was identified **correctly** by the sequencer. Depending on how we decide to express it, we can make the following calculations:

If we want the probability of error (E), we take:

$$ E = 10 ^{-\left(\frac{Q}{10}\right)} $$

And conversely, if we want to express this as the estimate of accuracy (A), we simply take

$$

\begin{eqnarray}

A &=& 1 - E \nonumber \

&=& 1 - 10 ^{-\left(\frac{Q}{10}\right)} \nonumber \

\end{eqnarray}

$$

Here is a table of how to interpret a range of Phred Quality Scores. It is largely adapted from the Wikipedia page for Phred Quality Score.

For many purposes, a Phred Score of 20 or above is acceptable, because this means that whatever it qualifies is 99% accurate, with a 1% chance of error.

Phred Quality Score | Error | Accuracy (1 - Error) |
---|---|---|

10 | 1/10 = 10% | 90% |

20 | 1/100 = 1% | 99% |

30 | 1/1000 = 0.1% | 99.9% |

40 | 1/10000 = 0.01% | 99.99% |

50 | 1/100000 = 0.001% | 99.999% |

60 | 1/1000000 = 0.0001% | 99.9999% |

And finally, here is a graphical representation of the Phred scores showing their relationship to accuracy and error probabilities.

The red line shows the error, and the blue line shows the accuracy. Of course, as error decreases, accuracy increases symmetrically.

Note: You can see that below Q20 (which is how we usually refer to a Phred score of 20), the curve is really steep, meaning that as the Phred score decreases, you lose confidence very rapidly. In contrast, above Q20, both of the graphs level out. This is why Q20 is a good cutoff score for many basic purposes.

Post edited by Geraldine_VdAuwera on

Tagged:

## Comments

Hi, I have a question about interpreting variant quality. Let say I have 8 reads aligned to position 100 with base phred quality as follow:

C (REF bases): 28, 29, 31, 31, 31, 32 (6 bases)

G (ALT bases): 33, 38 (2 bases)

Mapping quality = 60 for all aligned reads.

Here is what GATK reports:

. 100 . C G 25.78 . AC=1;AF=0.500;AN=2;BaseQRankSum=1.589;DP=8;Dels=0.00;FS=0.000;HaplotypeScore=0.0000;MLEAC=1;MLEAF=0.500;MQ=60.00;MQ0=0;MQRankSum=0.000;QD=3.22;ReadPosRankSum=-0.572 GT:AD:DP:GQ:PL 0/1:6,2:8:54:54,0,184

I tried to calculate Wilcoxon-rank sum test for 2 groups of base quality and got a p-value of 0.0301490435667 (assume that alpha is 0.05). Is that the reason why GATK call that ALT base? If so, how GATK get a quite high quality like that (25.78, i.e, accuracy is more than 99%) while there are only 2 ALT aligned bases among 8 aligned bases? How GATK calculate variant quality in this case?

Thank you for your help.

@namsyvo The calling algorithm is more complex; see the HaplotypeCaller docs for details.

Why in a previous post about SNP hard filtering you mention marking as 'FAIL' SNPs with a strand bias phred-scaled p-value > 60 (thus, marking as 'PASS' SNPs with a strand bias score < 60)? Should it not be the opposite if the confidence decrease with the phred score?

@user31888

Hi,

I hope this document and this document will answer your question.

-Sheila