Scalable Bayesian inversion with Poisson data