X-Git-Url: http://source.jalview.org/gitweb/?a=blobdiff_plain;f=src%2Fjalview%2Fschemes%2FBlosum62ColourScheme.java;h=02f9b3e3d9f55bf5159699a8da8eb145931a3ce6;hb=3d0101179759ef157b088ea135423cd909512d9f;hp=0f04168574b927a93093c03ea4408846f0cafdc0;hpb=a79f9e113c51c032070c670e45ce3eb464691166;p=jalview.git diff --git a/src/jalview/schemes/Blosum62ColourScheme.java b/src/jalview/schemes/Blosum62ColourScheme.java index 0f04168..02f9b3e 100755 --- a/src/jalview/schemes/Blosum62ColourScheme.java +++ b/src/jalview/schemes/Blosum62ColourScheme.java @@ -1,89 +1,119 @@ /* - * Jalview - A Sequence Alignment Editor and Viewer (Version 2.7) - * Copyright (C) 2011 J Procter, AM Waterhouse, J Engelhardt, LM Lui, G Barton, M Clamp, S Searle + * 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. - * + * 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 . + * 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.schemes; -import jalview.analysis.AAFrequency; +import jalview.analysis.scoremodels.ScoreModels; +import jalview.api.analysis.PairwiseScoreModelI; +import jalview.datamodel.AnnotatedCollectionI; +import jalview.datamodel.SequenceCollectionI; +import jalview.datamodel.SequenceI; +import jalview.util.Comparison; import java.awt.Color; -import jalview.datamodel.SequenceI; +import java.util.Map; + public class Blosum62ColourScheme extends ResidueColourScheme { + private static final Color LIGHT_BLUE = new Color(204, 204, 255); + + private static final Color DARK_BLUE = new Color(154, 154, 255); + public Blosum62ColourScheme() { 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 || j >= consensus.length || consensus[j] == null - || (threshold != 0 && !aboveThreshold(res, j))) + if (Comparison.isGap(res) || consensusResidue == null) { return Color.white; } - Color currentColour; + Color colour; - if (!jalview.util.Comparison.isGap(res)) + if (consensusResidue.indexOf(res) > -1) + { + colour = DARK_BLUE; + } + else { - String max = (String) consensus[j].get(AAFrequency.MAXRESIDUE); + float score = 0; - if (max.indexOf(res) > -1) + for (char consensus : consensusResidue.toCharArray()) { - currentColour = new Color(154, 154, 255); + score += sm.getPairwiseScore(consensus, res); } - else + + if (score > 0) { - 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 = new Color(204, 204, 255); - } - else - { - currentColour = Color.white; - } + colour = LIGHT_BLUE; } - - if (conservationColouring) + else { - currentColour = applyConservation(currentColour, j); + colour = Color.white; } } - else - { - return Color.white; - } + return colour; + } + + @Override + public boolean isPeptideSpecific() + { + return true; + } + + @Override + public String getSchemeName() + { + return JalviewColourScheme.Blosum62.toString(); + } - return currentColour; + @Override + public boolean isSimple() + { + return false; } }