powered by
plotNegBinomFit(x, breaks = c(-0.5:15.5, 1e+08), bar.col=rainbow(2), addLegend=TRUE, legend.names=c("data","background"), ...)
barplot