Question: Genotype dosage from VCF genotype likelihoods
0
gravatar for lkw222
4.6 years ago by
lkw22230
lkw22230 wrote:

Is it possible to create a file of genotype dosages from a VCF using the phred-scale likelihoods (PL) given for each possible genotype? For example the genotype call:

1 497 . C T 1882.89 PASS GT:AD:DP:GQ:PL 0/0:3,0:3:0:0,0,45

where the PL for 0/0 is 0 and for 0/1 is also 0 would be assigned a dosage of 0.5 because there is an exact uncertainty between whether the genotype should be homozygous reference (which is given a value of 0) and heterozygous (which is given a value of 1).

rna-seq genotype dosage vcf • 1.7k views
ADD COMMENTlink modified 4.3 years ago by Biostar ♦♦ 20 • written 4.6 years ago by lkw22230
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: 1287 users visited in the last hour