Stochastic virus dynamics modeled by a system of stochastic differential equations with Beddington-DeAngelis functional response and driven by white noise is investigated. The global existence of positive solutions and the existence of stationary distribution are proved. Upper and lower bounds of the pathwise and asymptotic moments for the positive solutions are sharply estimated. The absorbing property in time average is shown and the moment Lyapunov exponents are proved to be nonpositive.