Analysis of distributional variation through multi-scale Beta-Binomial modeling