Question: how much positions overlapped for each pair of the pairs obtained by findoverlap.
0
3.5 years ago by
vinod.acear30
India
vinod.acear30 wrote:

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.

modified 3.5 years ago by Hervé Pagès ♦♦ 13k • written 3.5 years ago by vinod.acear30
Answer: how much positions overlapped for each pair of the pairs obtained by findover
3
3.5 years ago by
Hervé Pagès ♦♦ 13k
United States
Hervé Pagès ♦♦ 13k wrote:

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.