Abstract
COVID-19 has become a world wide pandemic since its first appearance at the end of the year 2019. Although some vaccines have already been announced, a new mutant version has been reported in UK. We certainly should be more careful and make further investigations to the virus spread and dynamics. This work investigates dynamics in Lotka-Volterra based Models of COVID-19. The proposed models involve fractional derivatives which provide more adequacy and realistic description of the natural phenomena arising from such models. Existence and boundedness of non-negative solution of the fractional model is proved. Local stability is also discussed based on Matignon's stability conditions. Numerical results show that the fractional parameter has effect on flattening the curves of the coexistence steady state. This interesting foundation might be used among the public health strategies to control the spread of COVID-19 and its mutated versions.