The usage of Monte Carlo neutrino event generators (MCνEGs) is a norm within the high-energy ν scattering community. The relevance of quasielastic energy regimes to ν oscillation experiments implies that accurate calculations of ν-nucleus cross sections in this regime will be a key contributor to reducing the systematic uncertainties affecting the extraction of oscillation parameters. In spite of this, many MCνEGs utilize highly phenomenological, parameterized models of quasielastic scattering cross sections. Moreover, a culture of validation of MCνEGs against prolific electron scattering data has been historically lacking. In this work, we implement new electron-nucleus cross sections obtained from nuclear ab initio approaches in GENIE, the primary MCνEG utilized by the FNAL community. In particular, we utilize results from Quantum Monte Carlo computational methods which solve the many-body nuclear problem in the Short-Time Approximation (STA), allowing consistent retention of two-nucleon dynamical effects that are crucial to explain available nuclear electromagnetic (electroweak) data over a wide range of energy and momentum transfers. This new implementation in GENIE is fully tested against the world quasielastic electromagnetic data, finding agreement with available data below ∼ 2 GeV of beam energy with the aid of a scaling function formalism. The STA is currently limited to study A ≤ 12 nuclei, however, its exclusive multibody identity components are exportable to other many-body computational techniques such as Auxiliary Field Diffusion Monte Carlo which can reach A ≤ 40 systems while continuing to realize the semifinal states contained within the STA's multinucleon dynamics. Together, these developments promise to make future experiments such as the Deep Underground Neutrino Experiment ever more accurate in their assessment of MCνEG systematics, ν properties, and potentially empower the discovery of physics beyond the Standard Model. 2. We discuss our new, holistic framework within GENIE based on calculated electromagnetic (electroweak) nuclear response functions and supplemented by interpolation schemes to derive double differential electron scattering cross sections;