Noticed it would be cool to run ShortRead::qa() on a ShortReadQ that was sampled from a BAM file. But readAligned seems to only support the old "Solexa" objects and character paths, so I can't just pass it a BamSampler. This would fit nicely in a certain paper on scalable computing in genomics. And probably would be useful in general.
Thanks, Michael [[alternative HTML version deleted]] _______________________________________________ Bioc-devel@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/bioc-devel