allow specifying a name for a sequence

[mussa.git] / alg / mussa.cpp

diff --git a/alg/mussa.cpp b/alg/mussa.cpp
index 2391b0478080863d196271df6d1204e04d175787..5acd58e6bdfeaf1a9081ee3d93c5351be79054ba 100644 (file)
--- a/alg/mussa.cpp
+++ b/alg/mussa.cpp
@@ -274,13 +274,17 @@ Mussa::sequences() const
 }
 
 void Mussa::load_sequence(fs::path seq_file, fs::path annot_file, 
-                          int fasta_index, int sub_seq_start, int sub_seq_end)
+                          int fasta_index, int sub_seq_start, int sub_seq_end,
+                          std::string *name)
 {
   boost::shared_ptr<Sequence> aseq(new Sequence);
   aseq->load_fasta(seq_file, fasta_index, sub_seq_start, sub_seq_end);
   if ( not annot_file.empty() ) {
     aseq->load_annot(annot_file, sub_seq_start, sub_seq_end);
   }
+  if (name != 0 and name->size() > 0 ) {
+    aseq->set_species(*name);
+  }
   the_seqs.push_back(aseq);
 }