The CME governs probability distribution P(n,t) of molecule count n. For ∅ →λ X →μ ∅, the steady state is Poisson(λ/μ). Gillespie simulation samples exact stochastic trajectories.