Substantial variability in single-cell RNA-seq data stems from diverse sources including gene expression heterogeneity, cell cycle dynamics, and technical noise.Many of the single-cell RNA transcriptome data show expression bursts derived from the mixture of transcriptional on/off status.We hereby develop hierarchical Bayesian method to solve such bimodal signal.The inferred transcriptional status by this method reveals the inherent correlation of transcriptions in a co-bursting fashion.The co-burst of genes in single cells brings insight of into the cell population under interrogation.