CrossMap icon indicating copy to clipboard operation
CrossMap copied to clipboard

Issue with liftover of vcf file from T2T to hg38

Open Taghrid-M opened this issue 1 year ago • 3 comments

Dear Crossmap Maintainer,

I am reaching out to report an issue I've encountered while using Crossmap. I was trying to liftover my VCF file from T2T to hg38, but found a significant number of variants were rejected during the process. The total number of rejected variants was around 1,855,261, which is quite a significant figure.

To perform the liftover, I've used both the chain file from github.com/marbl/CHM13 and the UCSC chain file. Unfortunately, the issue persisted with both chain files.

Here is the command line: CrossMap.py vcf $chain_file $input_vcf $GRCh38_REF $output_vcf

I was wondering if you could help me understand why this is happening, and if there's a possible solution or workaround. Perhaps I'm missing something, or there could be an issue with the chain files I've used?

Thank you in advance.

Taghrid

Taghrid-M avatar Jun 06 '23 18:06 Taghrid-M