Scripts for 16S amplicon analysis