liftover hg19 to hg38 fragmented
0
0
Entering edit mode
10 months ago
heureuse • 0

Good morning, I performed a liftover from hg19 to hg38. The result is fragmented. I had 91 lines in hg19 and I end up with 7003 in hg38. I have a script to do the merge in 1 nucleotide steps but I'm not sure if that's what to do. Has anyone been in this situation before?

Thank you

hg38 chunk hg19 liftover • 464 views
ADD COMMENT
0
Entering edit mode

I have a script to do the merge in 1 nucleotide steps

it's not clear. What's the input ? what are the duplicated lines ? what's the script ? don't you fall in multiple-mapped areas ?

ADD REPLY
0
Entering edit mode

I use crossmap.py to do my liftover. The input is a bed file hg19 and the chain file (hg19 to hg38). The lines are not duplicated but fragmented. For example for one entry in hg19, I have 6 entry in the lifted file

ADD REPLY

Login before adding your answer.

Traffic: 1764 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6