makelookupcommand.h

   1 //
   2 //  makelookupcommand.h
   3 //  Mothur
   4 //
   5 //  Created by SarahsWork on 5/14/13.
   6 //  Copyright (c) 2013 Schloss Lab. All rights reserved.
   7 //
   8
   9 #ifndef Mothur_makelookupcommand_h
  10 #define Mothur_makelookupcommand_h
  11
  12 #include "command.hpp"
  13 #include "sequence.hpp"
  14
  15 /**************************************************************************************************/
  16
  17 class MakeLookupCommand : public Command {
  18 public:
  19     MakeLookupCommand(string);
  20     MakeLookupCommand();
  21     ~MakeLookupCommand(){}
  22
  23     vector<string> setParameters();
  24     string getCommandName()                     { return "make.lookup";                 }
  25     string getCommandCategory()         { return "Sequence Processing";         }
  26
  27     string getOutputPattern(string);
  28         string getHelpString();
  29     string getCitation() { return "Quince, C., A. Lanzén, T. P. Curtis, R. J. Davenport, N. Hall, I. M. Head, L. F. Read, and W. T. Sloan. 2009. Accurate determination of microbial diversity from 454 pyrosequencing data. Nat Methods 6:639-41. http://www.mothur.org/wiki/Make.lookup"; }
  30     string getDescription()             { return "Creates a lookup file for use with shhh.flows using user-supplied mock community data and flow grams"; }
  31
  32     int execute();
  33     void help() { m->mothurOut(getHelpString()); }
  34
  35 private:
  36     bool abort;
  37     string outputDir, flowFileName, errorFileName, flowOrder, refFastaFileName, barcodeSequence, keySequence;
  38     vector<string> outputNames;
  39     int thresholdCount;
  40
  41     vector<double> convertSeqToFlow(string sequence, string order);
  42     int alignFlowGrams(vector<double>& flowgram, vector<double>& refFlow, double gapOpening, vector<vector<double> > penaltyMatrix, string flowOrder);
  43     int regress(vector<double>& data, int N);
  44 };
  45
  46 /**************************************************************************************************/
  47
  48
  49
  50
  51 #endif