This study analyzes the fractional order SIRC epidemic model under stochastic fractional differential equations in the Caputo sense. This article describes Salmonella infection in animal herds. For the deterministic system, we explain solution positivity and boundness. We also prove that the fractional stochastic solution exists and is unique. Other criteria considered include non-negativity of solutions, local and global stability analyses, Hyers-Ulam stability analysis, and sensitivity analysis for the deterministic system. Furthermore, according to the truncated Ito-Taylor expansion, we apply a numerical method to solve the stochastic fractional SIRC epidemic model, namely the Milstein method, to solve the stochastic fractional SIRC epidemic model. A comparison of the approximation solution and the corresponding deterministic model for different sample paths shows the efficiency of the numerical method. In addition, graphs and error tables provide insight into numerical experiments' results. The stochastic nature of the model allows random fluctuations in Salmonella infection spread. By incorporating uncertainty into the model, we gain a more realistic understanding of the epidemic dynamics and can better evaluate the effectiveness of control measures. Additionally, the numerical method used to solve the stochastic fractional SIRC epidemic model provides valuable insights into the variability of the results. This enhances our ability to make informed decisions about managing and preventing bacterial infections in animal herds.