COUNTSEGREGATINGSITES - Counts segregating sites
This function counts segregating sites in alignment(aln).
Syntax: [s_site,v_site]=countsegregatingsites(aln)
Inputs:
aln - Alignment structure
Outputs:
s_site - Number of segregating sites
v_site - Number of valid sites = total length-gaps
m_num - Number of total mutation
sn_num - Number of singleton mutation site
sm_num - Number of singleton mutations
sm_numv - Vector of numbers of singleton mutations for each allele