X-Git-Url: https://git.donarmstrong.com/?a=blobdiff_plain;f=sffinfocommand.h;h=1909e2d43e6d83b3d48ef06922b2d369a0e104d5;hb=250e3b11b1c9c1e1ad458ab6c7e71ac2e67e11d9;hp=d8eb853dae6892711a9004ae9ad670032f47228e;hpb=3e2465c16d187247ce3befd29811c2d5dfc15ee8;p=mothur.git

diff --git a/sffinfocommand.h b/sffinfocommand.h
index d8eb853..1909e2d 100644
--- a/sffinfocommand.h
+++ b/sffinfocommand.h
@@ -11,13 +11,13 @@
  */
 
 #include "command.hpp"
-
+#include "groupmap.h"
 
 /**********************************************************/
 struct CommonHeader {
 	unsigned int magicNumber;
 	string version;
-	unsigned long int indexOffset;
+	unsigned long long indexOffset;
 	unsigned int indexLength;
 	unsigned int numReads;
 	unsigned short headerLength;
@@ -64,26 +64,42 @@ class SffInfoCommand : public Command {
 public:
 	SffInfoCommand(string);
 	SffInfoCommand();
-	~SffInfoCommand();
-	vector<string> getRequiredParameters();
-	vector<string> getValidParameters();
-	vector<string> getRequiredFiles();
-	map<string, vector<string> > getOutputFiles() { return outputTypes; }
-	int execute();
-	void help();
+	~SffInfoCommand(){}
+	
+	vector<string> setParameters();
+	string getCommandName()			{ return "sffinfo";					}
+	string getCommandCategory()		{ return "Sequence Processing";		}
+	
+	string getHelpString();	
+    string getOutputPattern(string);	
+	string getCitation() { return "http://www.mothur.org/wiki/Sffinfo"; }
+	string getDescription()		{ return "extract sequences reads from a .sff file"; }
+
+	int execute(); 
+	void help() { m->mothurOut(getHelpString()); }	
 	
 private:
-	string sffFilename, outputDir, accnosName;
-	vector<string> filenames, outputNames, accnosFileNames;
-	bool abort, fasta, qual, trim, flow, sfftxt, hasAccnos;
+	string sffFilename, sfftxtFilename, outputDir, accnosName, currentFileName, oligosfile, noMatchFile, groupfile;
+	vector<string> filenames, outputNames, accnosFileNames, oligosFileNames, groupFileNames;
+	bool abort, fasta, qual, trim, flow, sfftxt, hasAccnos, hasOligos, hasGroup;
+	int mycount, split, numFPrimers, numLinkers, numSpacers, pdiffs, bdiffs, ldiffs, sdiffs, tdiffs, numNoMatch;
 	set<string> seqNames;
-	map<string, vector<string> > outputTypes;
-	
-	int extractSffInfo(string, string);
+    map<string, int> barcodes;
+    map<string, int> primers;
+    GroupMap* groupMap;
+    vector<string> linker, spacer, primerNameVector, barcodeNameVector, revPrimer;
+    vector<vector<int> > numSplitReads;
+    vector<vector<string> > filehandles;
+    vector<vector<string> > filehandlesHeaders;
+    
+	//extract sff file functions
+	int extractSffInfo(string, string, string);
 	int readCommonHeader(ifstream&, CommonHeader&);
 	int readHeader(ifstream&, Header&);
-	int readSeqData(ifstream&, seqRead&, int, int);
+	bool readSeqData(ifstream&, seqRead&, int, Header&);
 	int decodeName(string&, string&, string&, string);
+    bool readOligos(string oligosFile);
+    bool readGroup(string oligosFile);
 	
 	int printCommonHeader(ofstream&, CommonHeader&); 
 	int printHeader(ofstream&, Header&);
@@ -92,7 +108,18 @@ private:
 	int printFastaSeqData(ofstream&, seqRead&, Header&);
 	int printQualSeqData(ofstream&, seqRead&, Header&);
 	int readAccnosFile(string);
-		
+	int parseSffTxt();
+	bool sanityCheck(Header&, seqRead&);
+    int adjustCommonHeader(CommonHeader);
+    int findGroup(Header header, seqRead read, int& barcode, int& primer);
+    int findGroup(Header header, seqRead read, int& barcode, int& primer, string);
+    string reverseOligo(string oligo);
+    
+	//parsesfftxt file functions
+	int parseHeaderLineToInt(ifstream&);
+	vector<unsigned short> parseHeaderLineToFloatVector(ifstream&, int);
+	vector<unsigned int> parseHeaderLineToIntVector(ifstream&, int);
+	string parseHeaderLineToString(ifstream&);
 };
 
 /**********************************************************/