unplaced, unlocalized, ... scaffolds from hg19 to GRCh38
0
2
Entering edit mode
8.2 years ago

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

Assembly • 3.4k views
ADD COMMENT
0
Entering edit mode

as far as I know these elements do not exist anymore in hg38.

ADD REPLY
0
Entering edit mode

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.

ADD REPLY
0
Entering edit mode

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.

ADD REPLY
0
Entering edit mode

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

ADD REPLY

Login before adding your answer.

Traffic: 1809 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