X-Git-Url: https://git.donarmstrong.com/?a=blobdiff_plain;f=getseqscommand.cpp;h=83fdf73e9af36fd9826f5a8c2efcc11655ac938e;hb=a0f1fca79d2ddfa7ad36b4485039c68b5704fe8d;hp=142023e7a40881c1c6397a761ba350132490b75c;hpb=af0a94ea8f02421b2b73e77e68753a2b4c37768e;p=mothur.git diff --git a/getseqscommand.cpp b/getseqscommand.cpp index 142023e..83fdf73 100644 --- a/getseqscommand.cpp +++ b/getseqscommand.cpp @@ -45,7 +45,7 @@ string GetSeqsCommand::getHelpString(){ helpString += "The get.seqs command reads an .accnos file and any of the following file types: fasta, name, group, count, list, taxonomy, quality or alignreport file.\n"; helpString += "It outputs a file containing only the sequences in the .accnos file.\n"; helpString += "The get.seqs command parameters are accnos, fasta, name, group, list, taxonomy, qfile, alignreport and dups. You must provide accnos unless you have a valid current accnos file, and at least one of the other parameters.\n"; - helpString += "The dups parameter allows you to add the entire line from a name file if you add any name from the line. default=false. \n"; + helpString += "The dups parameter allows you to add the entire line from a name file if you add any name from the line. default=true. \n"; helpString += "The get.seqs command should be in the following format: get.seqs(accnos=yourAccnos, fasta=yourFasta).\n"; helpString += "Example get.seqs(accnos=amazon.accnos, fasta=amazon.fasta).\n"; helpString += "Note: No spaces between parameter labels (i.e. fasta), '=' and parameters (i.e.yourFasta).\n"; @@ -430,9 +430,11 @@ int GetSeqsCommand::readFasta(){ if (!dups) {//adjust name if needed map::iterator it = uniqueMap.find(name); - if (it != uniqueMap.end()) { name = it->second; } + if (it != uniqueMap.end()) { currSeq.setName(it->second); } } + name = currSeq.getName(); + if (name != "") { //if this name is in the accnos file if (names.count(name) != 0) { @@ -710,7 +712,7 @@ int GetSeqsCommand::readName(){ if (m->control_pressed) { in.close(); out.close(); m->mothurRemove(outputFileName); return 0; } - in >> firstCol; + in >> firstCol; m->gobble(in); in >> secondCol; string hold = ""; @@ -734,6 +736,7 @@ int GetSeqsCommand::readName(){ selectedCount += parsedNames.size(); if (m->debug) { sanity["name"].insert(firstCol); } }else { + selectedCount += validSecond.size(); //if the name in the first column is in the set then print it and any other names in second column also in set @@ -752,6 +755,7 @@ int GetSeqsCommand::readName(){ //make first name in set you come to first column and then add the remaining names to second column }else { + //you want part of this row if (validSecond.size() != 0) {