Protein complexes often represent an ensemble of different assemblies with distinct functions and regulation. This increased complexity is enabled by the variety of protein diversification mechanisms that exist at every step of the protein biosynthesis pathway, such as alternative splicing and post transcriptional and translational modifications. The resulting variation in subunits can generate compositionally distinct protein assemblies. These different forms of a single protein complex may comprise functional variances that enable response and adaptation to varying cellular conditions. Despite the biological importance of this layer of complexity, relatively little is known about the compositional heterogeneity of protein complexes, mostly due to technical barriers of studying such closely related species. Here, we show that native mass spectrometry (MS) offers a way to unravel this inherent heterogeneity of protein assemblies. Our approach relies on the advanced Orbitrap mass spectrometer capable of multistage MS analysis across all levels of protein organization. Specifically, we have implemented a two-step fragmentation process in the inject flatapole device, which was converted to a linear ion trap, and can now probe the intact protein complex assembly, through its constituent subunits, to the primary sequence of each protein. We demonstrate our approach on the yeast homotetrameric FBP1 complex, the rate-limiting enzyme in gluconeogenesis. We show that the complex responds differently to changes in growth conditions by tuning phosphorylation dynamics. Our methodology deciphers, on a single instrument and in a single measurement, the stoichiometry, kinetics, and exact position of modifications, contributing to the exposure of the multilevel diversity of protein complexes.