A wide class of one-dimensional continuous-time discrete-space interacting particle systems are shown to be Pfaffian point processes at fixed times with kernels characterised by the solutions to associated two-dimensional ODEs. The models comprise instantaneously coalescing or annihilating random walks with fully spatially inhomogeneous jump rates and deterministic initial conditions, including additional pairwise immigration or branching in the pure interaction regimes. We formulate convergence of Pfaffian point processes via their kernels, enabling investigation of diffusive scaling limits, which boils down uniform convergence of lattice approximations to two-dimensional PDEs. Convergence to continuum point processes is developed for a subset of the discrete models. Finally, in the case of annihilating random walks with pairwise immigration we extend the picture to multiple times, establishing the extended Pfaffian property for the temporal process.