I asked this question following a discussion in my Mathematical Methods course and didn't get a satisfactory answer.
If we have an infinite dimensional Hilbert space, why do we need an orthonormal basis? For finite dimensional spaces, it's fine to just have a linearly independent set that spans the space. Why the extra conditions when we move up to infinite dimensions?
When we're working with Hilbert spaces, we have the extra structure given by the norm: continuous linear maps are a very important subclass of such. And on a Hilbert space our notion of isomorphism is a linear surjection $T: H_1 \to H_2$ that preserves the inner product.
The main value of a basis for a vector space is that we can determine the entire linear map from where the basis maps alone. If we restricting to continuous linear maps, a Schauder basis is all we need to determine the entire linear map; because if $v = \sum_{n=0}^\infty a_n e_n$ (where $e_n$ is a Schauder basis), then - by continuity and linearity - $T(v) = \sum_{n=0}^\infty a_n T(e_n)$.
So a Schauder basis is enough for any Banach space (since the entire above discussion works there, too!). Unfortunately, Banach spaces don't necessarily have a Schauder basis - we need the Hilbert space structure to prove they exist.
We've also got even more structure, and it would be nice if our basis took that into account. So suppose we want to check if a map $T: H_1 \to H_2$ is orthogonal; that is, if $\langle v_1, v_2\rangle = \langle T(v_1), T(v_2)\rangle$. Checking this is easier again if we have an orthonormal basis: you only need to check if $\langle T(e_i), T(e_j)\rangle = 0$ for all $i \neq j$, and $\langle T(e_i), T(e_i)\rangle = 1$ for all $i$. So - since we're often interested in orthogonal operators, say - it's really nice to have an orthonormal basis, since it reduces our workload.
Here's a nice corollary of the existence of orthonormal bases of a Hilbert space.
Proof: I'll assume the bases are countable for clarity; the same argument works with any cardinality basis. Let $H_1$ have orthonormal basis $e_1, \dots$ and let $H_2$ have orthonormal basis $e'_1, \dots$. Define a map $H_1 \to H_2$ by $$T\left(\sum_{n=0}^\infty a_n e_n\right) = \sum_{n=0}^\infty a_n e'_n.$$ Because these are orthonormal bases, this is a continuous invertible linear map with continuous inverse; and because they're, again, orthonormal bases, $\langle v, w\rangle = \langle T(v), T(w)\rangle$. So this is an isomorphism of Hilbert spaces, as desired.
Orthonormal bases have other nice properties (some of which are mentioned in the comments), so we often like to pick one when we have Hilbert spaces in the first place. They make things easier to work with. But, of course, they're not the only Schauder bases. But if we're picking a basis - we might as well pick a nice one.