Thanks Ryan; this is very helpful. It's too bad that there isn't a GLM/Bayesian approach (beyond pairwise comparisons) that can model the uncertainty from multi-mapped reads. -- output of sessionInfo(): -- Sent via the guest posting facility at bioconductor.org.