I'm trying to understand, as a self learner, the covariant formulation of Electromagnetism. In particular I've been stuck for a while on the Bianchi identity. As I've come to understand, when we introduce the field strength tensor $F_{\mu\nu}$ we can write the Maxwell's equation in an explicitly covariant form using observable quantities as the aforementioned tensor (because of its Gauge invariance). Maxwell's equation thus reduce to one single 4-tensor equation (Gauss's + Maxwell-Ampere's laws), plus an identity (Gauss's law for magnetism and Faraday-Maxwell's law): I don't understand how to derive the identity in terms of the field strength tensor. This is my reasoning up until now.
I will use latin indices that range from 1 to 3 and greek ones that range from 0 to 3. We can say that:
$$ F_{ij} = B_{ij} = \varepsilon_{ijk} B^k \Longrightarrow \varepsilon^{ijk} F_{ij}= \varepsilon_{ijk}\varepsilon^{ijk} B^k \Longrightarrow B^k = \frac{1}{3!} \varepsilon^{ijk} F_{ij}$$
and
$$ F_{0i} = -\frac{E_i}{c} \Longrightarrow E_i = -cF_{0i} $$
Let us consider a flat space-time, then Gauss's law for magnetism can be written as:
$$ \nabla \cdot \vec{B} = \partial_k B^k = 0 \Longrightarrow \partial_k \frac{1}{3!} \varepsilon^{ijk} F_{ij} = 0 \Longrightarrow \varepsilon^{ijk} \partial_k F_{ij} = 0 $$
and Faraday-Maxwell's law can be written as:
$$ \nabla \times \vec{E} + \frac{\partial \vec{B}}{\partial t} = \vec{0} \Longrightarrow \varepsilon^{kij} \partial_i E_j + c \partial_0 B^k = 0 \Longrightarrow \varepsilon^{kij} \partial_i \left( -c F_{0j} \right) + c \partial_0 \frac{1}{3!} \varepsilon^{ijk} F_{ij} = 0\Longrightarrow $$
$$ - \varepsilon^{kij} \partial_i F_{0j} + \frac{1}{3!} \varepsilon^{ijk} \partial_0 F_{ij} = 0$$
Now, if I'm not mistaken, $\varepsilon^{0 \beta \gamma \delta} = \varepsilon^{0ijk} = \varepsilon^{ijk}$, therefore I can sobsitute the three dimensional Levi-Civita symbol with its four dimensional variant:
$$ \varepsilon^{0ijk} \partial_k F_{ij} = 0 $$
$$ -\varepsilon^{0kij} \partial_i F_{0j} + \frac{1}{3!} \varepsilon^{0ijk} \partial_0 F_{ij} = 0$$
From here I don't know how to continue. The latter of the two equations really seems to involve some kind of Levi-Civita symbol, because of the swapping of the $0$ index between the 4-gradient component and the field tensor component, but the $\frac{1}{3!}$ is in the way. Moreover, even if I simplified this equation I wouldn't know how to unify the two to obtain the final identity in terms of the field strength tensor:
$$ \epsilon^{\alpha \beta \gamma \delta} \partial_{\beta} F_{\gamma \delta} = 0 $$
Thank you in advance for any advice!