rtracklayer: how to view coordinates that failed liftover?
1
1
Entering edit mode
chang02_23 ▴ 30
@chang02_23-7435
Last seen 6.6 years ago
United States

I am ilfting over SNV calls from hg38 to hg19 using rtracklayer. The liftover result has less number of lines that my input. How do i figure out which SNV calls failed liftover?

 

rtracklayer • 1.2k views
ADD COMMENT
3
Entering edit mode
@michael-lawrence-3846
Last seen 3.0 years ago
United States

liftOver() yields a GRangesList, with one element per input range. I'm not sure what you mean by "lines" but I guess you are unlisting the result and writing it to a file. Instead of unlisting, you can just look at the lengths() of the return value to determine the degree of mapping for each input range.

ADD COMMENT
0
Entering edit mode

Aha! I unlist the result so I couldn't figure out how to map each input range. Thank you.

ADD REPLY

Login before adding your answer.

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