This manuscript is devoted to investigating a fractional-order mathematical model of COVID-19. The corresponding derivative is taken in Caputo sense with power-law of fractional order μ and fractal dimension χ. We give some detailed analysis on the existence and uniqueness of the solution to the proposed problem. Furthermore, some results regarding basic reproduction number and stability are given. For the proposed theoretical analysis, we use fixed point theory while for numerical analysis fractional Adams–Bashforth iterative techniques are utilized. Using our numerical scheme is verified by using some real values of the parameters to plot the approximate solution to the considered model. Graphical presentations corresponding to different values of fractional order and fractal dimensions are given. Moreover, we provide some information regarding the real data of Saudi Arabia from 1 March 2020 till 22 April 2021, then calculated the fatality rates by utilizing the SPSS, Eviews and Expert Modeler procedure. We also built forecasts of infection for the period 23 April 2021 to 30 May 2021, with 95% confidence.