[BioC] "islands" and number of reads in the island in package chipseq/ShortRead
Nora Rieber
n.rieber at dkfz-heidelberg.de
Tue Dec 1 13:17:25 CET 2009
Dear all,
I'm looking for a fast way of computing how many reads map into a
certain location (e.g. between position 300 and 400 on chromosome 14).
Is there any function to compute this in the packages chipseq or
ShortRead? After reading the vignette for chipseq, I thought
viewSums(islands) computes the number of reads, but it seems to compute
the sum of the coverage.
I've written my own code to do this but it is pretty slow and I wonder
if there's a better way.
Thanks,
Nora
More information about the Bioconductor
mailing list