Merge branch 'develop' into update_212_Dec_merge_with_21125_chamges
[jalview.git] / test / jalview / renderer / seqfeatures / FeatureRendererTest.java
index 723f3b8..e1dac7d 100644 (file)
@@ -22,34 +22,47 @@ package jalview.renderer.seqfeatures;
 
 import static org.testng.Assert.assertEquals;
 import static org.testng.Assert.assertFalse;
+import static org.testng.Assert.assertNotNull;
 import static org.testng.Assert.assertNull;
+import static org.testng.Assert.assertSame;
 import static org.testng.Assert.assertTrue;
 
+import java.awt.Color;
+import java.util.ArrayList;
+import java.util.Arrays;
+import java.util.HashMap;
+import java.util.List;
+import java.util.Map;
+
+import org.testng.annotations.AfterMethod;
+import org.testng.annotations.Test;
+
+import jalview.analysis.GeneticCodes;
 import jalview.api.AlignViewportI;
 import jalview.api.FeatureColourI;
+import jalview.bin.Jalview;
+import jalview.datamodel.MappedFeatures;
 import jalview.datamodel.SequenceFeature;
 import jalview.datamodel.SequenceI;
 import jalview.datamodel.features.FeatureMatcher;
 import jalview.datamodel.features.FeatureMatcherSet;
 import jalview.datamodel.features.FeatureMatcherSetI;
 import jalview.gui.AlignFrame;
+import jalview.gui.AlignViewport;
+import jalview.gui.Desktop;
 import jalview.io.DataSourceType;
 import jalview.io.FileLoader;
 import jalview.schemes.FeatureColour;
 import jalview.util.matcher.Condition;
 import jalview.viewmodel.seqfeatures.FeatureRendererModel.FeatureSettingsBean;
 
-import java.awt.Color;
-import java.util.ArrayList;
-import java.util.Arrays;
-import java.util.HashMap;
-import java.util.List;
-import java.util.Map;
-
-import org.testng.annotations.Test;
-
 public class FeatureRendererTest
 {
+  @AfterMethod(alwaysRun = true)
+  public void tearDown()
+  {
+    Desktop.getInstance().closeAll_actionPerformed(null);
+  }
 
   @Test(groups = "Functional")
   public void testFindAllFeatures()
@@ -456,7 +469,8 @@ public class FeatureRendererTest
     gc = new FeatureColour(Color.red, Color.yellow, Color.red, Color.green,
             5f, 11f);
     fr.getFeatureColours().put("Cath", gc);
-    gc.setAutoScaled(false); // this does little other than save a checkbox setting!
+    gc.setAutoScaled(false); // this does little other than save a checkbox
+                             // setting!
     assertEquals(fr.getColour(sf2), new Color(255, 213, 0));
 
     /*
@@ -604,4 +618,69 @@ public class FeatureRendererTest
     filter.and(FeatureMatcher.byScore(Condition.LE, "0.4"));
     assertFalse(fr.isVisible(sf)); // feature doesn't match filter
   }
+
+  @Test(groups = "Functional")
+  public void testFindComplementFeaturesAtResidue()
+  {
+    Jalview.main(
+            new String[]
+            { "-nonews", "-props", "test/jalview/testProps.jvprops" });
+
+    // codons for MCWHSE
+    String cdsSeq = ">cds\nATGtgtTGGcacTCAgaa";
+    AlignFrame af = new FileLoader().LoadFileWaitTillLoaded(cdsSeq,
+            DataSourceType.PASTE);
+    af.showTranslation_actionPerformed(
+            GeneticCodes.getInstance().getStandardCodeTable());
+    af.closeMenuItem_actionPerformed(true);
+
+    /*
+     * find the complement frames (ugly)
+     */
+    AlignFrame[] frames = Desktop.getAlignFrames();
+    assertEquals(frames.length, 2);
+    AlignViewport av1 = frames[0].getViewport();
+    AlignViewport av2 = frames[1].getViewport();
+    AlignViewport cds = av1.getAlignment().isNucleotide() ? av1 : av2;
+    AlignViewport peptide = cds == av1 ? av2 : av1;
+    assertNotNull(cds);
+    assertNotNull(peptide);
+
+    /*
+     * add features to CDS at first codon, positions 2-3
+     */
+    SequenceI seq1 = cds.getAlignment().getSequenceAt(0);
+    SequenceFeature sf1 = new SequenceFeature("sequence_variant", "G,GT", 2,
+            2, "ensembl");
+    seq1.addSequenceFeature(sf1);
+    SequenceFeature sf2 = new SequenceFeature("sequence_variant", "C, CA",
+            3, 3, "ensembl");
+    seq1.addSequenceFeature(sf2);
+
+    /*
+     * 'find' mapped features from the peptide position
+     * - first with CDS features _not_ shown on peptide alignment
+     */
+    SequenceI seq2 = peptide.getAlignment().getSequenceAt(0);
+    FeatureRenderer frC = new FeatureRenderer(cds);
+    frC.featuresAdded();
+    MappedFeatures mf = frC.findComplementFeaturesAtResidue(seq2, 1);
+    assertNotNull(mf);
+    assertEquals(mf.features.size(), 2);
+    assertSame(mf.features.get(0), sf1);
+    assertSame(mf.features.get(1), sf2);
+
+    /*
+     * add exon feature and verify it is only returned once for a
+     * peptide position, even though it is on all 3 codon positions
+     */
+    SequenceFeature sf3 = new SequenceFeature("exon", "exon1", 4, 12,
+            "ensembl");
+    seq1.addSequenceFeature(sf3);
+    frC.featuresAdded();
+    mf = frC.findComplementFeaturesAtResidue(seq2, 3);
+    assertNotNull(mf);
+    assertEquals(mf.features.size(), 1);
+    assertSame(mf.features.get(0), sf3);
+  }
 }