I am interested in comparing alpha diversity indices among samples at a set OTU level. Right now, the best way to collate this data is to read.otu(list,group,label) and summary.single(label). This produces a separate .summary file for each unique sample in .groups. To collate the data, I am currently concatenating the files and tracking which line is associated with which file. It would be swell if there was a way to collate all of this within mothur via a command like summary.single(group,label) which would list alpha diversity by group for a given label much like the .shared or .relabund files.