+ /**
+ * Tests retrieval of one entry from EMBLCDS. Test is dependent on
+ * availability of network and the EMBLCDS service.
+ *
+ * @throws Exception
+ */
+ @Test(groups = { "External" })
+ public void testEmblCDSUniprotProductRecovery() throws Exception
+ {
+ String retrievalId = "AAH29712";
+ DbSourceProxy embl = new SequenceFetcher().getSourceProxy(
+ DBRefSource.EMBLCDS).get(0);
+ assertNotNull("Couldn't find the EMBL retrieval client", embl);
+ verifyProteinNucleotideXref(retrievalId, embl);
+ }
+
+ /**
+ * Helper method to perform database retrieval and verification of results.
+ *
+ * @param retrievalId
+ * @param embl
+ * @throws Exception
+ */
+ private void verifyProteinNucleotideXref(String retrievalId,
+ DbSourceProxy embl) throws Exception
+ {
+ AlignmentI alsq = embl.getSequenceRecords(retrievalId);
+ assertNotNull("Couldn't find the EMBL record " + retrievalId, alsq);
+ assertEquals("Didn't retrieve right number of records", 1,
+ alsq.getHeight());
+ SequenceI seq = alsq.getSequenceAt(0);
+ assertEquals("Wrong sequence name", embl.getDbSource() + "|"
+ + retrievalId, seq.getName());
+ SequenceFeature[] sfs = seq.getSequenceFeatures();
+ assertNotNull("Sequence features missing", sfs);
+ assertTrue(
+ "Feature not CDS",
+ FeatureProperties.isCodingFeature(embl.getDbSource(),
+ sfs[0].getType()));
+ assertEquals(embl.getDbSource(), sfs[0].getFeatureGroup());
+ DBRefEntry[] dr = DBRefUtils.selectRefs(seq.getDBRefs(),
+ new String[] { DBRefSource.UNIPROT });
+ assertNotNull(dr);
+ assertEquals("Expected a single Uniprot cross reference", 1, dr.length);
+ assertEquals("Expected cross reference map to be one amino acid", dr[0]
+ .getMap().getMappedWidth(), 1);
+ assertEquals("Expected local reference map to be 3 nucleotides", dr[0]
+ .getMap().getWidth(), 3);
+ AlignmentI sprods = new CrossRef(alsq.getSequencesArray(), alsq)
+ .findXrefSequences(dr[0].getSource(), true);
+ assertNotNull(
+ "Couldn't recover cross reference sequence from dataset. Was it ever added ?",
+ sprods);
+ assertEquals("Didn't xref right number of records", 1,
+ sprods.getHeight());
+ SequenceI proteinSeq = sprods.getSequenceAt(0);
+ assertEquals(proteinSeq.getSequenceAsString(), dr[0].getMap().getTo()
+ .getSequenceAsString());
+ assertEquals(dr[0].getSource() + "|" + dr[0].getAccessionId(),
+ proteinSeq.getName());
+ }