X-Git-Url: http://source.jalview.org/gitweb/?a=blobdiff_plain;f=src%2Fjalview%2Fviewmodel%2Fseqfeatures%2FFeatureRendererModel.java;h=f0901908e2d41634ee1c8f6a71b85daaf695ac30;hb=7a0d503181fe41452120a8a02ca63476392aa08c;hp=821f7f47d0661d31a34604026ce285a16a965791;hpb=3740241b6d0dfa109dc81847afe58f17497c39e8;p=jalview.git diff --git a/src/jalview/viewmodel/seqfeatures/FeatureRendererModel.java b/src/jalview/viewmodel/seqfeatures/FeatureRendererModel.java index 821f7f4..f090190 100644 --- a/src/jalview/viewmodel/seqfeatures/FeatureRendererModel.java +++ b/src/jalview/viewmodel/seqfeatures/FeatureRendererModel.java @@ -25,6 +25,7 @@ import jalview.api.FeatureColourI; import jalview.api.FeaturesDisplayedI; import jalview.datamodel.AlignedCodonFrame; import jalview.datamodel.AlignmentI; +import jalview.datamodel.MappedFeatures; import jalview.datamodel.Mapping; import jalview.datamodel.SearchResultMatchI; import jalview.datamodel.SearchResults; @@ -1016,8 +1017,7 @@ public abstract class FeatureRendererModel * displayed, and feature group is null or the empty string * or marked for display */ - Set visibleFeatures = getFeaturesDisplayed() - .getVisibleFeatures(); + List visibleFeatures = getDisplayedFeatureTypes(); String[] visibleTypes = visibleFeatures .toArray(new String[visibleFeatures.size()]); List features = sequence.getFeatures().findFeatures( @@ -1156,18 +1156,27 @@ public abstract class FeatureRendererModel } /** - * Answers a (possibly empty) list of features in this alignment at a position - * (or range) which is mappable from the given sequence residue position in a - * mapped alignment. + * Answers a bean containing a mapping, and a list of features in this + * alignment at a position (or range) which is mappable from the given + * sequence residue position in a mapped alignment. Features are returned in + * render order of feature type (on top last), with order within feature type + * undefined. If no features or mapping are found, answers null. * * @param sequence * @param pos * @return */ - public List findComplementFeaturesAtResidue(SequenceI sequence, int pos) + public MappedFeatures findComplementFeaturesAtResidue(SequenceI sequence, + int pos) { SequenceI ds = sequence.getDatasetSequence(); - List result = new ArrayList<>(); + if (ds == null) + { + ds = sequence; + } + final char residue = ds.getCharAt(pos - ds.getStart()); + + List found = new ArrayList<>(); List mappings = this.av.getAlignment() .getCodonFrame(sequence); @@ -1175,9 +1184,12 @@ public abstract class FeatureRendererModel * todo: direct lookup of CDS for peptide and vice-versa; for now, * have to search through an unordered list of mappings for a candidate */ + Mapping mapping = null; + SequenceI mapFrom = null; + for (AlignedCodonFrame acf : mappings) { - Mapping mapping = acf.getMappingForSequence(sequence, true); + mapping = acf.getMappingForSequence(sequence, true); if (mapping == null || mapping.getMap().getFromRatio() == mapping .getMap().getToRatio()) { @@ -1189,19 +1201,52 @@ public abstract class FeatureRendererModel { int fromRes = match.getStart(); int toRes = match.getEnd(); + mapFrom = match.getSequence(); List fs = findFeaturesAtResidue( match.getSequence(), fromRes, toRes); for (SequenceFeature sf : fs) { - if (!result.contains(sf)) + if (!found.contains(sf)) { - result.add(sf); + found.add(sf); + } + } + } + + /* + * just take the first mapped features we find + */ + if (!found.isEmpty()) + { + break; + } + } + if (found.isEmpty()) + { + return null; + } + + /* + * sort by renderorder, inefficiently + */ + List result = new ArrayList<>(); + for (String type : renderOrder) + { + for (SequenceFeature sf : found) + { + if (type.equals(sf.getType())) + { + result.add(sf); + if (result.size() == found.size()) + { + return new MappedFeatures(mapping, mapFrom, pos, residue, + result); } } } } - return result; + return new MappedFeatures(mapping, mapFrom, pos, residue, result); } }