This is a test version of Biostars. For the public version, visit https://www.biostars.org.
Allele flip in GATK

I have this variant (rs17863783) from Visscher et al. 2013 , where they reported alleles as A,C for this variant.

enter image description here

However, when I checked this in my VCF, it has G,T. I was wondering if this is due to strand bias or allele flip in GATK and should be considered the same as A,C. Can someone please help me understand this?

 chr2       233693631       chr2:233693631:G:T      G       T
gatk vcf
• 1,749 views • link
updated by benformatics benformatics • written by MAPK2 • MAPK2

Hi! Doing a quick search it looks like the genotype of that SNP is G,T. Have you checked the Visscher et al. 2013 publication to see how they carried out the genotyping? Perhaps they employed a methodology based on complementary sequence? I'd take that direction; I doubt strand bias has anything to do with this :).

• 0 views • link
written by iraun iraun

iraun Thank you. That seems a valid explanation. This is how they genotyped "DNA samples were genotyped using a custom 96-plex Illumina Veracode GoldenGate SNP genotyping assay according to manufacturer’s instructions (Illumina, San Diego, CA). This assay included an additional 63 non-study SNPs used for quality control purposes only". My guess is that GATK uses StrandOddsRatio or Fisher's Exact Test to evaluate strand bias and therefore the genotype by GATK is more valid.

• 0 views • link
written by MAPK2 • MAPK2

You would need to check the original BAM but it could be that they called the variant on the opposite strand.

• 0 views • link
written by benformatics benformatics

0 answers

No answers yet.

Log in to answer this question.

More posts like this