X-Git-Url: https://git.donarmstrong.com/?a=blobdiff_plain;f=trimseqscommand.h;h=93e8c776918ec9fa2d6f661d7c178380a4f47d66;hb=6e81846c8e5b2614f6b06643a9f558fb0e6669fa;hp=03ee1f5745222a1e3fb154c22b4099d26e77d029;hpb=0470f6d037aacb3563c3f7010708120a4a67d4e6;p=mothur.git

diff --git a/trimseqscommand.h b/trimseqscommand.h
index 03ee1f5..93e8c77 100644
--- a/trimseqscommand.h
+++ b/trimseqscommand.h
@@ -12,8 +12,8 @@
 
 #include "mothur.h"
 #include "command.hpp"
-#include "globaldata.hpp"
 #include "sequence.hpp"
+#include "qualityscores.h"
 
 class TrimSeqsCommand : public Command {
 public:
@@ -23,24 +23,45 @@ public:
 	void help();
 	
 private:
-	void getOligos();
-	bool stripBarcode(Sequence&, string&);
-	bool stripForward(Sequence&);
-	bool stripReverse(Sequence&);
+
+	struct linePair {
+		unsigned long int start;
+		unsigned long int end;
+		linePair(unsigned long int i, unsigned long int j) : start(i), end(j) {}
+	};
+
+	void getOligos(vector<string>&, vector<string>&);
+	int stripBarcode(Sequence&, QualityScores&, int&);
+	int stripForward(Sequence&, QualityScores&, int&);
+	bool stripReverse(Sequence&, QualityScores&);
 	bool cullLength(Sequence&);
 	bool cullHomoP(Sequence&);
 	bool cullAmbigs(Sequence&);
-	
-	GlobalData* globaldata;
-	OptionParser* parser;
-	map<string, string> parameters;
-	map<string, string>::iterator it;
+	bool compareDNASeq(string, string);
+	int countDiffs(string, string);
+
 	bool abort;
-	string fastafile;
-	bool oligos, flip;
-	int numFPrimers, numRPrimers, maxAmbig, maxHomoP, minLength, maxLength;
-	vector<string> forPrimer, revPrimer;
-	map<string, string> barcodes;
+	string fastaFile, oligoFile, qFileName, outputDir;
+	
+	bool flip, allFiles, qtrim;
+	int numFPrimers, numRPrimers, maxAmbig, maxHomoP, minLength, maxLength, processors, tdiffs, bdiffs, pdiffs, comboStarts;
+	int qWindowSize, qWindowStep;
+	double qRollAverage, qThreshold, qWindowAverage, qAverage;
+	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, string, string, vector<string>, vector<string>, linePair*, linePair*);	
+	int createProcessesCreateTrim(string, string, string, string, string, string, string, vector<string>, vector<string>);
+	int setLines(string, string, vector<unsigned long int>&, vector<unsigned long int>&);
+	
 };
 
 #endif