#include "fastamap.h"
#include "groupmap.h"
-class GlobalData;
-
typedef list<PCell>::iterator MatData;
+typedef map<int, float> SeqMap;
class GetOTURepCommand : public Command {
-
+
public:
- GetOTURepCommand(string);
+ GetOTURepCommand(string);
~GetOTURepCommand();
int execute();
- void help();
-
+ void help();
+
private:
GlobalData* globaldata;
SparseMatrix* matrix;
ofstream out;
ifstream in, inNames;
bool groupError;
-
+
bool abort, allLines;
set<int> lines; //hold lines to be used
set<string> labels; //holds labels to be used
map<string, int> nameToIndex; //maps sequence name to index in sparsematrix
-
+
+ vector<SeqMap> seqVec; // contains maps with sequence index and distance
+ // for all distances related to a certain sequence
+
+
void readNamesFile();
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
+ string findRep(int, string&, ListVector*, int&); // returns the name of the "representative" sequence of given bin,
+ // fills a string containing the groups in that bin if a groupfile is given,
+ // and returns the number of sequences in the given bin
};