helpString += "The dist.seqs command should be in the following format: \n";
helpString += "dist.seqs(fasta=yourFastaFile, calc=yourCalc, countends=yourEnds, cutoff= yourCutOff, processors=yourProcessors) \n";
helpString += "Example dist.seqs(fasta=amazon.fasta, calc=eachgap, countends=F, cutoff= 2.0, processors=3).\n";
- helpString += "Note: No spaces between parameter labels (i.e. calc), '=' and parameters (i.e.yourCalc).\n\n";
+ helpString += "Note: No spaces between parameter labels (i.e. calc), '=' and parameters (i.e.yourCalc).\n";
return helpString;
}
catch(exception& e) {
//allow user to run help
if(option == "help") { help(); abort = true; calledHelp = true; }
+ else if(option == "citation") { citation(); abort = true; calledHelp = true;}
else {
vector<string> myArray = setParameters();
m->openInputFile(fastafile, inFASTA);
alignDB = SequenceDB(inFASTA);
inFASTA.close();
+ m->setFastaFile(fastafile);
}
oldfastafile = validParameter.validFile(parameters, "oldfasta", true);
column = validParameter.validFile(parameters, "column", true);
if (column == "not found") { column = ""; }
else if (column == "not open") { abort = true; }
+ else { m->setColumnFile(column); }
//if the user changes the output directory command factory will send this info to us in the output parameter
outputDir = validParameter.validFile(parameters, "outputdir", false); if (outputDir == "not found"){
int numSeqs = alignDB.getNumSeqs();
cutoff += 0.005;
+ if (!alignDB.sameLength()) { m->mothurOut("[ERROR]: your sequences are not the same length, aborting."); m->mothurOutEndLine(); return 0; }
+
string outputFile;
if (output == "lt") { //does the user want lower triangle phylip formatted file