GATK-Allele frequency
1
3
Entering edit mode
13 months ago

Hi Guys,

I am running GATK on bam file for variant calling. In the output file, I noticed that the Allele frequency is computed as 0.5 and 1.00. What may be the reason for this? Is it calculated correctly?

VCF Allele GATK frequency • 1.5k views
ADD COMMENT
0
Entering edit mode

The information you provide is insufficient to get the right answer

is your VCF file multi-sample?

your variants are biallelic?

ADD REPLY
4
Entering edit mode
13 months ago

do you have only one sample ? the genotype is either heterozygous AN=2 AC=1 then AF=0.5 , or homozygous on ALT allele: AN=2, AC=2 then AF=1.0.

ADD COMMENT
0
Entering edit mode

Yes, there is only one sample.

ADD REPLY
0
Entering edit mode

Hi,

How this allele frequency is calculated? I did not get it. Can you just brief me and send me a link to read about it. ?

Thanks

ADD REPLY
0
Entering edit mode

Calculation of VAF (variant allele frequency)

freq(a) =( sum(samples_with_geno_aa x 2) + sum(samples_with_geno_Aa)) / (samples x 2)

ADD REPLY
0
Entering edit mode

Hi,

Thank you so much. I want to confirm whether Allele frequencies were calculated correctly for one sample. Also, If I have more than one samples, how to proceed for variant calling. I mean we have to calculate separately for each sample ??

Thanks

ADD REPLY
0
Entering edit mode

Hi Pierre,

I am new to this analysis. Can you please guide me how AF is coming as 0.500 and 1.000? I will really appreciate your support.

Thanks

ADD REPLY
0
Entering edit mode

I already answered.

ADD REPLY
0
Entering edit mode

Hi Pierre,

Thank you so much. The link you provided is very helpful which includes two aspects as Joint Variant Calling (from bam) and Joint Genotyping (from vcf). I found command for Joint genotyping https://biohpc.cornell.edu/lab/doc/Variant_exercise.pdf.

However, I could not find any command for variant calling from multiple bam files. Can you please help in this regard?

ADD REPLY
0
Entering edit mode

If you are looking to add VAF values, consider using:

bcftools +fill-tags in.bcf -Ob -o out.bcf -- -t FORMAT/VAF

See bcftools plugin docs here:

https://samtools.github.io/bcftools/howtos/plugin.fill-tags.html

And further discussion of MAF and VAF here:

Calculation of VAF (variant allele frequency)

ADD REPLY

Login before adding your answer.

Traffic: 1484 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6