* Copyright 2009 Schloss Lab UMASS Amherst. All rights reserved.
*
*/
+
+ /* The get.oturep command outputs a .fastarep file for each distance you specify, selecting one OTU representative for each bin. */
#include "command.hpp"
#include "globaldata.hpp"
#include "sparsematrix.hpp"
#include "listvector.hpp"
#include "inputdata.h"
-#include "readmatrix.hpp"
+#include "readotu.h"
#include "fastamap.h"
-
+#include "groupmap.h"
class GlobalData;
SparseMatrix* matrix;
ListVector* list;
ListVector* listOfNames;
- ReadMatrix* read;
+ ReadOTUFile* read;
InputData* input;
FastaMap* fasta;
- string filename, fastafile, namesfile;
+ GroupMap* groupMap;
+ string filename, fastafile, namesfile, groupfile;
ofstream out;
ifstream in, inNames;
+ bool groupError;
map<string, int> nameToIndex; //maps sequence name to index in sparsematrix
map<string, int>::iterator it3;
void readNamesFile();
- string FindRep(int); // returns name of "representative" sequence of given bin.
+ int process(ListVector*);
+ string FindRep(int, string&, ListVector*); // returns name of "representative" sequence of given bin. //and fill a string containing the groups in that bin if a groupfile is given
};