@@ -292,7 +292,7 @@ phylogeny::Tree* msa::io::readAlignments_and_buildTree(SequenceDB *database, Opt
292292
293293
294294 std::cerr << " Read " << T->allNodes .size () << " MSA files in " << alnReadTime.count () / 1e6 << " ms.\n " ;
295- std::cerr << " ============================\n " ;
295+ // std::cerr << "============================\n";
296296 return T;
297297}
298298
@@ -304,7 +304,7 @@ void msa::io::readBackboneAlignment(Tree* T, SequenceDB *database, Option *optio
304304 T->root ->seqsIncluded .push_back (subtreeIdx);
305305 auto alnReadEnd = std::chrono::high_resolution_clock::now ();
306306 std::chrono::nanoseconds alnReadTime = alnReadEnd - alnReadStart;
307- std::cerr << " Loaded backbone alignment in " << alnReadTime.count () / 1e6 << " ms. (Count: " << T->root ->alnNum << " , Length: " << T->root ->msaFreq .size () << " )\n " ;
307+ std::cerr << " Backbone alignment read in " << alnReadTime.count () / 1000000 << " ms. (Count: " << T->root ->alnNum << " , Length: " << T->root ->msaFreq .size () << " )\n " ;
308308 return ;
309309}
310310
@@ -357,7 +357,7 @@ int msa::io::update_and_writeAlignment(SequenceDB* database, Option* option, std
357357 bool nochange = false ;
358358 char gapType = (option->alnMode == PLACE_WO_TREE) ? ' .' : ' -' ;
359359 stringPairVec seqs_before, seqs_after;
360- if (option->alnMode == 0 && option->compressed ) fileName += " .gz" ;
360+ if (option->alnMode == DEFAULT_ALN && option->compressed ) fileName += " .gz" ;
361361 fs::path p (fileName);
362362 std::string filename = p.stem ().string ();
363363 std::string finalAlnFileName = option->tempDir + " /" + filename + " .final.aln" ;
@@ -431,6 +431,7 @@ int msa::io::update_and_writeAlignment(SequenceDB* database, Option* option, std
431431 seqs_after.clear ();
432432 seqs_before.clear ();
433433 }
434+ if (option->alnMode == PLACE_WO_TREE) std::cerr << " Final Alignment Length: " << database->subtreeAln [subtreeIdx].size () << ' \n ' ;
434435 database->subtreeAln [subtreeIdx].clear ();
435436 if (nochange) {
436437 std::string command = " cp " + fileName + " " + finalAlnFileName;
0 commit comments