Based on epidemiological parameter data, probe into a stochastically perturbed dominant variant of the COVID-19 pandemic model.

During the COVID-19 pandemic, the SARS-CoV-2 gene mutation has been rapidly emerging and spreading all over the world. Experts worldwide regularly monitor genetic mutations and variants through genome-sequence-based surveillance, laboratory testing, outbreak investigation, and epidemiological probing. Clinical pathologists and medical laboratory scientists prefer developing or endorsing COVID-19 vaccines with a broader immune response involving various antibodies and cells to protect against mutations or new variants. Randomness plays an enormous role in pathology and epidemiology. Hence, based on epidemiological parameter data, we construct and probe a stochastically perturbed dominant variant of the coronavirus epidemic model with three nonlinear saturated incidence rates. We reveal the existence of a unique global positive solution to the constructed stochastic COVID-19 model. The Lyapunov function method is used to determine the presence of a stationary distribution of positive solutions. We derive sufficient conditions for the coronavirus to be eradicated. Eventually, numerical simulations validate the effectiveness of our theoretical outcomes.