load and display a motif list

[mussa.git] / alg / test / test_mussa.cpp

diff --git a/alg/test/test_mussa.cpp b/alg/test/test_mussa.cpp
index 6d9036fdb05cb41e76a0c790ae73a836ed51904c..1ca85a4144fed65fc055b99fd4d2e75c7224cdf9 100644 (file)
--- a/alg/test/test_mussa.cpp
+++ b/alg/test/test_mussa.cpp
@@ -1,9 +1,13 @@
 #include <boost/test/auto_unit_test.hpp>
 
 #include <string>
+#include <sstream>
+#include <vector>
 
 #include "alg/mussa.hpp"
 
+using namespace std;
+
 //! can we initialize a mussa object?
 BOOST_AUTO_TEST_CASE( mussa_simple )
 {
@@ -106,3 +110,24 @@ BOOST_AUTO_TEST_CASE( mussa_load_full_path )
 
   BOOST_CHECK( m1.size() > 0);
 }
+
+BOOST_AUTO_TEST_CASE( mussa_load_motif )
+{
+  string data = "AAGG 1.0 1.0 0.0\n"
+                "GGTT 0.0 0.1 1.0\n"
+                "ZXY 2 1.9 0\n";
+
+  istringstream test_istream(data);
+
+  Mussa m1;
+  m1.add_a_seq("AAAAGGGGTTTT");
+  m1.add_a_seq("GGGCCCCTTGGTT");
+  m1.load_motifs(test_istream);
+
+  for (vector<Sequence>::const_iterator seq_i = m1.sequences().begin();
+       seq_i != m1.sequences().end();
+       ++seq_i)
+  {
+    BOOST_CHECK( seq_i->motifs().size() > 0 );
+  }
+}