3 import static org.testng.Assert.assertEquals;
4 import static org.testng.Assert.assertNotNull;
5 import static org.testng.Assert.assertTrue;
6 import static org.testng.Assert.fail;
8 import jalview.bin.Jalview;
9 import jalview.datamodel.AlignmentI;
10 import jalview.datamodel.HiddenMarkovModel;
11 import jalview.datamodel.SequenceI;
12 import jalview.gui.AlignFrame;
13 import jalview.gui.Desktop;
14 import jalview.io.HMMFile;
15 import jalview.ws.params.ArgumentI;
17 import java.io.IOException;
18 import java.net.MalformedURLException;
19 import java.util.ArrayList;
21 import org.testng.annotations.AfterClass;
22 import org.testng.annotations.BeforeClass;
23 import org.testng.annotations.Test;
25 public class HMMERTest {
29 @BeforeClass(alwaysRun = true)
30 public void setUpBeforeClass() throws Exception
34 { "-noquestionnaire", "-nonews", "-props",
35 "test/jalview/hmmer/testProps.jvprops", "-open",
36 "examples/uniref50.fa" });
37 // FastaFile file = null;
38 // file = new FastaFile(
39 // new FileParse("examples/uniref50.fa", DataSourceType.FILE));
40 // SequenceI[] seqs = file.getSeqsAsArray();
41 // AlignmentI al = new Alignment(seqs);
42 // frame = new AlignFrame(al, 150, 20);
43 frame = Desktop.getAlignFrames()[0];
46 @AfterClass(alwaysRun = true)
47 public static void tearDownAfterClass() throws Exception
49 Desktop.instance.closeAll_actionPerformed(null);
53 * Test with a dependency on locally installed hmmbuild binaries
55 * @throws MalformedURLException
58 @Test(groups = "External")
59 public void testHMMBuildThenHMMAlign()
60 throws MalformedURLException, IOException
63 * run hmmbuild - note the side-effect of selecting the HMM
64 * sequence that gets added to the alignment
67 HiddenMarkovModel hmm = frame.getSelectedHMM();
71 * now run hmmalign - with respect to the select HMM profile
76 public void testHMMBuild()
78 HMMBuild builder = new HMMBuild(frame,
79 new ArrayList<ArgumentI>());
82 SequenceI seq = frame.getViewport().getAlignment().getSequenceAt(0);
83 HiddenMarkovModel hmm = seq.getHMM();
86 assertEquals(hmm.getLength(), 148);
87 assertEquals(hmm.getAlphabetType(), "amino");
88 assertEquals(hmm.getName(), "Alignment_HMM");
89 assertEquals(hmm.getProperty(HMMFile.EFF_NUMBER_OF_SEQUENCES),
91 assertEquals(hmm.getConsensusAtAlignColumn(15), 's');
94 public void testHMMAlign()
96 HMMAlign thread = new HMMAlign(frame,
97 new ArrayList<ArgumentI>());
100 AlignFrame[] alignFrames = Desktop.getAlignFrames();
101 if (alignFrames == null)
103 fail("No align frame loaded");
107 * now have the original align frame, and another for realigned sequences
109 assertEquals(alignFrames.length, 2);
110 AlignmentI original = alignFrames[0].getViewport().getAlignment();
111 assertNotNull(original);
112 AlignmentI realigned = alignFrames[1].getViewport().getAlignment();
113 assertNotNull(realigned);
114 assertNotNull(original.getHmmConsensus());
115 assertNotNull(realigned.getHmmConsensus());
117 SequenceI ferCapan = original.findName("FER_CAPAN");
118 assertTrue(ferCapan.getSequenceAsString().startsWith("MA------SVSAT"));
120 SequenceI ferCapanRealigned = realigned.findName("FER_CAPAN");
121 assertTrue(ferCapanRealigned.getSequenceAsString()
122 .startsWith("-------m-A----SVSAT"));