work on sfs
site frequency spectra