Metagenomic amplicon analysis with stats