pop_tajimas_d function

Estimate Tajima's D for the whole genome