X-Git-Url: http://source.jalview.org/gitweb/?a=blobdiff_plain;f=test%2Fjalview%2Fdatamodel%2FHiddenMarkovModelTest.java;h=65241b4d0c483b1e2df35f9a9611c300f37566d8;hb=f452bcb5ddecb4378e952447031ddb076bdf8268;hp=7978da273c6cab041713caf7e93cfd62f4cb88c9;hpb=64144b92481279aa0aaf11ba05febd65f73cfa30;p=jalview.git

diff --git a/test/jalview/datamodel/HiddenMarkovModelTest.java b/test/jalview/datamodel/HiddenMarkovModelTest.java
index 7978da2..65241b4 100644
--- a/test/jalview/datamodel/HiddenMarkovModelTest.java
+++ b/test/jalview/datamodel/HiddenMarkovModelTest.java
@@ -18,6 +18,8 @@ public class HiddenMarkovModelTest {
 
   HiddenMarkovModel hmm;
 
+  HiddenMarkovModel alignmentHmm;
+
   @BeforeClass(alwaysRun = true)
   public void setUp() throws MalformedURLException, IOException
   {
@@ -27,6 +29,12 @@ public class HiddenMarkovModelTest {
     HMMFile file = new HMMFile(new FileParse(
             "test/jalview/io/test_PKinase_hmm.txt", DataSourceType.FILE));
     hmm = file.getHMM();
+
+    // used to check if consensus sequence is automatically aligned with alignment
+    HMMFile alignmentTest = new HMMFile(
+            new FileParse("test/jalview/io/HMMAlignmentTestHMM.hmm",
+                    DataSourceType.FILE));
+    alignmentHmm = alignmentTest.getHMM();
   }
 
   @Test(groups = "Functional")
@@ -115,6 +123,14 @@ public class HiddenMarkovModelTest {
     assertEquals(subStr, "yelleklGsG");
     subStr = seq.getSequenceAsString().substring(150, 161);
     assertEquals(subStr, "-dllk------");
+
+    // test to see if consensus sequence maps to alignment correctly
+    // see HMMAlignmentTest.sto for corresponding alignment file
+    SequenceI seq2 = alignmentHmm.getConsensusSequence();
+    assertEquals(seq2.getCharAt(0), '-');
+    assertEquals(seq2.getCharAt(7), '-');
+    assertEquals(seq2.getCharAt(8), 's');
+
   }
 
   /**