Question: Liftover errors with CrossMap and Picard
0
gravatar for rshoobs
21 months ago by
rshoobs10
rshoobs10 wrote:

Hi all,

I am attempting to use CrossMap to convert my dosage files output from the sanger imputation service from hg19/37 to hg38. Here is the command I am using

python3 ~/software/CrossMap-0.3.1/bin/CrossMap.py vcf ~/Data/liftover/hg19ToHg38.over.chain.gz 22_withMAF.vcf.gz ~/Data/GenCode/GRCh38.primary_assembly.fa test_output_chr22_hg38.vcf

and I receive the following error part way through

KeyError: "sequence 'b'chr22_KI270879v1_alt'' not present"

I get similar errors when I attempt to run other chromosomes. Has anyone encountered this problem before?

Additional note: I have also attempted to LiftOver using the picard LiftoverVcf tool. When I do this happens

100.0000% of variants were not successfully lifted over and written to the output

Any help that can be given towards either tool would be appreciated!

snp rna-seq software error • 1.0k views
ADD COMMENTlink written 21 months ago by rshoobs10

I've been trying further and find that it works with a bed file that I've created from my vcf. Only issue is that It doesn't preserve any mapping so i can't coordinate any genotypes!

ADD REPLYlink written 21 months ago by rshoobs10

I have ever met the problem: sh: liftOver.hg19Tohg38: command not found, but have never met the problem you mentioned.

ADD REPLYlink modified 9 months ago • written 9 months ago by Shicheng Guo8.3k
Please log in to add an answer.

Help
Access

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