In this paper, we first prove the existence and uniqueness of the solution to a variable-order Caputo–Fabrizio fractional stochastic differential equation driven by a multiplicative white noise, which describes random phenomena with non-local effects and non-singular kernels. The Euler–Maruyama scheme is extended to develop the Euler–Maruyama method, and the strong convergence of the proposed method is demonstrated. The main difference between our work and the existing literature is the fact that our assumptions on the nonlinear external forces are those of one-sided Lipschitz conditions on both the drift and the nonlinear intensity of the noise as well as the proofs of the higher integrability of the solution and the approximating sequence. Finally, to validate the numerical approach, current results from the numerical implementation are presented to test the efficiency of the scheme used in order to substantiate the theoretical analysis.