The stated mean and variance identify as a gamma distribution with shape and scale . A Poisson-gamma mixture is negative binomial, with
When , this is exactly the constant-shape variance model fitted by the negative binomial regression, so the most appropriate printed p-value is from model3.