Question: Convert VCF dosage file into dosage only file
gravatar for dp0b
3.0 years ago by
dp0b50 wrote:


Is it possibly to only output dosages from a generated vcf dosage file from minimac?

Input format is

##FORMAT=<ID=DS,Number=1,Type=Float,Description="Estimated Alternate Allele Dosage : [P(0/1)+2*P(1/1)]">
##INFO=<ID=AF,Number=1,Type=Float,Description="Estimated Alternate Allele Frequency">
##INFO=<ID=MAF,Number=1,Type=Float,Description="Estimated Minor Allele Frequency">
##INFO=<ID=R2,Number=1,Type=Float,Description="Estimated Imputation Accuracy">
##INFO=<ID=ER2,Number=1,Type=Float,Description="Empirical (Leave-One-Out) R-square (available only for genotyped variants)">
#CHROM  POS ID  REF ALT QUAL    FILTER  INFO    FORMAT  1242658141  1364665948  1242658615
29  11  Chr29:11    A   G   .   PASS    .   GT:DS   0|0:0.193   0|0:0.193   0|0:0.193

Is there an option in vcftools to output the dosages for each genotype seperatly, without the genotype call?

minimac dosage vcf • 3.4k views
ADD COMMENTlink modified 3.0 years ago by trausch1.5k • written 3.0 years ago by dp0b50
gravatar for trausch
3.0 years ago by
trausch1.5k wrote:

bcftools query -f "%CHROM\t%POS[\t%DS]" minimac3.dose.vcf.gz

ADD COMMENTlink written 3.0 years ago by trausch1.5k

Missing an \n at the end :):

bcftools query -f "%CHROM\t%POS[\t%DS]\n" minimac3.dose.vcf.gz
ADD REPLYlink written 2.8 years ago by endrebak830
Please log in to add an answer.


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