X-Git-Url: https://git.donarmstrong.com/?p=mothur.git;a=blobdiff_plain;f=sffinfocommand.h;h=5ec6e72239d0772548ce10f584ebccb435230d25;hp=903a589cd777b89b7e508e48d1173d550210f864;hb=050a3ff02473a3d4c0980964e1a9ebe52e55d6b8;hpb=ae57e166b2ed7b475ec3f466106bd76fabadd063 diff --git a/sffinfocommand.h b/sffinfocommand.h index 903a589..5ec6e72 100644 --- a/sffinfocommand.h +++ b/sffinfocommand.h @@ -69,7 +69,9 @@ public: vector 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"; } @@ -77,18 +79,24 @@ public: void help() { m->mothurOut(getHelpString()); } private: - string sffFilename, sfftxtFilename, outputDir, accnosName; - vector filenames, outputNames, accnosFileNames; - bool abort, fasta, qual, trim, flow, sfftxt, hasAccnos; - int mycount; + string sffFilename, sfftxtFilename, outputDir, accnosName, currentFileName, oligosfile, noMatchFile; + vector filenames, outputNames, accnosFileNames, oligosFileNames; + bool abort, fasta, qual, trim, flow, sfftxt, hasAccnos, hasOligos; + int mycount, split, numFPrimers, numLinkers, numSpacers, pdiffs, bdiffs, ldiffs, sdiffs, tdiffs; set seqNames; - + map barcodes; + map primers; + vector linker, spacer, primerNameVector, barcodeNameVector, revPrimer; + vector > numSplitReads; + vector > filehandles, filehandlesHeaders; + //extract sff file functions - int extractSffInfo(string, string); + int extractSffInfo(string, string, string); int readCommonHeader(ifstream&, CommonHeader&); - int readHeader(ifstream&, Header&); - int readSeqData(ifstream&, seqRead&, int, int); + //int readHeader(ifstream&, Header&); + int readSeqData(ifstream&, seqRead&, int, Header&); int decodeName(string&, string&, string&, string); + bool readOligos(string oligosFile); int printCommonHeader(ofstream&, CommonHeader&); int printHeader(ofstream&, Header&); @@ -98,7 +106,11 @@ private: 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); + string reverseOligo(string oligo); + //parsesfftxt file functions int parseHeaderLineToInt(ifstream&); vector parseHeaderLineToFloatVector(ifstream&, int);