TL;DR- General relativity (GR) is based upon general coordinate invariance; this means that physics is invariant under a general coordinate transformation (GCT). This invariance implies the contracted Bianchi identities ($\nabla_\nu G^{\mu \nu} =0$), which in turn give us constraints on the equations of motion. We needed some constraints anyway since we allowed for arbitrary GCT in the beginning. That Einstein's field equations do not determine $g_{\mu \nu}$ uniquely, but only up to $4$ arbitrary coordinate transformations, is explained by the contracted Bianchi identities.
Ignoring factors of $16$ and $\pi$ throughout the answer, the Einstein-Hilbert action is
$$S = \int_Vd^4x \sqrt{-g} R,$$
where $V$ is the spacetime region of integration. Taking the variation of $S$ (with respect to $g^{\mu \nu}$) gives
$$\delta S = \int_V d^4x \sqrt{-g} G_{\mu \nu} \delta g^{\mu \nu} \tag{1}.$$
Under arbitrary variations of the metric, $\delta g^{\mu \nu}$, the least action principle $\delta S=0$ gives us Einstein's equations of motion in vacuum: $G_{\mu \nu}=0$. (You can repeat the procedure by adding a matter Lagrangian). So $G_{\mu \nu}$ is the source-free part of the equations of motion for the metric field $g_{\mu \nu}$. In $3+1$ dimensions, $G_{\mu \nu}$ is the unique tensor (apart from the metric tensor itself) that is constructed from $g_{\mu \nu}$ and its first and second derivatives, is symmetric in its two indices and is divergence-free (Lovelock's theorem). In higher dimensions, $G_{\mu \nu}$ is no longer unique, if you allow non-linear functions of 2nd derivatives of the metric. But if you allow only linear functions of 2nd derivatives of the metric, $G_{\mu \nu}$ remains unique.
But without hurrying up to obtain the equations of motion, you can gain useful information just from the form of the action variation in $(1)$.
GR must be invariant under GCT, so $S$ must be invariant under GCT $\Rightarrow \delta S = 0$ under GCT: $x \rightarrow x'$. What could this tell us? It is sufficient to consider infinitesimal GCT. So, suppose
$$x^\mu \rightarrow x'^\mu = x^\mu + \epsilon^\mu,$$
where $\epsilon^\mu$ is arbitrary inside $V$ but constrained to vanish on the boundary of $V$: the hypersurface $\partial V$.
Under this infinitesimal GCT, evaluate the variation of the metric tensor field,
$$\delta g_{\mu \nu} = g'_{\mu \nu}(x) - g_{\mu \nu}(x),$$
up to first order in $\epsilon$, ignoring $O(\epsilon^2)$ and higher terms. Raising the indices gives you
$$\delta g^{\mu \nu} = \nabla^\mu \epsilon^\nu + \nabla^\nu \epsilon^\mu. \tag{2}$$
Note the symmetry of indices $\mu, \nu$. Substitute $(2)$ in $(1)$ to get
$$\delta S = \int_V d^4x \sqrt{-g} \ G^{\mu \nu} \nabla_\nu \epsilon_\mu.$$
Integrating by parts and using Gauss' theorem gives
$$\delta S = -\int_V d^4x \sqrt{-g} \ (\nabla_\nu G^{\mu \nu}) \epsilon_\mu + \oint_{\partial V} d \Sigma_\nu G^{\mu \nu} \epsilon_\mu.$$
Recalling that $\epsilon_\mu$ is arbitrary in $V$ and that $\epsilon_\mu = 0$ on $\partial V$, we see that asking for action invariance, $\delta S=0$, under GCT implies
$$\nabla_\nu G^{\mu \nu} = 0, \tag{3}$$
which are the contracted Bianchi identities. So we see that (3) can be seen as a consequence of imposing GCT on GR.
Einstein's equations $G_{\mu \nu} = T_{\mu \nu}$ seem like they imply that there are $10$ equations for $10$ unknowns in $g_{\mu \nu}$. But that's not the full story. Since we are allowed to do GCT, Einstein's equations do not uniquely determine $g_{\mu \nu}$, but only up to $4$ arbitrary coordinate transformations. The $4$ equations in $\nabla_\nu G^{\mu \nu} =0$ provide the missing link. Explicitly rewriting the Bianchi identities in terms of the timelike components of $G^{\mu \nu}$, we see
$$\partial_t G^{t \nu} = \text{some expression of $G$ and $\Gamma$ that contains at most 2nd derivatives of $g_{\alpha \beta}$},$$
so that $G^{t \nu}$ on the LHS must contain at most 1st derivatives of $g_{\alpha \beta}$. But specifying the metric and its first derivative accounts for specifying the initial conditions, and it does not really account for dynamical equations of motion. We therefore are able to get $4$ constraint equations, and $10-4=6$ truly dynamical equations of motion.