Let $M$ be a $n$-dimensional manifold and let
$$\star : \Lambda^p(M) \longrightarrow \Lambda^{n-p} : \beta \longmapsto \star \beta$$
be the Hodge star operator, such that
$$ ( \alpha, \beta ) = \int_M \alpha \wedge \star \beta$$
for all $\alpha \in \Lambda^p$. Here, $\big( \boldsymbol{\cdot}, \boldsymbol{\cdot} \big)$ denotes the inner product, defined as
$$ (\alpha, \beta) = \frac{1}{p!} \int_M \alpha_{i_1 \ldots i_p} \; \beta_{j_1 \ldots j_p} \; g^{i_1 j_1} \cdots g^{i_p j_p} \; dV$$
Show that
$$\star \big( dx^{i_1} \wedge \ldots \wedge dx^{i_p} \big) = \frac{\sqrt{\det g}}{(n-p)!} \cdot g^{i_1 j_1} \cdots g^{i_p j_p} \cdot \varepsilon_{j_1 \ldots j_n} \cdot dx^{j_{p+1}} \wedge \ldots \wedge dx^{j_n}$$
Now, I suppose this is a pretty standard problem, when dealing with differential geometry -- however, I neither managed to solve it myself nor did I find a complete proof of that proposition. Most of the proofs I found were either missing some vital step, or did not bother to expand the proof as detailed as I would need it to understand.
Can anyone help me by providing a source for a proof or can contribute a proof itself? I would really appreciate it.
P.S.: I found the answer to How to derive the coordinate expression of the Hodge dual? useful, but I don't understand some of the steps. I copied the proof and marked the steps I don't understand with a question mark (?)
Since both sides are linear in $\omega$, it suffices to check on a basis. Let $\omega := \mathrm{d}x^{1}\wedge...\wedge \mathrm{d}x^{k}$. Then, one immediately notices that $\eta \overset{?}{=} \varepsilon^{i_1 ... i_k} \omega$. Observe that \begin{align*} ⟨\omega,\eta⟩ \ vol \enspace &= \enspace \varepsilon^{i_1 ... i_k} \ ⟨\omega,\omega⟩ \sqrt{g} \ \mathrm{d} x^1 \wedge ... \wedge \mathrm{d}x^n\\ &\overset{?}{=} \enspace \sqrt{g} \ \varepsilon^{i_1 ... i_k} \ g^{1 j_1} ... g^{k j_k} \ \varepsilon_{j_1 ... j_k} \mathrm{d} x^1 \wedge ... \wedge \mathrm{d}x^n\\ &= \enspace \sqrt{g} \ \varepsilon^{i_1 ... i_k} \ \mathrm{d} x^1 \wedge ... \wedge \mathrm{d}x^n\\ &= \enspace \sqrt{g} \ \Big( \frac{1}{(n-k)!} g^{i_1 j_1} ... g^{i_k j_k} \ \varepsilon_{j_1 ... j_n} \ \varepsilon^{j_{k+1} ... j_n} \Big) \ \mathrm{d} x^1 \wedge ... \wedge \mathrm{d}x^n\\ &\overset{?}{=} \enspace \mathrm{d} x^1 \wedge ... \wedge \mathrm{d} x^k \wedge \Big( \frac{\sqrt{g}}{(n-k)!} g^{i_1 j_1} ... g^{i_k j_k} \ \varepsilon_{j_1 ... j_n} \Big) \ \mathrm{d} x^{j_{k+1}} \wedge ... \wedge \mathrm{d}x^{j_n}\\ &= \enspace \omega \wedge \widetilde{\eta} \end{align*}