X-Git-Url: http://source.jalview.org/gitweb/?a=blobdiff_plain;f=test%2Fjalview%2Fext%2Fensembl%2FEnsemblCdnaTest.java;h=9e9d9a422673456b7cb7e4bfe35c7810ace63501;hb=3c8a25936a2d805e7e3d7ab82f83b13135406d18;hp=6611e057bc5064938723fdee5678bb97cd0e436b;hpb=136c0793b90b72b928c4d77dc109dd5c644e00d3;p=jalview.git diff --git a/test/jalview/ext/ensembl/EnsemblCdnaTest.java b/test/jalview/ext/ensembl/EnsemblCdnaTest.java index 6611e05..9e9d9a4 100644 --- a/test/jalview/ext/ensembl/EnsemblCdnaTest.java +++ b/test/jalview/ext/ensembl/EnsemblCdnaTest.java @@ -25,6 +25,7 @@ import static org.testng.AssertJUnit.assertFalse; import static org.testng.AssertJUnit.assertNull; import static org.testng.AssertJUnit.assertTrue; +import jalview.datamodel.Sequence; import jalview.datamodel.SequenceDummy; import jalview.datamodel.SequenceFeature; import jalview.datamodel.SequenceI; @@ -78,19 +79,19 @@ public class EnsemblCdnaTest // exon at (start+10000) length 501 SequenceFeature sf = new SequenceFeature("exon", "", 20000, 20500, 0f, null); - sf.setValue("Parent", "transcript:" + transcriptId); + sf.setValue("Parent", transcriptId); sf.setStrand("-"); genomic.addSequenceFeature(sf); // exon (sub-type) at (start + exon_variant) length 101 sf = new SequenceFeature("coding_exon", "", 10500, 10600, 0f, null); - sf.setValue("Parent", "transcript:" + transcriptId); + sf.setValue("Parent", transcriptId); sf.setStrand("-"); genomic.addSequenceFeature(sf); // exon belonging to a different transcript doesn't count sf = new SequenceFeature("exon", "", 11500, 12600, 0f, null); - sf.setValue("Parent", "transcript:anotherOne"); + sf.setValue("Parent", "anotherOne"); genomic.addSequenceFeature(sf); // transcript feature doesn't count @@ -133,19 +134,19 @@ public class EnsemblCdnaTest // exon at (start+10000) length 501 SequenceFeature sf = new SequenceFeature("exon", "", 20000, 20500, 0f, null); - sf.setValue("Parent", "transcript:" + transcriptId); + sf.setValue("Parent", transcriptId); sf.setStrand("+"); genomic.addSequenceFeature(sf); // exon (sub-type) at (start + exon_variant) length 101 sf = new SequenceFeature("coding_exon", "", 10500, 10600, 0f, null); - sf.setValue("Parent", "transcript:" + transcriptId); + sf.setValue("Parent", transcriptId); sf.setStrand("+"); genomic.addSequenceFeature(sf); // exon belonging to a different transcript doesn't count sf = new SequenceFeature("exon", "", 11500, 12600, 0f, null); - sf.setValue("Parent", "transcript:anotherOne"); + sf.setValue("Parent", "anotherOne"); genomic.addSequenceFeature(sf); // transcript feature doesn't count @@ -225,11 +226,14 @@ public class EnsemblCdnaTest assertTrue(testee.retainFeature(sf, accId)); // other feature with desired parent is retained - sf.setValue("Parent", "transcript:" + accId); + sf.setValue("Parent", accId); assertTrue(testee.retainFeature(sf, accId)); + // test is not case-sensitive + assertTrue(testee.retainFeature(sf, accId.toLowerCase())); + // feature with wrong parent is not retained - sf.setValue("Parent", "transcript:XYZ"); + sf.setValue("Parent", "XYZ"); assertFalse(testee.retainFeature(sf, accId)); } @@ -238,37 +242,51 @@ public class EnsemblCdnaTest * accession id as parent */ @Test(groups = "Functional") - public void testIdentifiesSequence() + public void testGetIdentifyingFeatures() { String accId = "ABC123"; - EnsemblCdna testee = new EnsemblCdna(); + SequenceI seq = new Sequence(accId, "MKLNFRQIE"); - // exon with no parent not valid - SequenceFeature sf = new SequenceFeature("exon", "", 1, 2, 0f, null); - assertFalse(testee.identifiesSequence(sf, accId)); + // exon with no parent: not valid + SequenceFeature sf1 = new SequenceFeature("exon", "", 1, 2, 0f, null); + seq.addSequenceFeature(sf1); - // exon with wrong parent not valid - sf.setValue("Parent", "transcript:XYZ"); - assertFalse(testee.identifiesSequence(sf, accId)); + // exon with wrong parent: not valid + SequenceFeature sf2 = new SequenceFeature("exon", "", 1, 2, 0f, null); + sf2.setValue("Parent", "XYZ"); + seq.addSequenceFeature(sf2); // exon with right parent is valid - sf.setValue("Parent", "transcript:" + accId); - assertTrue(testee.identifiesSequence(sf, accId)); + SequenceFeature sf3 = new SequenceFeature("exon", "", 1, 2, 0f, null); + sf3.setValue("Parent", accId); + seq.addSequenceFeature(sf3); // exon sub-type with right parent is valid - sf = new SequenceFeature("coding_exon", "", 1, 2, 0f, null); - sf.setValue("Parent", "transcript:" + accId); - assertTrue(testee.identifiesSequence(sf, accId)); + SequenceFeature sf4 = new SequenceFeature("coding_exon", "", 1, 2, 0f, + null); + sf4.setValue("Parent", accId); + seq.addSequenceFeature(sf4); // transcript not valid: - sf = new SequenceFeature("transcript", "", 1, 2, 0f, null); - sf.setValue("Parent", "transcript:" + accId); - assertFalse(testee.identifiesSequence(sf, accId)); + SequenceFeature sf5 = new SequenceFeature("transcript", "", 1, 2, 0f, + null); + sf5.setValue("Parent", accId); + seq.addSequenceFeature(sf5); // CDS not valid: - sf = new SequenceFeature("CDS", "", 1, 2, 0f, null); - sf.setValue("Parent", "transcript:" + accId); - assertFalse(testee.identifiesSequence(sf, accId)); + SequenceFeature sf6 = new SequenceFeature("transcript", "", 1, 2, 0f, + null); + sf6.setValue("Parent", accId); + seq.addSequenceFeature(sf6); + + List sfs = new EnsemblCdna() + .getIdentifyingFeatures(seq, accId); + assertFalse(sfs.contains(sf1)); + assertFalse(sfs.contains(sf2)); + assertTrue(sfs.contains(sf3)); + assertTrue(sfs.contains(sf4)); + assertFalse(sfs.contains(sf5)); + assertFalse(sfs.contains(sf6)); } @Test(groups = "Functional")