Ref and Alt switch after imputation
0
0
Entering edit mode
8.2 years ago
sxl919 ▴ 10

I don't know if anyone have experienced this problem before. After having my original vcf from GATK, I impute my vcf using Shapeit2 and minimac3. However, I got something strange.

Original vcf:

2 25467408 . C T 107.66 PASS .
GT:AB:AD:DP:FT:GQ:MLPSAC:MLPSAF:MQ0:PL 0/0:.:192,0:44:PASS:99:0:0.00:0:0,105,1575

After imputed vcf (Shapeit2)

2 25467408 2:25467408 T C .
PASS;GENOTYPED_ONLY AF=0.99926;MAF=0.00074 GT:DS 1|1:2.000
1|1:2.000 1|1:2.000

As you can see the Ref and Alt for some variants have switch, making it from no call all the way to calls everywhere. Why does this happen? Any way to fix it?

Thank you

vcf imputation Ref Alt • 2.0k views
ADD COMMENT
0
Entering edit mode

For hg38, at that position, ref is T.

ADD REPLY
0
Entering edit mode

This can occur depending on which genotyping platform was used, because these genotype different SNPs.

Also read this to gain some enlightenment about the 'reference' genomes that we use: A: Alternate nucleotide is more frequent than reference nucleotide. OMG I'm dizzy.

ADD REPLY

Login before adding your answer.

Traffic: 2088 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