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";
+ "GGTT 0.0 0.1 1.0\n";
istringstream test_istream(data);
m1.append_sequence("GGGCCCCTTCCAATT");
m1.load_motifs(test_istream);
+ BOOST_CHECK_EQUAL( m1.motifs().size(), 2);
for (Mussa::vector_sequence_type::const_iterator seq_i = m1.sequences().begin();
seq_i != m1.sequences().end();
++seq_i)
}
}
+BOOST_AUTO_TEST_CASE( mussa_load_broken_motif )
+{
+ string data = "AAGG 1.0 1.0 0.0\n"
+ "GGTT 0.0 0.1 1.0 1.0\n"
+ "ZZCTA 0.1 0.0 1.0\n";
+
+ istringstream test_istream(data);
+
+ Mussa m1;
+ m1.append_sequence("AAAAGGGGTTTT");
+ m1.append_sequence("GGGCCCCTTCCAATT");
+ BOOST_CHECK_THROW(m1.load_motifs(test_istream), motif_load_error);
+
+ BOOST_CHECK_EQUAL( m1.motifs().size(), 0);
+}
+
BOOST_AUTO_TEST_CASE( mussa_named_motif )
{
string data = "CCAATT cat 0.1 0.2 0.3\n";