Service Notice: Due to the blizzard currently hammering the US Northeast, the Broad is shut down and the GATK forum will be mostly unattended while we hunker down and sip hot cocoa with marshmallows. Assuming the power stays on and we're able to dig ourselves out of the snow when it's all over, normal service should resume Wednesday or Thursday.

Why is ABHet not appearing in the vcf file's annotations?

ro6ertro6ert Posts: 4Member

We are using the unified genotyper to create a vcf file. In the vcf file we see the info tags;

INFO=<ID=ABHet,Number=1,Type=Float,Description="Allele Balance for hets (ref/(ref+alt))">

INFO=<ID=ABHom,Number=1,Type=Float,Description="Allele Balance for homs (A/(A+O))">

but when we look at the exported annotations

ABHom=1.00;AC=2;AF=1.00;AN=2;DB;DP=2;Dels=0.00;FS=0.000;HaplotypeScore=0.0000;MLEAC=2;MLEAF=1.00;MQ=40.00;MQ0=0;QD=12.21

we only see ABHom not ABHet.

Best Answer

Answers

Sign In or Register to comment.