Two consecutive symmetry transformation generated via Poisson brackets

I spent some time double checking this, and actually, unlike my comment indicated, the resolution to your issue has nothing to do with integrations by parts, but rather with a minus sign you seem to have missed somewhere in your calculation. Show that there is indeed a missing minus sign, let me first reproduce the argument for identifying the charge algebra with (a possible central extension of) the Killing vector algebra.

To begin, let me make a point on notation. There is a well-defined notion of a differential on the phase space, and so the notation $delta_xi F$ where $xi$ is one of your Killing vectors (not actually important that it's a Killing vector, only that it generates a symmetry) and $F$ is any function on the phase space is really a shorthand for the contraction of the phase space vector field associated to $xi$ contracted on the phase space differential of $F$. This is, almost by definition, precisely equivalent to the action of the spacetime Lie derivative of $F$ with respect to $xi$. So with this in mind, we know, by a matter of definition, that these $delta_f$ and $delta_g$ operators satisfy the same Lie bracket algebra as the Killing vectors we started with. The task then is to show that the Noether charges also obey the same algebra.

Hence, I will assume that I have some collection of vector fields $xi_a$ over spacetime where $a$ indexes which vector field I'm talking about. These vector fields should satisfy some algebra $$ [xi_a,xi_b]=f^{c}_{ ab}xi_c. $$ I denote the Noether charge associated to each of these by $Q_a$. Letting $F$ be any function on the phase space, it then follows that $$ [delta_{xi_a},delta_{xi_b}]F=delta_{xi_a}delta_{xi_b}F-delta_{xi_b}delta_{xi_a}F={{F,Q_b},Q_a}-{{F,Q_a},Q_b}. $$ Using the Jacobi identity written in the form (and note that we haven't missed a minus sign because all these terms can be put on the same side of the equality with cyclic orders $a,b,F$) $$ {{F,Q_b},Q_a}-{{F,Q_a},Q_b}=-{F,{Q_a,Q_b}} $$ we immediately obtain $$ [delta_{xi_a},delta_{xi_b}]F=-{F,{Q_a,Q_b}} $$ so the action of the commutator of the Killing fields is minus the action of the Poisson bracket of the Noether charges. If we were to put the structure constants I defined earlier into here, we would find that the algebra of Noether charges is isomorphic to that of the Killing vectors, up to possible central extension, but with negated structure constants.

There's an important point to be made here by this minus sign. Vector fields are naturally defined to act from the left on functions. However if we look at the action of the Noether charges via the Poisson bracket, ${cdot,Q_a}$, this is actually defining an action from the right. Hence the difference in the minus sign.

So, now that we have taken the time to find this minus sign and point out that the Poisson bracket is actually defining a right action, how does this answer your question? Well let's do the calculation you're asking about: We know that $$ delta_fdelta_gphi=f^nupartial_nu(g^mupartial_muphi). $$ Let's do the calculation in terms of the Poisson brackets: $$ {{phi,Q_g},Q_f}={g^mupartial_muphi,Q_f}=g^mupartial_mu{phi,Q_f}=g^mupartial_mu(f^nupartial_nuphi), $$ so we see that the Poisson bracket calculation actually gives the correct result, but ordered differently. This is what I mean by the bracket defining the right action as opposed to a right action. The ordering we supposed should give the result we wanted was based on thinking about actions from the left, not right. It's easy to see from here that noting the bracket gives a right action will immediately tell us that the algebra should have negated structure constants.