I want to derive the Dirac massless equation in curved spacetime from the action. I have the symmetric form of the Dirac action: $$S = \frac{1}{2} \int \bigg[i\bar{\psi} \gamma^\mu D_\mu \psi - i D_\mu \bar{\psi} \gamma^\mu \psi \bigg] \sqrt{-g} d^4 x$$ where $D_\mu \bar{\psi}= \partial_\mu + \frac{i}{4} \omega_{\mu a b} \sigma^{ab} \bar{\psi}$ ($\omega$ is the spin connection and $\sigma^{ab} = \frac{i}{2} [\gamma^a, \gamma^b]$) and I proceed to do $$\frac{\delta S}{\delta \bar{\psi}} = \frac{1}{2} \int \bigg[i \gamma^\mu D_\mu \psi - i \frac{\delta }{\delta \bar{\psi}} (D_\mu \bar{\psi}) \gamma^\mu \psi \bigg] \sqrt{-g} d^4 x$$ If I'm not mistaken, $\delta$ and $D_\mu$ commute since I'm varying the matter fields over the background of a fixed metric and so $$\frac{\delta S}{\delta \bar{\psi}} = \frac{1}{2} \int \bigg[i \gamma^\mu D_\mu \psi - i D_\mu \frac{\delta \bar{\psi}}{\delta \bar{\psi}} \gamma^\mu \psi \bigg] \sqrt{-g} \, d^4 x \overbrace{=}^{\color{red}{!}} \frac{1}{2} \int \bigg[i \gamma^\mu D_\mu \psi - i D_\mu \gamma^\mu \psi \bigg] \sqrt{-g} d^4 x $$ If this is correct, I'm not sure how I can proceed. The Dirac equation is supposed to be $$i \gamma^\mu D_\mu \psi = 0$$ and thus the only way for that to happen is IF $$- i D_\mu \gamma^\mu \psi = i \gamma^\mu D_\mu \psi$$ but this doesn't seem to be correct. What am I misunderstanding here?
EDIT: $\frac{\delta \bar{\psi}}{\delta \bar{\psi}} \neq 1$, as mentioned in the comments. Therefore, we have $D_\mu (something) * \gamma^\mu \psi$. This means that the next term (with the red exclamation mark) in my computations is wrong. By doing partial integration, we get rid of a surface term and we change the derivative to $\gamma^\mu \psi$. Therefore, the correct relation is
$$\frac{\delta S}{\delta \bar{\psi}} = \frac{1}{2} \int \bigg[i \gamma^\mu D_\mu \psi + i D_\mu (\gamma^\mu \psi) \bigg] \sqrt{-g} d^4 x $$
where indeed the sign has changed. Now, as mentioned in the answer, $D_\mu \gamma^\nu = 0$ and thus
$$D_\mu (\gamma^\mu \psi) = (D_\mu \gamma^\mu) \psi + \gamma^\mu D_\mu \psi = \gamma^\mu D_\mu \psi$$
effectively changing the order of $D_\mu$ and $\gamma^\mu$. This means that
$$\frac{\delta S}{\delta \bar{\psi}} = \int \bigg[i \gamma^\mu D_\mu \psi \bigg] \sqrt{-g} d^4 x $$
and thus we get the desired Dirac equation in curved spacetime.