Entering edit mode
7.5 years ago
musaevdamir
•
0
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