consensusseqscommand.h

   1 #ifndef CONSENSUSSEQSCOMMAND_H
   2 #define CONSENSUSSEQSCOMMAND_H
   3 //test
   4 /*
   5  *  consensusseqscommand.h
   6  *  Mothur
   7  *
   8  *  Created by westcott on 11/23/10.
   9  *  Copyright 2010 Schloss Lab. All rights reserved.
  10  *
  11  */
  12
  13
  14 #include "command.hpp"
  15 #include "listvector.hpp"
  16
  17 class ConsensusSeqsCommand : public Command {
  18 public:
  19         ConsensusSeqsCommand(string);
  20         ConsensusSeqsCommand();
  21         ~ConsensusSeqsCommand(){}
  22
  23         vector<string> setParameters();
  24         string getCommandName()                 { return "consensus.seqs";              }
  25         string getCommandCategory()             { return "Sequence Processing"; }
  26         string getOutputFileNameTag(string, string);
  27         string getHelpString();
  28         string getCitation() { return "http://www.mothur.org/wiki/Consensus.seqs"; }
  29         string getDescription()         { return "create a consensus sequence for each OTU or for a fasta file"; }
  30
  31
  32         int execute();
  33         void help() { m->mothurOut(getHelpString()); }
  34
  35 private:
  36
  37         bool abort, allLines;
  38         string fastafile, listfile, namefile, label, outputDir;
  39         set<string> labels;
  40         vector<string> outputNames;
  41         map<string, string> fastaMap;
  42         map<string, string> nameMap;
  43         map<string, string> nameFileMap;
  44         int cutoff;
  45
  46         int readFasta();
  47         int readNames();
  48         int processList(ListVector*&);
  49         string getConsSeq(string, ofstream&, string&, int);
  50         char getBase(vector<int>, int);
  51 };
  52
  53 #endif
  54
  55
  56
  57