X-Git-Url: http://source.jalview.org/gitweb/?a=blobdiff_plain;f=test%2Fjalview%2Frenderer%2Fseqfeatures%2FFeatureRendererTest.java;h=a3d0a7ca4c8a449660d10daa71deb5693e117cad;hb=57738a1f3c19b1c3a00bd3ac5108f8cd0af32f99;hp=a9e3754615f2873f93187fb877c3a767f8f5d72d;hpb=ef7d9bed59341ba8a035dc4578cd9c62b655286b;p=jalview.git diff --git a/test/jalview/renderer/seqfeatures/FeatureRendererTest.java b/test/jalview/renderer/seqfeatures/FeatureRendererTest.java index a9e3754..a3d0a7c 100644 --- a/test/jalview/renderer/seqfeatures/FeatureRendererTest.java +++ b/test/jalview/renderer/seqfeatures/FeatureRendererTest.java @@ -22,17 +22,24 @@ 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 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; @@ -285,8 +292,8 @@ public class FeatureRendererTest * give "Type3" features a graduated colour scheme * - first with no threshold */ - FeatureColourI gc = new FeatureColour(Color.yellow, Color.red, null, 0f, - 10f); + FeatureColourI gc = new FeatureColour(Color.green, Color.yellow, + Color.red, null, 0f, 10f); fr.getFeatureColours().put("Type3", gc); features = fr.findFeaturesAtColumn(seq, 8); assertTrue(features.contains(sf4)); @@ -428,8 +435,8 @@ public class FeatureRendererTest * graduated colour by score, no threshold, no score * */ - FeatureColourI gc = new FeatureColour(Color.yellow, Color.red, - Color.green, 1f, 11f); + FeatureColourI gc = new FeatureColour(Color.red, Color.yellow, + Color.red, Color.green, 1f, 11f); fr.getFeatureColours().put("Cath", gc); assertEquals(fr.getColour(sf1), Color.green); @@ -453,9 +460,11 @@ public class FeatureRendererTest * threshold is min-max; now score 6 is 1/6 of the way from 5 to 11 * or from yellow(255, 255, 0) to red(255, 0, 0) */ - gc = new FeatureColour(Color.yellow, Color.red, Color.green, 5f, 11f); + 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)); /* @@ -476,7 +485,8 @@ public class FeatureRendererTest * colour by feature attribute value * first with no value held */ - gc = new FeatureColour(Color.yellow, Color.red, Color.green, 1f, 11f); + gc = new FeatureColour(Color.red, Color.yellow, Color.red, Color.green, + 1f, 11f); fr.getFeatureColours().put("Cath", gc); gc.setAttributeName("AF"); assertEquals(fr.getColour(sf2), Color.green); @@ -535,4 +545,136 @@ public class FeatureRendererTest csqData.put("Feature", "ENST01234"); assertEquals(fr.getColour(sf2), expected); } + + @Test(groups = "Functional") + public void testIsVisible() + { + String seqData = ">s1\nMLQGIFPRS\n"; + AlignFrame af = new FileLoader().LoadFileWaitTillLoaded(seqData, + DataSourceType.PASTE); + AlignViewportI av = af.getViewport(); + FeatureRenderer fr = new FeatureRenderer(av); + SequenceI seq = av.getAlignment().getSequenceAt(0); + SequenceFeature sf = new SequenceFeature("METAL", "Desc", 10, 10, 1f, + "Group"); + sf.setValue("AC", "11"); + sf.setValue("CLIN_SIG", "Likely Pathogenic"); + seq.addSequenceFeature(sf); + + assertFalse(fr.isVisible(null)); + + /* + * initial state FeatureRenderer hasn't 'found' feature + * and so its feature type has not yet been set visible + */ + assertFalse(fr.getDisplayedFeatureCols().containsKey("METAL")); + assertFalse(fr.isVisible(sf)); + + fr.findAllFeatures(true); + assertTrue(fr.isVisible(sf)); + + /* + * feature group not visible + */ + fr.setGroupVisibility("Group", false); + assertFalse(fr.isVisible(sf)); + fr.setGroupVisibility("Group", true); + assertTrue(fr.isVisible(sf)); + + /* + * feature score outwith colour threshold (score > 2) + */ + FeatureColourI fc = new FeatureColour(null, Color.white, Color.black, + Color.white, 0, 10); + fc.setAboveThreshold(true); + fc.setThreshold(2f); + fr.setColour("METAL", fc); + assertFalse(fr.isVisible(sf)); // score 1 is not above threshold 2 + fc.setBelowThreshold(true); + assertTrue(fr.isVisible(sf)); // score 1 is below threshold 2 + + /* + * colour with threshold on attribute AC (value is 11) + */ + fc.setAttributeName("AC"); + assertFalse(fr.isVisible(sf)); // value 11 is not below threshold 2 + fc.setAboveThreshold(true); + assertTrue(fr.isVisible(sf)); // value 11 is above threshold 2 + + fc.setAttributeName("AF"); // attribute AF is absent in sf + assertTrue(fr.isVisible(sf)); // feature is not excluded by threshold + + FeatureMatcherSetI filter = new FeatureMatcherSet(); + filter.and(FeatureMatcher.byAttribute(Condition.Contains, "pathogenic", + "CLIN_SIG")); + fr.setFeatureFilter("METAL", filter); + assertTrue(fr.isVisible(sf)); // feature matches filter + 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); + } }