Question: why and how multi allele gets reported during variant calling in vcf
0
gravatar for genie66
2.1 years ago by
genie6620
United States
genie6620 wrote:

This might be a very basic question for many here. With the basic understanding of inheritance, eventhough there is a possibility of multiple genotypes due to multi alleles, the resulting genotype can only have two alleles(paternal and maternal) that way after variant calling, an allele at one position can be homozygous or heterozygous. So there can be max of two alleles, but why do we see multiple alleles at a given position in single sample VCF. I am trying to understand the science behind this. Please help out. Thanks!

Eg:

chr5    127640782   .   AG  A,AA    .   .   .   GT:AD:DP    1/2:0,28,409:437
ADD COMMENTlink modified 2.1 years ago • written 2.1 years ago by genie6620

the resulting genotype can only have two alleles(paternal and maternal)

for one sample, yes. If there is more than one sample in the VCF for the same variant: sample 1 is A/T, sample2 is T/G ...

ADD REPLYlink written 2.1 years ago by Pierre Lindenbaum131k
2

To add to what Pierre said, even in situations where there is a single sample, multi-alleles appear due to the fact that NGS data has such high error rates. It is difficult for the in silico tools to deal with this messy data in a way that is faithful.

Indels should dominate the multi-allelic regions, too, as indels are inherently difficult to call for a variant caller.

ADD REPLYlink written 2.1 years ago by Kevin Blighe66k

Please don't ask the same question on multiple boards at the same time. If you already have a good answer on one board, it's a waste of people time to repeat it here.

ADD REPLYlink written 2.1 years ago by swbarnes28.9k

Some times, both maternal and paternal are alternate alleles if reference is not exact reference.

ADD REPLYlink written 2.1 years ago by cpad011214k
Please log in to add an answer.

Help
Access

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