X-Git-Url: https://git.donarmstrong.com/?a=blobdiff_plain;f=trimseqscommand.h;h=2d0eabc1aa5a1504d23d1fc146f525504f3213f9;hb=ca9ac1d80c62f57270b0dcd49410ebe08a8aecd6;hp=2f8fbd7c499559bf65f9f5393982bd00dddfbcc3;hpb=d9b668f68b99f92ecdc71dd8cd363cb4e27107f9;p=mothur.git

diff --git a/trimseqscommand.h b/trimseqscommand.h
index 2f8fbd7..2d0eabc 100644
--- a/trimseqscommand.h
+++ b/trimseqscommand.h
@@ -20,13 +20,15 @@ class TrimSeqsCommand : public Command {
 public:
 	TrimSeqsCommand(string);
 	TrimSeqsCommand();
-	~TrimSeqsCommand();
-	vector<string> getRequiredParameters();
-	vector<string> getValidParameters();
-	vector<string> getRequiredFiles();
-	map<string, vector<string> > getOutputFiles() { return outputTypes; }
-	int execute();
-	void help();
+	~TrimSeqsCommand(){}
+	
+	vector<string> setParameters();
+	string getCommandName()			{ return "trim.seqs";	}
+	string getCommandCategory()		{ return "Sequence Processing";		}
+	string getHelpString();	
+	
+	int execute(); 
+	void help() { m->mothurOut(getHelpString()); }	
 	
 private:
 	
@@ -37,24 +39,27 @@ private:
 		unsigned long int end;
 		linePair(unsigned long int i, unsigned long int j) : start(i), end(j) {}
 	};
-
-	void getOligos(vector<string>&, vector<string>&);
+	
+	void getOligos(vector<vector<string> >&, vector<vector<string> >&);
 	int stripBarcode(Sequence&, QualityScores&, int&);
 	int stripForward(Sequence&, QualityScores&, int&);
 	bool stripReverse(Sequence&, QualityScores&);
+	
+	bool keepFirstTrim(Sequence&, QualityScores&);
+	bool removeLastTrim(Sequence&, QualityScores&);
+
 	bool cullLength(Sequence&);
 	bool cullHomoP(Sequence&);
 	bool cullAmbigs(Sequence&);
 	bool compareDNASeq(string, string);
 	int countDiffs(string, string);
-	map<string, vector<string> > outputTypes;
 
 	bool abort;
 	string fastaFile, oligoFile, qFileName, groupfile, outputDir;
 	
 	bool flip, allFiles, qtrim;
 	int numFPrimers, numRPrimers, maxAmbig, maxHomoP, minLength, maxLength, processors, tdiffs, bdiffs, pdiffs, comboStarts;
-	int qWindowSize, qWindowStep;
+	int qWindowSize, qWindowStep, keepFirst, removeLast;
 	double qRollAverage, qThreshold, qWindowAverage, qAverage;
 	vector<string> revPrimer, outputNames;
 	set<string> filesToRemove;
@@ -63,15 +68,17 @@ private:
 	map<string, int> primers;
 	map<string, int> combos;
 	map<string, int> groupToIndex;
-	
+	vector<string> primerNameVector;	//needed here?
+	vector<string> barcodeNameVector;	//needed here?
+	map<string, int> groupCounts;  
+
 	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 driverCreateTrim(string, string, string, string, string, string, string, vector<vector<string> >, vector<vector<string> >, linePair*, linePair*);	
+	int createProcessesCreateTrim(string, string, string, string, string, string, string, vector<vector<string> >, vector<vector<string> >);
 	int setLines(string, string, vector<unsigned long int>&, vector<unsigned long int>&);
-	
 };
 
 #endif