Entering edit mode
3 months ago
Nadia
•
0
Hello, I have a VCF file with GT:HDS (haploid dosage) values, but I need to generate GT:DS (diploid dosage), since this is the most common format in vcf files and used commonly for association analyses.
Is there any way, command or plugin to compute or convert the GT:DS field from GT:HDS?
The VCF FORMAT column is GT:HDS. A sample line looks like: chr22 10510303 ... GT:HDS 0|0:0,0 0|0:0.001,0 ...
VCF Header Snippet:
FORMAT=<ID=HDS,Number=2,Type=Float,Description="Allele Probs for two haplotypes=Haploid Alternate Allele Dosage">
Thank you in advance for your guidance.
OK this is really a Minimac4 / Michigan Imputation Server concept, not a general VCF thing