I have a dataset consisting in COI sequences from several catshark populations with different sample size (7 < N < 24). In order to obtain comparable estimates of genetic diversity, I would like to employ a rarefaction algorithm. Can anyone suggest a rarefaction method similar to those employed with microsatellites and if there is any software performing such analysis on mtDNA data?