Significance Cardiolipins are a unique class of phospholipids in mitochondrial membranes that are crucial for cellular bioenergetics as they stabilize respiratory chain complexes. In contrast to most other phospholipids, cardiolipins are substituted with four, rather than only two fatty acyl side chains. Consequently, this opens up a vast number of different theoretically possible molecular lipid species. Experimentally assessing the molecular diversity of cardiolipin species is analytically challenging. In this study we successfully combine tandem mass spectrometry with a mathematical structural modeling approach, to achieve the comprehensive characterization of complex biological cardiolipin compositions.