The Feynman propagator is given by the expectation value of two time-ordered (scalar) field operators (evaluated in the vacuum): $$ G_\mathrm{F}(x,y) \equiv \langle 0 | \mathcal{T}\big( \hat{\phi}(x) \hat{\phi}(y) \big) | 0 \rangle $$
This is a Green's function for the Klein-Gordon operator in the sense that $$ \bigg( - \left( \tfrac{\partial}{\partial x^0} \right)^2 + \boldsymbol{\nabla}_{\mathbf{x}}^2 - m^2\bigg) G_\mathrm{F}(x,y) = - \delta^{(4)}(x-y) $$ See for example, Equation (6.2.17) in Weinberg's Volume 1.
Normally, mathematicians define a generic Green's function $G$ for an operator $\hat{L}$ as one satisfying; $$ \hat{L} G(x;y) = + \delta(x-y) $$ Note the difference in minus sign.
My question is, why do physicists have an extra minus sign for the Green's function they always use? Is it just a convention that caught on, or is there a practical reason?