It seems that type = c("equal", "within") is not implemented for findOverlaps for GenomicRanges, despite being part of the IRanges interface. To the casual eye, this should be extremely easy to fix (just extend the type argument), but there may be some problems with the .findOverlaps.circle. Kasper