+
+ //add headers
+ ofstream r;
+ m->openOutputFile(rare, r);
+ r << rareHeader << endl;
+ r.close();
+ m->appendFiles(rare+".temp", rare);
+ m->mothurRemove(rare+".temp");
+
+ ofstream a;
+ m->openOutputFile(abund, a);
+ a << abundHeader << endl;
+ a.close();
+ m->appendFiles(abund+".temp", abund);
+ m->mothurRemove(abund+".temp");