Fast and accurate estimation of multidimensional site frequency spectra from low-coverage high-throughput sequencing data.

The new implementation allows for accurate estimation of population genetic parameters from arbitrarily large, low-coverage datasets, thus facilitating cost-effective sequencing experiments in model and non-model organisms.