This fits my purpose best.

>  gr[countOverlaps(gr, gr) <= 1L]
GRanges with 3 ranges and 0 metadata columns:
      seqnames     ranges strand
         <Rle>  <IRanges>  <Rle>
  [1]     chr1 [  4,  10]      *
  [2]     chr1 [ 12,  19]      *
  [3]     chr1 [105, 200]      *
  ---
  seqlengths:
   chr1
     NA

Could you recommend some further help for countOverlaps. Just because I did
not understand the syntax.

Thanks
Hermann

	[[alternative HTML version deleted]]

