This is a test version of Biostars. For the public version, visit https://www.biostars.org.
Annotate VCF to include dosage

Is there an easy way to annotate a VCF to include dosages? I have a VCF with phred-scaled likelihoods (PL), but no dosages (DS).

I know that there is a bcftools plugin to compute dosages based off of the PL tag, but it isn't clear how to incorporate these output into the original VCF w/o having to write my own custom scripts.

Is the bcftools annotate able to handle the output from the plugin or is there another way to annotate variant calls with the DS info?

vcf bcftools dosages
• 2,576 views • link
updated by zx8754 zx8754 • written by Nicholas Mancuso Nicholas Mancuso

1 answer

I ended up writing my own c++ program to do this quickly as possible. It is fairly bare-bones, but shouldn't be too hard to extend it to properly take care of other user's edge cases.

https://github.com/quattro/DS2VCF

• 0 views • link
written by Nicholas Mancuso Nicholas Mancuso

Log in to answer this question.

More posts like this