src/jalview/io/AnnotationFile.java

   1 /*
   2  * Jalview - A Sequence Alignment Editor and Viewer ($$Version-Rel$$)
   3  * Copyright (C) $$Year-Rel$$ The Jalview Authors
   4  *
   5  * This file is part of Jalview.
   6  *
   7  * Jalview is free software: you can redistribute it and/or
   8  * modify it under the terms of the GNU General Public License
   9  * as published by the Free Software Foundation, either version 3
  10  * of the License, or (at your option) any later version.
  11  *
  12  * Jalview is distributed in the hope that it will be useful, but
  13  * WITHOUT ANY WARRANTY; without even the implied warranty
  14  * of MERCHANTABILITY or FITNESS FOR A PARTICULAR
  15  * PURPOSE.  See the GNU General Public License for more details.
  16  *
  17  * You should have received a copy of the GNU General Public License
  18  * along with Jalview.  If not, see <http://www.gnu.org/licenses/>.
  19  * The Jalview Authors are detailed in the 'AUTHORS' file.
  20  */
  21 package jalview.io;
  22
  23 import jalview.analysis.Conservation;
  24 import jalview.api.AlignViewportI;
  25 import jalview.datamodel.AlignmentAnnotation;
  26 import jalview.datamodel.AlignmentI;
  27 import jalview.datamodel.Annotation;
  28 import jalview.datamodel.ColumnSelection;
  29 import jalview.datamodel.GraphLine;
  30 import jalview.datamodel.HiddenSequences;
  31 import jalview.datamodel.SequenceGroup;
  32 import jalview.datamodel.SequenceI;
  33 import jalview.jsdev.GenericFileAdapter;
  34 import jalview.schemes.ColourSchemeI;
  35 import jalview.schemes.ColourSchemeProperty;
  36 import jalview.schemes.ResidueProperties;
  37 import jalview.schemes.UserColourScheme;
  38 import jalview.util.Comparison;
  39 import jalview.util.Format;
  40
  41 import java.io.BufferedReader;
  42 import java.io.InputStreamReader;
  43 import java.io.StringReader;
  44 import java.util.ArrayList;
  45 import java.util.BitSet;
  46 import java.util.Enumeration;
  47 import java.util.Hashtable;
  48 import java.util.List;
  49 import java.util.Map;
  50 import java.util.StringTokenizer;
  51 import java.util.Vector;
  52
  53 public class AnnotationFile
  54 {
  55   public AnnotationFile()
  56   {
  57     init();
  58   }
  59
  60   /**
  61    * character used to write newlines
  62    */
  63   protected String newline = System.getProperty("line.separator");
  64
  65   /**
  66    * set new line string and reset the output buffer
  67    *
  68    * @param nl
  69    */
  70   public void setNewlineString(String nl)
  71   {
  72     newline = nl;
  73     init();
  74   }
  75
  76   public String getNewlineString()
  77   {
  78     return newline;
  79   }
  80
  81   StringBuffer text;
  82
  83   private void init()
  84   {
  85     text = new StringBuffer("JALVIEW_ANNOTATION" + newline + "# Created: "
  86             + new java.util.Date() + newline + newline);
  87     refSeq = null;
  88     refSeqId = null;
  89   }
  90
  91   /**
  92    * convenience method for pre-2.8.3 annotation files which have no view,
  93    * hidden columns or hidden row keywords.
  94    *
  95    * @param annotations
  96    * @param list
  97    * @param properties
  98    * @return annotation file as a string.
  99    */
 100   public String printAnnotations(AlignmentAnnotation[] annotations,
 101           List<SequenceGroup> list, Hashtable properties)
 102   {
 103     return printAnnotations(annotations, list, properties, null, null, null);
 104
 105   }
 106
 107   /**
 108    * hold all the information about a particular view definition read from or
 109    * written out in an annotations file.
 110    */
 111   public class ViewDef
 112   {
 113     public String viewname;
 114
 115     public HiddenSequences hidseqs;
 116
 117     public ColumnSelection hiddencols;
 118
 119     public Vector visibleGroups;
 120
 121     public Hashtable hiddenRepSeqs;
 122
 123     public ViewDef(String viewname, HiddenSequences hidseqs,
 124             ColumnSelection hiddencols, Hashtable hiddenRepSeqs)
 125     {
 126       this.viewname = viewname;
 127       this.hidseqs = hidseqs;
 128       this.hiddencols = hiddencols;
 129       this.hiddenRepSeqs = hiddenRepSeqs;
 130     }
 131   }
 132
 133   /**
 134    * Prepare an annotation file given a set of annotations, groups, alignment
 135    * properties and views.
 136    *
 137    * @param annotations
 138    * @param list
 139    * @param properties
 140    * @param views
 141    * @return annotation file
 142    */
 143   public String printAnnotations(AlignmentAnnotation[] annotations,
 144           List<SequenceGroup> list, Hashtable properties,
 145           ColumnSelection cs, AlignmentI al, ViewDef view)
 146   {
 147     if (view != null)
 148     {
 149       if (view.viewname != null)
 150       {
 151         text.append("VIEW_DEF\t" + view.viewname + "\n");
 152       }
 153       if (list == null)
 154       {
 155         list = view.visibleGroups;
 156       }
 157       if (cs == null)
 158       {
 159         cs = view.hiddencols;
 160       }
 161       if (al == null)
 162       {
 163         // add hidden rep sequences.
 164       }
 165     }
 166     // first target - store and restore all settings for a view.
 167     if (al != null && al.hasSeqrep())
 168     {
 169       text.append("VIEW_SETREF\t" + al.getSeqrep().getName() + "\n");
 170     }
 171     if (cs != null && cs.hasHiddenColumns())
 172     {
 173       text.append("VIEW_HIDECOLS\t");
 174       List<int[]> hc = cs.getHiddenColumns();
 175       boolean comma = false;
 176       for (int[] r : hc)
 177       {
 178         if (!comma)
 179         {
 180           comma = true;
 181         }
 182         else
 183         {
 184           text.append(",");
 185         }
 186         text.append("" + r[0]);
 187         text.append("-");
 188         text.append("" + r[1]);
 189       }
 190       text.append("\n");
 191     }
 192     // TODO: allow efficient recovery of annotation data shown in several
 193     // different views
 194     if (annotations != null)
 195     {
 196       boolean oneColour = true;
 197       AlignmentAnnotation row;
 198       String comma;
 199       SequenceI refSeq = null;
 200       SequenceGroup refGroup = null;
 201
 202       StringBuffer colours = new StringBuffer();
 203       StringBuffer graphLine = new StringBuffer();
 204       StringBuffer rowprops = new StringBuffer();
 205       Hashtable<Integer, String> graphGroup = new Hashtable<Integer, String>();
 206       Hashtable<Integer, Object[]> graphGroup_refs = new Hashtable<Integer, Object[]>();
 207       BitSet graphGroupSeen = new BitSet();
 208
 209       java.awt.Color color;
 210
 211       for (int i = 0; i < annotations.length; i++)
 212       {
 213         row = annotations[i];
 214
 215         if (!row.visible
 216                 && !row.hasScore()
 217                 && !(row.graphGroup > -1 && graphGroupSeen
 218                         .get(row.graphGroup)))
 219         {
 220           continue;
 221         }
 222
 223         color = null;
 224         oneColour = true;
 225
 226         // mark any sequence references for the row
 227         writeSequence_Ref(refSeq, row.sequenceRef);
 228         refSeq = row.sequenceRef;
 229         // mark any group references for the row
 230         writeGroup_Ref(refGroup, row.groupRef);
 231         refGroup = row.groupRef;
 232
 233         boolean hasGlyphs = row.hasIcons, hasLabels = row.hasText, hasValues = row.hasScore, hasText = false;
 234         // lookahead to check what the annotation row object actually contains.
 235         for (int j = 0; row.annotations != null
 236                 && j < row.annotations.length
 237                 && (!hasGlyphs || !hasLabels || !hasValues); j++)
 238         {
 239           if (row.annotations[j] != null)
 240           {
 241             hasLabels |= (row.annotations[j].displayCharacter != null
 242                     && row.annotations[j].displayCharacter.length() > 0 && !row.annotations[j].displayCharacter
 243                     .equals(" "));
 244             hasGlyphs |= (row.annotations[j].secondaryStructure != 0 && row.annotations[j].secondaryStructure != ' ');
 245             hasValues |= (!Float.isNaN(row.annotations[j].value)); // NaNs can't  // BH
 246             // be
 247             // rendered..
 248             hasText |= (row.annotations[j].description != null && row.annotations[j].description
 249                     .length() > 0);
 250           }
 251         }
 252
 253         if (row.graph == AlignmentAnnotation.NO_GRAPH)
 254         {
 255           text.append("NO_GRAPH\t");
 256           hasValues = false; // only secondary structure
 257           // hasLabels = false; // and annotation description string.
 258         }
 259         else
 260         {
 261           if (row.graph == AlignmentAnnotation.BAR_GRAPH)
 262           {
 263             text.append("BAR_GRAPH\t");
 264             hasGlyphs = false; // no secondary structure
 265
 266           }
 267           else if (row.graph == AlignmentAnnotation.LINE_GRAPH)
 268           {
 269             hasGlyphs = false; // no secondary structure
 270             text.append("LINE_GRAPH\t");
 271           }
 272
 273           if (row.getThreshold() != null)
 274           {
 275             graphLine.append("GRAPHLINE\t");
 276             graphLine.append(row.label);
 277             graphLine.append("\t");
 278             graphLine.append(row.getThreshold().value);
 279             graphLine.append("\t");
 280             graphLine.append(row.getThreshold().label);
 281             graphLine.append("\t");
 282             graphLine
 283                     .append(Format.getHexString(row.getThreshold().colour));
 284             graphLine.append(newline);
 285           }
 286
 287           if (row.graphGroup > -1)
 288           {
 289             graphGroupSeen.set(row.graphGroup);
 290             Integer key = new Integer(row.graphGroup);
 291             if (graphGroup.containsKey(key))
 292             {
 293               graphGroup.put(key, graphGroup.get(key) + "\t" + row.label);
 294
 295             }
 296             else
 297             {
 298               graphGroup_refs.put(key, new Object[]
 299               { refSeq, refGroup });
 300               graphGroup.put(key, row.label);
 301             }
 302           }
 303         }
 304
 305         text.append(row.label + "\t");
 306         if (row.description != null)
 307         {
 308           text.append(row.description + "\t");
 309         }
 310         for (int j = 0; row.annotations != null
 311                 && j < row.annotations.length; j++)
 312         {
 313           if (refSeq != null
 314                   && Comparison.isGap(refSeq.getCharAt(j)))
 315           {
 316             continue;
 317           }
 318
 319           if (row.annotations[j] != null)
 320           {
 321             comma = "";
 322             if (hasGlyphs) // could be also hasGlyphs || ...
 323             {
 324
 325               text.append(comma);
 326               if (row.annotations[j].secondaryStructure != ' ')
 327               {
 328                 // only write out the field if its not whitespace.
 329                 text.append(row.annotations[j].secondaryStructure);
 330               }
 331               comma = ",";
 332             }
 333             if (hasValues)
 334             {
 335               if (!Float.isNaN(row.annotations[j].value))
 336               {
 337                 text.append(comma + row.annotations[j].value);
 338               }
 339               else
 340               {
 341                 System.err.println("Skipping NaN - not valid value.");
 342                 text.append(comma + 0f);// row.annotations[j].value);
 343               }
 344               comma = ",";
 345             }
 346             if (hasLabels)
 347             {
 348               // TODO: labels are emitted after values for bar graphs.
 349               if // empty labels are allowed, so
 350               (row.annotations[j].displayCharacter != null
 351                       && row.annotations[j].displayCharacter.length() > 0
 352                       && !row.annotations[j].displayCharacter.equals(" "))
 353               {
 354                 text.append(comma + row.annotations[j].displayCharacter);
 355                 comma = ",";
 356               }
 357             }
 358             if (hasText)
 359             {
 360               if (row.annotations[j].description != null
 361                       && row.annotations[j].description.length() > 0
 362                       && !row.annotations[j].description
 363                               .equals(row.annotations[j].displayCharacter))
 364               {
 365                 text.append(comma + row.annotations[j].description);
 366                 comma = ",";
 367               }
 368             }
 369             if (color != null && !color.equals(row.annotations[j].colour))
 370             {
 371               oneColour = false;
 372             }
 373
 374             color = row.annotations[j].colour;
 375
 376             if (row.annotations[j].colour != null
 377                     && row.annotations[j].colour != java.awt.Color.black)
 378             {
 379               text.append(comma
 380                       + "["
 381                       + Format
 382                               .getHexString(row.annotations[j].colour)
 383                       + "]");
 384               comma = ",";
 385             }
 386           }
 387           text.append("|");
 388         }
 389
 390         if (row.hasScore())
 391         {
 392           text.append("\t" + row.score);
 393         }
 394
 395         text.append(newline);
 396
 397         if (color != null && color != java.awt.Color.black && oneColour)
 398         {
 399           colours.append("COLOUR\t");
 400           colours.append(row.label);
 401           colours.append("\t");
 402           colours.append(Format.getHexString(color));
 403           colours.append(newline);
 404         }
 405         if (row.scaleColLabel || row.showAllColLabels
 406                 || row.centreColLabels)
 407         {
 408           rowprops.append("ROWPROPERTIES\t");
 409           rowprops.append(row.label);
 410           rowprops.append("\tscaletofit=");
 411           rowprops.append(row.scaleColLabel);
 412           rowprops.append("\tshowalllabs=");
 413           rowprops.append(row.showAllColLabels);
 414           rowprops.append("\tcentrelabs=");
 415           rowprops.append(row.centreColLabels);
 416           rowprops.append(newline);
 417         }
 418         if (graphLine.length() > 0)
 419         {
 420           text.append(graphLine.toString());
 421           graphLine.setLength(0);
 422         }
 423       }
 424
 425       text.append(newline);
 426
 427       text.append(colours.toString());
 428       if (graphGroup.size() > 0)
 429       {
 430         SequenceI oldRefSeq = refSeq;
 431         SequenceGroup oldRefGroup = refGroup;
 432         for (Map.Entry<Integer, String> combine_statement : graphGroup
 433                 .entrySet())
 434         {
 435           Object[] seqRefAndGroup = graphGroup_refs.get(combine_statement
 436                   .getKey());
 437
 438           writeSequence_Ref(refSeq, (SequenceI) seqRefAndGroup[0]);
 439           refSeq = (SequenceI) seqRefAndGroup[0];
 440
 441           writeGroup_Ref(refGroup, (SequenceGroup) seqRefAndGroup[1]);
 442           refGroup = (SequenceGroup) seqRefAndGroup[1];
 443           text.append("COMBINE\t");
 444           text.append(combine_statement.getValue());
 445           text.append(newline);
 446         }
 447         writeSequence_Ref(refSeq, oldRefSeq);
 448         refSeq = oldRefSeq;
 449
 450         writeGroup_Ref(refGroup, oldRefGroup);
 451         refGroup = oldRefGroup;
 452       }
 453       text.append(rowprops.toString());
 454     }
 455
 456     if (list != null)
 457     {
 458       printGroups(list);
 459     }
 460
 461     if (properties != null)
 462     {
 463       text.append(newline);
 464       text.append(newline);
 465       text.append("ALIGNMENT");
 466       Enumeration en = properties.keys();
 467       while (en.hasMoreElements())
 468       {
 469         String key = en.nextElement().toString();
 470         text.append("\t");
 471         text.append(key);
 472         text.append("=");
 473         text.append(properties.get(key));
 474       }
 475       // TODO: output alignment visualization settings here if required
 476       // iterate through one or more views, defining, marking columns and rows as visible/hidden, and emmitting view properties.
 477       // View specific annotation is
 478     }
 479
 480     return text.toString();
 481   }
 482
 483   private Object writeGroup_Ref(SequenceGroup refGroup,
 484           SequenceGroup next_refGroup)
 485   {
 486     if (next_refGroup == null)
 487     {
 488
 489       if (refGroup != null)
 490       {
 491         text.append(newline);
 492         text.append("GROUP_REF\t");
 493         text.append("ALIGNMENT");
 494         text.append(newline);
 495       }
 496       return true;
 497     }
 498     else
 499     {
 500       if (refGroup == null || refGroup != next_refGroup)
 501       {
 502         text.append(newline);
 503         text.append("GROUP_REF\t");
 504         text.append(next_refGroup.getName());
 505         text.append(newline);
 506         return true;
 507       }
 508     }
 509     return false;
 510   }
 511
 512   private boolean writeSequence_Ref(SequenceI refSeq, SequenceI next_refSeq)
 513   {
 514
 515     if (next_refSeq == null)
 516     {
 517       if (refSeq != null)
 518       {
 519         text.append(newline);
 520         text.append("SEQUENCE_REF\t");
 521         text.append("ALIGNMENT");
 522         text.append(newline);
 523         return true;
 524       }
 525     }
 526     else
 527     {
 528       if (refSeq == null || refSeq != next_refSeq)
 529       {
 530         text.append(newline);
 531         text.append("SEQUENCE_REF\t");
 532         text.append(next_refSeq.getName());
 533         text.append(newline);
 534         return true;
 535       }
 536     }
 537     return false;
 538   }
 539
 540   public void printGroups(List<SequenceGroup> list)
 541   {
 542     SequenceI seqrep = null;
 543     for (SequenceGroup sg : list)
 544     {
 545       if (!sg.hasSeqrep())
 546       {
 547         text.append("SEQUENCE_GROUP\t" + sg.getName() + "\t"
 548                 + (sg.getStartRes() + 1) + "\t" + (sg.getEndRes() + 1)
 549                 + "\t" + "-1\t");
 550         seqrep = null;
 551       }
 552       else
 553       {
 554         seqrep = sg.getSeqrep();
 555         text.append("SEQUENCE_REF\t");
 556         text.append(seqrep.getName());
 557         text.append(newline);
 558         text.append("SEQUENCE_GROUP\t");
 559         text.append(sg.getName());
 560         text.append("\t");
 561         text.append("" + (seqrep.findPosition(sg.getStartRes())));
 562         text.append("\t");
 563         text.append((seqrep.findPosition(sg.getEndRes())));
 564         text.append("\t");
 565         text.append("-1\t");
 566       }
 567       for (int s = 0; s < sg.getSize(); s++)
 568       {
 569         text.append(sg.getSequenceAt(s).getName());
 570         text.append("\t");
 571       }
 572       text.append(newline);
 573       text.append("PROPERTIES\t");
 574       text.append(sg.getName());
 575       text.append("\t");
 576
 577       if (sg.getDescription() != null)
 578       {
 579         text.append("description=");
 580         text.append(sg.getDescription());
 581         text.append("\t");
 582       }
 583       if (sg.cs != null)
 584       {
 585         text.append("colour=");
 586         text.append(ColourSchemeProperty.getColourName(sg.cs));
 587         text.append("\t");
 588         if (sg.cs.getThreshold() != 0)
 589         {
 590           text.append("pidThreshold=");
 591           text.append(sg.cs.getThreshold());
 592         }
 593         if (sg.cs.conservationApplied())
 594         {
 595           text.append("consThreshold=");
 596           text.append(sg.cs.getConservationInc());
 597           text.append("\t");
 598         }
 599       }
 600       text.append("outlineColour=");
 601       text.append(Format.getHexString(sg.getOutlineColour()));
 602       text.append("\t");
 603
 604       text.append("displayBoxes=");
 605       text.append(sg.getDisplayBoxes());
 606       text.append("\t");
 607       text.append("displayText=");
 608       text.append(sg.getDisplayText());
 609       text.append("\t");
 610       text.append("colourText=");
 611       text.append(sg.getColourText());
 612       text.append("\t");
 613       text.append("showUnconserved=");
 614       text.append(sg.getShowNonconserved());
 615       text.append("\t");
 616       if (sg.textColour != java.awt.Color.black)
 617       {
 618         text.append("textCol1=");
 619         text.append(Format.getHexString(sg.textColour));
 620         text.append("\t");
 621       }
 622       if (sg.textColour2 != java.awt.Color.white)
 623       {
 624         text.append("textCol2=");
 625         text.append(Format.getHexString(sg.textColour2));
 626         text.append("\t");
 627       }
 628       if (sg.thresholdTextColour != 0)
 629       {
 630         text.append("textColThreshold=");
 631         text.append(sg.thresholdTextColour);
 632         text.append("\t");
 633       }
 634       if (sg.idColour != null)
 635       {
 636         text.append("idColour=");
 637         text.append(Format.getHexString(sg.idColour));
 638         text.append("\t");
 639       }
 640       if (sg.isHidereps())
 641       {
 642         text.append("hide=true\t");
 643       }
 644       if (sg.isHideCols())
 645       {
 646         text.append("hidecols=true\t");
 647       }
 648       if (seqrep != null)
 649       {
 650         // terminate the last line and clear the sequence ref for the group
 651         text.append(newline);
 652         text.append("SEQUENCE_REF");
 653       }
 654       text.append(newline);
 655       text.append(newline);
 656
 657     }
 658   }
 659
 660   SequenceI refSeq = null;
 661
 662   String refSeqId = null;
 663
 664   public boolean annotateAlignmentView(AlignViewportI viewport,
 665           String file, String protocol)
 666   {
 667     ColumnSelection colSel = viewport.getColumnSelection();
 668     if (colSel == null)
 669     {
 670       colSel = new ColumnSelection();
 671     }
 672     boolean rslt = readAnnotationFile(viewport.getAlignment(), colSel,
 673             file, protocol);
 674     if (rslt
 675             && (colSel.hasSelectedColumns() || colSel.hasHiddenColumns()))
 676     {
 677       viewport.setColumnSelection(colSel);
 678     }
 679
 680     return rslt;
 681   }
 682   public boolean readAnnotationFile(AlignmentI al, String file,
 683           String protocol)
 684   {
 685     return readAnnotationFile(al, null, file, protocol);
 686   }
 687
 688   public boolean readAnnotationFile(AlignmentI al, ColumnSelection colSel,
 689           String file, String protocol)
 690   {
 691     BufferedReader in = null;
 692     try
 693     {
 694       if (protocol.equals(AppletFormatAdapter.FILE))
 695       {
 696         in = GenericFileAdapter.getReader(file, false);
 697       }
 698       else if (protocol.equals(AppletFormatAdapter.URL))
 699       {
 700         in = GenericFileAdapter.getReader(file, true);
 701       }
 702       else if (protocol.equals(AppletFormatAdapter.PASTE))
 703       {
 704         in = new BufferedReader(new StringReader(file));
 705       }
 706       else if (protocol.equals(AppletFormatAdapter.CLASSLOADER))
 707       {
 708         java.io.InputStream is = getClass().getResourceAsStream("/" + file);
 709         if (is != null)
 710         {
 711           in = new BufferedReader(new InputStreamReader(is));
 712         }
 713       }
 714       if (in != null)
 715       {
 716         return parseAnnotationFrom(al, colSel, in);
 717       }
 718
 719     } catch (Exception ex)
 720     {
 721       ex.printStackTrace();
 722       System.out.println("Problem reading annotation file: " + ex);
 723       if (nlinesread > 0)
 724       {
 725         System.out.println("Last read line " + nlinesread + ": '"
 726                 + lastread + "' (first 80 chars) ...");
 727       }
 728       return false;
 729     }
 730     return false;
 731   }
 732
 733   long nlinesread = 0;
 734
 735   String lastread = "";
 736
 737   private static String GRAPHLINE = "GRAPHLINE", COMBINE = "COMBINE";
 738
 739   public boolean parseAnnotationFrom(AlignmentI al, ColumnSelection colSel,
 740           BufferedReader in)
 741           throws Exception
 742   {
 743     nlinesread = 0;
 744     ArrayList<Object[]> combineAnnotation_calls = new ArrayList<Object[]>();
 745     ArrayList<Object[]> deferredAnnotation_calls = new ArrayList<Object[]>();
 746     boolean modified = false;
 747     String groupRef = null;
 748     Hashtable groupRefRows = new Hashtable();
 749
 750     Hashtable autoAnnots = new Hashtable();
 751     {
 752       String line, label, description, token;
 753       int graphStyle, index;
 754       int refSeqIndex = 1;
 755       int existingAnnotations = 0;
 756       // when true - will add new rows regardless of whether they are duplicate
 757       // auto-annotation like consensus or conservation graphs
 758       boolean overrideAutoAnnot = false;
 759       if (al.getAlignmentAnnotation() != null)
 760       {
 761         existingAnnotations = al.getAlignmentAnnotation().length;
 762         if (existingAnnotations > 0)
 763         {
 764           AlignmentAnnotation[] aa = al.getAlignmentAnnotation();
 765           for (int aai = 0; aai < aa.length; aai++)
 766           {
 767             if (aa[aai].autoCalculated)
 768             {
 769               // make a note of the name and description
 770               autoAnnots.put(
 771                       autoAnnotsKey(aa[aai], aa[aai].sequenceRef,
 772                               (aa[aai].groupRef == null ? null
 773                                       : aa[aai].groupRef.getName())),
 774                       new Integer(1));
 775             }
 776           }
 777         }
 778       }
 779
 780       int alWidth = al.getWidth();
 781
 782       StringTokenizer st;
 783       Annotation[] annotations;
 784       AlignmentAnnotation annotation = null;
 785
 786       // First confirm this is an Annotation file
 787       boolean jvAnnotationFile = false;
 788       while ((line = in.readLine()) != null)
 789       {
 790         nlinesread++;
 791         lastread = new String(line);
 792         if (line.indexOf("#") == 0)
 793         {
 794           continue;
 795         }
 796
 797         if (line.indexOf("JALVIEW_ANNOTATION") > -1)
 798         {
 799           jvAnnotationFile = true;
 800           break;
 801         }
 802       }
 803
 804       if (!jvAnnotationFile)
 805       {
 806         in.close();
 807         return false;
 808       }
 809
 810       while ((line = in.readLine()) != null)
 811       {
 812         nlinesread++;
 813         lastread = new String(line);
 814         if (line.indexOf("#") == 0
 815                 || line.indexOf("JALVIEW_ANNOTATION") > -1
 816                 || line.length() == 0)
 817         {
 818           continue;
 819         }
 820
 821         st = new StringTokenizer(line, "\t");
 822         token = st.nextToken();
 823         if (token.equalsIgnoreCase("COLOUR"))
 824         {
 825           // TODO: use graduated colour def'n here too
 826           colourAnnotations(al, st.nextToken(), st.nextToken());
 827           modified = true;
 828           continue;
 829         }
 830
 831         else if (token.equalsIgnoreCase(COMBINE))
 832         {
 833           // keep a record of current state and resolve groupRef at end
 834           combineAnnotation_calls.add(new Object[]
 835           { st, refSeq, groupRef });
 836           modified = true;
 837           continue;
 838         }
 839         else if (token.equalsIgnoreCase("ROWPROPERTIES"))
 840         {
 841           addRowProperties(al, st);
 842           modified = true;
 843           continue;
 844         }
 845         else if (token.equalsIgnoreCase(GRAPHLINE))
 846         {
 847           // resolve at end
 848           deferredAnnotation_calls.add(new Object[]
 849           { GRAPHLINE, st, refSeq, groupRef });
 850           modified = true;
 851           continue;
 852         }
 853
 854         else if (token.equalsIgnoreCase("SEQUENCE_REF"))
 855         {
 856           if (st.hasMoreTokens())
 857           {
 858             refSeq = al.findName(refSeqId = st.nextToken());
 859             if (refSeq == null)
 860             {
 861               refSeqId = null;
 862             }
 863             try
 864             {
 865               refSeqIndex = Integer.parseInt(st.nextToken());
 866               if (refSeqIndex < 1)
 867               {
 868                 refSeqIndex = 1;
 869                 System.out
 870                         .println("WARNING: SEQUENCE_REF index must be > 0 in AnnotationFile");
 871               }
 872             } catch (Exception ex)
 873             {
 874               refSeqIndex = 1;
 875             }
 876           }
 877           else
 878           {
 879             refSeq = null;
 880             refSeqId = null;
 881           }
 882           continue;
 883         }
 884         else if (token.equalsIgnoreCase("GROUP_REF"))
 885         {
 886           // Group references could be forward or backwards, so they are
 887           // resolved after the whole file is read in
 888           groupRef = null;
 889           if (st.hasMoreTokens())
 890           {
 891             groupRef = st.nextToken();
 892             if (groupRef.length() < 1)
 893             {
 894               groupRef = null; // empty string
 895             }
 896             else
 897             {
 898               if (groupRefRows.get(groupRef) == null)
 899               {
 900                 groupRefRows.put(groupRef, new Vector());
 901               }
 902             }
 903           }
 904           continue;
 905         }
 906         else if (token.equalsIgnoreCase("SEQUENCE_GROUP"))
 907         {
 908           addGroup(al, st);
 909           modified = true;
 910           continue;
 911         }
 912
 913         else if (token.equalsIgnoreCase("PROPERTIES"))
 914         {
 915           addProperties(al, st);
 916           modified = true;
 917           continue;
 918         }
 919
 920         else if (token.equalsIgnoreCase("BELOW_ALIGNMENT"))
 921         {
 922           setBelowAlignment(al, st);
 923           modified = true;
 924           continue;
 925         }
 926         else if (token.equalsIgnoreCase("ALIGNMENT"))
 927         {
 928           addAlignmentDetails(al, st);
 929           modified = true;
 930           continue;
 931         }
 932         // else if (token.equalsIgnoreCase("VIEW_DEF"))
 933         // {
 934         // addOrSetView(al,st);
 935         // modified = true;
 936         // continue;
 937         // }
 938         else if (token.equalsIgnoreCase("VIEW_SETREF"))
 939         {
 940           if (refSeq != null)
 941           {
 942             al.setSeqrep(refSeq);
 943           }
 944           modified = true;
 945           continue;
 946         }
 947         else if (token.equalsIgnoreCase("VIEW_HIDECOLS"))
 948         {
 949           if (st.hasMoreTokens())
 950           {
 951             if (colSel == null)
 952             {
 953               colSel = new ColumnSelection();
 954             }
 955             parseHideCols(colSel, st.nextToken());
 956           }
 957           modified = true;
 958           continue;
 959         }
 960         else if (token.equalsIgnoreCase("HIDE_INSERTIONS"))
 961         {
 962           SequenceI sr = refSeq == null ? al.getSeqrep() : refSeq;
 963           if (sr == null)
 964           {
 965             sr = al.getSequenceAt(0);
 966           }
 967           if (sr != null)
 968           {
 969             if (colSel == null)
 970             {
 971               System.err
 972                       .println("Cannot process HIDE_INSERTIONS without an alignment view: Ignoring line: "
 973                               + line);
 974             }
 975             else
 976             {
 977               // consider deferring this till after the file has been parsed ?
 978               colSel.hideInsertionsFor(sr);
 979             }
 980           }
 981           modified = true;
 982           continue;
 983         }
 984
 985         // Parse out the annotation row
 986         graphStyle = AlignmentAnnotation.getGraphValueFromString(token);
 987         label = st.nextToken();
 988
 989         index = 0;
 990         annotations = new Annotation[alWidth];
 991         description = null;
 992         float score = Float.NaN;
 993
 994         if (st.hasMoreTokens())
 995         {
 996           line = st.nextToken();
 997
 998           if (line.indexOf("|") == -1)
 999           {
1000             description = line;
1001             if (st.hasMoreTokens())
1002             {
1003               line = st.nextToken();
1004             }
1005           }
1006
1007           if (st.hasMoreTokens())
1008           {
1009             // This must be the score
1010             score = Float.valueOf(st.nextToken()).floatValue();
1011           }
1012
1013           st = new StringTokenizer(line, "|", true);
1014
1015           boolean emptyColumn = true;
1016           boolean onlyOneElement = (st.countTokens() == 1);
1017
1018           while (st.hasMoreElements() && index < alWidth)
1019           {
1020             token = st.nextToken().trim();
1021
1022             if (onlyOneElement)
1023             {
1024               try
1025               {
1026                 score = Float.valueOf(token).floatValue();
1027                 break;
1028               } catch (NumberFormatException ex)
1029               {
1030               }
1031             }
1032
1033             if (token.equals("|"))
1034             {
1035               if (emptyColumn)
1036               {
1037                 index++;
1038               }
1039
1040               emptyColumn = true;
1041             }
1042             else
1043             {
1044               annotations[index++] = parseAnnotation(token, graphStyle);
1045               emptyColumn = false;
1046             }
1047           }
1048
1049         }
1050
1051         annotation = new AlignmentAnnotation(label, description,
1052                 (index == 0) ? null : annotations, 0, 0, graphStyle);
1053
1054         annotation.score = score;
1055         if (!overrideAutoAnnot
1056                 && autoAnnots.containsKey(autoAnnotsKey(annotation, refSeq,
1057                         groupRef)))
1058         {
1059           // skip - we've already got an automatic annotation of this type.
1060           continue;
1061         }
1062         // otherwise add it!
1063         if (refSeq != null)
1064         {
1065
1066           annotation.belowAlignment = false;
1067           // make a copy of refSeq so we can find other matches in the alignment
1068           SequenceI referedSeq = refSeq;
1069           do
1070           {
1071             // copy before we do any mapping business.
1072             // TODO: verify that undo/redo with 1:many sequence associated
1073             // annotations can be undone correctly
1074             AlignmentAnnotation ann = new AlignmentAnnotation(annotation);
1075             annotation
1076                     .createSequenceMapping(referedSeq, refSeqIndex, false);
1077             annotation.adjustForAlignment();
1078             referedSeq.addAlignmentAnnotation(annotation);
1079             al.addAnnotation(annotation);
1080             al.setAnnotationIndex(annotation,
1081                     al.getAlignmentAnnotation().length
1082                             - existingAnnotations - 1);
1083             if (groupRef != null)
1084             {
1085               ((Vector) groupRefRows.get(groupRef)).addElement(annotation);
1086             }
1087             // and recover our virgin copy to use again if necessary.
1088             annotation = ann;
1089
1090           } while (refSeqId != null
1091                   && (referedSeq = al.findName(referedSeq, refSeqId, true)) != null);
1092         }
1093         else
1094         {
1095           al.addAnnotation(annotation);
1096           al.setAnnotationIndex(annotation,
1097                   al.getAlignmentAnnotation().length - existingAnnotations
1098                           - 1);
1099           if (groupRef != null)
1100           {
1101             ((Vector) groupRefRows.get(groupRef)).addElement(annotation);
1102           }
1103         }
1104         // and set modification flag
1105         modified = true;
1106       }
1107       // Resolve the groupRefs
1108       Hashtable<String, SequenceGroup> groupRefLookup = new Hashtable<String, SequenceGroup>();
1109       Enumeration en = groupRefRows.keys();
1110
1111       while (en.hasMoreElements())
1112       {
1113         groupRef = (String) en.nextElement();
1114         boolean matched = false;
1115         // Resolve group: TODO: add a getGroupByName method to alignments
1116         for (SequenceGroup theGroup : al.getGroups())
1117         {
1118           if (theGroup.getName().equals(groupRef))
1119           {
1120             if (matched)
1121             {
1122               // TODO: specify and implement duplication of alignment annotation
1123               // for multiple group references.
1124               System.err
1125                       .println("Ignoring 1:many group reference mappings for group name '"
1126                               + groupRef + "'");
1127             }
1128             else
1129             {
1130               matched = true;
1131               Vector rowset = (Vector) groupRefRows.get(groupRef);
1132               groupRefLookup.put(groupRef, theGroup);
1133               if (rowset != null && rowset.size() > 0)
1134               {
1135                 AlignmentAnnotation alan = null;
1136                 for (int elm = 0, elmSize = rowset.size(); elm < elmSize; elm++)
1137                 {
1138                   alan = (AlignmentAnnotation) rowset.elementAt(elm);
1139                   alan.groupRef = theGroup;
1140                 }
1141               }
1142             }
1143           }
1144         }
1145         ((Vector) groupRefRows.get(groupRef)).removeAllElements();
1146       }
1147       // process any deferred attribute settings for each context
1148       for (Object[] _deferred_args : deferredAnnotation_calls)
1149       {
1150         if (_deferred_args[0] == GRAPHLINE)
1151         {
1152           addLine(al,
1153                   (StringTokenizer) _deferred_args[1], // st
1154                   (SequenceI) _deferred_args[2], // refSeq
1155                   (_deferred_args[3] == null) ? null : groupRefLookup
1156                           .get(_deferred_args[3]) // the reference
1157                                                            // group, or null
1158           );
1159         }
1160       }
1161
1162       // finally, combine all the annotation rows within each context.
1163       /**
1164        * number of combine statements in this annotation file. Used to create
1165        * new groups for combined annotation graphs without disturbing existing
1166        * ones
1167        */
1168       int combinecount = 0;
1169       for (Object[] _combine_args : combineAnnotation_calls)
1170       {
1171         combineAnnotations(al,
1172                 ++combinecount,
1173                 (StringTokenizer) _combine_args[0], // st
1174                 (SequenceI) _combine_args[1], // refSeq
1175                 (_combine_args[2] == null) ? null : groupRefLookup
1176                         .get(_combine_args[2]) // the reference group,
1177                                                         // or null
1178         );
1179       }
1180     }
1181     return modified;
1182   }
1183
1184   private void parseHideCols(ColumnSelection colSel, String nextToken)
1185   {
1186     StringTokenizer inval = new StringTokenizer(nextToken,",");
1187     while (inval.hasMoreTokens())
1188     {
1189       String range = inval.nextToken().trim();
1190       int from, to = range.indexOf("-");
1191       if (to == -1)
1192       {
1193         from = to = Integer.parseInt(range);
1194         if (from >= 0)
1195         {
1196           colSel.hideColumns(from, to);
1197         }
1198       }
1199       else
1200       {
1201         from = Integer.parseInt(range.substring(0, to));
1202         if (to < range.length() - 1)
1203         {
1204           to = Integer.parseInt(range.substring(to + 1));
1205         }
1206         else
1207         {
1208           to = from;
1209         }
1210         if (from > 0 && to >= from)
1211         {
1212           colSel.hideColumns(from, to);
1213         }
1214       }
1215     }
1216   }
1217
1218   private Object autoAnnotsKey(AlignmentAnnotation annotation,
1219           SequenceI refSeq, String groupRef)
1220   {
1221     return annotation.graph + "\t" + annotation.label + "\t"
1222             + annotation.description + "\t"
1223             + (refSeq != null ? refSeq.getDisplayId(true) : "");
1224   }
1225
1226   Annotation parseAnnotation(String string, int graphStyle)
1227   {
1228     boolean hasSymbols = (graphStyle == AlignmentAnnotation.NO_GRAPH); // don't
1229     // do the
1230     // glyph
1231     // test
1232     // if we
1233     // don't
1234     // want
1235     // secondary
1236     // structure
1237     String desc = null, displayChar = null;
1238     char ss = ' '; // secondaryStructure
1239     float value = 0;
1240     boolean parsedValue = false, dcset = false;
1241
1242     // find colour here
1243     java.awt.Color colour = null;
1244     int i = string.indexOf("[");
1245     int j = string.indexOf("]");
1246     if (i > -1 && j > -1)
1247     {
1248       UserColourScheme ucs = new UserColourScheme();
1249
1250       colour = ucs.getColourFromString(string.substring(i + 1, j));
1251       if (i > 0 && string.charAt(i - 1) == ',')
1252       {
1253         // clip the preceding comma as well
1254         i--;
1255       }
1256       string = string.substring(0, i) + string.substring(j + 1);
1257     }
1258
1259     StringTokenizer st = new StringTokenizer(string, ",", true);
1260     String token;
1261     boolean seenContent = false;
1262     int pass = 0;
1263     while (st.hasMoreTokens())
1264     {
1265       pass++;
1266       token = st.nextToken().trim();
1267       if (token.equals(","))
1268       {
1269         if (!seenContent && parsedValue && !dcset)
1270         {
1271           // allow the value below the bar/line to be empty
1272           dcset = true;
1273           displayChar = " ";
1274         }
1275         seenContent = false;
1276         continue;
1277       }
1278       else
1279       {
1280         seenContent = true;
1281       }
1282
1283       if (!parsedValue)
1284       {
1285         try
1286         {
1287           displayChar = token;
1288           // foo
1289           value = new Float(token).floatValue();
1290           parsedValue = true;
1291           continue;
1292         } catch (NumberFormatException ex)
1293         {
1294         }
1295       }
1296       else
1297       {
1298         if (token.length() == 1)
1299         {
1300           displayChar = token;
1301         }
1302       }
1303       if (hasSymbols
1304               && (token.equals("H") || token.equals("E")
1305                       || token.equals("S") || token.equals(" ")))
1306       {
1307         // Either this character represents a helix or sheet
1308         // or an integer which can be displayed
1309         ss = token.charAt(0);
1310         if (displayChar.equals(token.substring(0, 1)))
1311         {
1312           displayChar = "";
1313         }
1314       }
1315       else if (desc == null || (parsedValue && pass > 2))
1316       {
1317         desc = token;
1318       }
1319
1320     }
1321     // if (!dcset && string.charAt(string.length() - 1) == ',')
1322     // {
1323     // displayChar = " "; // empty display char symbol.
1324     // }
1325     if (displayChar != null && desc != null && desc.length() == 1)
1326     {
1327       if (displayChar.length() > 1)
1328       {
1329         // switch desc and displayChar - legacy support
1330         String tmp = displayChar;
1331         displayChar = desc;
1332         desc = tmp;
1333       }
1334       else
1335       {
1336         if (displayChar.equals(desc))
1337         {
1338           // duplicate label - hangover from the 'robust parser' above
1339           desc = null;
1340         }
1341       }
1342     }
1343     Annotation anot = new Annotation(displayChar, desc, ss, value);
1344
1345     anot.colour = colour;
1346
1347     return anot;
1348   }
1349
1350   void colourAnnotations(AlignmentI al, String label, String colour)
1351   {
1352     UserColourScheme ucs = new UserColourScheme(colour);
1353     Annotation[] annotations;
1354     for (int i = 0; i < al.getAlignmentAnnotation().length; i++)
1355     {
1356       if (al.getAlignmentAnnotation()[i].label.equalsIgnoreCase(label))
1357       {
1358         annotations = al.getAlignmentAnnotation()[i].annotations;
1359         for (int j = 0; j < annotations.length; j++)
1360         {
1361           if (annotations[j] != null)
1362           {
1363             annotations[j].colour = ucs.findColour('A');
1364           }
1365         }
1366       }
1367     }
1368   }
1369
1370   void combineAnnotations(AlignmentI al, int combineCount,
1371           StringTokenizer st, SequenceI seqRef, SequenceGroup groupRef)
1372   {
1373     String group = st.nextToken();
1374     // First make sure we are not overwriting the graphIndex
1375     int graphGroup = 0;
1376     if (al.getAlignmentAnnotation() != null)
1377     {
1378       for (int i = 0; i < al.getAlignmentAnnotation().length; i++)
1379       {
1380         AlignmentAnnotation aa = al.getAlignmentAnnotation()[i];
1381
1382         if (aa.graphGroup > graphGroup)
1383         {
1384           // try to number graphGroups in order of occurence.
1385           graphGroup = aa.graphGroup + 1;
1386         }
1387         if (aa.sequenceRef == seqRef && aa.groupRef == groupRef
1388                 && aa.label.equalsIgnoreCase(group))
1389         {
1390           if (aa.graphGroup > -1)
1391           {
1392             graphGroup = aa.graphGroup;
1393           }
1394           else
1395           {
1396             if (graphGroup <= combineCount)
1397             {
1398               graphGroup = combineCount + 1;
1399             }
1400             aa.graphGroup = graphGroup;
1401           }
1402           break;
1403         }
1404       }
1405
1406       // Now update groups
1407       while (st.hasMoreTokens())
1408       {
1409         group = st.nextToken();
1410         for (int i = 0; i < al.getAlignmentAnnotation().length; i++)
1411         {
1412           AlignmentAnnotation aa = al.getAlignmentAnnotation()[i];
1413           if (aa.sequenceRef == seqRef && aa.groupRef == groupRef
1414                   && aa.label.equalsIgnoreCase(group))
1415           {
1416             aa.graphGroup = graphGroup;
1417             break;
1418           }
1419         }
1420       }
1421     }
1422     else
1423     {
1424       System.err
1425               .println("Couldn't combine annotations. None are added to alignment yet!");
1426     }
1427   }
1428
1429   void addLine(AlignmentI al, StringTokenizer st, SequenceI seqRef,
1430           SequenceGroup groupRef)
1431   {
1432     String group = st.nextToken();
1433     AlignmentAnnotation annotation = null, alannot[] = al
1434             .getAlignmentAnnotation();
1435     float value = new Float(st.nextToken()).floatValue();
1436     String label = st.hasMoreTokens() ? st.nextToken() : null;
1437     java.awt.Color colour = null;
1438     if (st.hasMoreTokens())
1439     {
1440       UserColourScheme ucs = new UserColourScheme(st.nextToken());
1441       colour = ucs.findColour('A');
1442     }
1443     if (alannot != null)
1444     {
1445       for (int i = 0; i < alannot.length; i++)
1446       {
1447         if (alannot[i].label.equalsIgnoreCase(group)
1448                 && (seqRef == null || alannot[i].sequenceRef == seqRef)
1449                 && (groupRef == null || alannot[i].groupRef == groupRef))
1450         {
1451           alannot[i].setThreshold(new GraphLine(value, label, colour));
1452         }
1453       }
1454     }
1455     if (annotation == null)
1456     {
1457       return;
1458     }
1459   }
1460
1461   void addGroup(AlignmentI al, StringTokenizer st)
1462   {
1463     SequenceGroup sg = new SequenceGroup();
1464     sg.setName(st.nextToken());
1465     String rng = "";
1466     try
1467     {
1468       rng = st.nextToken();
1469       if (rng.length() > 0 && !rng.startsWith("*"))
1470       {
1471         sg.setStartRes(Integer.parseInt(rng) - 1);
1472       }
1473       else
1474       {
1475         sg.setStartRes(0);
1476       }
1477       rng = st.nextToken();
1478       if (rng.length() > 0 && !rng.startsWith("*"))
1479       {
1480         sg.setEndRes(Integer.parseInt(rng) - 1);
1481       }
1482       else
1483       {
1484         sg.setEndRes(al.getWidth() - 1);
1485       }
1486     } catch (Exception e)
1487     {
1488       System.err
1489               .println("Couldn't parse Group Start or End Field as '*' or a valid column or sequence index: '"
1490                       + rng + "' - assuming alignment width for group.");
1491       // assume group is full width
1492       sg.setStartRes(0);
1493       sg.setEndRes(al.getWidth() - 1);
1494     }
1495
1496     String index = st.nextToken();
1497     if (index.equals("-1"))
1498     {
1499       while (st.hasMoreElements())
1500       {
1501         sg.addSequence(al.findName(st.nextToken()), false);
1502       }
1503     }
1504     else
1505     {
1506       StringTokenizer st2 = new StringTokenizer(index, ",");
1507
1508       while (st2.hasMoreTokens())
1509       {
1510         String tmp = st2.nextToken();
1511         if (tmp.equals("*"))
1512         {
1513           for (int i = 0; i < al.getHeight(); i++)
1514           {
1515             sg.addSequence(al.getSequenceAt(i), false);
1516           }
1517         }
1518         else if (tmp.indexOf("-") >= 0)
1519         {
1520           StringTokenizer st3 = new StringTokenizer(tmp, "-");
1521
1522           int start = (Integer.parseInt(st3.nextToken()));
1523           int end = (Integer.parseInt(st3.nextToken()));
1524
1525           if (end > start)
1526           {
1527             for (int i = start; i <= end; i++)
1528             {
1529               sg.addSequence(al.getSequenceAt(i - 1), false);
1530             }
1531           }
1532         }
1533         else
1534         {
1535           sg.addSequence(al.getSequenceAt(Integer.parseInt(tmp) - 1), false);
1536         }
1537       }
1538     }
1539
1540     if (refSeq != null)
1541     {
1542       sg.setStartRes(refSeq.findIndex(sg.getStartRes() + 1) - 1);
1543       sg.setEndRes(refSeq.findIndex(sg.getEndRes() + 1) - 1);
1544       sg.setSeqrep(refSeq);
1545     }
1546
1547     if (sg.getSize() > 0)
1548     {
1549       al.addGroup(sg);
1550     }
1551   }
1552
1553   void addRowProperties(AlignmentI al, StringTokenizer st)
1554   {
1555     String label = st.nextToken(), keyValue, key, value;
1556     boolean scaletofit = false, centerlab = false, showalllabs = false;
1557     while (st.hasMoreTokens())
1558     {
1559       keyValue = st.nextToken();
1560       key = keyValue.substring(0, keyValue.indexOf("="));
1561       value = keyValue.substring(keyValue.indexOf("=") + 1);
1562       if (key.equalsIgnoreCase("scaletofit"))
1563       {
1564         scaletofit = Boolean.valueOf(value).booleanValue();
1565       }
1566       if (key.equalsIgnoreCase("showalllabs"))
1567       {
1568         showalllabs = Boolean.valueOf(value).booleanValue();
1569       }
1570       if (key.equalsIgnoreCase("centrelabs"))
1571       {
1572         centerlab = Boolean.valueOf(value).booleanValue();
1573       }
1574       AlignmentAnnotation[] alr = al.getAlignmentAnnotation();
1575       if (alr != null)
1576       {
1577         for (int i = 0; i < alr.length; i++)
1578         {
1579           if (alr[i].label.equalsIgnoreCase(label))
1580           {
1581             alr[i].centreColLabels = centerlab;
1582             alr[i].scaleColLabel = scaletofit;
1583             alr[i].showAllColLabels = showalllabs;
1584           }
1585         }
1586       }
1587     }
1588   }
1589
1590   void addProperties(AlignmentI al, StringTokenizer st)
1591   {
1592
1593     // So far we have only added groups to the annotationHash,
1594     // the idea is in the future properties can be added to
1595     // alignments, other annotations etc
1596     if (al.getGroups() == null)
1597     {
1598       return;
1599     }
1600
1601     String name = st.nextToken();
1602     SequenceGroup sg = null;
1603     for (SequenceGroup _sg : al.getGroups())
1604     {
1605       if ((sg = _sg).getName().equals(name))
1606       {
1607         break;
1608       }
1609       else
1610       {
1611         sg = null;
1612       }
1613     }
1614
1615     if (sg != null)
1616     {
1617       String keyValue, key, value;
1618       ColourSchemeI def = sg.cs;
1619       sg.cs = null;
1620       while (st.hasMoreTokens())
1621       {
1622         keyValue = st.nextToken();
1623         key = keyValue.substring(0, keyValue.indexOf("="));
1624         value = keyValue.substring(keyValue.indexOf("=") + 1);
1625
1626         if (key.equalsIgnoreCase("description"))
1627         {
1628           sg.setDescription(value);
1629         }
1630         else if (key.equalsIgnoreCase("colour"))
1631         {
1632           sg.cs = ColourSchemeProperty.getColour(al, value);
1633         }
1634         else if (key.equalsIgnoreCase("pidThreshold"))
1635         {
1636           sg.cs.setThreshold(Integer.parseInt(value), true);
1637
1638         }
1639         else if (key.equalsIgnoreCase("consThreshold"))
1640         {
1641           sg.cs.setConservationInc(Integer.parseInt(value));
1642           Conservation c = new Conservation("Group",
1643                   ResidueProperties.propHash, 3, sg.getSequences(null),
1644                   sg.getStartRes(), sg.getEndRes() + 1);
1645
1646           c.calculate();
1647           c.verdict(false, 25); // TODO: refer to conservation percent threshold
1648
1649           sg.cs.setConservation(c);
1650
1651         }
1652         else if (key.equalsIgnoreCase("outlineColour"))
1653         {
1654           sg.setOutlineColour(new UserColourScheme(value).findColour('A'));
1655         }
1656         else if (key.equalsIgnoreCase("displayBoxes"))
1657         {
1658           sg.setDisplayBoxes(Boolean.valueOf(value).booleanValue());
1659         }
1660         else if (key.equalsIgnoreCase("showUnconserved"))
1661         {
1662           sg.setShowNonconserved(Boolean.valueOf(value).booleanValue());
1663         }
1664         else if (key.equalsIgnoreCase("displayText"))
1665         {
1666           sg.setDisplayText(Boolean.valueOf(value).booleanValue());
1667         }
1668         else if (key.equalsIgnoreCase("colourText"))
1669         {
1670           sg.setColourText(Boolean.valueOf(value).booleanValue());
1671         }
1672         else if (key.equalsIgnoreCase("textCol1"))
1673         {
1674           sg.textColour = new UserColourScheme(value).findColour('A');
1675         }
1676         else if (key.equalsIgnoreCase("textCol2"))
1677         {
1678           sg.textColour2 = new UserColourScheme(value).findColour('A');
1679         }
1680         else if (key.equalsIgnoreCase("textColThreshold"))
1681         {
1682           sg.thresholdTextColour = Integer.parseInt(value);
1683         }
1684         else if (key.equalsIgnoreCase("idColour"))
1685         {
1686           // consider warning if colour doesn't resolve to a real colour
1687           sg.setIdColour((def = new UserColourScheme(value))
1688                   .findColour('A'));
1689         }
1690         else if (key.equalsIgnoreCase("hide"))
1691         {
1692           // see bug https://mantis.lifesci.dundee.ac.uk/view.php?id=25847
1693           sg.setHidereps(true);
1694         }
1695         else if (key.equalsIgnoreCase("hidecols"))
1696         {
1697           // see bug https://mantis.lifesci.dundee.ac.uk/view.php?id=25847
1698           sg.setHideCols(true);
1699         }
1700         sg.recalcConservation();
1701       }
1702       if (sg.cs == null)
1703       {
1704         sg.cs = def;
1705       }
1706     }
1707   }
1708
1709   void setBelowAlignment(AlignmentI al, StringTokenizer st)
1710   {
1711     String token;
1712     AlignmentAnnotation aa, ala[] = al.getAlignmentAnnotation();
1713     if (ala == null)
1714     {
1715       System.err
1716               .print("Warning - no annotation to set below for sequence associated annotation:");
1717     }
1718     while (st.hasMoreTokens())
1719     {
1720       token = st.nextToken();
1721       if (ala == null)
1722       {
1723         System.err.print(" " + token);
1724       }
1725       else
1726       {
1727         for (int i = 0; i < al.getAlignmentAnnotation().length; i++)
1728         {
1729           aa = al.getAlignmentAnnotation()[i];
1730           if (aa.sequenceRef == refSeq && aa.label.equals(token))
1731           {
1732             aa.belowAlignment = true;
1733           }
1734         }
1735       }
1736     }
1737     if (ala == null)
1738     {
1739       System.err.print("\n");
1740     }
1741   }
1742
1743   void addAlignmentDetails(AlignmentI al, StringTokenizer st)
1744   {
1745     String keyValue, key, value;
1746     while (st.hasMoreTokens())
1747     {
1748       keyValue = st.nextToken();
1749       key = keyValue.substring(0, keyValue.indexOf("="));
1750       value = keyValue.substring(keyValue.indexOf("=") + 1);
1751       al.setProperty(key, value);
1752     }
1753   }
1754
1755   /**
1756    * Write annotations as a CSV file of the form 'label, value, value, ...' for
1757    * each row.
1758    *
1759    * @param annotations
1760    * @return CSV file as a string.
1761    */
1762   public String printCSVAnnotations(AlignmentAnnotation[] annotations)
1763   {
1764     StringBuffer sp = new StringBuffer();
1765     for (int i = 0; i < annotations.length; i++)
1766     {
1767       String atos = annotations[i].toString();
1768       int p = 0;
1769       do
1770       {
1771         int cp = atos.indexOf("\n", p);
1772         sp.append(annotations[i].label);
1773         sp.append(",");
1774         if (cp > p)
1775         {
1776           sp.append(atos.substring(p, cp + 1));
1777         }
1778         else
1779         {
1780           sp.append(atos.substring(p));
1781           sp.append(newline);
1782         }
1783         p = cp + 1;
1784       } while (p > 0);
1785     }
1786     return sp.toString();
1787   }
1788
1789   public String printAnnotationsForView(AlignViewportI viewport)
1790   {
1791     return printAnnotations(viewport.isShowAnnotation() ? viewport
1792             .getAlignment().getAlignmentAnnotation() : null, viewport
1793             .getAlignment().getGroups(), viewport.getAlignment()
1794             .getProperties(), viewport.getColumnSelection(),
1795             viewport.getAlignment(), null);
1796   }
1797
1798   public String printAnnotationsForAlignment(AlignmentI al)
1799   {
1800     return printAnnotations(al.getAlignmentAnnotation(), al.getGroups(),
1801             al.getProperties(), null, al, null);
1802   }
1803 }