src/jalview/appletgui/AlignViewport.java

   1 /*\r
   2  * Jalview - A Sequence Alignment Editor and Viewer\r
   3  * Copyright (C) 2006 AM Waterhouse, J Procter, G Barton, M Clamp, S Searle\r
   4  *\r
   5  * This program is free software; you can redistribute it and/or\r
   6  * modify it under the terms of the GNU General Public License\r
   7  * as published by the Free Software Foundation; either version 2\r
   8  * of the License, or (at your option) any later version.\r
   9  *\r
  10  * This program is distributed in the hope that it will be useful,\r
  11  * but WITHOUT ANY WARRANTY; without even the implied warranty of\r
  12  * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the\r
  13  * GNU General Public License for more details.\r
  14  *\r
  15  * You should have received a copy of the GNU General Public License\r
  16  * along with this program; if not, write to the Free Software\r
  17  * Foundation, Inc., 51 Franklin Street, Fifth Floor, Boston, MA  02110-1301, USA\r
  18  */\r
  19 \r
  20 package jalview.appletgui;\r
  21 \r
  22 import java.util.*;\r
  23 \r
  24 import java.awt.*;\r
  25 \r
  26 import jalview.analysis.*;\r
  27 import jalview.bin.*;\r
  28 import jalview.datamodel.*;\r
  29 import jalview.schemes.*;\r
  30 \r
  31 public class AlignViewport\r
  32 {\r
  33   int startRes;\r
  34   int endRes;\r
  35 \r
  36   int startSeq;\r
  37   int endSeq;\r
  38 \r
  39 \r
  40   boolean cursorMode = false;\r
  41 \r
  42   boolean showJVSuffix = true;\r
  43   boolean showText = true;\r
  44   boolean showColourText = false;\r
  45   boolean showBoxes = true;\r
  46   boolean wrapAlignment = false;\r
  47   boolean renderGaps = true;\r
  48   boolean showSequenceFeatures = false;\r
  49   boolean showAnnotation = true;\r
  50   boolean showConservation = true;\r
  51   boolean showQuality = true;\r
  52   boolean showConsensus = true;\r
  53 \r
  54   boolean colourAppliesToAllGroups = true;\r
  55   ColourSchemeI globalColourScheme = null;\r
  56   boolean conservationColourSelected = false;\r
  57   boolean abovePIDThreshold = false;\r
  58 \r
  59   SequenceGroup selectionGroup;\r
  60 \r
  61   int charHeight;\r
  62   int charWidth;\r
  63   int wrappedWidth;\r
  64 \r
  65   Font font = new Font("SansSerif", Font.PLAIN, 10);\r
  66   boolean validCharWidth = true;\r
  67   AlignmentI alignment;\r
  68 \r
  69   ColumnSelection colSel = new ColumnSelection();\r
  70 \r
  71   int threshold;\r
  72   int increment;\r
  73 \r
  74   NJTree currentTree = null;\r
  75 \r
  76   boolean scaleAboveWrapped = true;\r
  77   boolean scaleLeftWrapped = true;\r
  78   boolean scaleRightWrapped = true;\r
  79 \r
  80   // The following vector holds the features which are\r
  81  // currently visible, in the correct order or rendering\r
  82   Hashtable featuresDisplayed;\r
  83 \r
  84   boolean hasHiddenColumns = false;\r
  85   boolean hasHiddenRows = false;\r
  86   boolean showHiddenMarkers = true;\r
  87 \r
  88 \r
  89   public Vector vconsensus;\r
  90   AlignmentAnnotation consensus;\r
  91   AlignmentAnnotation conservation;\r
  92   AlignmentAnnotation quality;\r
  93 \r
  94   boolean autocalculateConsensus = true;\r
  95 \r
  96   public int ConsPercGaps = 25; // JBPNote : This should be a scalable property!\r
  97 \r
  98   private java.beans.PropertyChangeSupport changeSupport = new java.beans.PropertyChangeSupport(this);\r
  99 \r
 100   boolean ignoreGapsInConsensusCalculation = false;\r
 101 \r
 102   jalview.bin.JalviewLite applet;\r
 103 \r
 104   boolean MAC = false;\r
 105 \r
 106   public AlignViewport(AlignmentI al, JalviewLite applet)\r
 107   {\r
 108     this.applet = applet;\r
 109     setAlignment(al);\r
 110     this.startRes = 0;\r
 111     this.endRes = al.getWidth() - 1;\r
 112     this.startSeq = 0;\r
 113     this.endSeq = al.getHeight() - 1;\r
 114     setFont(font);\r
 115 \r
 116     if(System.getProperty("os.name").startsWith("Mac"))\r
 117       MAC = true;\r
 118 \r
 119     if (applet != null)\r
 120     {\r
 121       String param = applet.getParameter("showFullId");\r
 122       if (param != null)\r
 123       {\r
 124         showJVSuffix = Boolean.valueOf(param).booleanValue();\r
 125       }\r
 126 \r
 127       param = applet.getParameter("showAnnotation");\r
 128       if (param != null)\r
 129       {\r
 130         showAnnotation = Boolean.valueOf(param).booleanValue();\r
 131       }\r
 132 \r
 133       param = applet.getParameter("showConservation");\r
 134       if (param != null)\r
 135       {\r
 136         showConservation = Boolean.valueOf(param).booleanValue();\r
 137       }\r
 138 \r
 139       param = applet.getParameter("showQuality");\r
 140       if (param != null)\r
 141       {\r
 142         showQuality = Boolean.valueOf(param).booleanValue();\r
 143       }\r
 144 \r
 145       param = applet.getParameter("showConsensus");\r
 146       if (param != null)\r
 147       {\r
 148         showConsensus = Boolean.valueOf(param).booleanValue();\r
 149       }\r
 150     }\r
 151     // We must set conservation and consensus before setting colour,\r
 152     // as Blosum and Clustal require this to be done\r
 153     updateConservation();\r
 154     updateConsensus();\r
 155 \r
 156 \r
 157     if (applet != null)\r
 158     {\r
 159       String colour = applet.getParameter("defaultColour");\r
 160 \r
 161       if(colour == null)\r
 162       {\r
 163         colour = applet.getParameter("userDefinedColour");\r
 164         if(colour !=null)\r
 165           colour = "User Defined";\r
 166       }\r
 167 \r
 168       if(colour != null)\r
 169       {\r
 170         globalColourScheme = ColourSchemeProperty.getColour(alignment, colour);\r
 171         if (globalColourScheme != null)\r
 172         {\r
 173           globalColourScheme.setConsensus(vconsensus);\r
 174         }\r
 175       }\r
 176 \r
 177       if(applet.getParameter("userDefinedColour")!=null)\r
 178       {\r
 179         ((UserColourScheme)globalColourScheme).parseAppletParameter(\r
 180             applet.getParameter("userDefinedColour"));\r
 181       }\r
 182 \r
 183 \r
 184     }\r
 185   }\r
 186 \r
 187   public void showSequenceFeatures(boolean b)\r
 188   {\r
 189     showSequenceFeatures = b;\r
 190   }\r
 191 \r
 192   public boolean getShowSequenceFeatures()\r
 193   {\r
 194     return showSequenceFeatures;\r
 195   }\r
 196 \r
 197 \r
 198   public void updateConservation()\r
 199   {\r
 200     if(alignment.isNucleotide())\r
 201           return;\r
 202 \r
 203     Conservation cons = new jalview.analysis.Conservation("All",\r
 204         jalview.schemes.ResidueProperties.propHash, 3,\r
 205         alignment.getSequences(), 0,\r
 206         alignment.getWidth() - 1);\r
 207     cons.calculate();\r
 208     cons.verdict(false, ConsPercGaps);\r
 209     cons.findQuality();\r
 210     int alWidth = alignment.getWidth();\r
 211     Annotation[] annotations = new Annotation[alWidth];\r
 212     Annotation[] qannotations = new Annotation[alWidth];\r
 213     String sequence = cons.getConsSequence().getSequence();\r
 214     float minR, minG, minB, maxR, maxG, maxB;\r
 215     minR = 0.3f;\r
 216     minG = 0.0f;\r
 217     minB = 0f;\r
 218     maxR = 1.0f - minR;\r
 219     maxG = 0.9f - minG;\r
 220     maxB = 0f - minB; // scalable range for colouring both Conservation and Quality\r
 221     float min = 0f;\r
 222     float max = 11f;\r
 223     float qmin = cons.qualityRange[0].floatValue();\r
 224     float qmax = cons.qualityRange[1].floatValue();\r
 225 \r
 226     for (int i = 0; i < alWidth; i++)\r
 227     {\r
 228       float value = 0;\r
 229       try\r
 230       {\r
 231         value = Integer.parseInt(sequence.charAt(i) + "");\r
 232       }\r
 233       catch (Exception ex)\r
 234       {\r
 235         if (sequence.charAt(i) == '*')\r
 236         {\r
 237           value = 11;\r
 238         }\r
 239         if (sequence.charAt(i) == '+')\r
 240         {\r
 241           value = 10;\r
 242         }\r
 243       }\r
 244       float vprop = value - min;\r
 245       vprop /= max;\r
 246 \r
 247       annotations[i] = new Annotation(sequence.charAt(i) + "",\r
 248                                       "", ' ', value,\r
 249                                       new Color(minR + maxR * vprop,\r
 250                                                 minG + maxG * vprop,\r
 251                                                 minB + maxB * vprop));\r
 252       // Quality calc\r
 253       value = ( (Double) cons.quality.elementAt(i)).floatValue();\r
 254       vprop = value - qmin;\r
 255       vprop /= qmax;\r
 256       qannotations[i] = new Annotation(" ",\r
 257                                        String.valueOf(value), ' ', value,\r
 258                                        new\r
 259                                        Color(minR + maxR * vprop,\r
 260                                              minG + maxG * vprop,\r
 261                                              minB + maxB * vprop));\r
 262     }\r
 263 \r
 264     if (conservation == null)\r
 265     {\r
 266       conservation = new AlignmentAnnotation("Conservation",\r
 267                                              "Conservation of total alignment less than " +\r
 268                                              ConsPercGaps + "% gaps",\r
 269                                              annotations,\r
 270                                              0f, // cons.qualityRange[0].floatValue(),\r
 271                                              11f, // cons.qualityRange[1].floatValue()\r
 272                                              AlignmentAnnotation.BAR_GRAPH);\r
 273       if (showConservation)\r
 274       {\r
 275         alignment.addAnnotation(conservation);\r
 276       }\r
 277       quality = new AlignmentAnnotation("Quality",\r
 278                                         "Alignment Quality based on Blosum62 scores",\r
 279                                         qannotations,\r
 280                                         cons.qualityRange[0].floatValue(),\r
 281                                         cons.qualityRange[1].floatValue(),\r
 282                                         AlignmentAnnotation.BAR_GRAPH);\r
 283       if (showQuality)\r
 284       {\r
 285         alignment.addAnnotation(quality);\r
 286       }\r
 287     }\r
 288     else\r
 289     {\r
 290       conservation.annotations = annotations;\r
 291       quality.annotations = qannotations;\r
 292       quality.graphMax = cons.qualityRange[1].floatValue();\r
 293     }\r
 294 \r
 295   }\r
 296 \r
 297   public void updateConsensus()\r
 298   {\r
 299     Annotation[] annotations = new Annotation[alignment.getWidth()];\r
 300 \r
 301     // this routine prevents vconsensus becoming a new object each time\r
 302     // consenus is calculated. Important for speed of Blosum62\r
 303     // and PID colouring of alignment\r
 304     if (vconsensus == null)\r
 305     {\r
 306       vconsensus = alignment.getAAFrequency();\r
 307     }\r
 308     else\r
 309     {\r
 310       Vector temp = alignment.getAAFrequency();\r
 311       vconsensus.removeAllElements();\r
 312       Enumeration e = temp.elements();\r
 313       while (e.hasMoreElements())\r
 314       {\r
 315         vconsensus.addElement(e.nextElement());\r
 316       }\r
 317     }\r
 318     Hashtable hash = null;\r
 319     for (int i = 0; i < alignment.getWidth(); i++)\r
 320     {\r
 321       hash = (Hashtable) vconsensus.elementAt(i);\r
 322       float value = 0;\r
 323       if(ignoreGapsInConsensusCalculation)\r
 324         value = ((Float)hash.get("pid_nogaps")).floatValue();\r
 325       else\r
 326         value = ((Float)hash.get("pid_gaps")).floatValue();\r
 327 \r
 328       String maxRes = hash.get("maxResidue").toString();\r
 329       String mouseOver = hash.get("maxResidue") + " ";\r
 330       if (maxRes.length() > 1)\r
 331       {\r
 332         mouseOver = "[" + maxRes + "] ";\r
 333         maxRes = "+";\r
 334       }\r
 335 \r
 336 \r
 337       mouseOver += (int) value + "%";\r
 338       annotations[i] = new Annotation(maxRes, mouseOver, ' ', value);\r
 339 \r
 340     }\r
 341 \r
 342     if (consensus == null)\r
 343     {\r
 344       consensus = new AlignmentAnnotation("Consensus",\r
 345                                           "PID", annotations, 0f, 100f, AlignmentAnnotation.BAR_GRAPH);\r
 346       if (showConsensus)\r
 347       {\r
 348         alignment.addAnnotation(consensus);\r
 349       }\r
 350     }\r
 351     else\r
 352     {\r
 353       consensus.annotations = annotations;\r
 354     }\r
 355 \r
 356     if(globalColourScheme!=null)\r
 357           globalColourScheme.setConsensus(vconsensus);\r
 358 \r
 359   }\r
 360 \r
 361   public SequenceGroup getSelectionGroup()\r
 362   {\r
 363     return selectionGroup;\r
 364   }\r
 365 \r
 366   public void setSelectionGroup(SequenceGroup sg)\r
 367   {\r
 368     selectionGroup = sg;\r
 369   }\r
 370 \r
 371   public boolean getConservationSelected()\r
 372   {\r
 373     return conservationColourSelected;\r
 374   }\r
 375 \r
 376   public void setConservationSelected(boolean b)\r
 377   {\r
 378     conservationColourSelected = b;\r
 379   }\r
 380 \r
 381   public boolean getAbovePIDThreshold()\r
 382   {\r
 383     return abovePIDThreshold;\r
 384   }\r
 385 \r
 386   public void setAbovePIDThreshold(boolean b)\r
 387   {\r
 388     abovePIDThreshold = b;\r
 389   }\r
 390 \r
 391   public int getStartRes()\r
 392   {\r
 393     return startRes;\r
 394   }\r
 395 \r
 396   public int getEndRes()\r
 397   {\r
 398     return endRes;\r
 399   }\r
 400 \r
 401   public int getStartSeq()\r
 402   {\r
 403     return startSeq;\r
 404   }\r
 405 \r
 406   public void setGlobalColourScheme(ColourSchemeI cs)\r
 407   {\r
 408     globalColourScheme = cs;\r
 409   }\r
 410 \r
 411   public ColourSchemeI getGlobalColourScheme()\r
 412   {\r
 413     return globalColourScheme;\r
 414   }\r
 415 \r
 416   public void setStartRes(int res)\r
 417   {\r
 418     this.startRes = res;\r
 419   }\r
 420 \r
 421   public void setStartSeq(int seq)\r
 422   {\r
 423     this.startSeq = seq;\r
 424   }\r
 425 \r
 426   public void setEndRes(int res)\r
 427   {\r
 428     if (res > alignment.getWidth() - 1)\r
 429     {\r
 430       // log.System.out.println(" Corrected res from " + res + " to maximum " + (alignment.getWidth()-1));\r
 431       res = alignment.getWidth() - 1;\r
 432     }\r
 433     if (res < 0)\r
 434     {\r
 435       res = 0;\r
 436     }\r
 437     this.endRes = res;\r
 438   }\r
 439 \r
 440   public void setEndSeq(int seq)\r
 441   {\r
 442     if (seq > alignment.getHeight())\r
 443     {\r
 444       seq = alignment.getHeight();\r
 445     }\r
 446     if (seq < 0)\r
 447     {\r
 448       seq = 0;\r
 449     }\r
 450     this.endSeq = seq;\r
 451   }\r
 452 \r
 453   public int getEndSeq()\r
 454   {\r
 455     return endSeq;\r
 456   }\r
 457 \r
 458   java.awt.Frame nullFrame;\r
 459   public void setFont(Font f)\r
 460   {\r
 461     font = f;\r
 462     if(nullFrame == null)\r
 463     {\r
 464       nullFrame = new java.awt.Frame();\r
 465       nullFrame.addNotify();\r
 466     }\r
 467 \r
 468     java.awt.FontMetrics fm = nullFrame.getGraphics().getFontMetrics(font);\r
 469     setCharHeight(fm.getHeight());\r
 470     charWidth = fm.charWidth('M');\r
 471   }\r
 472 \r
 473   public Font getFont()\r
 474   {\r
 475     return font;\r
 476   }\r
 477 \r
 478   public int getCharWidth()\r
 479   {\r
 480     return charWidth;\r
 481   }\r
 482 \r
 483   public void setCharHeight(int h)\r
 484   {\r
 485     this.charHeight = h;\r
 486   }\r
 487 \r
 488   public int getCharHeight()\r
 489   {\r
 490     return charHeight;\r
 491   }\r
 492 \r
 493   public void setWrappedWidth(int w)\r
 494   {\r
 495     this.wrappedWidth = w;\r
 496   }\r
 497 \r
 498   public int getwrappedWidth()\r
 499   {\r
 500     return wrappedWidth;\r
 501   }\r
 502 \r
 503   public AlignmentI getAlignment()\r
 504   {\r
 505     return alignment;\r
 506   }\r
 507 \r
 508   public void setAlignment(AlignmentI align)\r
 509   {\r
 510     this.alignment = align;\r
 511   }\r
 512 \r
 513   public void setWrapAlignment(boolean state)\r
 514   {\r
 515     wrapAlignment = state;\r
 516   }\r
 517 \r
 518   public void setShowText(boolean state)\r
 519   {\r
 520     showText = state;\r
 521   }\r
 522 \r
 523   public void setRenderGaps(boolean state)\r
 524   {\r
 525     renderGaps = state;\r
 526   }\r
 527 \r
 528   public boolean getColourText()\r
 529   {\r
 530     return showColourText;\r
 531   }\r
 532 \r
 533   public void setColourText(boolean state)\r
 534   {\r
 535     showColourText = state;\r
 536   }\r
 537 \r
 538   public void setShowBoxes(boolean state)\r
 539   {\r
 540     showBoxes = state;\r
 541   }\r
 542 \r
 543   public boolean getWrapAlignment()\r
 544   {\r
 545     return wrapAlignment;\r
 546   }\r
 547 \r
 548   public boolean getShowText()\r
 549   {\r
 550     return showText;\r
 551   }\r
 552 \r
 553   public boolean getShowBoxes()\r
 554   {\r
 555     return showBoxes;\r
 556   }\r
 557 \r
 558   public char getGapCharacter()\r
 559   {\r
 560     return getAlignment().getGapCharacter();\r
 561   }\r
 562 \r
 563   public void setGapCharacter(char gap)\r
 564   {\r
 565     if (getAlignment() != null)\r
 566     {\r
 567       getAlignment().setGapCharacter(gap);\r
 568     }\r
 569   }\r
 570 \r
 571   public void setThreshold(int thresh)\r
 572   {\r
 573     threshold = thresh;\r
 574   }\r
 575 \r
 576   public int getThreshold()\r
 577   {\r
 578     return threshold;\r
 579   }\r
 580 \r
 581   public void setIncrement(int inc)\r
 582   {\r
 583     increment = inc;\r
 584   }\r
 585 \r
 586   public int getIncrement()\r
 587   {\r
 588     return increment;\r
 589   }\r
 590 \r
 591   public void setHiddenColumns(ColumnSelection colsel)\r
 592   {\r
 593     this.colSel = colsel;\r
 594     if(colSel.getHiddenColumns()!=null)\r
 595       hasHiddenColumns = true;\r
 596   }\r
 597 \r
 598   public ColumnSelection getColumnSelection()\r
 599   {\r
 600     return colSel;\r
 601   }\r
 602 \r
 603   public void resetSeqLimits(int height)\r
 604   {\r
 605     setEndSeq(height / getCharHeight());\r
 606   }\r
 607 \r
 608   public void setCurrentTree(NJTree tree)\r
 609   {\r
 610     currentTree = tree;\r
 611   }\r
 612 \r
 613   public NJTree getCurrentTree()\r
 614   {\r
 615     return currentTree;\r
 616   }\r
 617 \r
 618   public void setColourAppliesToAllGroups(boolean b)\r
 619   {\r
 620     colourAppliesToAllGroups = b;\r
 621   }\r
 622 \r
 623   public boolean getColourAppliesToAllGroups()\r
 624   {\r
 625     return colourAppliesToAllGroups;\r
 626   }\r
 627 \r
 628   public boolean getShowJVSuffix()\r
 629   {\r
 630     return showJVSuffix;\r
 631   }\r
 632 \r
 633   public void setShowJVSuffix(boolean b)\r
 634   {\r
 635     showJVSuffix = b;\r
 636   }\r
 637 \r
 638   public boolean getShowAnnotation()\r
 639   {\r
 640     return showAnnotation;\r
 641   }\r
 642 \r
 643   public void setShowAnnotation(boolean b)\r
 644   {\r
 645     showAnnotation = b;\r
 646   }\r
 647 \r
 648   public boolean getScaleAboveWrapped()\r
 649   {\r
 650     return scaleAboveWrapped;\r
 651   }\r
 652 \r
 653   public boolean getScaleLeftWrapped()\r
 654   {\r
 655     return scaleLeftWrapped;\r
 656   }\r
 657 \r
 658   public boolean getScaleRightWrapped()\r
 659   {\r
 660     return scaleRightWrapped;\r
 661   }\r
 662 \r
 663   public void setScaleAboveWrapped(boolean b)\r
 664   {\r
 665     scaleAboveWrapped = b;\r
 666   }\r
 667 \r
 668   public void setScaleLeftWrapped(boolean b)\r
 669   {\r
 670     scaleLeftWrapped = b;\r
 671   }\r
 672 \r
 673   public void setScaleRightWrapped(boolean b)\r
 674   {\r
 675     scaleRightWrapped = b;\r
 676   }\r
 677 \r
 678   public void setIgnoreGapsConsensus(boolean b)\r
 679   {\r
 680     ignoreGapsInConsensusCalculation = b;\r
 681     updateConsensus();\r
 682     if (globalColourScheme!=null)\r
 683     {\r
 684       globalColourScheme.setThreshold(globalColourScheme.getThreshold(),\r
 685           ignoreGapsInConsensusCalculation);\r
 686 \r
 687     }\r
 688   }\r
 689 \r
 690   /**\r
 691    * Property change listener for changes in alignment\r
 692    *\r
 693    * @param listener DOCUMENT ME!\r
 694    */\r
 695   public void addPropertyChangeListener(\r
 696       java.beans.PropertyChangeListener listener)\r
 697   {\r
 698       changeSupport.addPropertyChangeListener(listener);\r
 699   }\r
 700 \r
 701   /**\r
 702    * DOCUMENT ME!\r
 703    *\r
 704    * @param listener DOCUMENT ME!\r
 705    */\r
 706   public void removePropertyChangeListener(\r
 707       java.beans.PropertyChangeListener listener)\r
 708   {\r
 709       changeSupport.removePropertyChangeListener(listener);\r
 710   }\r
 711 \r
 712   /**\r
 713    * Property change listener for changes in alignment\r
 714    *\r
 715    * @param prop DOCUMENT ME!\r
 716    * @param oldvalue DOCUMENT ME!\r
 717    * @param newvalue DOCUMENT ME!\r
 718    */\r
 719   public void firePropertyChange(String prop, Object oldvalue, Object newvalue)\r
 720   {\r
 721       changeSupport.firePropertyChange(prop, oldvalue, newvalue);\r
 722   }\r
 723 \r
 724 \r
 725 \r
 726   public boolean getIgnoreGapsConsensus()\r
 727   {\r
 728     return ignoreGapsInConsensusCalculation;\r
 729   }\r
 730   public void hideSelectedColumns()\r
 731   {\r
 732     if (colSel.size() < 1)\r
 733       return;\r
 734 \r
 735     colSel.hideSelectedColumns();\r
 736     setSelectionGroup(null);\r
 737 \r
 738     hasHiddenColumns = true;\r
 739   }\r
 740 \r
 741   public void invertColumnSelection()\r
 742   {\r
 743     int column;\r
 744     for (int i = 0; i < alignment.getWidth(); i++)\r
 745     {\r
 746       column = i;\r
 747 \r
 748       if (colSel.contains(column))\r
 749         colSel.removeElement(column);\r
 750       else\r
 751         colSel.addElement(column);\r
 752 \r
 753     }\r
 754   }\r
 755 \r
 756 \r
 757   public void hideColumns(int start, int end)\r
 758   {\r
 759     if(start==end)\r
 760       colSel.hideColumns(start);\r
 761     else\r
 762       colSel.hideColumns(start, end);\r
 763 \r
 764     hasHiddenColumns = true;\r
 765   }\r
 766 \r
 767   public void hideSequence(SequenceI seq)\r
 768   {\r
 769     if(seq!=null)\r
 770     {\r
 771       alignment.getHiddenSequences().hideSequence(seq);\r
 772       hasHiddenRows = true;\r
 773       firePropertyChange("alignment", null, alignment.getSequences());\r
 774     }\r
 775   }\r
 776 \r
 777   public void hideAllSelectedSeqs()\r
 778   {\r
 779     if (selectionGroup == null)\r
 780       return;\r
 781 \r
 782     SequenceI[] seqs = selectionGroup.getSequencesInOrder(alignment);\r
 783 \r
 784     for (int i = 0; i < seqs.length; i++)\r
 785     {\r
 786       alignment.getHiddenSequences().hideSequence(seqs[i]);\r
 787     }\r
 788     firePropertyChange("alignment", null, alignment.getSequences());\r
 789     hasHiddenRows = true;\r
 790     setSelectionGroup(null);\r
 791   }\r
 792 \r
 793   public void showColumn(int col)\r
 794   {\r
 795     colSel.revealHiddenColumns(col);\r
 796     if(colSel.getHiddenColumns()==null)\r
 797       hasHiddenColumns = false;\r
 798   }\r
 799 \r
 800   public void showAllHiddenColumns()\r
 801   {\r
 802     colSel.revealAllHiddenColumns();\r
 803     hasHiddenColumns = false;\r
 804   }\r
 805 \r
 806   public void showAllHiddenSeqs()\r
 807   {\r
 808     if(alignment.getHiddenSequences().getSize()>0)\r
 809     {\r
 810       if(selectionGroup==null)\r
 811       {\r
 812         selectionGroup = new SequenceGroup();\r
 813         selectionGroup.setEndRes(alignment.getWidth()-1);\r
 814       }\r
 815       Vector tmp = alignment.getHiddenSequences().showAll();\r
 816       for(int t=0; t<tmp.size(); t++)\r
 817       {\r
 818         selectionGroup.addSequence(\r
 819             (SequenceI)tmp.elementAt(t), false\r
 820             );\r
 821       }\r
 822       firePropertyChange("alignment", null, alignment.getSequences());\r
 823       hasHiddenRows = false;\r
 824     }\r
 825   }\r
 826 \r
 827   public int adjustForHiddenSeqs(int alignmentIndex)\r
 828   {\r
 829     return alignment.getHiddenSequences().adjustForHiddenSeqs(alignmentIndex);\r
 830   }\r
 831 \r
 832   /**\r
 833    * This method returns the a new SequenceI [] with\r
 834    * the selection sequence and start and end points adjusted\r
 835    * @return String[]\r
 836    */\r
 837   public SequenceI[] getSelectionAsNewSequence()\r
 838   {\r
 839     SequenceI[] sequences;\r
 840 \r
 841     if (selectionGroup == null)\r
 842       sequences = alignment.getSequencesArray();\r
 843     else\r
 844       sequences = selectionGroup.getSelectionAsNewSequences(alignment);\r
 845 \r
 846     return sequences;\r
 847   }\r
 848 \r
 849   /**\r
 850    * This method returns the visible alignment as text, as\r
 851    * seen on the GUI, ie if columns are hidden they will not\r
 852    * be returned in the result.\r
 853    * Use this for calculating trees, PCA, redundancy etc on views\r
 854    * which contain hidden columns.\r
 855    * @return String[]\r
 856    */\r
 857   public jalview.datamodel.CigarArray getViewAsCigars(boolean selectedRegionOnly)\r
 858   {\r
 859     CigarArray selection=null;\r
 860     SequenceI [] seqs= null;\r
 861     int i, iSize;\r
 862     int start = 0, end = 0;\r
 863     if(selectedRegionOnly && selectionGroup!=null)\r
 864     {\r
 865       iSize = selectionGroup.getSize(false);\r
 866       seqs = selectionGroup.getSequencesInOrder(alignment);\r
 867       start = selectionGroup.getStartRes();\r
 868       end = selectionGroup.getEndRes(); // inclusive for start and end in SeqCigar constructor\r
 869     }\r
 870     else\r
 871     {\r
 872       iSize = alignment.getHeight();\r
 873       seqs = alignment.getSequencesArray();\r
 874       end = alignment.getWidth()-1;\r
 875     }\r
 876     SeqCigar[] selseqs = new SeqCigar[iSize];\r
 877     for(i=0; i<iSize; i++)\r
 878     {\r
 879       selseqs[i] = new SeqCigar(seqs[i], start, end);\r
 880     }\r
 881     selection=new CigarArray(selseqs);\r
 882     // now construct the CigarArray operations\r
 883     if (hasHiddenColumns) {\r
 884       Vector regions = colSel.getHiddenColumns();\r
 885       int [] region;\r
 886       int hideStart, hideEnd;\r
 887       int last=start;\r
 888       for (int j = 0; last<end & j < regions.size(); j++)\r
 889       {\r
 890         region = (int[]) regions.elementAt(j);\r
 891         hideStart = region[0];\r
 892         hideEnd = region[1];\r
 893         // edit hidden regions to selection range\r
 894         if(hideStart<last) {\r
 895           if (hideEnd > last)\r
 896           {\r
 897             hideStart = last;\r
 898           } else\r
 899             continue;\r
 900         }\r
 901 \r
 902         if (hideStart>end)\r
 903           break;\r
 904 \r
 905         if (hideEnd>end)\r
 906           hideEnd=end;\r
 907 \r
 908         if (hideStart>hideEnd)\r
 909           break;\r
 910         /**\r
 911          * form operations...\r
 912          */\r
 913         if (last<hideStart)\r
 914           selection.addOperation(CigarArray.M, hideStart-last);\r
 915         selection.addOperation(CigarArray.D, 1+hideEnd-hideStart);\r
 916         last = hideEnd+1;\r
 917       }\r
 918       // Final match if necessary.\r
 919       if (last<end)\r
 920         selection.addOperation(CigarArray.M, end-last);\r
 921     } else {\r
 922       selection.addOperation(CigarArray.M, end-start);\r
 923     }\r
 924     return selection;\r
 925   }\r
 926   /**\r
 927    * return a compact representation of the current alignment selection to\r
 928    * pass to an analysis function\r
 929    * @param selectedOnly boolean true to just return the selected view\r
 930    * @return AlignmentView\r
 931    */\r
 932   jalview.datamodel.AlignmentView getAlignmentView(boolean selectedOnly) {\r
 933     // JBPNote:\r
 934     // this is here because the AlignmentView constructor modifies the CigarArray\r
 935     // object. Refactoring of Cigar and alignment view representation should\r
 936     // be done to remove redundancy.\r
 937     CigarArray aligview = getViewAsCigars(selectedOnly);\r
 938     if (aligview!=null)\r
 939       return new AlignmentView(aligview);\r
 940     return null;\r
 941   }\r
 942   /**\r
 943    * This method returns the visible alignment as text, as\r
 944    * seen on the GUI, ie if columns are hidden they will not\r
 945    * be returned in the result.\r
 946    * Use this for calculating trees, PCA, redundancy etc on views\r
 947    * which contain hidden columns.\r
 948    * @return String[]\r
 949    */\r
 950   public String [] getViewAsString(boolean selectedRegionOnly)\r
 951   {\r
 952     String [] selection = null;\r
 953     SequenceI [] seqs= null;\r
 954     int i, iSize;\r
 955     int start = 0, end = 0;\r
 956     if(selectedRegionOnly && selectionGroup!=null)\r
 957     {\r
 958       iSize = selectionGroup.getSize(false);\r
 959       seqs = selectionGroup.getSequencesInOrder(alignment);\r
 960       start = selectionGroup.getStartRes();\r
 961       end = selectionGroup.getEndRes()+1;\r
 962     }\r
 963     else\r
 964     {\r
 965       iSize = alignment.getHeight();\r
 966       seqs = alignment.getSequencesArray();\r
 967       end = alignment.getWidth();\r
 968     }\r
 969 \r
 970     selection = new String[iSize];\r
 971 \r
 972 \r
 973     for(i=0; i<iSize; i++)\r
 974     {\r
 975       if (hasHiddenColumns)\r
 976       {\r
 977            StringBuffer visibleSeq = new StringBuffer();\r
 978            Vector regions = colSel.getHiddenColumns();\r
 979 \r
 980            int blockStart = start, blockEnd=end;\r
 981            int [] region;\r
 982            int hideStart, hideEnd;\r
 983 \r
 984            for (int j = 0; j < regions.size(); j++)\r
 985            {\r
 986              region = (int[]) regions.elementAt(j);\r
 987              hideStart = region[0];\r
 988              hideEnd = region[1];\r
 989 \r
 990              if(hideStart < start)\r
 991              {\r
 992                continue;\r
 993              }\r
 994 \r
 995              blockStart = Math.min(blockStart, hideEnd+1);\r
 996              blockEnd = Math.min(blockEnd, hideStart);\r
 997 \r
 998              if(blockStart>blockEnd)\r
 999              {\r
1000                 break;\r
1001              }\r
1002 \r
1003 \r
1004              visibleSeq.append(seqs[i].getSequence(blockStart, blockEnd));\r
1005 \r
1006              blockStart = hideEnd+1;\r
1007              blockEnd = end;\r
1008            }\r
1009 \r
1010            if(end>blockStart)\r
1011              visibleSeq.append(seqs[i].getSequence(blockStart, end));\r
1012 \r
1013            selection[i] = visibleSeq.toString();\r
1014       }\r
1015       else\r
1016       {\r
1017         selection[i] = seqs[i].getSequence(start, end);\r
1018       }\r
1019     }\r
1020 \r
1021     return selection;\r
1022   }\r
1023 \r
1024   public boolean getShowHiddenMarkers()\r
1025   {\r
1026     return showHiddenMarkers;\r
1027   }\r
1028 \r
1029   public void setShowHiddenMarkers(boolean show)\r
1030   {\r
1031     showHiddenMarkers = show;\r
1032   }\r
1033 \r
1034 \r
1035 }\r