In this manuscript, a qualitative analysis of the mathematical model of novel coronavirus ($\mathrm{COVID}$-19) involving anew devised fractal-fractional operator in the Caputo sense having the fractional-order $\mathsf{q}$ and the fractal dimension $\mathsf{p}$ is considered. The concerned model is composed of eight compartments: susceptible, exposed, infected, super-spreaders, asymptomatic, hospitalized, recovery and fatality. When, choosing the fractal order one we obtain fractional order, and when choosing the fractional order one a fractal system is obtained. Considering both the operators together we present a model with fractal-fractional. Under the new derivative the existence and uniqueness of the solution for considered model are proved using Schaefer’s and Banach type fixed point approaches. Additionally, with the help of nonlinear functional analysis, the condition for Ulam’s type of stability of the solution to the considered model is established. For numerical simulation of proposed model, a fractional type of two-step Lagrange polynomial known as fractional Adams-Bashforth $\left(\mathrm{AB}\right)$ method is applied to simulate the results. At last, the results are tested with real data from $\mathrm{COVID}$-19 outbreak in Wuhan City, Hubei Province of China from 4 January to 9 March 2020, taken from a source (Ndaïrou, 2020). The Numerical results are presented in terms of graphs for different fractional-order $\mathsf{q}$ and fractal dimensions $\mathsf{p}$ to describe the transmission dynamics of disease infection.