VCF failed liftover, incorrect alternate allele
0
0
Entering edit mode
7.5 years ago

Hi all,

I lifted over hg38 vcf to hg19 one. However, some positions were corrupted - for some genes hg38 had sense strand while hg19 had anti-sense strand. So, after liftover the reference base was changed (e.g. if in chr1:12345 you had base A in hg38, it changed to T in hg19), but the alternate allele was not changed accordingly (e.g. if the alternate allele was C in hg38, it stayed C in hg19 instead of being changed to G).

The command below was run:

python CrossMap.py vcf /chain_file/hg38ToHg19.over.chain /projects2/vcfs/Mito_4208hg38.vcf /projects/human_reference/ucsc.hg19.fasta /projects2/vcfs/Mito_4208Hg19.vcf

Do you guys know what might have gone wrong?

Thanks, Damir

genome next-gen alignment sequence • 1.5k views
ADD COMMENT

Login before adding your answer.

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