X-Git-Url: https://git.donarmstrong.com/?a=blobdiff_plain;f=trimseqscommand.cpp;h=4b90cda9355112d746999e2a90f9752f5befdf38;hb=e946299389ed49b4afc4262a4e8070ee3d804769;hp=709202f7d690c17e8a6edb7f5a7812047aec9177;hpb=8d341208b0cd0c7a4ff8b77c6b164e06f2f39610;p=mothur.git diff --git a/trimseqscommand.cpp b/trimseqscommand.cpp index 709202f..4b90cda 100644 --- a/trimseqscommand.cpp +++ b/trimseqscommand.cpp @@ -684,15 +684,11 @@ int TrimSeqsCommand::driverCreateTrim(string filename, string qFileName, string //create reoriented primer and barcode pairs map rpairedPrimers, rpairedBarcodes; for (map::iterator it = pairedPrimers.begin(); it != pairedPrimers.end(); it++) { - cout << "primer " << (it->second).forward << '\t' << (it->second).reverse << '\t' << primerNameVector[it->first] << endl; - cout << "rprimer " << trimOligos->reverseOligo((it->second).reverse) << '\t' << (trimOligos->reverseOligo((it->second).forward)) << endl; - oligosPair tempPair(reverseOligo((it->second).reverse), (reverseOligo((it->second).forward))); //reversePrimer, rc ForwardPrimer + oligosPair tempPair(reverseOligo((it->second).reverse), (reverseOligo((it->second).forward))); //reversePrimer, rc ForwardPrimer rpairedPrimers[it->first] = tempPair; } for (map::iterator it = pairedBarcodes.begin(); it != pairedBarcodes.end(); it++) { - cout << "barcode " << (it->second).forward << '\t' << (it->second).reverse << '\t' << barcodeNameVector[it->first] << endl; - cout << "rbarcode " << trimOligos->reverseOligo((it->second).reverse) << '\t' << (trimOligos->reverseOligo((it->second).forward)) << endl; - oligosPair tempPair(reverseOligo((it->second).reverse), (reverseOligo((it->second).forward))); //reverseBarcode, rc ForwardBarcode + oligosPair tempPair(reverseOligo((it->second).reverse), (reverseOligo((it->second).forward))); //reverseBarcode, rc ForwardBarcode rpairedBarcodes[it->first] = tempPair; } rtrimOligos = new TrimOligos(pdiffs, bdiffs, 0, 0, rpairedPrimers, rpairedBarcodes); numBarcodes = rpairedBarcodes.size(); @@ -1561,6 +1557,8 @@ bool TrimSeqsCommand::getOligos(vector >& fastaFileNames, vector< } oligosPair newPrimer(oligo, roligo); + + if (m->debug) { m->mothurOut("[DEBUG]: primer pair " + newPrimer.forward + " " + newPrimer.reverse + ", and group = " + group + ".\n"); } //check for repeat barcodes string tempPair = oligo+roligo; @@ -1604,7 +1602,8 @@ bool TrimSeqsCommand::getOligos(vector >& fastaFileNames, vector< if(reverseBarcode[i] == 'U') { reverseBarcode[i] = 'T'; } } - oligosPair newPair(oligo, reverseOligo(reverseBarcode)); + reverseBarcode = reverseOligo(reverseBarcode); + oligosPair newPair(oligo, reverseBarcode); if (m->debug) { m->mothurOut("[DEBUG]: barcode pair " + newPair.forward + " " + newPair.reverse + ", and group = " + group + ".\n"); }