X-Git-Url: http://source.jalview.org/gitweb/?a=blobdiff_plain;f=src%2Fjalview%2Fext%2Fensembl%2FEnsemblCdna.java;h=f60125b94b97bdef9237ed8a57e017a463680d52;hb=ef9282b464dc189faf9ce40a4b7420a204266668;hp=e4eb8730529226e7e7d339ce7fc31cd7f298faa9;hpb=f9b80711054b61e8c2257488a1637e15616cb9c9;p=jalview.git diff --git a/src/jalview/ext/ensembl/EnsemblCdna.java b/src/jalview/ext/ensembl/EnsemblCdna.java index e4eb873..f60125b 100644 --- a/src/jalview/ext/ensembl/EnsemblCdna.java +++ b/src/jalview/ext/ensembl/EnsemblCdna.java @@ -1,14 +1,21 @@ package jalview.ext.ensembl; import jalview.datamodel.SequenceFeature; -import jalview.io.gff.SequenceOntology; +import jalview.io.gff.SequenceOntologyFactory; +import jalview.io.gff.SequenceOntologyI; + +import java.util.List; import com.stevesoft.pat.Regex; public class EnsemblCdna extends EnsemblSeqProxy { + // TODO modify to accept other species e.g. ENSMUSPnnn + private static final Regex ACCESSION_REGEX = new Regex( + "(ENST|ENSG|CCDS)[0-9.]{3,}$"); + /* - * fetch exon features on genomic sequence (to identify the cdnaregions) + * fetch exon features on genomic sequence (to identify the cdna regions) * and cds and variation features (to retain) */ private static final EnsemblFeatureType[] FEATURES_TO_FETCH = { @@ -35,7 +42,7 @@ public class EnsemblCdna extends EnsemblSeqProxy @Override public Regex getAccessionValidator() { - return new Regex("((ENST|ENSG|CCDS)[0-9.]{3,})"); + return ACCESSION_REGEX; } @Override @@ -45,20 +52,17 @@ public class EnsemblCdna extends EnsemblSeqProxy } /** - * Answers true unless the feature type is 'exon' (or a sub-type of exon in - * the Sequence Ontology). Exon features are only retrieved in order to - * identify the exon sequence loci, and are redundant information on the exon - * sequence itself. + * Answers true unless the feature type is 'transcript' (or a sub-type in the + * Sequence Ontology). */ @Override protected boolean retainFeature(SequenceFeature sf, String accessionId) { - if (SequenceOntology.getInstance().isA(sf.getType(), - SequenceOntology.EXON)) + if (isTranscript(sf.getType())) { return false; } - return super.retainFeature(sf, accessionId); + return featureMayBelong(sf, accessionId); } /** @@ -69,8 +73,8 @@ public class EnsemblCdna extends EnsemblSeqProxy @Override protected boolean identifiesSequence(SequenceFeature sf, String accId) { - if (SequenceOntology.getInstance().isA(sf.getType(), - SequenceOntology.EXON)) + if (SequenceOntologyFactory.getInstance().isA(sf.getType(), + SequenceOntologyI.EXON)) { String parentFeature = (String) sf.getValue(PARENT); if (("transcript:" + accId).equals(parentFeature)) @@ -81,4 +85,12 @@ public class EnsemblCdna extends EnsemblSeqProxy return false; } + @Override + protected List getCrossReferenceDatabases() + { + return super.getCrossReferenceDatabases(); + // 30/01/16 also found Vega_transcript, OTTT, ENS_LRG_transcript, UCSC, + // HGNC_trans_name, RefSeq_mRNA, RefSeq_mRNA_predicted + } + }