Small suggestion: Add group column header to nmds axes file

It would be nice if mothur wrote the “group” column header in nmds() axes file output. This way, I could plot it with the same code as pcoa() axes file output. Easy enough to work around, obviously.