X-Git-Url: http://source.jalview.org/gitweb/?a=blobdiff_plain;f=src%2Fjalview%2Fdatamodel%2FAlignment.java;h=cec2abcabe44a367208caa009544dc0ca7153ca9;hb=77b038f418aea84563953ad94f839592317ea500;hp=8b159a63de117095db79ad5bbb7d6c22e784d88c;hpb=1d6bdfa4ce1d88c055a257e542de6b8542144d91;p=jalview.git diff --git a/src/jalview/datamodel/Alignment.java b/src/jalview/datamodel/Alignment.java index 8b159a6..cec2abc 100755 --- a/src/jalview/datamodel/Alignment.java +++ b/src/jalview/datamodel/Alignment.java @@ -1,6 +1,6 @@ /* * Jalview - A Sequence Alignment Editor and Viewer - * Copyright (C) 2006 AM Waterhouse, J Procter, G Barton, M Clamp, S Searle + * Copyright (C) 2007 AM Waterhouse, J Procter, G Barton, M Clamp, S Searle * * This program is free software; you can redistribute it and/or * modify it under the terms of the GNU General Public License @@ -18,644 +18,713 @@ */ package jalview.datamodel; -import jalview.analysis.*; - import java.util.*; +import jalview.analysis.*; + /** Data structure to hold and manipulate a multiple sequence alignment */ -public class Alignment implements AlignmentI +public class Alignment + implements AlignmentI { - protected Alignment dataset; - protected Vector sequences; - protected Vector groups = new Vector(); - protected char gapCharacter = '-'; - protected int type = NUCLEOTIDE; - public static final int PROTEIN = 0; - public static final int NUCLEOTIDE = 1; + protected Alignment dataset; + protected Vector sequences; + protected Vector groups = new Vector(); + protected char gapCharacter = '-'; + protected int type = NUCLEOTIDE; + public static final int PROTEIN = 0; + public static final int NUCLEOTIDE = 1; - /** DOCUMENT ME!! */ - public AlignmentAnnotation[] annotations; + /** DOCUMENT ME!! */ + public AlignmentAnnotation[] annotations; - HiddenSequences hiddenSequences = new HiddenSequences(this); + HiddenSequences hiddenSequences = new HiddenSequences(this); - private void initAlignment(SequenceI[] seqs) { - int i=0; + public Hashtable alignmentProperties; - if( jalview.util.Comparison.isNucleotide(seqs)) - type = NUCLEOTIDE; - else - type = PROTEIN; - - sequences = new Vector(); - - for (i = 0; i < seqs.length; i++) - { - sequences.addElement(seqs[i]); - } + private void initAlignment(SequenceI[] seqs) + { + int i = 0; - } - /** Make an alignment from an array of Sequences. - * - * @param sequences - */ - public Alignment(SequenceI[] seqs) + if (jalview.util.Comparison.isNucleotide(seqs)) { - initAlignment(seqs); + type = NUCLEOTIDE; } - /** - * Make a new alignment from an array of SeqCigars - * @param seqs SeqCigar[] - */ - public Alignment(SeqCigar[] alseqs) { - SequenceI[] seqs = SeqCigar.createAlignmentSequences(alseqs, gapCharacter, new ColumnSelection(), null); - initAlignment(seqs); - } - /** - * Make a new alignment from an CigarArray - * JBPNote - can only do this when compactAlignment does not contain hidden regions. - * JBPNote - must also check that compactAlignment resolves to a set of SeqCigars - or construct them appropriately. - * @param compactAlignment CigarArray - */ - public static AlignmentI createAlignment(CigarArray compactAlignment) { - throw new Error("Alignment(CigarArray) not yet implemented"); - // this(compactAlignment.refCigars); + else + { + type = PROTEIN; } - /** - * DOCUMENT ME! - * - * @return DOCUMENT ME! - */ - public Vector getSequences() + sequences = new Vector(); + + for (i = 0; i < seqs.length; i++) { - return sequences; + sequences.addElement(seqs[i]); } - public SequenceI [] getSequencesArray() + } + + /** Make an alignment from an array of Sequences. + * + * @param sequences + */ + public Alignment(SequenceI[] seqs) + { + initAlignment(seqs); + } + + /** + * Make a new alignment from an array of SeqCigars + * @param seqs SeqCigar[] + */ + public Alignment(SeqCigar[] alseqs) + { + SequenceI[] seqs = SeqCigar.createAlignmentSequences(alseqs, gapCharacter, + new ColumnSelection(), null); + initAlignment(seqs); + } + + /** + * Make a new alignment from an CigarArray + * JBPNote - can only do this when compactAlignment does not contain hidden regions. + * JBPNote - must also check that compactAlignment resolves to a set of SeqCigars - or construct them appropriately. + * @param compactAlignment CigarArray + */ + public static AlignmentI createAlignment(CigarArray compactAlignment) + { + throw new Error("Alignment(CigarArray) not yet implemented"); + // this(compactAlignment.refCigars); + } + + /** + * DOCUMENT ME! + * + * @return DOCUMENT ME! + */ + public Vector getSequences() + { + return sequences; + } + + public SequenceI[] getSequencesArray() + { + SequenceI[] reply = new SequenceI[sequences.size()]; + for (int i = 0; i < sequences.size(); i++) { - SequenceI [] reply = new SequenceI[sequences.size()]; - for(int i=0; i -1 && i < getHeight()) { - deleteSequence(findIndex(s)); + sequences.removeElementAt(i); + hiddenSequences.adjustHeightSequenceDeleted(i); } + } - /** - * DOCUMENT ME! - * - * @param i DOCUMENT ME! - */ - public void deleteSequence(int i) + /** */ + public SequenceGroup findGroup(SequenceI s) + { + for (int i = 0; i < this.groups.size(); i++) { - if(i>-1 && i 0) + { + int i, iSize = sg.getSize(); + for (i = 0; i < iSize; i++) { - SequenceGroup sg = (SequenceGroup) groups.elementAt(i); - if(sg==null || sg.getSequences(null)==null) - { - this.deleteGroup(sg); - gSize--; - continue; - } - - if (sg.getSequences(null).contains(s)) - { - temp.addElement(sg); - } + if (!sequences.contains(sg.getSequenceAt(i))) + { + sg.deleteSequence(sg.getSequenceAt(i), false); + iSize--; + i--; + } } - SequenceGroup[] ret = new SequenceGroup[temp.size()]; - - for (int i = 0; i < temp.size(); i++) + if (sg.getSize() < 1) { - ret[i] = (SequenceGroup) temp.elementAt(i); + return; } + } - return ret; + groups.addElement(sg); } + } + /** + * DOCUMENT ME! + */ + public void deleteAllGroups() + { + groups.removeAllElements(); + } - - /** */ - public void addGroup(SequenceGroup sg) + /** */ + public void deleteGroup(SequenceGroup g) + { + if (groups.contains(g)) { - if (!groups.contains(sg)) - { - if(hiddenSequences.getSize()>0) - { - int i, iSize = sg.getSize(); - for (i = 0; i < iSize; i++) - { - if (!sequences.contains(sg.getSequenceAt(i))) - { - sg.deleteSequence(sg.getSequenceAt(i), false); - iSize--; - i--; - } - } - - if (sg.getSize() < 1) - return; - } - - groups.addElement(sg); - } + groups.removeElement(g); } + } - /** - * DOCUMENT ME! - */ - public void deleteAllGroups() + /** */ + public SequenceI findName(String name) + { + int i = 0; + + while (i < sequences.size()) { - groups.removeAllElements(); + if (getSequenceAt(i).getName().equals(name)) + { + return getSequenceAt(i); + } + + i++; } - /** */ - public void deleteGroup(SequenceGroup g) + return null; + } + + public SequenceI[] findSequenceMatch(String name) + { + Vector matches = new Vector(); + int i = 0; + + while (i < sequences.size()) { - if (groups.contains(g)) - { - groups.removeElement(g); - } + if (getSequenceAt(i).getName().equals(name)) + { + matches.addElement(getSequenceAt(i)); + } + i++; } - /** */ - public SequenceI findName(String name) + SequenceI[] result = new SequenceI[matches.size()]; + for (i = 0; i < result.length; i++) { - int i = 0; + result[i] = (SequenceI) matches.elementAt(i); + } - while (i < sequences.size()) - { - if (getSequenceAt(i).getName().equals(name)) - { - return getSequenceAt(i); - } + return result; - i++; - } + } - return null; - } + /** */ + public int findIndex(SequenceI s) + { + int i = 0; - public SequenceI [] findSequenceMatch(String name) + while (i < sequences.size()) { - Vector matches = new Vector(); - int i = 0; - - while (i < sequences.size()) + if (s == getSequenceAt(i)) { - if (getSequenceAt(i).getName().equals(name)) - { - matches.addElement(getSequenceAt(i)); - } - i++; + return i; } - SequenceI [] result = new SequenceI[matches.size()]; - for(i=0; i maxLength) + { + maxLength = getSequenceAt(i).getLength(); + } + } - i++; - } + return maxLength; + } - return -1; - } + /** + * DOCUMENT ME! + * + * @param gc DOCUMENT ME! + */ + public void setGapCharacter(char gc) + { + gapCharacter = gc; - /** - * DOCUMENT ME! - * - * @return DOCUMENT ME! - */ - public int getHeight() + for (int i = 0; i < sequences.size(); i++) { - return sequences.size(); + Sequence seq = (Sequence) sequences.elementAt(i); + seq.setSequence(seq.getSequenceAsString() + .replace('.', gc) + .replace('-', gc) + .replace(' ', gc) + ); } + } - /** - * DOCUMENT ME! - * - * @return DOCUMENT ME! - */ - public int getWidth() - { - int maxLength = -1; + /** + * DOCUMENT ME! + * + * @return DOCUMENT ME! + */ + public char getGapCharacter() + { + return gapCharacter; + } - for (int i = 0; i < sequences.size(); i++) - { - if (getSequenceAt(i).getLength() > maxLength) - { - maxLength = getSequenceAt(i).getLength(); - } - } + /** + * DOCUMENT ME! + * + * @return DOCUMENT ME! + */ + public boolean isAligned() + { + int width = getWidth(); - return maxLength; + for (int i = 0; i < sequences.size(); i++) + { + if (getSequenceAt(i).getLength() != width) + { + return false; + } } + return true; + } + /* (non-Javadoc) + * @see jalview.datamodel.AlignmentI#deleteAnnotation(jalview.datamodel.AlignmentAnnotation) + */ + public boolean deleteAnnotation(AlignmentAnnotation aa) + { + int aSize = 1; - /** - * DOCUMENT ME! - * - * @param gc DOCUMENT ME! - */ - public void setGapCharacter(char gc) + if (annotations != null) { - gapCharacter = gc; - - for (int i = 0; i < sequences.size(); i++) - { - Sequence seq = (Sequence) sequences.elementAt(i); - seq.setSequence(seq.getSequenceAsString() - .replace('.', gc) - .replace('-', gc) - .replace(' ', gc) - ); - } + aSize = annotations.length; } - /** - * DOCUMENT ME! - * - * @return DOCUMENT ME! - */ - public char getGapCharacter() + if (aSize < 1) { - return gapCharacter; + return false; } + AlignmentAnnotation[] temp = new AlignmentAnnotation[aSize - 1]; - /** - * DOCUMENT ME! - * - * @return DOCUMENT ME! - */ - public boolean isAligned() - { - int width = getWidth(); + boolean swap=false; + int tIndex = 0; - for (int i = 0; i < sequences.size(); i++) - { - if (getSequenceAt(i).getLength() != width) - { - return false; - } - } + for (int i = 0; tIndex 1) - { - for (i = 0; i < (aSize-1); i++) - { - temp[i] = annotations[i]; - } - } - - annotations = temp; + if (aSize > 1) + { + for (i = 0; i < (aSize - 1); i++) + { + temp[i] = annotations[i]; + } } - public void setAnnotationIndex(AlignmentAnnotation aa, int index) + annotations = temp; + } + + public void setAnnotationIndex(AlignmentAnnotation aa, int index) + { + if (aa == null || annotations == null || annotations.length - 1 < index) { - if(aa==null || annotations==null || annotations.length-1 maxLength; j--) { - current = getSequenceAt(i); - for (int j = current.getLength(); j > maxLength; j--) + if (j > maxLength && !jalview.util.Comparison.isGap( + current.getCharAt(j))) { - if (j > maxLength && !jalview.util.Comparison.isGap( - current.getCharAt(j))) - { - maxLength = j; - break; - } + maxLength = j; + break; } } + } - maxLength++; + maxLength++; - int cLength; - for (int i = 0; i < sequences.size(); - i++) - { - current = getSequenceAt(i); - cLength = current.getLength(); + int cLength; + for (int i = 0; i < sequences.size(); + i++) + { + current = getSequenceAt(i); + cLength = current.getLength(); - if (cLength < maxLength) - { - current.insertCharAt(cLength, - maxLength-cLength, gapCharacter); - modified=true; - } - else if(current.getLength() > maxLength) - { - current.deleteChars(maxLength, current.getLength()); - } + if (cLength < maxLength) + { + current.insertCharAt(cLength, + maxLength - cLength, gapCharacter); + modified = true; + } + else if (current.getLength() > maxLength) + { + current.deleteChars(maxLength, current.getLength()); } - return modified; } + return modified; + } - public HiddenSequences getHiddenSequences() - { - return hiddenSequences; - } + public HiddenSequences getHiddenSequences() + { + return hiddenSequences; + } public CigarArray getCompactAlignment() { SeqCigar alseqs[] = new SeqCigar[sequences.size()]; - for (int i=0; i