how much positions overlapped for each pair of the pairs obtained by findoverlap.
1
0
Entering edit mode
vinod.acear ▴ 50
@vinodacear-8884
Last seen 21 months ago
India

I have got overlapping pairs list using findoverlap function on two granges, i want to know  how much positions overlapped for each  pair of  the pairs obtained by findoverlap.

3
Entering edit mode
@herve-pages-1542
Last seen 21 hours ago
Seattle, WA, United States

Hi,

So you have a Hits object that you got by calling findOverlaps() on your GRanges objects query and subject, e.g. by doing something like:

hits <- findOverlaps(query, subject, ...)

Now you can get the corresponding overlap as a range for each pair in hits with:

overlap <- pintersect(query[queryHits(hits)],
subject[subjectHits(hits)])

and the number of positions in each overlap with:

width(overlap)

Cheers,

H.

0
Entering edit mode

Thanks herve