Question: VCF failed liftover, incorrect alternate allele
gravatar for musaevdamir
2.5 years ago by
musaevdamir0 wrote:

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 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

ADD COMMENTlink written 2.5 years ago by musaevdamir0
Please log in to add an answer.


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