Freebayes variant calling understanding
0
0
Entering edit mode
2.4 years ago
User000 ▴ 580

Hello,

I am using freebayes to do variant calling in my polyploid species (with ploidy 2 anyway). Before filtering I was looking at the bam file alignment and a vcf file. Could you, please explain what does alleles 2/2, 0/2, 1/2 means? 0/0 is reference allele, 1/1 alternative and 0/1 heterozygote? Because if I filter for biallelic alleles, 2/2, 0/2, 1/2 will be removed, although in bam file I can see a very nice high quality SNP which in vcf results as 2/2.

Moreover, in the situation like this, even if there is only one reference allele with a high mapping quality freebayes would result in a 0/1 call. Any suggestions/explanations? image here

freebayes SNP • 1.3k views
1
Entering edit mode

In some loci you have more than just two alleles and for those loci the caller use notations like 2,3, etc. This variants are multiallelic and were filtered out when you tried to leave only biallelic ones.

0
Entering edit mode

that sounds reasonable, but when I have a look at the loci in bam file, I can see only reference allele and alternative allele, something like this: image here

0
Entering edit mode

Ah I see, so you mean the allele is different from other genotypes!

0
Entering edit mode

Try to check this position in different samples, perhaps you'll find several alternative alleles for that position.

1
Entering edit mode

Please paste the actual entries from the VCF file. I believe that Denis is correct in that these are likely multi-allelic sites. You can split these via the following command:

bcftools norm -m-any


0
Entering edit mode

thank you Kevin, I still haven't found a solution for the second part of my question, but if you want to close the thread, go ahead.

0
Entering edit mode

Oh no, I meant 'close off' in the non-literal sense, i.e., find a solution to the problem. Can you possibly paste here these records from your VCF?

Traffic: 1415 users visited in the last hour
FAQ
API
Stats

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