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

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
• 1,217 views • link
written by heureuse • heureuse

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 ?

• 0 views • link
written by Pierre Lindenbaum Pierre Lindenbaum

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

• 0 views • link
written by heureuse • heureuse

0 answers

No answers yet.

Log in to answer this question.

More posts like this