X-Git-Url: http://source.jalview.org/gitweb/?a=blobdiff_plain;f=src%2Fjalview%2Fext%2Fensembl%2FEnsemblCds.java;h=ec5780fb05c92867cff5819c938a0ab8e1e58ab9;hb=4b3e12d9f3c2e2ad8a0317160881cef075690e97;hp=897371df7dab4f2bca0fbceb0d20da296c91b366;hpb=b03ec66ae6238b44bd20d2403d1157cadc5f0e01;p=jalview.git diff --git a/src/jalview/ext/ensembl/EnsemblCds.java b/src/jalview/ext/ensembl/EnsemblCds.java index 897371d..ec5780f 100644 --- a/src/jalview/ext/ensembl/EnsemblCds.java +++ b/src/jalview/ext/ensembl/EnsemblCds.java @@ -1,16 +1,21 @@ package jalview.ext.ensembl; import jalview.datamodel.SequenceFeature; -import jalview.io.gff.SequenceOntology; +import jalview.datamodel.SequenceI; +import jalview.io.gff.SequenceOntologyFactory; +import jalview.io.gff.SequenceOntologyI; + +import java.util.List; public class EnsemblCds extends EnsemblSeqProxy { /* * fetch cds features on genomic sequence (to identify the CDS regions) - * and variation features (to retain) + * and exon and variation features (to retain for display) */ private static final EnsemblFeatureType[] FEATURES_TO_FETCH = { - EnsemblFeatureType.cds, EnsemblFeatureType.variation }; + EnsemblFeatureType.cds, EnsemblFeatureType.exon, + EnsemblFeatureType.variation }; /** * Constructor @@ -45,9 +50,14 @@ public class EnsemblCds extends EnsemblSeqProxy * itself. */ @Override - protected boolean retainFeature(String type) + protected boolean retainFeature(SequenceFeature sf, String accessionId) { - return !SequenceOntology.getInstance().isA(type, SequenceOntology.CDS); + if (SequenceOntologyFactory.getInstance().isA(sf.getType(), + SequenceOntologyI.CDS)) + { + return false; + } + return featureMayBelong(sf, accessionId); } /** @@ -58,10 +68,10 @@ public class EnsemblCds extends EnsemblSeqProxy @Override protected boolean identifiesSequence(SequenceFeature sf, String accId) { - if (SequenceOntology.getInstance().isA(sf.getType(), - SequenceOntology.CDS)) + if (SequenceOntologyFactory.getInstance().isA(sf.getType(), + SequenceOntologyI.CDS)) { - String parentFeature = (String) sf.getValue("Parent"); + String parentFeature = (String) sf.getValue(PARENT); if (("transcript:" + accId).equals(parentFeature)) { return true; @@ -70,4 +80,18 @@ public class EnsemblCds extends EnsemblSeqProxy return false; } + /** + * Overrides this method to trivially return a range which is the whole of the + * nucleotide sequence. This is both faster than scanning for CDS features, + * 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) + { + int len = dnaSeq.getLength(); + ranges.add(new int[] { 1, len }); + return len; + } + }