X-Git-Url: http://source.jalview.org/gitweb/?a=blobdiff_plain;f=src%2Fjalview%2Fext%2Fensembl%2FEnsemblCds.java;h=8f13d99df94f71654fb11655176bc43dd2ddcfed;hb=3c8a25936a2d805e7e3d7ab82f83b13135406d18;hp=22c0a067013615c027a6c6faf6194c76a52ad3cc;hpb=56b5f4d5ca50971a34c9284bbb4b0507f7ba8a71;p=jalview.git diff --git a/src/jalview/ext/ensembl/EnsemblCds.java b/src/jalview/ext/ensembl/EnsemblCds.java index 22c0a06..8f13d99 100644 --- a/src/jalview/ext/ensembl/EnsemblCds.java +++ b/src/jalview/ext/ensembl/EnsemblCds.java @@ -1,11 +1,43 @@ +/* + * Jalview - A Sequence Alignment Editor and Viewer ($$Version-Rel$$) + * Copyright (C) $$Year-Rel$$ The Jalview Authors + * + * This file is part of Jalview. + * + * Jalview is free software: you can redistribute it and/or + * modify it under the terms of the GNU General Public License + * as published by the Free Software Foundation, either version 3 + * of the License, or (at your option) any later version. + * + * Jalview is distributed in the hope that it will be useful, but + * WITHOUT ANY WARRANTY; without even the implied warranty + * of MERCHANTABILITY or FITNESS FOR A PARTICULAR + * PURPOSE. See the GNU General Public License for more details. + * + * You should have received a copy of the GNU General Public License + * along with Jalview. If not, see . + * The Jalview Authors are detailed in the 'AUTHORS' file. + */ package jalview.ext.ensembl; import jalview.datamodel.SequenceFeature; import jalview.datamodel.SequenceI; -import jalview.io.gff.SequenceOntology; +import jalview.io.gff.SequenceOntologyFactory; +import jalview.io.gff.SequenceOntologyI; +import java.util.ArrayList; import java.util.List; +/** + * A client for direct fetching of CDS sequences from Ensembl (i.e. that part of + * the genomic sequence that is translated to protein) + * + * TODO: not currently used as CDS sequences are computed from CDS features on + * transcripts - delete this class? + * + * @author gmcarstairs + * + */ public class EnsemblCds extends EnsemblSeqProxy { /* @@ -17,13 +49,23 @@ public class EnsemblCds extends EnsemblSeqProxy EnsemblFeatureType.variation }; /** - * Constructor + * Default constructor (to use rest.ensembl.org) */ public EnsemblCds() { super(); } + /** + * Constructor given the target domain to fetch data from + * + * @param d + */ + public EnsemblCds(String d) + { + super(d); + } + @Override public String getDbName() { @@ -51,8 +93,8 @@ public class EnsemblCds extends EnsemblSeqProxy @Override protected boolean retainFeature(SequenceFeature sf, String accessionId) { - if (SequenceOntology.getInstance().isA(sf.getType(), - SequenceOntology.CDS)) + if (SequenceOntologyFactory.getInstance().isA(sf.getType(), + SequenceOntologyI.CDS)) { return false; } @@ -60,23 +102,26 @@ public class EnsemblCds extends EnsemblSeqProxy } /** - * Answers true if the sequence feature type is 'CDS' (or a subtype of CDS in - * the Sequence Ontology), and the Parent of the feature is the transcript we - * are retrieving + * Answers a list of sequence features (if any) whose type is 'CDS' (or a + * subtype of CDS in the Sequence Ontology), and whose Parent is the + * transcript we are retrieving */ @Override - protected boolean identifiesSequence(SequenceFeature sf, String accId) + protected List getIdentifyingFeatures(SequenceI seq, + String accId) { - if (SequenceOntology.getInstance().isA(sf.getType(), - SequenceOntology.CDS)) + List result = new ArrayList<>(); + List sfs = seq.getFeatures() + .getFeaturesByOntology(SequenceOntologyI.CDS); + for (SequenceFeature sf : sfs) { String parentFeature = (String) sf.getValue(PARENT); - if (("transcript:" + accId).equals(parentFeature)) + if ( accId.equals(parentFeature)) { - return true; + result.add(sf); } } - return false; + return result; } /** @@ -85,12 +130,12 @@ public class EnsemblCds extends EnsemblSeqProxy * and also means we don't need to keep CDS features on CDS sequence (where * they are redundant information). */ - @Override - protected int getCdsRanges(SequenceI dnaSeq, List ranges) + protected List getCdsRanges(SequenceI dnaSeq) { int len = dnaSeq.getLength(); + List ranges = new ArrayList<>(); ranges.add(new int[] { 1, len }); - return len; + return ranges; } }