The bounds of the integral have no dependence on any of the variables, and hence we may move the differential operator into the integrand,
δδη(z)∫d4ySF(z−y)η(y)=∫d4ySF(z−y)δ(4)(z−y)
Evaluating the integral using the standard delta distribution identity, we obtain your result, namely SF(z−z). In this case, the final answer does not pick a minus sign, even though η is Grassmann-valued. See Peskin and Schroeder's text on QFT for a summary of Berezin/Grassmann integration.
This post imported from StackExchange Physics at 2014-04-13 14:30 (UCT), posted by SE-user JamalS