If we let $E$ and $F$ be Hilbert spaces and $A: E \rightarrow F$ a bounded linear operator between them, then $A$ is a Hilbert-Schmidt operator if there exists a Hilbert basis for $E$ (i.e. a complete, orthonormal sequence in $E$), $(e_n)_{n=1}^\infty$, such that $$\displaystyle\sum_{n=1}^\infty ||A(e_n)||_F^2<\infty.$$ What I'd like to show is that if $(f_n)_{n=1}^\infty$ is any other complete orthonormal sequence in $E$, then $$\displaystyle\sum_{n=1}^\infty ||A(e_n)||_F^2 = \displaystyle\sum_{n=1}^\infty ||A(f_n)||_F^2.$$ It seems like the easiest way of doing this is using the fact that we can represent each $e_n$ and $f_n$ in terms of each other, since for any $x\in E$, $x = \sum_{n=1}^\infty \langle x, e_n \rangle e_n = \sum_{n=1}^\infty \langle x, f_n \rangle f_n$, and we may allow that $x$ to be any of the $e_n$ or $f_n$.
But I'm getting a bit lost in the details. Is this a sensible way of proving the statement, and if so, am I just a tweak away from finishing the proof?