“…The telegraph model is theoretically extensively analysed, and it is known that it admits a Beta-Poisson distribution at steady-state ( Peccoud and Ycart 1995 , Kepler and Elston 2001 , Raj et al 2006 , Shahrezaei and Swain 2008a , Kim and Marioni 2013 ). At the bursty limit of transcription, the solution of the telegraph model can be approximated as a negative-binomial distribution characterized by the burst size and burst frequency ( Raj et al 2006 , Shahrezaei and Swain 2008a , Kumar et al 2015 , Amrhein et al 2019 , Thomas 2020 ). Moreover, the negative binomial (NB) distribution is a versatile over-dispersed distribution that is commonly used in bulk and scRNA-seq studies to model gene expression capturing both biological and technical dispersion ( Anders and Huber 2012 , Love et al 2014 , Tang et al 2020 , Svensson 2020 ).…”