Hi, I'm trying to understand how does UnifiedGenotyper usually calculate the GQ of a SNP.
I understood the model to calculate the Likelihood of all the genotypes (AA,AC, AG,AT,CC,CG,CT,GG,GT,TT).
Once all the likelihood have been calculated, correct me if I'm wrong, the likelihood of the best genotype is normilized to 1 and all the other likelihoods are normilezed according to that scale.
So the PL field in the VCF should be the Phred-scaled values of the normalized values.
But is not clear to me how do they finally calculate the GQ value.
What values do they use to calculate that quality (normalized or phred scaled)? And what's the right formula?
I've tried to debug the code but it ends to be really tricky.
I really hope that anyone of you know the answer and would help me. I would be really thankfull for that.
I understood the model to calculate the Likelihood of all the genotypes (AA,AC, AG,AT,CC,CG,CT,GG,GT,TT).
Once all the likelihood have been calculated, correct me if I'm wrong, the likelihood of the best genotype is normilized to 1 and all the other likelihoods are normilezed according to that scale.
So the PL field in the VCF should be the Phred-scaled values of the normalized values.
But is not clear to me how do they finally calculate the GQ value.
What values do they use to calculate that quality (normalized or phred scaled)? And what's the right formula?
I've tried to debug the code but it ends to be really tricky.
I really hope that anyone of you know the answer and would help me. I would be really thankfull for that.
Comment