Question: Crossmap liftover: Chromosome changed between hg38 and hg19???
gravatar for skim
2.5 years ago by
skim30 wrote:

I have done a liftover on a hg38 bed file to a hg19 bed file on Crossmap(Link) and found out that multiple queries were changed from chromosome 22 to chromosome 14. Should I just filter this or did I do something wrong?

Thank you very much.

Data: hg38

chr22   15690509    15690536    POTEH   28.8238274624   -
chr22   15698714    15698741    POTEH   69.8367597081   -
chr22   15698718    15698745    POTEH   63.510218446    +
chr22   15698719    15698746    POTEH   64.7710373827   +
chr22   15698720    15698747    POTEH   71.7040268333   -
chr22   15698720    15698747    POTEH   64.4982928966   +


chr14   19553737    19553764    POTEH   28.8238274624   -
chr14   19562046    19562073    POTEH   69.8367597081   -
chr14   19562052    19562079    POTEH   71.7040268333   -

Some of these were already filtered based on 27bp length. All of the data are 27bp. Chain file is this:

sequencing liftover next-gen • 1.8k views
ADD COMMENTlink modified 2.5 years ago by Devon Ryan94k • written 2.5 years ago by skim30
gravatar for Devon Ryan
2.5 years ago by
Devon Ryan94k
Freiburg, Germany
Devon Ryan94k wrote:

You didn't do anything wrong, if you zgrep chr22 hg38ToHg19.over.chain.gz | grep chr14 you'll see there's a bunch of sequence originally in chr14 that was moved to chr22. If you go over to the gencode website and scroll through the GRCh37 assembly issues, you'll see a LOT of them were fixed with the GRCh38 release.

ADD COMMENTlink written 2.5 years ago by Devon Ryan94k
Please log in to add an answer.


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