  We study various statistics related to the eigenvalues and eigenfunctions of
random Hamiltonians in the localized regime. Consider a random Hamiltonian at
an energy $E$ in the localized phase. Assume the density of states function is
not too flat near $E$. Restrict it to some large cube $\Lambda$. Consider now
$I_\Lambda$, a small energy interval centered at $E$ that asymptotically
contains infintely many eigenvalues when the volume of the cube $\Lambda$ grows
to infinity. We prove that, with probability one in the large volume limit, the
eigenvalues of the random Hamiltonian restricted to the cube inside the
interval are given by independent identically distributed random variables, up
to an error of size an arbitrary power of the volume of the cube. As a
consequence, we derive * uniform Poisson behavior of the locally unfolded
eigenvalues, * a.s. Poisson behavior of the joint distibutions of the unfolded
energies and unfolded localization centers in a large range of scales. * the
distribution of the unfolded level spacings, locally and globally, * the
distribution of the unfolded localization centers, locally and globally.
