This is to be able to run and evaluate our software (old and new) on empirical data. The software is for estimating parameters under a population genetic split model and simultaneously performing goodness of fit of these models. The input data consists of whole genome data from different individuals, either as vcfs or as bams.