Pošalji zapis e-poštom: Single-gene negative binomial regression models for RNA-Seq data with higher-order asymptotic inference