Spectral statistics of quantum systems have been studied in detail using the nearest neighbour level spacings, which for generic chaotic systems follows random matrix theory predictions. In this work, the probability density of the and spacings from a given level are introduced. Analytical predictions are derived using a matrix model. The density is generalized to the spacing density, which allows for investigating long-range correlations. For larger the probability density of spacings is well described by a Gaussian. Using these spacings we propose the ratio of the to the as an alternative to the ratio of successive spacings. For a Poissonian spectrum the density of the ratio is flat, whereas for the three Gaussian ensembles repulsion at small values is found. The ordered spacing statistics and their ratio are numerically studied for the integrable circle billiard, the chaotic cardioid billiard, the standard map and the zeroes of the Riemann zeta function. Very good agreement with the predictions is found.