In this article, we present a simple bond-based energy correction method for reaction energies. The correction contributions for each bond are estimated using a linear regression approach that learns the differences between a high (G4MP2) and low level (DFT, PBE) method. The regression coefficients are derived from the reactions derived within the QM9 dataset. Even though we applied this method to reactions where the bond change between reactants and products involves only C, O, and N, the method can be easily extended to other types of bonds given the availability of an appropriate dataset. The method relies on the assumption that only a small part of a molecule changes in an elementary step of a chemical reaction. Therefore, the major contribution to the error in DFT calculated reaction energies comes from the changing bonds. It is worth mentioning that a bond-based reaction energy correction method, specific to the C --O bond, was previously carried out for the CO 2 reduction reaction. Our approach generalizes it to encompass a broad variety of molecular reactions. As the model is based on linear regression, the resulting coefficients can simply be added to the corresponding DFT calculated reaction energies. We demonstrate the applicability of the method with three DFT functionals of the different rung, as well as both single and multi-step reactions with molecules much beyond the size of the molecules included in the QM9 dataset.