SARS-CoV-2 is a highly virulent and deadly RNA virus causing the Covid-19 pandemic and several deaths across the world. The pandemic is so fast that any concrete theory of sudden widespread of this disease is still not known. In this work, we studied and analyzed a large number of publicly available SARS-CoV-2 genomes across the world using the multifractal approach. The mutation events in the isolates obey the Markov process and exhibit very high mutational rates, which occur in six specific genes and highest in orf1ab gene, leading to virulent nature. f ( $\alpha$ ) analysis indicated that the isolates are highly asymmetric (left-skewed), revealing the richness of complexity and dominance by large fluctuations in genome structure organization. The values of Hq and Dq are found to be significantly large, showing heterogeneous genome structure self-organization, strong positive correlation in organizing the isolates and quite sensitive to fluctuations in and around it. We then present multiple-isolates hosts-virus interaction models and derived Price equation for the model. The phase plane analysis of the model showed asymptotic stability type bifurcation. The competition among the mutant isolates drives the trade-off of the dominant mutant isolates, otherwise confined to the present hosts. ### Competing Interest Statement The authors have declared no competing interest.