The purpose of this paper is to propose and investigate a stochastic SIQR epidemic model with saturated incidence rate. Firstly, we give some conditions to guarantee the stochastic SIQR epidemic model has a unique global positive solution. Then we verify that the disease in this model will die out exponentially if R s 0 < 1, while the disease will be persistent in the mean if R s 0 > 1. Moreover, by constructing suitable Lyapunov functions, we establish some sufficient conditions for the existence of an ergodic stationary distribution for the model. Finally, we provide some numerical simulations to illustrate the analytical results.