We introduce a global analysis of relativistic heavy ion collisions using Trajectum of a significantly higher precision and including a new option to vary the normalization of the centrality estimator. We use the posterior distribution of our parameters to generate a set of high-statistics samples that allows us to make precise predictions including statistical and systematic uncertainties estimated from the model parameter distribution. The results are systematically compared with experiment whereby we also include many observables not included in the global analysis. This includes in particular (extremely) ultracentral anisotropic flow and mean transverse momentum, whereby we find satisfactory agreement with experiment where data are available. Lastly, we compute spectra and anisotropic flow for oxygen-oxygen collisions performed at the BNL Relativistic Heavy Ion Collider (RHIC) and to be performed at the the CERN Large Hadron Collider (LHC) and comment on how these collisions may inform us on properties of the quark-gluon plasma.