I have calculated genomewide Tajima's D using VCFtools --Tajima option and got a positive Tajima value (diploid organism, butterfly species). This could either mean balancing selection or population contraction. Is there a way to distinguish between the two? I have also done pairwise Fst and did find outlier loci in my populations. Can I use that as a proxy for selection along with positive Tajima's D?