added list.seqs command

[mothur.git] / blastalign.cpp

diff --git a/blastalign.cpp b/blastalign.cpp
index 0a1350fac0c7fa655d5d04585156c28445c4d0d5..0645fb031f0d2a4e98b4caba5695f30ecd555474 100644 (file)
--- a/blastalign.cpp
+++ b/blastalign.cpp
@@ -11,7 +11,6 @@
  *
  */
 
-using namespace std;
 
 #include "alignment.hpp"
 #include "blastalign.hpp"
@@ -23,6 +22,10 @@ BlastAlignment::BlastAlignment(float go, float ge, float m, float mm) :
                        match(m),                               //      This is the score to award for two nucleotides matching (match >= 0)
                        mismatch(mm)                    //      This is the penalty to assess for a mismatch (mismatch <= 0)
 {
+       globaldata = GlobalData::getInstance();
+       path = globaldata->argv;
+       path = path.substr(0, (path.find_last_of('m')));
+       
        gapOpen = abs(go);                              //      This is the penalty to assess for opening a gap (gapOpen >= 0)
        gapExtend = abs(ge);                            //      This is the penalty to assess for extending a gap (gapExtend >= 0)
                
@@ -54,7 +57,7 @@ void BlastAlignment::align(string seqA, string seqB){ //Use blastn to align the
        
        //      The blastCommand assumes that we have DNA sequences (blastn) and that they are fairly similar (-e 0.001) and
        //      that we don't want to apply any kind of complexity filtering (-F F)
-       string blastCommand = "~/Pipeline/src/cpp/production/blast/bin/bl2seq -p blastn -i " + candidateFileName + " -j " + templateFileName + " -e 0.0001 -F F -o " + blastFileName + " -W 11";
+       string blastCommand = path + "blast/bin/bl2seq -p blastn -i " + candidateFileName + " -j " + templateFileName + " -e 0.0001 -F F -o " + blastFileName + " -W 11";
        blastCommand += " -r " + toString(match) + " -q " + toString(mismatch);
        blastCommand += " -G " + toString(gapOpen) + " -E " + toString(gapExtend);
        
@@ -79,10 +82,10 @@ void BlastAlignment::setPairwiseSeqs(){     //      This method call assigns the blast ge
        
        string candidateName, templateName;
        
-       while(d=blastFile.get() != '='){};
+       while(d=blastFile.get() != '='){}
        blastFile >> candidateName;                                     //      Get the candidate sequence name from flatfile
        
-       while(d=blastFile.get() != '('){};
+       while(d=blastFile.get() != '('){}
        blastFile >> candidateLength;                           //      Get the candidate sequence length from flatfile
        
        while(d=blastFile.get()){
@@ -99,9 +102,9 @@ void BlastAlignment::setPairwiseSeqs(){      //      This method call assigns the blast ge
                        pairwiseLength = 0;
                        
 //                     string dummy;
-//                     while(dummy != "query:"){       cout << dummy << endl;blastFile >> dummy;       }
+//                     while(dummy != "query:"){       mothurOut(dummy, ""); mothurOutEndLine(); blastFile >> dummy;   }
 //                     blastFile >> seqBend;
-//                     cout << seqBend << endl;
+//                     mothurOut(toString(seqBend), ""); mothurOutEndLine();
 //                     for(int i=0;i<seqBend;i++){
 //                             seqAaln += 'Z';
 //                             seqBaln += 'X';
@@ -111,10 +114,10 @@ void BlastAlignment::setPairwiseSeqs(){   //      This method call assigns the blast ge
                }
        }
        
-       while(d=blastFile.get() != '='){};
+       while(d=blastFile.get() != '='){}
        blastFile >> templateLength;                            //      Get the template sequence length from flatfile
                
-       while(d=blastFile.get() != 'Q'){};                      //      Suck up everything else until we get to the start of the alignment
+       while(d=blastFile.get() != 'Q'){}                       //      Suck up everything else until we get to the start of the alignment
        int queryStart, sbjctStart, queryEnd, sbjctEnd;
        string queryLabel, sbjctLabel, query, sbjct;
 
@@ -151,6 +154,7 @@ void BlastAlignment::setPairwiseSeqs(){     //      This method call assigns the blast ge
                seqAaln += 'Z';                                                 //      again need ot pad the sequences so that they extend to the length
                seqBaln += 'X';                                                 //      of the template sequence
        }
+       blastFile.close();
 }
 
 //**************************************************************************************************/