Varying $\partial_\lambda\phi\,\partial^\lambda\phi$ wrt the metric tensor $g_{\mu\nu}$ in two different ways gives me different results. Obviously I'm doing something wrong. Where am I going wrong?

Method 1: \begin{align} (\delta g_{\mu\nu})\,\partial^\mu\phi\,\partial^\nu\phi \end{align}

Method 2: \begin{align} &(\delta g^{\mu\nu})\,\partial_\mu\phi\,\partial_\nu\phi \nonumber \\ =&(-g^{\mu\rho}g^{\nu\sigma}\delta g_{\rho\sigma})\,\partial_\mu\phi\,\partial_\nu\phi \quad(\because \delta g^{\mu\nu}=-g^{\mu\rho}g^{\nu\sigma}\delta g_{\rho\sigma} \,\,\text{as can be checked by varying the identity}\,\, g^{\mu\lambda}g_{\lambda\nu}=\delta^\mu_\nu) \nonumber\\ =&-(\delta g_{\rho\sigma})\,\partial^\rho\phi\,\partial^\sigma\phi \end{align} The second result differs from the first one by a minus sign. What's going wrong?