X-Git-Url: http://source.jalview.org/gitweb/?p=jalview.git;a=blobdiff_plain;f=test%2Fjalview%2Fio%2FSequenceAnnotationReportTest.java;h=bb5d1c4f1a47716ac71ec292a466cdc721917aac;hp=7e00caac8566cb3658feb09bdafbb91994050679;hb=5d003f492184b893aa4c2f437c6021626b573c45;hpb=5a631296dd1dcc1df7b50487a647c27333696c74 diff --git a/test/jalview/io/SequenceAnnotationReportTest.java b/test/jalview/io/SequenceAnnotationReportTest.java index 7e00caa..bb5d1c4 100644 --- a/test/jalview/io/SequenceAnnotationReportTest.java +++ b/test/jalview/io/SequenceAnnotationReportTest.java @@ -33,6 +33,8 @@ import org.testng.annotations.Test; import jalview.api.FeatureColourI; import jalview.datamodel.DBRefEntry; +import jalview.datamodel.MappedFeatures; +import jalview.datamodel.Mapping; import jalview.datamodel.Sequence; import jalview.datamodel.SequenceFeature; import jalview.datamodel.SequenceI; @@ -40,6 +42,7 @@ import jalview.gui.JvOptionPane; import jalview.io.gff.GffConstants; import jalview.renderer.seqfeatures.FeatureRenderer; import jalview.schemes.FeatureColour; +import jalview.util.MapList; import jalview.viewmodel.seqfeatures.FeatureRendererModel; import junit.extensions.PA; @@ -59,8 +62,8 @@ public class SequenceAnnotationReportTest SequenceAnnotationReport sar = new SequenceAnnotationReport(false); StringBuilder sb = new StringBuilder(); sb.append("123456"); - SequenceFeature sf = new SequenceFeature("disulfide bond", "desc", 1, - 3, 1.2f, "group"); + SequenceFeature sf = new SequenceFeature("disulfide bond", "desc", 1, 3, + 1.2f, "group"); // residuePos == 2 does not match start or end of feature, nothing done: sar.appendFeature(sb, 2, null, sf, null, 0); @@ -102,8 +105,12 @@ public class SequenceAnnotationReportTest sfl.add(sf); sfl.add(sf); sfl.add(sf); - int n = sar.appendFeatures(sb, 1, sfl, - new FeatureRenderer(null), 200); // text should terminate before 200 characters + int n = sar.appendFeatures(sb, 1, sfl, new FeatureRenderer(null), 200); // text + // should + // terminate + // before + // 200 + // characters String s = sb.toString(); assertTrue(s.length() < 200); assertEquals(n, 7); // should be 7 features left over @@ -284,19 +291,20 @@ public class SequenceAnnotationReportTest /* * positional features are ignored */ - seq.addSequenceFeature(new SequenceFeature("Domain", "Ferredoxin", 5, - 10, 1f, null)); + seq.addSequenceFeature( + new SequenceFeature("Domain", "Ferredoxin", 5, 10, 1f, null)); sar.createSequenceAnnotationReport(sb, seq, true, true, null); - assertEquals("SeqDesc", sb.toString()); + assertEquals("SeqDesc\n" + "\n" + "", sb.toString()); /* * non-positional feature */ - seq.addSequenceFeature(new SequenceFeature("Type1", "Nonpos", 0, 0, 1f, - null)); + seq.addSequenceFeature( + new SequenceFeature("Type1", "Nonpos", 0, 0, 1f, null)); sb.setLength(0); sar.createSequenceAnnotationReport(sb, seq, true, true, null); - String expected = "SeqDesc
Type1 ; Nonpos Score=1.0
"; + String expected = "SeqDesc\n" + "\n" + + "
Type1 ; Nonpos Score=1.0
"; assertEquals(expected, sb.toString()); /* @@ -304,7 +312,7 @@ public class SequenceAnnotationReportTest */ sb.setLength(0); sar.createSequenceAnnotationReport(sb, seq, true, false, null); - assertEquals("SeqDesc", sb.toString()); + assertEquals("SeqDesc\n\n", sb.toString()); /* * add non-pos feature with score inside min-max range for feature type @@ -312,8 +320,8 @@ public class SequenceAnnotationReportTest * score is only appended for positional features so ignored here! * minMax are not recorded for non-positional features */ - seq.addSequenceFeature(new SequenceFeature("Metal", "Desc", 0, 0, 5f, - null)); + seq.addSequenceFeature( + new SequenceFeature("Metal", "Desc", 0, 0, 5f, null)); FeatureRendererModel fr = new FeatureRenderer(null); Map minmax = fr.getMinMax(); @@ -321,9 +329,10 @@ public class SequenceAnnotationReportTest sb.setLength(0); sar.createSequenceAnnotationReport(sb, seq, true, true, fr); - expected = "SeqDesc
Metal ; Desc
Type1 ; Nonpos
"; + expected = "SeqDesc\n" + "\n" + + "
Metal ; Desc
Type1 ; Nonpos
"; assertEquals(expected, sb.toString()); - + /* * 'linkonly' features are ignored; this is obsolete, as linkonly * is only set by DasSequenceFetcher, and DAS is history @@ -340,13 +349,14 @@ public class SequenceAnnotationReportTest * 'clinical_significance' attribute is only included in description * when used for feature colouring */ - SequenceFeature sf2 = new SequenceFeature("Variant", "Havana", 0, 0, - 5f, null); + SequenceFeature sf2 = new SequenceFeature("Variant", "Havana", 0, 0, 5f, + null); sf2.setValue(GffConstants.CLINICAL_SIGNIFICANCE, "benign"); seq.addSequenceFeature(sf2); sb.setLength(0); sar.createSequenceAnnotationReport(sb, seq, true, true, fr); - expected = "SeqDesc
Metal ; Desc
Type1 ; Nonpos
Variant ; Havana
"; + expected = "SeqDesc\n" + "\n" + + "
Metal ; Desc
Type1 ; Nonpos
Variant ; Havana
"; assertEquals(expected, sb.toString()); /* @@ -367,13 +377,15 @@ public class SequenceAnnotationReportTest fc.setAttributeName("clinical_significance"); fr.setColour("Variant", fc); sar.createSequenceAnnotationReport(sb, seq, true, true, fr); - expected = "SeqDesc
UNIPROT P30419
PDB 3iu1
Metal ; Desc
" - + "Type1 ; Nonpos
Variant ; Havana; clinical_significance=benign
"; + expected = "SeqDesc\n" + "
\n" + "UNIPROT P30419
\n" + + "PDB 3iu1\n" + + "
Metal ; Desc
Type1 ; Nonpos
Variant ; Havana; clinical_significance=benign
"; assertEquals(expected, sb.toString()); // with showNonPositionalFeatures = false sb.setLength(0); sar.createSequenceAnnotationReport(sb, seq, true, false, fr); - expected = "SeqDesc
UNIPROT P30419
PDB 3iu1
"; + expected = "SeqDesc\n" + "
\n" + "UNIPROT P30419
\n" + + "PDB 3iu1\n" + "
"; assertEquals(expected, sb.toString()); /* @@ -383,7 +395,8 @@ public class SequenceAnnotationReportTest sf2.setDescription( "This is a very long description which should be truncated"); sar.createSequenceAnnotationReport(sb, seq, false, true, fr); - expected = "SeqDesc
Metal ; Desc
Type1 ; Nonpos
Variant ; This is a very long description which sh...; clinical_significance=benign
"; + expected = "SeqDesc\n" + "\n" + + "
Metal ; Desc
Type1 ; Nonpos
Variant ; This is a very long description which sh...; clinical_significance=benign
"; assertEquals(expected, sb.toString()); // see other tests for treatment of status and html @@ -399,7 +412,7 @@ public class SequenceAnnotationReportTest { SequenceAnnotationReport sar = new SequenceAnnotationReport(false); StringBuilder sb = new StringBuilder(); - + SequenceI seq = new Sequence("s1", "ABC"); int maxSources = (int) PA.getValue(sar, "MAX_SOURCES"); @@ -407,20 +420,80 @@ public class SequenceAnnotationReportTest { seq.addDBRef(new DBRefEntry("PDB" + i, "0", "3iu1")); } - + int maxRefs = (int) PA.getValue(sar, "MAX_REFS_PER_SOURCE"); for (int i = 0; i <= maxRefs; i++) { seq.addDBRef(new DBRefEntry("Uniprot", "0", "P3041" + i)); } - + sar.createSequenceAnnotationReport(sb, seq, true, true, null, true); String report = sb.toString(); - assertTrue(report - .startsWith( - "
UNIPROT P30410, P30411, P30412, P30413,...
PDB0 3iu1")); - assertTrue(report - .endsWith( - "
PDB7 3iu1
PDB8,...
(Output Sequence Details to list all database references)
")); + assertTrue(report.startsWith("\n" + "
\n" + "UNIPROT P30410,\n" + + " P30411,\n" + " P30412,\n" + " P30413,...
\n" + + "PDB0 3iu1
\n" + "PDB1 3iu1
")); + assertTrue(report.endsWith("PDB5 3iu1
\n" + "PDB6 3iu1
\n" + + "PDB7 3iu1
\n" + "PDB8,...
\n" + + "(Output Sequence Details to list all database references)\n" + + "
")); + } + + /** + * Test adding a linked feature to the tooltip + */ + @Test(groups = "Functional") + public void testAppendFeature_virtualFeature() + { + /* + * map CDS to peptide sequence + */ + SequenceI cds = new Sequence("Cds/101-121", "CCTttgAGAtttCAAatgGAT"); + SequenceI peptide = new Sequence("Peptide/8-14", "PLRFQMD"); + MapList map = new MapList(new int[] { 101, 118 }, new int[] { 8, 13 }, + 3, 1); + Mapping mapping = new Mapping(peptide, map); + + /* + * assume variant feature found at CDS position 106 G>C + */ + List features = new ArrayList<>(); + // vary ttg (Leu) to ttc (Phe) + SequenceFeature sf = new SequenceFeature("variant", "G,C", 106, 106, + Float.NaN, null); + features.add(sf); + MappedFeatures mf = new MappedFeatures(mapping, cds, 9, 'L', features); + + StringBuilder sb = new StringBuilder(); + SequenceAnnotationReport sar = new SequenceAnnotationReport(false); + sar.appendFeature(sb, 1, null, sf, mf, 0); + + /* + * linked feature shown in tooltip in protein coordinates + */ + assertEquals("variant 9; G,C", sb.toString()); + + /* + * adding "alleles" attribute to variant allows peptide consequence + * to be calculated and added to the tooltip + */ + sf.setValue("alleles", "G,C"); + sb = new StringBuilder(); + sar.appendFeature(sb, 1, null, sf, mf, 0); + assertEquals("variant 9; G,C p.Leu9Phe", sb.toString()); + + /* + * now a virtual peptide feature on CDS + * feature at 11-12 on peptide maps to 110-115 on CDS + * here we test for tooltip at 113 (t) + */ + SequenceFeature sf2 = new SequenceFeature("metal", "Fe", 11, 12, 2.3f, + "Uniprot"); + features.clear(); + features.add(sf2); + mapping = new Mapping(peptide, map); + mf = new MappedFeatures(mapping, peptide, 113, 't', features); + sb = new StringBuilder(); + sar.appendFeature(sb, 1, null, sf2, mf, 0); + assertEquals("metal 110 115; Fe Score=2.3", sb.toString()); } }