[BioC] Computing the coverage on a bamfile with a huge amount of seqlevels

Stefanie Tauber stefanie.tauber at univie.ac.at
Thu Oct 17 16:37:11 CEST 2013

What is the best way to compute the coverage on a ReadGAlignments object
when I have many, many (about 50 000) seqlevels?

Just computing coverage(object) is terribly slow due to the huge amount of seqlevels.

At the moment, I am subsetting the ReadGAlignments object for each seqlevel,
 I omit the redundant seqlevels and calculate the coverage.
This is fast if you do it for a few seqlevels.
But if one really would like to compute the coverage for each seqlevel individually,
this really takes a long time..

I would appreciate any comments!


More information about the Bioconductor mailing list