We apply our four-state 70-vibration vibronic-coupling model for the properties of the photosynthetic special-pair radical cation to: (1) interpret the observed correlations between the midpoint potential and the distribution of spin density between the two bacteriochlorophylls for 30 mutants of Rhodobacter sphaeroides, (2) interpret the observed average intervalence hole-transfer absorption energies as a function of spin density for six mutants, and (3) simulate the recently obtained intervalence electroabsorption Stark spectrum of the wild-type reaction center. While three new parameters describing the location of the sites of mutation with respect to the special pair are required to describe the midpoint-potential data, a priori predictions are made for the transition energies and the Stark spectrum. In general, excellent predictions are made of the observed quantities, with deviations being typically of the order of twice the experimental uncertainties. A unified description of many chemical and spectroscopic properties of the bacterial reaction center is thus provided. Central to the analysis is the assumption that the perturbations made to the reaction center, either via mutations of protein residues or by application of an external electric field, act only to independently modify the oxidation potentials of the two halves of the special pair and hence the redox asymmetry Eo. While this appears to be a good approximation, clear evidence is presented that effects of mutation can be more extensive than what is allowed for. A thorough set of analytical equations describing the observed properties is obtained using the Born-Oppenheimer adiabatic approximation. These equations are generally appropriate for intervalence charge-transfer problems and include, for the first time, full treatment of both symmetric and antisymmetric vibrational motions. The limits of validity of the adiabatic approach to the full nonadiabatic problem are obtained.