powered by
baseQualityHist(object, xlab="Quality score", ylab="Number of bases", col="firebrick1", breaks=40, ...)