hap.div function

Haplotype Diversity