Test-drive the GATK tools and Best Practices pipelines on Terra


Check out this blog post to learn how you can get started with GATK and try out the pipelines in preconfigured workspaces (with a user-friendly interface!) without having to install anything.

CombineVariants resets AF to 0.5 and 1.0

Dear GATK Team,

I have been using CombineVariants to merge SNPs.vcf and Indels.vcf of the same sample into a single vcf file.

In the SNP.vcf and Indels.vcf files, I have changed the value for AF tag in INFO filed by calculating the 'actual frequency' from the AD tag.

So, when I merged the snp and indel vcf files, I expected the CombineVariants output vcf file to retain the 'actual frequency' that I calculated, but the tool CombineVariants has reset the values to 0.5 and 1.0

In the following example, the AF was set to 0.500 by HaplotypeCaller, I have used the AD tag values 50,29 to calculate the 'actual frequency' as 29/(50+29) = 0.367. And then I have changed the value of AF to 0.367.

VCF record with original AF value:
1 881627 rs2272757 G A 457.77 . AC=1;AF=0.500;AN=2;BaseQRankSum=-1.582e+00;ClippingRankSum=-4.930e-01;DB;DP=79;ExcessHet=3.0103;FS=0.885;MLEAC=1;MLEAF=0.500;MQ=59.83;MQRankSum=1.97;QD=5.79;ReadPosRankSum=2.70;SOR=0.906 GT:AD:DP:GQ:PL 0/1:50,29:79:99:486,0,1069

VCF record after re-calculating actual AF value:
1 881627 rs2272757 G A 457.77 . AC=1;AF=0.367;AN=2;BaseQRankSum=-1.582e+00;ClippingRankSum=-4.930e-01;DB;DP=79;ExcessHet=3.0103;FS=0.885;MLEAC=1;MLEAF=0.500;MQ=59.83;MQRankSum=1.97;QD=5.79;ReadPosRankSum=2.70;SOR=0.906 GT:AD:DP:GQ:PL 0/1:50,29:79:99:486,0,1069

VCF record after merging files with CombineVariants :
1 881627 rs2272757 G A 457.77 . AC=1;AF=0.500;AN=2;BaseQRankSum=-1.582e+00;ClippingRankSum=-4.930e-01;DB;DP=79;ExcessHet=3.0103;FS=0.885;MLEAC=1;MLEAF=0.500;MQ=59.83;MQRankSum=1.97;QD=5.79;ReadPosRankSum=2.70;SOR=0.906 GT:AD:DP:GQ:PL 0/1:50,29:79:99:486,0,1069

Is this the expected behavior of tool 'CombineVariants' that it resets AF value to 0.5 and 1.0 ?

Is there an option in CombineVariants that I can use to make this tool not reset the AF values ?

Thanks,
mglclinical

Issue · Github
by Sheila

Issue Number
1776
State
closed
Last Updated
Assignee
Array
Milestone
Array
Closed By
chandrans

Best Answer

Answers

  • mglclinicalmglclinical USAMember

    @Sheila , than you for pointing me to VariantAnnotator

Sign In or Register to comment.