Motivated by the problem of uncertainty quantification in self-organization, we study a spatially extended Sel’kov–Strogatz model of glycolysis. A variety of coexisting patterns induced by the Turing instability is studied in the parametric zones where the original local model without diffusion exhibits stable equilibria or self-oscillations. A phenomenon of the suppression of homogeneous self-oscillations by diffusion with formation of stationary non-homogeneous patterns-attractors is revealed. To quantify the uncertainty in the number and modality of patterns-attractors and to perform an advanced parametric analysis, we use the spectral coefficients technique and Shannon entropy.