make sure subsequences have a species name

[mussa.git] / alg / sequence.cpp

diff --git a/alg/sequence.cpp b/alg/sequence.cpp
index 7523ab679f3278ab9c2a9446dfb3497b8f1949c9..65ff57d2e326c42bea38346fabcdc2966d34e6a3 100644 (file)
--- a/alg/sequence.cpp
+++ b/alg/sequence.cpp
@@ -457,7 +457,7 @@ Sequence::subseq(int start, int count) const
   }
   Sequence new_seq(std::string::substr(start, count));
   new_seq.set_fasta_header(get_fasta_header());
-  //new_seq.set_species(get_species());
+  new_seq.set_species(get_species());
 
   new_seq.motif_list = motif_list;
   // attempt to copy & reannotate position based annotations