A hybrid Lattice Boltzmann(LB)-finite difference(FD) numerical scheme for the simulation of reacting flows at low Mach numbers is presented. The FD solver is used to model the energy and species fields while the LB model computes the flow field. The proposed LB solver is a modified version of the classical iso-thermal weakly compressible LB scheme with the hydrodynamic pressure as its zeroth-order moment instead of density, recovering the well-known low Mach number approximation for thermo-compressible flows. The proposed numerical approach is used to model a variety of test-cases ranging from 1-D to 3-D configurations, involving premixed and non-premixed flames. In all cases results obtained by this solver are shown to agree very well with reference data. (C) 2019 The Combustion Institute. Published by Elsevier Inc. All rights reserved.