This is a test version of Biostars. For the public version, visit https://www.biostars.org.
liftOver VCF from hg38 to hg19

Hi,

I am trying to lift the New York Genome Centre 1000 Genomes 30x on GRCh38 to GRCh37. As this file is a vcf I am trying to use picard's liftOverVCF tool. However, using the UCSC hg38 to hg19 chain file and the above mentioned vcf I get the following output:

INFO    2022-10-27 13:37:40 LiftoverVcf Loading up the target reference genome.
INFO    2022-10-27 13:37:41 LiftoverVcf Lifting variants over and sorting (not yet writing the output file.)
ERROR   2022-10-27 13:37:41 LiftoverVcf Encountered a contig, chr2 that is not part of the target reference.

I have seen multiple instances of this error online but I cannot find a single solution! I don't know if it is an issue with the chain file I am using being not exactly compatible to the New York Genome Centre release as it includes decoys and HLA, or whether it is a problem with the fact that chr21 in build 38 is just referred to as 21 in build 37. I tried to rename the header according to the hg19 reference, however, the latest version of bcftools reheader doesn't seem to support this function.

genetics genome liftover

0 answers

No answers yet.

Log in to answer this question.