powered by
gcContentHist(object, xlab="GC content", ylab="Number of sequences", col="firebrick1", breaks=50, ...)