Stochastic modelling and analysis of metabolic heterogeneity in single cells