We consider eigenvector statistics of large symmetric random matrices. When the matrix entries are sampled from independent Gaussian random variables, eigenvectors are uniformly distributed on the sphere and numerous properties can be computed exactly. In particular, we can bound their extremal coordinates with high probability. There has been an extensive amount of work on generalizing such a result, known as delocalization, to more general entry distributions. After giving a brief overview of the previous results going in this direction, we present an optimal delocalization result for matrices with sub-exponential entries for all eigenvectors. The proof is based on the dynamical method introduced by Erdos-Yau, an analysis of high moments of eigenvectors as well as new level repulsion estimates which will be presented during the talk. This is based on a joint work with P. Lopatto.