Probabilistic modeling of bifurcations in single-cell gene expression data using a Bayesian mixture of factor analyzers