Hi,
I need to convert hg19 coordinates to hg38 (or opposite, doesn't matter). I tried to use liftOver from uscc, but without results.
My coordinate in hg38 2:230442937-230450255
In hg19: chr2:231307651-231314970
And after conversion in liftOver the results is: chr2:231307652-231314970
Where I made a mistake? I understand that that the one base is due to the different database, but why is difference in length?
And how is the easiest way to convert this?
Best, d.
1 answer
A chromosomal rearrangement could give you regions of different lengths.
One good way to run liftOver is to do it in both directions, find what maps in both directions, and use those intervals:
From
hg19tohg38, taking inputAand generating interval setA'liftOver(A) -> A`From set
A', go back tohg19via a second liftOver fromhg38tohg19, to get setA''liftOver(A`) -> A``Compare sets
AandA''via set operations, to determine which regions are identical inhg19space:bedopsElementOfExact(A, A``) -> A```Run a final liftOver on
A''', fromhg19tohg38:liftOver(A```) -> A````
Log in to answer this question.