]> git.donarmstrong.com Git - mothur.git/blobdiff - distancecommand.cpp
paralellized the indicator command
[mothur.git] / distancecommand.cpp
index 34206a9460bf21eac543ee3d91031f8fe4af96ad..4bdeaaeb5d9b388ac56fc2553f57dc10c0d598d9 100644 (file)
@@ -152,6 +152,7 @@ DistanceCommand::DistanceCommand(string option) {
                                m->openInputFile(fastafile, inFASTA);
                                alignDB = SequenceDB(inFASTA); 
                                inFASTA.close();
+                               m->setFastaFile(fastafile);
                        }
                        
                        oldfastafile = validParameter.validFile(parameters, "oldfasta", true);
@@ -161,6 +162,7 @@ DistanceCommand::DistanceCommand(string option) {
                        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"){  
@@ -247,6 +249,8 @@ int DistanceCommand::execute(){
                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