We report an analytical bond energy from bond orders and populations (BEBOP) model that provides intramolecular bond energy decompositions for chemical insight into the thermochemistry of molecules. The implementation reported here employs a minimum basis set Mulliken population analysis on well-conditioned Hartree-Fock orbitals to decompose total electronic energies into physically interpretable contributions. The model's parametrization scheme is based on atom-specific parameters for hybridization and atom pair-specific parameters for short-range repulsion and extended Hückel-type bond energy term fitted to reproduce CBS-QB3 thermochemistry data. The current implementation is suitable for molecules involving H, Li, Be, B, C, N, O, and F atoms, and it can be used to analyze intramolecular bond energies of molecular structures at optimized stationary points found from other computational methods. This first-generation model brings the computational cost of a Hartree-Fock calculation using a large triple-ζ basis set, and its atomization energies are comparable to those from widely used hybrid Kohn-Sham density functional theory (DFT, as benchmarked to 109 species from the G2/97 test set and an additional 83 reference species). This model should be useful for the community by interpreting overall ab initio molecular energies in terms of physically insightful bond energy contributions, e.g., bond dissociation energies, resonance energies, molecular strain energies, and qualitative energetic contributions to the activation barrier in chemical reaction mechanisms. This work reports a critical benchmarking of this method as well as discussions of its strengths and weaknesses compared to hybrid DFT (i.e., B3LYP, M062X, PBE0, and APF methods), and other cost-effective approximate Hamiltonian semiempirical quantum methods (i.e., AM1, PM6, PM7, and DFTB3).