started work on sffinfo command. fixed bug across all paralellized commands if the...

[mothur.git] / trimseqscommand.h

diff --git a/trimseqscommand.h b/trimseqscommand.h
index ee79d777a8d11ffaa2e279e69463314f9c3b8938..ab692d3ee0dbd3ce77ce158ab9a0e82aca7719e6 100644 (file)
--- a/trimseqscommand.h
+++ b/trimseqscommand.h
@@ -22,25 +22,45 @@ public:
        void help();
        
 private:
-       void getOligos(vector<ofstream*>&, vector<ofstream*>&);
+
+       struct linePair {
+               unsigned long int start;
+               int num;
+               linePair(unsigned long int i, int j) : start(i), num(j) {}
+       };
+
+       void getOligos(vector<string>&);
        bool stripQualThreshold(Sequence&, ifstream&);
        bool cullQualAverage(Sequence&, ifstream&);
-       bool stripBarcode(Sequence&, int&);
-       bool stripForward(Sequence&);
+       int stripBarcode(Sequence&, int&);
+       int stripForward(Sequence&, int&);
        bool stripReverse(Sequence&);
        bool cullLength(Sequence&);
        bool cullHomoP(Sequence&);
        bool cullAmbigs(Sequence&);
        bool compareDNASeq(string, string);
+       int countDiffs(string, string);//, int, int&, int);
 
        bool abort;
-       string fastaFile, oligoFile, qFileName;
+       string fastaFile, oligoFile, qFileName, outputDir;
        
-       bool flip, allFiles;
-       int numFPrimers, numRPrimers, maxAmbig, maxHomoP, minLength, maxLength, qThreshold, qAverage;
-       vector<string> forPrimer, revPrimer;
+       bool flip, allFiles, qtrim;
+       int numFPrimers, numRPrimers, maxAmbig, maxHomoP, minLength, maxLength, qThreshold, qAverage, processors, tdiffs, bdiffs, pdiffs, comboStarts;
+       vector<string> revPrimer, outputNames;
+       set<string> filesToRemove;
        map<string, int> barcodes;
        vector<string> groupVector;
+       map<string, int> primers;
+       map<string, int> combos;
+       
+       vector<int> processIDS;   //processid
+       vector<linePair*> lines;
+       vector<linePair*> qLines;
+       
+       int driverCreateTrim(string, string, string, string, string, vector<string>, linePair*, linePair*);     
+       int createProcessesCreateTrim(string, string, string, string, string, vector<string>);
+       int setLines(string, vector<linePair*>&);
+       
 };
 
 #endif