Adeno-associated viruses (AAVs) are small, non-enveloped, and have a T=1 icosahedral capsid. They belong to the Parvoviridae, genus Dependoparvovirus. Interest in AAVs has grown over recent years as they have emerged as promising gene therapy vectors. The AAV capsid, encapsulating the transgene, consists, in total, of 60 subunits made up from three distinct viral proteins (VPs) originating from the same cap gene (VP1, VP2, and VP3), which vary only in their N-terminus. While all three VPs play a crucial and specific role in cell-entry and transduction, their exact stoichiometry and organization in AAV capsids has, despite the availability of several high-resolution structures remained elusive. Here we obtained a set of native mass spectra of intact AAV capsids (Mw ≈ 3.8 MDa) that display both highly resolved regions and regions wherein interferences occur. Through spectrum simulation we resolved and elucidated this spectral complexity, allowing us to directly assess the VP stoichiometries in a panel of serotypes from different production platforms. The data reveals an extremely heterogeneous population of capsids of variable composition. The relative abundance for each of the hundreds of co-occurring capsid compositions is accurately described by a model based upon stochastic assembly from a mixed pool of expressed VP1, VP2, and VP3. We show that even the single-most abundant VP stoichiometry represents only a few percent of the total AAV population. We estimate that virtually every AAV capsid in a particular preparation has a unique composition and arrangement, i.e. no particle is identical. The systematic scoring of the stochastic assembly model against experimental high-resolution native MS data offers a sensitive and accurate new method to characterize these exceptionally heterogeneous gene-delivery vectors.