In this article, we study a mathematical system modelling the dynamic of the collective behaviour of oxygen-driven swimming bacteria in an aquatic fluid flowing in a two dimensional bounded domain under stochastic perturbation. This model can be seen as a stochastic version of Chemotaxis–Navier–Stokes model. We prove the existence of a unique (probabilistic) strong solution. In addition, we establish some properties of the strong solution. More precisely, we prove that the unique solution is positive and satisfies the mass conservation property and an energy inequality.