I am trying to better understand Noether identities, i.e. relations between equations of motion in the presence of gauge symmetries for the example of the relativistic point particle.
Formally, a Noether identity arises because the variation of an action leading to an equation of motion requires identical manipulations as the variation under a gauge symmetry. Thus, if an action is invariant under a gauge symmetry, we have the relation
$$ \delta_g S[\phi] = \int d^nx\,\frac{\delta L}{\delta \phi}\, \delta_g\phi = 0,\tag{1} $$ where $\phi$ is a generic field, $\delta_g$ the variation under a gauge symmetry and $\delta L/\delta \phi$ is the equation of motion for the field $\phi$. Substituting the variation and isolating the (continuous) variation parameter then results in the constraint.
However, I cannot make sense of it applying this logic to the action of the relativistic point particle
$$
S[x] = m \int d\tau \sqrt{-\dot{x}_\mu \dot{x}^\mu}.\tag{2}
$$
The gauge symmetry here is parametrization invariance $\tau \to \tau'$ or infinitesimally
$$
\delta \tau = \varepsilon(\tau),\tag{3}
$$
for $\varepsilon$ an arbitrary (continous and monotonous) function. Calculating the variation on the coordinates gives
$$
\delta x^\mu = - \dot{x}^\mu \varepsilon.\tag{4}
$$
Correspondingly, the Noether identity reads
$$
0 = \delta S[x] = m \int d\tau\,\frac{d}{d\tau} \left( \frac{\dot{x}_\mu}{\sqrt{- \dot{x}^2}} \right) \delta x^\mu \\
= -m \int d\tau\,\frac{d}{d\tau} \left( \frac{\dot{x}_\mu}{\sqrt{- \dot{x}^2}} \right) \dot{x}^\mu \varepsilon,\tag{5}
$$
i.e.
$$
\frac{d}{d\tau} \left( \frac{\dot{x}_\mu}{\sqrt{- \dot{x}^2}} \right) \dot{x}^\mu = 0.\tag{6}
$$
However, I find the result hard to interpret as the relation seems to be completely arbitrary.
Can anyone maybe expand on the interpretation of this result?