Gravar-mail: Covariate-dependent negative binomial factor analysis of RNA sequencing data