<p>See <a href="https://gatkforums.broadinstitute.org/gatk/discussion/2403/gatks-quality-scores-on-vcf-files" rel="nofollow noreferrer">this post</a> on the GATK forums.</p>

<p>It is normal to see values above 255.  The QUAL score is the posterior probability that all samples in your callset are homozygous reference. I quote Sheila's response from the thread posted above:</p>

<blockquote>
  <p>The QUAL score is the Phred-scaled posterior of AC = 0. We use the AC priors and the PLs to get the likelihood of the data given each AC, then use those to get the posterior probability for each AC. From there, the calculation is 1 - Pr{AC > 0}.</p>
</blockquote>

<p>Please see the GATK post in Methods and Algorithms for the <a href="https://software.broadinstitute.org/gatk/documentation/article.php?id=7258" rel="nofollow noreferrer">nitty gritty details</a>. Beware, Bayesian statistics abound.</p>

CC BY-SA 4.0
2018-07-02T18:43:41.017Z
2018-07-02T20:21:16.957Z
2018-07-02T20:21:16.957Z
Kohl Kinning
<p>GATK's <code>HaplotypeCaller</code> and <code>freebayes</code> are two more variant caller that output QUAL values greater than 255. The reason for that is, that each read have its part on calculation that value. This results in the higher the read depth at this position the higher the QUAL value.</p>

<p>To get something useful out of this value, the usual thing that is done, is to calculate the quality per depth (QD) by just divide the QUAL value by the read depth.</p>

<p>Note whenever you try a different variant caller, you have to find out whether the QUAL value is read depth dependend or not.</p>

CC BY-SA 4.0
2019-02-07T08:51:06.240Z
2019-02-07T08:51:06.240Z
finswimmer
<p>I am working on ION Torrent S5 XL and I am facing to an issue with VCF output generated by the variant caller of thermofisher. It stating on VCF 4.1 that Phred Quality Score is on a range of 0 to 255. In my VCF I have values over than 255. Do you know if there is a way to transform these values?</p>

CC BY-SA 4.0
2018-06-28T15:03:28.217Z
Let me rephrase to see if I understood: your output exceeds the max value explained by its manual? Could you post an example? My only guess is that either the structure is wrong or there is a bug, or the program is not well documented.
llrs
CC BY-SA 4.0
2018-07-02T11:58:22.927Z
Hi Llopis, here's an example: chr17 7577104 . AGGA AAAA,AAGA,AGA,AGAA,AGGG,AGTA 35144.8 PASS. 35144.8 is QUAL Score into the VCF
user979974
CC BY-SA 4.0
2018-07-02T12:30:22.720Z
Please can you [edit] the question to include it in the right formatting? Also could you post the command you used to get this VCF ?
llrs
<vcf>
High Phred Quality score VCF file
# High Phred Quality score VCF file

(Asked by: username_0 on 2018-06-28)

I am working on ION Torrent S5 XL and I am facing to an issue with VCF output generated by the variant caller of thermofisher. It stating on VCF 4.1 that Phred Quality Score is on a range of 0 to 255. In my VCF I have values over than 255. Do you know if there is a way to transform these values?

**Tags:** vcf
> **Comments:**
> **username_1**: Let me rephrase to see if I understood: your output exceeds the max value explained by its manual? Could you post an example? My only guess is that either the structure is wrong or there is a bug, or the program is not well documented.  
> **username_0**: Hi Llopis, here's an example: chr17 7577104 . AGGA AAAA,AAGA,AGA,AGAA,AGGG,AGTA 35144.8 PASS. 35144.8 is QUAL Score into the VCF  
> **username_1**: Please can you [edit] the question to include it in the right formatting? Also could you post the command you used to get this VCF ?  

---

## 2 answers

### Answer by username_2 on 2018-07-02. Score: 2

See [this post](https://gatkforums.broadinstitute.org/gatk/discussion/2403/gatks-quality-scores-on-vcf-files) on the GATK forums.

It is normal to see values above 255. The QUAL score is the posterior probability that all samples in your callset are homozygous reference. I quote Sheila's response from the thread posted above:

> The QUAL score is the Phred-scaled posterior of AC = 0. We use the AC priors and the PLs to get the likelihood of the data given each AC, then use those to get the posterior probability for each AC. From there, the calculation is 1 - Pr{AC > 0}.

Please see the GATK post in Methods and Algorithms for the [nitty gritty details](https://software.broadinstitute.org/gatk/documentation/article.php?id=7258). Beware, Bayesian statistics abound.


---

### Answer by username_3 on 2019-02-07. Score: 0

GATK's `HaplotypeCaller` and `freebayes` are two more variant caller that output QUAL values greater than 255. The reason for that is, that each read have its part on calculation that value. This results in the higher the read depth at this position the higher the QUAL value.

To get something useful out of this value, the usual thing that is done, is to calculate the quality per depth (QD) by just divide the QUAL value by the read depth.

Note whenever you try a different variant caller, you have to find out whether the QUAL value is read depth dependend or not.


---

