Question: Quality Scores In A Multi Sample Vcf File
0
gravatar for Kssr
6.5 years ago by
Kssr110
Kssr110 wrote:

How is error probability estimated while calculating phred quality scores? I have multi samples in my vcf file and I see many SNP's having a quality score of '999'. Is there any threshold value after which the SNP's are assigned a quality score of '999'.

ADD COMMENTlink modified 6.5 years ago by Istvan Albert ♦♦ 81k • written 6.5 years ago by Kssr110
0
gravatar for Istvan Albert
6.5 years ago by
Istvan Albert ♦♦ 81k
University Park, USA
Istvan Albert ♦♦ 81k wrote:

The limits are probably hard coded into the software - a 999 value is probably unrealistic anyway. If not from mathematical perspective, there may be good reasons for a value to be this but from biological or experimental perspective. The underlying error rates of the experimental protocol or sequencing platform most likely far exceed the reported value that simply relies on the results that are most certainly biased above those values.

ADD COMMENTlink written 6.5 years ago by Istvan Albert ♦♦ 81k

I agree with the point that '999' value might be hard coded in to the software.I am not clear about what you have tried to explain.

I am using samtools mpileup output followed by bcftools to do the variant calling.I am assuming bcftools incorporates this value in to the vcf file.I see quality scores ranging from some min value up to 300 and rest all are placed in 999 category. Wondering whether these snp qualities couldn't be defined?In general, not sure how to interpret this quality score.

Any help appreciated.

ADD REPLYlink written 6.5 years ago by Kssr110

my guess would be that under a certain minimum it generates a value that corresponds to the probability value of 0 That is 1E-99

ADD REPLYlink modified 6.5 years ago • written 6.5 years ago by Istvan Albert ♦♦ 81k
Please log in to add an answer.

Help
Access

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 705 users visited in the last hour