clusterfragmentscommand.h

   1 #ifndef CLUSTERFRAGMENTSCOMMAND_H
   2 #define CLUSTERFRAGMENTSCOMMAND_H
   3
   4 /*
   5  *  clusterfragmentscommand.h
   6  *  Mothur
   7  *
   8  *  Created by westcott on 9/23/10.
   9  *  Copyright 2010 Schloss Lab. All rights reserved.
  10  *
  11  */
  12
  13
  14 #include "command.hpp"
  15 #include "sequence.hpp"
  16
  17 /************************************************************/
  18 struct seqRNode {
  19         int numIdentical;
  20         int length;
  21         Sequence seq;
  22         string names;
  23         bool active;
  24         seqRNode() {}
  25         seqRNode(int n, Sequence s, string nm, int l) : numIdentical(n), seq(s), names(nm), active(1), length(l) {}
  26         ~seqRNode() {}
  27 };
  28 /************************************************************/
  29
  30 class ClusterFragmentsCommand : public Command {
  31
  32 public:
  33         ClusterFragmentsCommand(string);
  34         ClusterFragmentsCommand();
  35         ~ClusterFragmentsCommand() {}
  36
  37         vector<string> setParameters();
  38         string getCommandName()                 { return "cluster.fragments";           }
  39         string getCommandCategory()             { return "Sequence Processing";         }
  40         string getHelpString();
  41
  42         int execute();
  43         void help() { m->mothurOut(getHelpString()); }
  44
  45 private:
  46         bool abort;
  47         string fastafile, namefile, outputDir;
  48         int diffs, percent;
  49         vector<seqRNode> alignSeqs;
  50         map<string, string> names; //represents the names file first column maps to second column
  51         map<string, int> sizes;  //this map a seq name to the number of identical seqs in the names file
  52         map<string, int>::iterator itSize;
  53         vector<string> outputNames;
  54
  55         int readFASTA();
  56         void readNameFile();
  57         void printData(string, string); //fasta filename, names file name
  58         bool isFragment(string, string);
  59
  60 };
  61
  62 /************************************************************/
  63
  64 #endif
  65