[BioC] countGenomicOverlaps output

Mete Civelek mcivelek at mednet.ucla.edu
Fri Jul 22 20:55:15 CEST 2011

Hi All,

I want countGenomicOverlaps to output a count of uniquely mapping reads
within a genomic feature. Will setting the resolution parameter to 'none'
allow countGenomicOverlaps to ignore reads which map to multiple locations
in the genome? If so, countGenomicOverlaps doesn't behave the way I expect
it to. I am using the Bioconductor GenomicRanges package version 1.4.6.


subj = GRangesList(feature1=GRanges(seq='1', IRanges(10,30), strand='+'))
qry = GRangesList(read1=GRanges(seq='1', IRanges(c(10,60,100),c(20,70,110)),
countGenomicOverlaps(qry, subj, resolution='none')

I would have expected the hit count to be 0 but instead it reports it as
1/3. Am I using this function correctly?



IMPORTANT WARNING:  This email (and any attachments) is ...{{dropped:9}}

More information about the Bioconductor mailing list