This is a famous result of Fatou and while the usual proofs use the properties of the Poisson kernel, there is a cool proof using just Fourier series.
Let $f(z)=\sum{a_nz^n}, |f(z)| \le 1, |z| <1$. Since $\sum{|a_n|^2r^{2n}}=\frac{1}{2\pi}\int_0^{2\pi}|f(re^{it}|^2dt \le 1$, it follows immediately $\sum{|a_n|^2} \le 1$, hence $f(t)=\sum{a_ne^{it}} \in L^2(dt)$ the Hilbert (!) space of square integrable functions on the unit circle with the normalized Lebesgue measure (or if you want the usual periodic real functions etc).
The hypothesis implies that $f(t)=0$ on an arc, hence on a set $E$ of nonzero measure on the unit circle - this is actually what we need and of course, the crucial fact that $f(t)$ has no negative index Fourier series terms as it comes from an analytic functions inside the unit disc. If we prove that $f=0$ a.e. on the unit circle, we get $a_n=0$ for all $n$ hence $f=0$ in the unit disc too
Assume $f$ not identically zero (on the circle) and assume wlog $a_0 \ne 0$ (as otherwise, we take $e^{-ikt}f(t)$ where $a_k \ne 0$ is the first nonzero coefficient). Consider the convex set $C_0=f(e^{it})(1+b_1e^{it}+...b_me^{imt})$ where $m \ge 1, b_k$ arbitrary and take its closure $C$ in $L^2$. This has a unique element $g$ of minimal norm by the usual basic facts about Hilbert spaces. We claim that $|g|$ is constant (a.e.) and since obviously $g$ vanishes where $f$ does by construction, so in particular on a set of positive measure, we get $|g|=0$, hence $g=0$ (a.e) while $g$ has constant Fourier term $a_0 \ne 0$ which is a contradiction.
So let's prove the claim about $g$. By construction $g+\alpha e^{int}g \in C, n \ge 1$ (here it is crucial that $f$ hence $C$ consists of functions with Fourier series that start with $a_0$ but have only nonnegative index terms as otherwise obviously the statement above is not true as we can introduce constant terms when we multiply a $e^{-int}$ term with $\alpha e^{int}$).
By minimality $||g+\alpha e^{int}g||^2=||g||^2(1+|\alpha|^2)+2\Re \alpha \frac{1}{2\pi}\int_0^{2\pi}{|g|^2e^{int}}$ has a minimum at $\alpha =0$ and obviously this implies $\int_0^{2\pi}{|g|^2e^{int}}dt=0, n \ge 1$ (otherwise as usual the $||g||^2|\alpha|^2$ being quadratic in small $\alpha$, will be overwhelmed by the linear in $\alpha$ integral term and by choosing arguments appropriately we can make it negative of course). Conjugating we get the same equality for $n=-1,-2,...$, hence $|g|^2$ is constant as all its Fourier terms except the constant one vanish. Done!