This work focuses on the low-temperature (5 K) photochemical (transient) hole-burned (HB) spectra within the P870 absorption band, and their theoretical analysis, for the (M)L214G mutant of the photosynthetic Rhodobacter sphaeroides bacterial reaction center (bRC). To provide insight into system-bath interactions of the bacteriochlorophyll a (BChl a) special pair, i.e., P870, in the mutated bRC, the optical line shape function for the P870 band is calculated numerically. On the basis of the modeling studies, we demonstrate that (M)L214G mutation leads to a heterogeneous population of bRCs with modified (increased) total electron-phonon coupling strength of the special pair BChl a and larger inhomogeneous broadening. Specifically, we show that after mutation in the (M)L214G bRC a large fraction (∼50%) of the bacteriopheophytin (HA) chromophores shifts red and the 800 nm absorption band broadens, while the remaining fraction of HA cofactors retains nearly the same site energy as HA in the wild-type bRC. Modeling using these two subpopulations allowed for fits of the absorption and nonresonant (transient) HB spectra of the mutant bRC in the charge neutral, oxidized, and charge-separated states using the Frenkel exciton Hamiltonian, providing new insight into the mutant's complex electronic structure. Although the average (M)L214G mutant quantum efficiency of P(+)QA(-) state formation seems to be altered in comparison with the wild-type bRC, the average electron transfer time (measured via resonant transient HB spectra within the P870 band) was not affected. Thus, mutation in the vicinity of the electron acceptor (HA) does not tune the charge separation dynamics. Finally, quenching of the (M)L214G mutant excited states by P(+) is addressed by persistent HB spectra burned within the B band in chemically oxidized samples.