In this paper we extend the deterministic model for the epidemics induced by virulent phages on bacteria in marine environment introduced by Beretta and Kuang [Math. Biosci. 149 (1998) 57], allowing random fluctuations around the positive equilibrium. The stochastic stability properties of the model are investigated both analytically and numerically suggesting that the deterministic model is robust with respect to stochastic perturbations.