Search
Question: Detect Sites with Minimum Gap in GRanges Object
0
7 months ago by
g.k0
g.k0 wrote:

Hello,

I have a GRange object and I am using the findOverlaps() function inorder to control the minOverlap and maxGap.

However I was wondering if there is a way to get the opposite? Is there a way to get the sites that have a  maxOverlap and the minGap?

 seqnames     ranges        strand   chr1              [106, 116]      +   chr1              [118, 120]      +   chr1              [124, 130]      +   chr2              [131, 140]      +

So if we chose a minGap of 2 , then the first 3 sites would be selected.

modified 7 months ago • written 7 months ago by g.k0

It would help to know more about your high-level use case. Why are you doing this?