minor edits in otuSize2(). ots.ids <- sapply(split(mcols(sites.gr)$newotuID,mcols(sites.gr)$readID),unique)