X-Git-Url: https://git.donarmstrong.com/?a=blobdiff_plain;f=sffinfocommand.cpp;h=c50255aeb2ae9c202074b78d89d460651cefdc43;hb=e10c72304ee071c0c40e0218a06d89dc4731cbc2;hp=d7ca230449fbebd3488a512c5fbd25d3159c1156;hpb=f509429e06e545bde69c97cacc0eb436775bd329;p=mothur.git diff --git a/sffinfocommand.cpp b/sffinfocommand.cpp index d7ca230..c50255a 100644 --- a/sffinfocommand.cpp +++ b/sffinfocommand.cpp @@ -1037,7 +1037,7 @@ int SffInfoCommand::readSeqData(ifstream& in, seqRead& read, int numFlowReads, H int SffInfoCommand::findGroup(Header header, seqRead read, int& barcode, int& primer) { try { //find group read belongs to - TrimOligos trimOligos(pdiffs, bdiffs, ldiffs, sdiffs, primers, barcodes, rbarcodes, revPrimer, linker, spacer); + TrimOligos trimOligos(pdiffs, bdiffs, ldiffs, sdiffs, primers, barcodes, revPrimer, linker, spacer); int success = 1; string trashCode = ""; @@ -1082,12 +1082,6 @@ int SffInfoCommand::findGroup(Header header, seqRead read, int& barcode, int& pr else{ currentSeqsDiffs += success; } } - if(rbarcodes.size() != 0){ - success = trimOligos.stripRBarcode(currSeq, currQual, barcode); - if(success > bdiffs) { trashCode += 'b'; } - else{ currentSeqsDiffs += success; } - } - if(numSpacers != 0){ success = trimOligos.stripSpacer(currSeq, currQual); if(success > sdiffs) { trashCode += 's'; } @@ -1675,36 +1669,13 @@ bool SffInfoCommand::readOligos(string oligoFile){ } else if(type == "BARCODE"){ inOligos >> group; - - //barcode lines can look like BARCODE atgcatgc groupName - for 454 seqs - //or BARCODE atgcatgc atgcatgc groupName - for illumina data that has forward and reverse info - string temp = ""; - while (!inOligos.eof()) { - char c = inOligos.get(); - if (c == 10 || c == 13){ break; } - else if (c == 32 || c == 9){;} //space or tab - else { temp += c; } - } - //then this is illumina data with 4 columns - if (temp != "") { - string reverseBarcode = reverseOligo(group); //reverse barcode - group = temp; - - //check for repeat barcodes - map::iterator itBar = rbarcodes.find(reverseBarcode); - if (itBar != rbarcodes.end()) { m->mothurOut("barcode " + reverseBarcode + " is in your oligos file already."); m->mothurOutEndLine(); } - - rbarcodes[reverseBarcode]=indexBarcode; - } - //check for repeat barcodes map::iterator itBar = barcodes.find(oligo); if (itBar != barcodes.end()) { m->mothurOut("barcode " + oligo + " is in your oligos file already."); m->mothurOutEndLine(); } barcodes[oligo]=indexBarcode; indexBarcode++; barcodeNameVector.push_back(group); - }else if(type == "LINKER"){ linker.push_back(oligo); }else if(type == "SPACER"){