to do hypothesis testing and see if two negative binomial distributions are different, I suggest generalized linear models in R, specifically glm.nb
