don't seqcomp sequences that are too small

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

diff --git a/alg/test/test_mussa.cpp b/alg/test/test_mussa.cpp
index 7eb301572777dc2e9037cfbe4410c3fe6ea416e9..4e0340df08123ebecb66ef2ff4e904f3f88866b3 100644 (file)
--- a/alg/test/test_mussa.cpp
+++ b/alg/test/test_mussa.cpp
@@ -475,6 +475,19 @@ BOOST_AUTO_TEST_CASE( three_way_local_alignment )
   }   
 }
 
+BOOST_AUTO_TEST_CASE( mussa_window_larger_than_sequence )
+{
+  string s0("AGCAGGG");
+  string s1("CAGCGGG");
+  
+  Mussa analysis;
+  analysis.append_sequence(s0);
+  analysis.append_sequence(s1);
+  analysis.set_threshold(23);
+  analysis.set_window(30);
+  BOOST_CHECK_THROW(analysis.analyze(), seqcomp_error);
+}
+
 BOOST_AUTO_TEST_CASE( subanalysis )
 {
   Sequence s1("AATGAAGATTTTAATGCTTTAATTTTGTTTTGTAAACTTCGAATTTCCAAAATTTGAAA");