liftOverPlink output 0 bytes
Hi, I have the plink bfiles for SNP genotyping map to build38 and I want to convert to build 37. Here is my command
$ python liftOverPlink.py --map SHARE_MESA_AA_chr9_position.map --out lifted --chain hg38ToHg19.over.chain.gz
Converting MAP file to UCSC BED file...
SUCC: map->bed succ
Lifting BED file...
Reading liftover chains
q end mismatch 12611573 vs 249240621 line 87 of hg38ToHg19.over.chain.gz
SUCC: liftBed succ
Converting lifted BED file back to MAP...
SUCC: bed->map succ
cleaning up BED files.
..
However, the output lifted.map is 0 bytes. I don't know where is wrong. Anyone has some suggestion for my problem? I am very appreciated your help!
• 1,906 views
•
link
updated
by
GenoMax
•
written
by
honghong0220 •
•
written
by
honghong0220 •
0 answers
No answers yet.
Log in to answer this question.
More posts like this
-
High ref mismatch rate after liftOver from 23andme hg19 to hg38
written by berndmann •
-
gatk LiftoverVCF finishes with no output VCF nor error.
written by jpablo •
-
Liftover and Then Update RSIDs
written by hi.there •
-
Tutorial: Automated dbSNP lookup by rsID position, plus genome build liftover
written by Kevin Blighe
-
plink Map change location from hg18 to hg19
written by enci.xue96 •
Try to confirm that your co-ordinates are indeed 'mappable' from hg38 to hg19 via the UCSC liftover tool: https://genome.ucsc.edu/cgi-bin/hgLiftOver
Many regions in hg38 did not exist in hg19.
Thanks Kevin for your suggestion. Already solve this. The downloaded chain file has problem.