Recent works have stressed the important role that random mutations have in the development of cancer phenotype. We challenge this current view by means of bioinformatic data analysis and computational modelling approaches. Not all the mutations are equally important for the development of metastasis. The survival of cancer cells from the primary tumour site to the secondary seeding sites depends on the occurrence of very few driver mutations promoting oncogenic cell behaviours and on the order with which these mutations occur. We introduce a model in the framework of Cellular Automata to investigate the effects of metabolic mutations and mutation order on cancer stemness and tumour cell migration in bone metastasised breast cancers. The metabolism of the cancer cell is a key factor in its proliferation rate. Bioinformatics analysis on a cancer mutation database shows that metabolism-modifying alterations constitute an important class of key cancer mutations. Our approach models three types of mutations: drivers, the order of which is relevant for the dynamics, metabolic which support cancer growth and are estimated from existing databases, and non--driver mutations. Our results provide a quantitative basis of how the order of driver mutations and the metabolic mutations in different cancer clones could impact proliferation of therapy-resistant clonal populations and patient survival. Further mathematical modelling of the order of mutations is presented in terms of operators. We believe our work is novel because it quantifies two important factors in cancer spreading models: the order of driver mutations and the effects of metabolic mutations.