annotation which can then be used to sort the alignment via the Sort
by→Score menu.</em> <br></li>
<li><strong>Translate as cDNA</strong> (not applet)<br><em>This option is visible for nucleotide alignments.
- Selecting this option shows the DNA's calculated protein product in a new window. Note that the
+ Selecting this option shows the DNA's calculated protein product in a new <a href="../features/splitView.html">split frame</a> window. Note that the
translation is not frame- or intron-aware; it simply translates all codons in each sequence, using the
standard <a href="../misc/geneticCode.html">genetic code</a> (any incomplete final codon is discarded).
You can perform this action on the whole alignment,
or selected rows, columns, or regions.</em> <br></li>
<li><strong>Get Cross-References</strong> (not applet)<br><em>This option is visible where sequences have cross-references to
other standard databases; for example, an EMBL entry may have cross-references to one or more UNIPROT entries.
- Select the database to view all cross-referenced sequences in a new window.</em> <br></li>
+ Select the database to view all cross-referenced sequences in a new <a href="../features/splitView.html">split frame</a> window.</em> <br></li>
<li><strong>Autocalculate Consensus</strong><br> <em>For
large alignments it can be useful to deselect "Autocalculate
Consensus" when editing. This prevents the sometimes lengthy
import jalview.datamodel.AlignmentI;
import jalview.datamodel.DBRefEntry;
import jalview.datamodel.DBRefSource;
+import jalview.datamodel.FeatureProperties;
+import jalview.datamodel.SequenceFeature;
+import jalview.datamodel.SequenceI;
import jalview.util.DBRefUtils;
import jalview.ws.SequenceFetcher;
{
}
+ /**
+ * Tests that standard protein database sources include Uniprot (as the first)
+ * and also PDB. (Additional sources are dependent on available of DAS
+ * services.)
+ */
@Test(groups ={ "Functional" })
public void testStandardProtDbs()
{
srces.addAll(srcesfordb);
}
}
- DbSourceProxy uniprot = null;
+
int i = 0;
+ int uniprotPos = -1;
+ int pdbPos = -1;
// append the selected sequence sources to the default dbs
for (DbSourceProxy s : srces)
{
- if (s.getDbSource().equalsIgnoreCase(DBRefSource.UNIPROT))
+ if (s instanceof jalview.ws.dbsources.Uniprot)
{
- i++;
+ uniprotPos = i;
}
-
- if (s instanceof jalview.ws.dbsources.Uniprot)
+ if (s instanceof jalview.ws.dbsources.Pdb)
{
- uniprot = s;
- break;
+ pdbPos = i;
}
+ i++;
}
assertTrue("Failed to find Uniprot source as first source amongst "
- + srces.size() + " sources (source was at position " + i + ")",
- uniprot != null && i < 2);
+ + srces.size() + " sources (source was at position "
+ + uniprotPos + ")", uniprotPos == 0);
+ assertTrue("Failed to find PDB source amongst " + srces.size()
+ + " sources", pdbPos >= 0);
}
+ /**
+ * Tests retrieval of one entry from EMBL. Test is dependent on availability
+ * of network and the EMBL service.
+ *
+ * @throws Exception
+ */
@Test(groups =
{ "External" })
public void testEmblUniprotProductRecovery() throws Exception
{
- String retrievalId = "CAA23748"; // "V00488";
+ String retrievalId = "V00488";
DbSourceProxy embl = new SequenceFetcher().getSourceProxy(DBRefSource.EMBL).get(0);
assertNotNull("Couldn't find the EMBL retrieval client", embl);
verifyProteinNucleotideXref(retrievalId, embl);
}
+ /**
+ * 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
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());
- DBRefEntry[] dr = DBRefUtils.selectRefs(alsq.getSequenceAt(0).getDBRef(), DBRefSource.PROTEINSEQ);
+ 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.getDBRef(), DBRefSource.PROTEINSEQ);
assertNotNull(dr);
assertEquals("Expected a single Uniprot cross reference", 1, dr.length);
- assertEquals("Expected cross refernce map to be one amino acid", dr[0]
+ assertEquals("Expected cross reference map to be one amino acid", dr[0]
.getMap().getMappedWidth(), 1);
- assertEquals("Expected local refernce map to be 3 nucleotides", dr[0]
+ assertEquals("Expected local reference map to be 3 nucleotides", dr[0]
.getMap().getWidth(), 3);
AlignmentI sprods = CrossRef.findXrefSequences(alsq.getSequencesArray(), true, dr[0].getSource(), alsq.getDataset());
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());
}
}