  Level statistics is a crucial tool in the exploration of localization
physics. The level spacing distribution of the disordered localized phase
follows Poisson statistics, and many studies naturally apply it to the
quasiperiodic localized phase. Here we analytically obtain the level spacing
distribution of the quasiperiodic localized phase, and find that it deviates
from Poisson statistics. Moreover, based on this level statistics, we derive
the ratio of adjacent gaps and find that for a single sample, it is a $\delta$
function, which is in excellent agreement with numerical studies. Additionally,
unlike disordered systems, in quasiperiodic systems, there are variations in
the level spacing distribution across different regions of the spectrum, and
increasing the size and increasing the sample are non-equivalent. Our findings
carry significant implications for the reevaluation of level statistics in
quasiperiodic systems and a profound understanding of the distinct effects of
quasiperiodic potentials and disorder induced localization.
