1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58
|
#ifndef PHYLOTYPECOMMAND_H
#define PHYLOTYPECOMMAND_H
/*
* phylotypecommand.h
* Mothur
*
* Created by westcott on 11/20/09.
* Copyright 2009 Schloss Lab. All rights reserved.
*
*/
#include "taxonomyequalizer.h"
#include "command.hpp"
/*************************************************************************/
class PhylotypeCommand : public Command {
public:
PhylotypeCommand(string);
~PhylotypeCommand(){}
vector<string> setParameters();
string getCommandName() { return "phylotype"; }
string getCommandCategory() { return "Clustering"; }
string getHelpString();
string getOutputPattern(string);
string getCitation() { return "http://www.mothur.org/wiki/Phylotype"; }
string getDescription() { return "cluster your sequences into OTUs based on their classifications"; }
int execute();
void help() { m->mothurOut(getHelpString()); }
private:
bool abort, allLines;
string taxonomyFileName, label, namefile, countfile;
set<string> labels; //holds labels to be used
int cutoff;
map<string, string> namemap;
vector<string> outputNames;
map<int, int> currentNodes;
map<int, int> parentNodes;
map<int, int>::iterator itCurrent;
};
/*************************************************************************/
#endif
|