Need help in IRanges

[Kumar10]

Hi
I am using IRanges for finding overlaps in genomic ranges. so I would like to know whether there is any function or any method by which i can determine that i will get a result of overlaps which are above the threshold level which i determine. Thank you in advance.

[Dario1984]

Hello,

You should try using the package GenomicRanges. Click here for its home page. It's much easier to handle chromosomal data with it. It's made by the same people who made IRanges. In this example I show you how to find all regions that have 2 or more reads.

library(GenomicRanges)
g <- GRanges(c("chr1", "chr1"), IRanges(start = c(5, 10), end = c(15, 20)), strand = '+')
g
[COLOR="Blue"]GRanges with 2 ranges and 0 elementMetadata values:
      seqnames    ranges strand
         <Rle> <IRanges>  <Rle>
  [1]     chr1  [ 5, 15]      +
  [2]     chr1  [10, 20]      +
  ---
  seqlengths:
   chr1
     NA[/COLOR]
slice(coverage(g), 2)[["chr1"]]
[COLOR="Blue"]Views on a 20-length Rle subject

views:
    start end width
[1]    10  15     6 [2 2 2 2 2 2]
[/COLOR]