Entering edit mode
8.7 years ago
tabone.olivier
▴
20
Dear all,
I transformed genome coordinates from hg19 to coordinates of GRCh38 using the function liftover()
from rtracklayer package in R. I used the hg19ToHg38.over.chain
file from UCSC. But it only contains chromosome from primary assembly (i.e., chr1, chr2, ..., chrY). It doesn't work for element coming from unplaced, unlocalized and alternative portions of chromosomes (for example chr6_cox_hap2
, or GL000191.1
, ...)
So my question is: Is there a way to find correspondence between coordinates of alternative chromosomes from hg19 and the new ones on GRCh38?
Thank you,
Olivier
as far as I know these elements do not exist anymore in hg38.
Thank you for your answer. Ok they don't exist but i am still surprised there is no correspondance between unlocalized scaffolds (known chromosome, but unknown position within chromosome) from hg19 and the "good" position within chromosome in hg38.
Was this resolved somewhere else besides this post? I am having the same issue but for lifting over variants ("LiftoverVcf Encountered a contig, 2_KI270773v1_alt that is not part of the target reference."). I am using the latest version of Picard, and the proper liftover file. Is the best way forward to just drop the hg19 variants that are on these contigs that are not part of hg38? Any help on this would be appreciated.
I am also very curious. But looking at all the time passed, I am assuming everyone moved on or silently found a solution? I have the same issue :-P