The emergence of efficient fragmentation methods such as electron capture dissociation (ECD) and electron transfer dissociation (ETD) provides the opportunity for detailed structural characterization of heavily covalently modified large peptides and small proteins such as intact histones. Even with effective gas phase ion isolation so that a single molecular precursor ion is selected, the MSMS spectrum of a heavily modified peptide may reveal the presence of a mixture of peptides with the same amino acid sequence and the same total number of posttranslational modification (PTM) moieties (same PTM composition) but with different PTM configurations or site-specific occupancy isoforms. Currently available data analysis methods depend on a deisotoping procedure, which becomes less effective when spectra (fragmentation patterns) contain many overlapping isotopic distributions. Peptide database search engines can only identify the most abundant PTM configuration (PTM arrangement on different residues) in such mixtures. To identify all the PTM configurations present in these mixtures and to estimate their relative abundances, we extended our fragment assignment by visual assistance program to search for ions representing all possible configurations, subjected to the total PTM composition constraint. This resulted in the identification of PTM configurations supported by unique fragment ions, and their relative abundances were estimated by use of a non-negative least squares procedure.