src/jalview/schemes/ClustalxColourScheme.java

   1 /* Jalview - a java multiple alignment editor\r
   2  * Copyright (C) 1998  Michele Clamp\r
   3  *\r
   4  * This program is free software; you can redistribute it and/or\r
   5  * modify it under the terms of the GNU General Public License\r
   6  * as published by the Free Software Foundation; either version 2\r
   7  * of the License, or (at your option) any later version.\r
   8  *\r
   9  * This program is distributed in the hope that it will be useful,\r
  10  * but WITHOUT ANY WARRANTY; without even the implied warranty of\r
  11  * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the\r
  12  * GNU General Public License for more details.\r
  13  *\r
  14  * You should have received a copy of the GNU General Public License\r
  15  * along with this program; if not, write to the Free Software\r
  16  * Foundation, Inc., 59 Temple Place - Suite 330, Boston, MA  02111-1307, USA.\r
  17  */\r
  18 \r
  19 package jalview.schemes;\r
  20 \r
  21 import jalview.datamodel.*;\r
  22 import java.util.*;\r
  23 import java.awt.*;\r
  24 \r
  25 public class ClustalxColourScheme extends ResidueColourScheme\r
  26 {\r
  27   Hashtable[] cons;\r
  28   int[][] cons2;\r
  29   ConsensusColour[] colours;\r
  30   ConsensusColour[] ResidueColour;\r
  31   int size;\r
  32     Consensus[] conses = new Consensus[32];\r
  33     public static Hashtable colhash = new Hashtable();\r
  34 \r
  35     {\r
  36         colhash.put("RED",new Color((float)0.9,(float)0.2,(float)0.1));\r
  37         colhash.put("BLUE",new Color((float)0.5,(float)0.7,(float)0.9));\r
  38         colhash.put("GREEN",new Color((float)0.1,(float)0.8,(float)0.1));\r
  39         colhash.put("ORANGE",new Color((float)0.9,(float)0.6,(float)0.3));\r
  40         colhash.put("CYAN",new Color((float)0.1,(float)0.7,(float)0.7));\r
  41         colhash.put("PINK",new Color((float)0.9,(float)0.5,(float)0.5));\r
  42         colhash.put("MAGENTA",new Color((float)0.8,(float)0.3,(float)0.8));\r
  43         colhash.put("YELLOW",new Color((float)0.8,(float)0.8,(float)0.0));\r
  44 \r
  45     }\r
  46   Vector colourTable = new Vector();\r
  47 \r
  48   public ClustalxColourScheme(Vector seqs, int maxWidth)\r
  49   {\r
  50     resetClustalX(seqs, maxWidth);\r
  51   }\r
  52 \r
  53   public void resetClustalX(Vector seqs, int maxWidth)\r
  54   {\r
  55     cons2 = new int[maxWidth][24];\r
  56     int start = 0;\r
  57 \r
  58     // Initialize the array\r
  59     for (int j=0;j<24;j++)\r
  60       for (int i=0; i < maxWidth;i++)\r
  61         cons2[i][j] = 0;\r
  62 \r
  63       int res,i, j = 0;\r
  64       String seq;\r
  65       while(j < seqs.size())\r
  66       {\r
  67         seq = ( (SequenceI) seqs.elementAt(j)).getSequence();\r
  68         int end_j = seq.length()-1;\r
  69         for (i = start; i <= end_j; i++)\r
  70         {\r
  71           if(seq.length()-1<i)\r
  72             res = 23;\r
  73           else\r
  74             res = ( (Integer) ResidueProperties.aaHash.get(seq.charAt(i) + "")).intValue();\r
  75           cons2[i][res]++;\r
  76         }\r
  77         j++;\r
  78       }\r
  79 \r
  80     this.size = seqs.size();\r
  81     makeColours();\r
  82   }\r
  83   public void makeColours() {\r
  84     conses[0] = new Consensus("WLVIMAFCYHP",60);\r
  85     conses[1] = new Consensus("WLVIMAFCYHP",80);\r
  86     conses[2] = new Consensus("ED",50);\r
  87     conses[3] = new Consensus("KR",60);\r
  88     conses[4] = new Consensus("G",50);\r
  89     conses[5] = new Consensus("N",50);\r
  90     conses[6] = new Consensus("QE",50);\r
  91     conses[7] = new Consensus("P",50);\r
  92     conses[8] = new Consensus("TS",50);\r
  93 \r
  94     conses[26] = new Consensus("A",85);\r
  95     conses[27] = new Consensus("C",85);\r
  96     conses[10] = new Consensus("E",85);\r
  97     conses[11] = new Consensus("F",85);\r
  98     conses[12] = new Consensus("G",85);\r
  99     conses[13] = new Consensus("H",85);\r
 100     conses[14] = new Consensus("I",85);\r
 101     conses[15] = new Consensus("L",85);\r
 102     conses[16] = new Consensus("M",85);\r
 103     conses[17] = new Consensus("N",85);\r
 104     conses[18] = new Consensus("P",85);\r
 105     conses[19] = new Consensus("Q",85);\r
 106     conses[20] = new Consensus("R",85);\r
 107     conses[21] = new Consensus("S",85);\r
 108     conses[22] = new Consensus("T",85);\r
 109     conses[23] = new Consensus("V",85);\r
 110     conses[24] = new Consensus("W",85);\r
 111     conses[25] = new Consensus("Y",85);\r
 112     conses[28] = new Consensus("K",85);\r
 113     conses[29] = new Consensus("D",85);\r
 114 \r
 115     conses[30] = new Consensus("G",0);\r
 116     conses[31] = new Consensus("P",0);\r
 117 \r
 118     // We now construct the colours\r
 119     colours = new ConsensusColour[11];\r
 120 \r
 121 \r
 122     Consensus[] tmp8 = new Consensus[1];\r
 123     tmp8[0] = conses[30];//G\r
 124     colours[7] =new ConsensusColour((Color)colhash.get("ORANGE"),tmp8);\r
 125 \r
 126     Consensus[] tmp9 = new Consensus[1];\r
 127     tmp9[0] = conses[31];//P\r
 128     colours[8] =new ConsensusColour((Color)colhash.get("YELLOW"),tmp9);\r
 129 \r
 130     Consensus[] tmp10 = new Consensus[1];\r
 131     tmp10[0] = conses[27];//C\r
 132     colours[9] =new ConsensusColour((Color)colhash.get("PINK"),tmp8);\r
 133 \r
 134     Consensus[] tmp1 = new Consensus[14];\r
 135     tmp1[0] = conses[0]; //%\r
 136     tmp1[1] = conses[1]; //#\r
 137     tmp1[2] = conses[26]; //A\r
 138     tmp1[3] = conses[27]; //C\r
 139     tmp1[4] = conses[11]; //F\r
 140     tmp1[5] = conses[13]; //H\r
 141     tmp1[6] = conses[14]; //I\r
 142     tmp1[7] = conses[15]; //L\r
 143     tmp1[8] = conses[16]; //M\r
 144     tmp1[9] = conses[23]; //V\r
 145     tmp1[10] = conses[24]; //W\r
 146     tmp1[11] = conses[25]; //Y\r
 147     tmp1[12] = conses[18]; //P\r
 148     tmp1[13] = conses[19]; //p\r
 149     colours[0] = new ConsensusColour((Color)colhash.get("BLUE"),tmp1);\r
 150 \r
 151     colours[10] = new ConsensusColour((Color)colhash.get("CYAN"),tmp1);\r
 152 \r
 153     Consensus[] tmp2 = new Consensus[5];\r
 154     tmp2[0] = conses[8];  //t\r
 155     tmp2[1] = conses[21]; //S\r
 156     tmp2[2] = conses[22]; //T\r
 157     tmp2[3] = conses[0];  //%\r
 158     tmp2[4] = conses[1];  //#\r
 159     colours[1] = new ConsensusColour((Color)colhash.get("GREEN"),tmp2);\r
 160 \r
 161     Consensus[] tmp3 = new Consensus[3];\r
 162 \r
 163     tmp3[0] = conses[17]; //N\r
 164     tmp3[1] = conses[29]; //D\r
 165     tmp3[2] = conses[5];  //n\r
 166     colours[2] = new ConsensusColour((Color)colhash.get("GREEN"),tmp3);\r
 167 \r
 168     Consensus[] tmp4 = new Consensus[6];\r
 169     tmp4[0] = conses[6];  // q = QE\r
 170     tmp4[1] = conses[19]; //Q\r
 171     tmp4[2] = conses[22]; //E\r
 172     tmp4[3] = conses[3];  //+\r
 173     tmp4[4] = conses[28]; //K\r
 174     tmp4[5] = conses[20]; //R\r
 175     colours[3] = new ConsensusColour((Color)colhash.get("GREEN"),tmp4);\r
 176     Consensus[] tmp5 = new Consensus[4];\r
 177     tmp5[0] = conses[3];  //+\r
 178     tmp5[1] = conses[28]; //K\r
 179     tmp5[2] = conses[20]; //R\r
 180     tmp5[3] = conses[19]; //Q\r
 181     colours[4] = new ConsensusColour((Color)colhash.get("RED"),tmp5);\r
 182     Consensus[] tmp6 = new Consensus[5];\r
 183     tmp6[0] = conses[3];  //-\r
 184     tmp6[1] = conses[29]; //D\r
 185     tmp6[2] = conses[10]; //E\r
 186     tmp6[3] = conses[6];  //q\r
 187     tmp6[4] = conses[19]; //Q\r
 188     colours[5] = new ConsensusColour((Color)colhash.get("MAGENTA"),tmp6);\r
 189 \r
 190     Consensus[] tmp7 = new Consensus[5];\r
 191     tmp7[0] = conses[3];  //-\r
 192     tmp7[1] = conses[29]; //D\r
 193     tmp7[2] = conses[10]; //E\r
 194     tmp7[3] = conses[17];  //N\r
 195     tmp7[4] = conses[2]; //DE\r
 196     colours[6] = new ConsensusColour((Color)colhash.get("MAGENTA"),tmp7);\r
 197 \r
 198     // Now attach the ConsensusColours to the residue letters\r
 199     ResidueColour = new ConsensusColour[20];\r
 200     ResidueColour[0] = colours[0]; // A\r
 201     ResidueColour[1] = colours[4]; // R\r
 202     ResidueColour[2] = colours[2]; // N\r
 203     ResidueColour[3] = colours[6]; // D\r
 204     ResidueColour[4] = colours[0]; // C\r
 205     ResidueColour[5] = colours[3]; // Q\r
 206     ResidueColour[6] = colours[5]; // E\r
 207     ResidueColour[7] = colours[7]; // G\r
 208     ResidueColour[8] = colours[10]; // H\r
 209     ResidueColour[9] = colours[0]; // I\r
 210     ResidueColour[10] = colours[0]; // L\r
 211     ResidueColour[11] = colours[4]; // K\r
 212     ResidueColour[12] = colours[0]; // M\r
 213     ResidueColour[13] = colours[0]; // F\r
 214     ResidueColour[14] = colours[8]; // P\r
 215     ResidueColour[15] = colours[1]; // S\r
 216     ResidueColour[16] = colours[1]; // T\r
 217     ResidueColour[17] = colours[0]; // W\r
 218     ResidueColour[18] = colours[10]; // Y\r
 219     ResidueColour[19] = colours[0]; // V\r
 220   }\r
 221 \r
 222   public Color findColour(String s)\r
 223   {\r
 224     return Color.pink;\r
 225   }\r
 226 \r
 227   public Color findColour(String s, int j)\r
 228    {\r
 229      if( threshold!=0 && !aboveThreshold(s,j))\r
 230       return Color.white;\r
 231 \r
 232 \r
 233      int i = ((Integer)ResidueProperties.aaHash.get(s)).intValue();\r
 234 \r
 235      Color c = Color.white;\r
 236 \r
 237      if(i>19)\r
 238        return c;\r
 239 \r
 240      for (int k=0; k < ResidueColour[i].conses.length ; k++)\r
 241        if (ResidueColour[i].conses[k].isConserved(cons2,j,size))\r
 242          c =  ResidueColour[i].c;\r
 243 \r
 244      if (i == 4)\r
 245     {\r
 246       if (conses[27].isConserved(cons2, j, size))\r
 247         c = (Color) colhash.get("PINK");\r
 248     }\r
 249 \r
 250 \r
 251      return c;\r
 252 \r
 253    }\r
 254 }\r
 255 \r
 256 class ConsensusColour {\r
 257 \r
 258   Consensus[] conses;\r
 259   Color c;\r
 260 \r
 261   public ConsensusColour(Color c,Consensus[] conses) {\r
 262     this.conses = conses;\r
 263     //    this.list = list;\r
 264     this.c = c;\r
 265   }\r
 266 }\r
 267 \r
 268 \r
 269 \r
 270 \r
 271 \r
 272 \r
 273 \r
 274 \r
 275 \r
 276 \r
 277 \r
 278 \r
 279 \r
 280 \r
 281 \r