X-Git-Url: http://source.jalview.org/gitweb/?a=blobdiff_plain;f=src%2Fjalview%2Fschemes%2FBlosum62ColourScheme.java;h=70f49109249aba5b43c54b9f8b7810dde2844973;hb=136c0793b90b72b928c4d77dc109dd5c644e00d3;hp=c47f17180bc2311648e4cf13693be66c197d6ca7;hpb=8677e6e34e291edc58c1da2fc9c958473754143f;p=jalview.git diff --git a/src/jalview/schemes/Blosum62ColourScheme.java b/src/jalview/schemes/Blosum62ColourScheme.java index c47f171..70f4910 100755 --- a/src/jalview/schemes/Blosum62ColourScheme.java +++ b/src/jalview/schemes/Blosum62ColourScheme.java @@ -20,6 +20,8 @@ */ package jalview.schemes; +import jalview.analysis.scoremodels.ScoreModels; +import jalview.api.analysis.PairwiseScoreModelI; import jalview.datamodel.AnnotatedCollectionI; import jalview.datamodel.SequenceCollectionI; import jalview.datamodel.SequenceI; @@ -38,73 +40,79 @@ public class Blosum62ColourScheme extends ResidueColourScheme super(); } + /** + * Returns a new instance of this colour scheme with which the given data may + * be coloured + */ @Override - public Color findColour(char res, int j, SequenceI seq) + public ColourSchemeI getInstance(AnnotatedCollectionI coll, + Map hrs) { + return new Blosum62ColourScheme(); + } + + @Override + public Color findColour(char res, int j, SequenceI seq, + String consensusResidue, float pid) + { + PairwiseScoreModelI sm = ScoreModels.getInstance().getBlosum62(); + + /* + * compare as upper case; note consensusResidue is + * always computed as uppercase + */ if ('a' <= res && res <= 'z') { - // TO UPPERCASE !!! res -= ('a' - 'A'); } - if (consensus == null || consensus.get(j) == null - || (threshold != 0 && !aboveThreshold(res, j))) + if (Comparison.isGap(res) || consensusResidue == null) { return Color.white; } - Color currentColour; + Color colour; - if (!Comparison.isGap(res)) + if (consensusResidue.indexOf(res) > -1) { - /* - * test if this is the consensus (or joint consensus) residue - */ - String max = consensus.get(j).getModalResidue(); + colour = DARK_BLUE; + } + else + { + float score = 0; - if (max.indexOf(res) > -1) + for (char consensus : consensusResidue.toCharArray()) { - currentColour = DARK_BLUE; + score += sm.getPairwiseScore(consensus, res); } - else - { - int c = 0; - int max_aa = 0; - int n = max.length(); - do - { - c += ResidueProperties.getBLOSUM62(max.charAt(max_aa), res); - } while (++max_aa < n); - - if (c > 0) - { - currentColour = LIGHT_BLUE; - } - else - { - currentColour = Color.white; - } + if (score > 0) + { + colour = LIGHT_BLUE; } - - if (conservationColouring) + else { - currentColour = applyConservation(currentColour, j); + colour = Color.white; } } - else - { - return Color.white; - } + return colour; + } - return currentColour; + @Override + public boolean isPeptideSpecific() + { + return true; + } + + @Override + public String getSchemeName() + { + return JalviewColourScheme.Blosum62.toString(); } @Override - public ColourSchemeI applyTo(AnnotatedCollectionI sg, - Map hiddenRepSequences) + public boolean isSimple() { - ColourSchemeI newcs = super.applyTo(sg, hiddenRepSequences); - return newcs; + return false; } }