coverage(), form Views() on the coverage for your genes and then viewSums().
On Fri, Jul 22, 2011 at 2:36 PM, Kunbin Qu <k...@genomichealth.com> wrote: > Hi, > > If I want to count the coverage base by base, not read by read, can I still > use countOverlaps? I have a human transcriptome done, and would like to > count the coverage for each gene based on the mapping. As there are some > reads mapped across the junctions (ie, one read is splitted into two > portions), or partially mapped into the introns, it would be good to count > the coverage by base pair, instead of by read number? Could anybody suggest > a way doing that? Thanks. > > -Kunbin > > > > ______________________________________________________________________ > The contents of this electronic message, including any attachments, are > intended only for the use of the individual or entity to which they are > addressed and may contain confidential information. If you are not the > intended recipient, you are hereby notified that any use, dissemination, > distribution, or copying of this message or any attachment is strictly > prohibited. If you have received this transmission in error, please send an > e-mail to postmas...@genomichealth.com and delete this message, along with > any attachments, from your computer. > [[alternative HTML version deleted]] > > _______________________________________________ > Bioc-sig-sequencing mailing list > Bioc-sig-sequencing@r-project.org > https://stat.ethz.ch/mailman/listinfo/bioc-sig-sequencing > [[alternative HTML version deleted]] _______________________________________________ Bioc-sig-sequencing mailing list Bioc-sig-sequencing@r-project.org https://stat.ethz.ch/mailman/listinfo/bioc-sig-sequencing