]> git.donarmstrong.com Git - mothur.git/commitdiff
changes for chimera slayer
authorwestcott <westcott>
Wed, 4 May 2011 15:25:08 +0000 (15:25 +0000)
committerwestcott <westcott>
Wed, 4 May 2011 15:25:08 +0000 (15:25 +0000)
blastdb.cpp
chimeraslayer.cpp

index cc85989dad80dd1b0f85997b3428006b7d353c5c..40855d0a262a49bf64488ec53f8f9b58b1beffea 100644 (file)
@@ -141,7 +141,7 @@ vector<int> BlastDB::findClosestMegaBlast(Sequence* seq, int n, int minPerID) {
                
                while(!m8FileHandle.eof()){
                        m8FileHandle >> dummy >> templateAccession >> searchScore >> numBases >> mismatch >> gap >> startQuery >> endQuery >> startRef >> endRef >> eScore >> score;
-                       //cout << dummy << '\t' << templateAccession << '\t' << searchScore << '\t';
+                       //cout << dummy << '\t' << templateAccession << '\t' << searchScore << '\t' << numBases << '\t' << mismatch << '\t' << gap << '\t' << startQuery << '\t' << endQuery << '\t' << startRef << '\t' << endRef << '\t' << eScore << '\t' << score << endl; 
                        
                        //get rest of junk in line
                        //while (!m8FileHandle.eof())   {       char c = m8FileHandle.get(); if (c == 10 || c == 13){   break;  }else{ cout << c; }     } //
index cd268d3d6f1a513e932940332bab8ed15c59ada1..5b1075902d6a038a4e863913bee7969c0ae52154 100644 (file)
@@ -1030,7 +1030,7 @@ vector<Sequence*> ChimeraSlayer::getBlastSeqs(Sequence* q, vector<Sequence*>& db
                
                vector<int> tempIndexesLeft = databaseLeft->findClosestMegaBlast(queryLeft, num+1, minSim);
                vector<int> tempIndexesRight = databaseLeft->findClosestMegaBlast(queryRight, num+1, minSim);
-               cout << q->getName() << '\t' << leftQuery << '\t' << "leftMatches = " << tempIndexesLeft.size() << '\t' << rightQuery   << " rightMatches = " << tempIndexesRight.size() << endl;
+               //cout << q->getName() << '\t' << leftQuery << '\t' << "leftMatches = " << tempIndexesLeft.size() << '\t' << rightQuery << " rightMatches = " << tempIndexesRight.size() << endl;
                vector<int> smaller;
                vector<int> larger;