X-Git-Url: http://source.jalview.org/gitweb/?a=blobdiff_plain;f=src%2Fjalview%2Fdatamodel%2FSequence.java;h=7df3eaf5c7c1dd18cc0be284cab4cf8634bf9991;hb=d587f1aa61946dc14f6f089cf1dc2a3116cfb773;hp=dee67d374f3d7e4111b645763fca0eb3e0f068c3;hpb=c431a65bbe794a0dd27f3f806d3c5cf37c53279c;p=jalview.git diff --git a/src/jalview/datamodel/Sequence.java b/src/jalview/datamodel/Sequence.java index dee67d3..7df3eaf 100755 --- a/src/jalview/datamodel/Sequence.java +++ b/src/jalview/datamodel/Sequence.java @@ -1,6 +1,6 @@ /* * Jalview - A Sequence Alignment Editor and Viewer -* Copyright (C) 2005 AM Waterhouse, J Procter, G Barton, M Clamp, S Searle +* Copyright (C) 2006 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 @@ -31,683 +31,665 @@ import java.util.*; */ public class Sequence implements SequenceI { - SequenceI datasetSequence; - String name; - private String sequence; - String description; - int start; - int end; - Color color = Color.white; - Vector pdbIds; - String vamsasId; - DBRefEntry [] dbrefs; - - /** This annotation is displayed below the alignment but the - * positions are tied to the residues of this sequence */ - Vector annotation; - - /** DOCUMENT ME!! */ - public SequenceFeature [] sequenceFeatures; - - /** This array holds hidden sequences - * of which this sequence is the representitive member of a group - */ - SequenceGroup hiddenSequences; - - /** - * Creates a new Sequence object. - * - * @param name DOCUMENT ME! - * @param sequence DOCUMENT ME! - * @param start DOCUMENT ME! - * @param end DOCUMENT ME! - */ - public Sequence(String name, String sequence, int start, int end) - { - this.name = name; - this.sequence = sequence; - this.start = start; - this.end = end; - - parseId(); - - checkValidRange(); - } - - com.stevesoft.pat.Regex limitrx = new com.stevesoft.pat.Regex( - "[/][0-9]{1,}[-][0-9]{1,}$"); - com.stevesoft.pat.Regex endrx = new com.stevesoft.pat.Regex( - "[0-9]{1,}$"); - - void parseId() - { - // Does sequence have the /start-end signiature? - if(limitrx.search(name)) - { - name = limitrx.left(); - endrx.search(limitrx.stringMatched()); - setStart( Integer.parseInt( limitrx.stringMatched().substring(1,endrx.matchedFrom()-1 ))); - setEnd( Integer.parseInt( endrx.stringMatched() )); - } - } - - void checkValidRange() - { - if (end < 1) + SequenceI datasetSequence; + String name; + private String sequence; + String description; + int start; + int end; + Vector pdbIds; + String vamsasId; + DBRefEntry[] dbrefs; + + /** This annotation is displayed below the alignment but the + * positions are tied to the residues of this sequence */ + Vector annotation; + + /** DOCUMENT ME!! */ + public SequenceFeature[] sequenceFeatures; + + /** This array holds hidden sequences + * of which this sequence is the representitive member of a group + */ + SequenceGroup hiddenSequences; + + /** + * Creates a new Sequence object. + * + * @param name DOCUMENT ME! + * @param sequence DOCUMENT ME! + * @param start DOCUMENT ME! + * @param end DOCUMENT ME! + */ + public Sequence(String name, String sequence, int start, int end) + { + this.name = name; + this.sequence = sequence; + this.start = start; + this.end = end; + + parseId(); + + checkValidRange(); + } + + com.stevesoft.pat.Regex limitrx = new com.stevesoft.pat.Regex( + "[/][0-9]{1,}[-][0-9]{1,}$"); + com.stevesoft.pat.Regex endrx = new com.stevesoft.pat.Regex( + "[0-9]{1,}$"); + + void parseId() + { + // Does sequence have the /start-end signiature? + if (limitrx.search(name)) + { + name = limitrx.left(); + endrx.search(limitrx.stringMatched()); + setStart(Integer.parseInt(limitrx.stringMatched().substring(1, + endrx.matchedFrom() - 1))); + setEnd(Integer.parseInt(endrx.stringMatched())); + } + } + + void checkValidRange() + { + if (end < 1) + { + int endRes = 0; + char ch; + for (int j = 0; j < sequence.length(); j++) { - int endRes = 0; - char ch; - for (int j = 0; j < sequence.length(); j++) + ch = sequence.charAt(j); + if (!jalview.util.Comparison.isGap( (ch))) { - ch = sequence.charAt(j); - if (!jalview.util.Comparison.isGap( (ch))) - { - endRes++; - } + endRes++; } - if (endRes > 0) - { - endRes += start - 1; - } - - this.end = endRes; - } - - } - - /** - * Creates a new Sequence object. - * - * @param name DOCUMENT ME! - * @param sequence DOCUMENT ME! - */ - public Sequence(String name, String sequence) - { - this(name, sequence, 1, -1); - } - - /** - * Creates a new Sequence object. - * - * @param seq DOCUMENT ME! - */ - public Sequence(SequenceI seq) - { - this(seq.getName(), seq.getSequence(), seq.getStart(), seq.getEnd()); - } - - /** - * DOCUMENT ME! - * - * @param v DOCUMENT ME! - */ - public void setSequenceFeatures(SequenceFeature [] features) - { - sequenceFeatures = features; - } - - public synchronized void addSequenceFeature(SequenceFeature sf) - { - if(sequenceFeatures==null) - { - sequenceFeatures = new SequenceFeature[0]; } - - for(int i=0; i 0) { - if(sequenceFeatures[i].equals(sf)) - { - return; - } + endRes += start - 1; } - SequenceFeature [] temp = new SequenceFeature[sequenceFeatures.length+1]; - System.arraycopy(sequenceFeatures, 0, temp, 0, sequenceFeatures.length); - temp[sequenceFeatures.length] = sf; - - sequenceFeatures = temp; + this.end = endRes; } + } - /** - * DOCUMENT ME! - * - * @return DOCUMENT ME! - */ - public SequenceFeature [] getSequenceFeatures() - { - return sequenceFeatures; - } - - public void addPDBId(PDBEntry entry) - { - if(pdbIds == null) - pdbIds = new Vector(); + /** + * Creates a new Sequence object. + * + * @param name DOCUMENT ME! + * @param sequence DOCUMENT ME! + */ + public Sequence(String name, String sequence) + { + this(name, sequence, 1, -1); + } - pdbIds.addElement(entry); - } + /** + * Creates a new Sequence object. + * + * @param seq DOCUMENT ME! + */ + public Sequence(SequenceI seq) + { + this(seq.getName(), seq.getSequence(), seq.getStart(), seq.getEnd()); + } - /** - * DOCUMENT ME! - * - * @param id DOCUMENT ME! - */ - public void setPDBId(Vector id) - { - pdbIds = id; - } + /** + * DOCUMENT ME! + * + * @param v DOCUMENT ME! + */ + public void setSequenceFeatures(SequenceFeature[] features) + { + sequenceFeatures = features; + } - /** - * DOCUMENT ME! - * - * @return DOCUMENT ME! - */ - public Vector getPDBId() + public synchronized void addSequenceFeature(SequenceFeature sf) + { + if (sequenceFeatures == null) { - return pdbIds; + sequenceFeatures = new SequenceFeature[0]; } - /** - * DOCUMENT ME! - * - * @return DOCUMENT ME! - */ - public String getDisplayId(boolean jvsuffix) + for (int i = 0; i < sequenceFeatures.length; i++) { - StringBuffer result = new StringBuffer(name); - if (jvsuffix) + if (sequenceFeatures[i].equals(sf)) { - result.append("/" + start + "-" + end); + return; } - - return result.toString(); } - /** - * DOCUMENT ME! - * - * @param name DOCUMENT ME! - */ - public void setName(String name) - { - this.name = name; - this.parseId(); - } + SequenceFeature[] temp = new SequenceFeature[sequenceFeatures.length + 1]; + System.arraycopy(sequenceFeatures, 0, temp, 0, sequenceFeatures.length); + temp[sequenceFeatures.length] = sf; - /** - * DOCUMENT ME! - * - * @return DOCUMENT ME! - */ - public String getName() - { - return this.name; - } + sequenceFeatures = temp; + } - /** - * DOCUMENT ME! - * - * @param start DOCUMENT ME! - */ - public void setStart(int start) - { - this.start = start; - } + public void deleteFeature(SequenceFeature sf) + { + if(sequenceFeatures==null) + return; - /** - * DOCUMENT ME! - * - * @return DOCUMENT ME! - */ - public int getStart() + int index=0; + for (index = 0; index < sequenceFeatures.length; index++) { - return this.start; - } - - /** - * DOCUMENT ME! - * - * @param end DOCUMENT ME! - */ - public void setEnd(int end) - { - this.end = end; - } - - /** - * DOCUMENT ME! - * - * @return DOCUMENT ME! - */ - public int getEnd() - { - return this.end; + if (sequenceFeatures[index].equals(sf)) + { + break; + } } - /** - * DOCUMENT ME! - * - * @return DOCUMENT ME! - */ - public int getLength() - { - return this.sequence.length(); - } - /** - * DOCUMENT ME! - * - * @param seq DOCUMENT ME! - */ - public void setSequence(String seq) - { - this.sequence = seq; - checkValidRange(); - } + if(index==sequenceFeatures.length) + return; - /** - * DOCUMENT ME! - * - * @return DOCUMENT ME! - */ - public String getSequence() + int sfLength = sequenceFeatures.length; + if(sfLength<2) { - return this.sequence; + sequenceFeatures = null; } - - /** - * DOCUMENT ME! - * - * @param start DOCUMENT ME! - * @param end DOCUMENT ME! - * - * @return DOCUMENT ME! - */ - public String getSequence(int start, int end) + else { - // JBPNote - left to user to pad the result here (TODO:Decide on this policy) - if (start >= sequence.length()) - { - return ""; - } + SequenceFeature[] temp = new SequenceFeature[sfLength-1]; + System.arraycopy(sequenceFeatures, 0, temp, 0, index); - if (end >= sequence.length()) - { - end = sequence.length(); - } + if(index= sequence.length()) + { + return ""; + } + + if (end >= sequence.length()) + { + end = sequence.length(); + } + + return this.sequence.substring(start, end); + } + + /** + * make a new Sequence object from start to end (including gaps) over this seqeunce + * @param start int + * @param end int + * @return SequenceI + */ + public SequenceI getSubSequence(int start, int end) + { + if (start < 0) + start = 0; + String seq = getSequence(start, end); + if (seq == "") + return null; + int nstart = findPosition(start); + int nend = findPosition(end) - 1; + // JBPNote - this is an incomplete copy. + SequenceI nseq = new Sequence(this.getName(), seq, nstart, nend); + nseq.setDescription(description); + nseq.setDatasetSequence(getDatasetSequence()); + return nseq; + } + + /** + * DOCUMENT ME! + * + * @param i DOCUMENT ME! + * + * @return DOCUMENT ME! + */ + public char getCharAt(int i) + { + if (i < sequence.length()) + { + return sequence.charAt(i); + } + else + { + return ' '; + } + } + + /** + * DOCUMENT ME! + * + * @param desc DOCUMENT ME! + */ + public void setDescription(String desc) + { + this.description = desc; + } + + /** + * DOCUMENT ME! + * + * @return DOCUMENT ME! + */ + public String getDescription() + { + return this.description; + } + + /** + * DOCUMENT ME! + * + * @param pos DOCUMENT ME! + * + * @return DOCUMENT ME! + */ + public int findIndex(int pos) + { + // returns the alignment position for a residue + int j = start; + int i = 0; + + while ( (i < sequence.length()) && (j <= end) && (j <= pos)) + { + if (!jalview.util.Comparison.isGap(sequence.charAt(i))) + { + j++; + } - /** - * DOCUMENT ME! - * - * @return DOCUMENT ME! - */ - public String getDescription() - { - return this.description; + i++; } - /** - * DOCUMENT ME! - * - * @param pos DOCUMENT ME! - * - * @return DOCUMENT ME! - */ - public int findIndex(int pos) + if ( (j == end) && (j < pos)) { - // returns the alignment position for a residue - int j = start; - int i = 0; - - while ((i < sequence.length()) && (j <= end) && (j <= pos)) - { - if (!jalview.util.Comparison.isGap(sequence.charAt(i))) - { - j++; - } - - i++; - } - - if ((j == end) && (j < pos)) - { - return end + 1; - } - else - { - return i; - } + return end + 1; } - - /** - * Returns the sequence position for an alignment position - * - * @param i column index in alignment (from 1) - * - * @return residue number for residue (left of and) nearest ith column - */ - public int findPosition(int i) + else { - int j = 0; - int pos = start; - int seqlen=sequence.length(); - while ((j < i) && (j < seqlen)) - { - if (!jalview.util.Comparison.isGap((sequence.charAt(j)))) - { - pos++; - } - - j++; - } - - return pos; + return i; } + } - /** - * DOCUMENT ME! - * - * @return DOCUMENT ME! - */ - public int[] gapMap() + /** + * Returns the sequence position for an alignment position + * + * @param i column index in alignment (from 1) + * + * @return residue number for residue (left of and) nearest ith column + */ + public int findPosition(int i) + { + int j = 0; + int pos = start; + int seqlen = sequence.length(); + while ( (j < i) && (j < seqlen)) { - // Returns an int array giving the position of each residue in the sequence in the alignment - String seq = jalview.analysis.AlignSeq.extractGaps(jalview.util.Comparison.GapChars, sequence); - int[] map = new int[seq.length()]; - int j = 0; - int p = 0; - - while (j < sequence.length()) - { - if (!jalview.util.Comparison.isGap(sequence.charAt(j))) - { - map[p++] = j; - } - - j++; - } + if (!jalview.util.Comparison.isGap( (sequence.charAt(j)))) + { + pos++; + } - return map; + j++; } - /** - * DOCUMENT ME! - * - * @param i DOCUMENT ME! - */ - public void deleteCharAt(int i) - { - if (i >= sequence.length()) - { - return; - } + return pos; + } - sequence = sequence.substring(0, i) + sequence.substring(i + 1); - } + /** + * Returns an int array where indices correspond to each residue in the sequence and the element value gives its position in the alignment + * + * @return int[SequenceI.getEnd()-SequenceI.getStart()+1] or null if no residues in SequenceI object + */ + public int[] gapMap() + { + String seq = jalview.analysis.AlignSeq.extractGaps(jalview.util.Comparison. + GapChars, sequence); + int[] map = new int[seq.length()]; + int j = 0; + int p = 0; - /** - * DOCUMENT ME! - * - * @param i DOCUMENT ME! - * @param j DOCUMENT ME! - */ - public void deleteChars(int i, int j) + while (j < sequence.length()) { - if (i >= sequence.length()) - { - return; - } + if (!jalview.util.Comparison.isGap(sequence.charAt(j))) + { + map[p++] = j; + } - if (j >= sequence.length()) - { - sequence = sequence.substring(0, i); - } - else - { - sequence = sequence.substring(0, i) + sequence.substring(j); - } + j++; } + return map; + } - /** - * DOCUMENT ME! - * - * @param i DOCUMENT ME! - * @param c DOCUMENT ME! - * @param chop DOCUMENT ME! - */ - public void insertCharAt(int i, char c) - { - String tmp = new String(sequence); - - if (i < sequence.length()) - { - sequence = tmp.substring(0, i) + String.valueOf(c) + - tmp.substring(i); - } - else - { - // JBPNote : padding char at end of sequence. We'll not get away with this when we insert residues, I bet! - char[] ch = new char[(1 + i) - sequence.length()]; - - for (int j = 0, k = ch.length; j < k; j++) - ch[j] = c; - - sequence = tmp + String.valueOf(ch); - } - } - - /** - * DOCUMENT ME! - * - * @param c DOCUMENT ME! - */ - public void setColor(Color c) + /** + * DOCUMENT ME! + * + * @param i DOCUMENT ME! + */ + public void deleteCharAt(int i) + { + if (i >= sequence.length()) { - this.color = c; + return; } - /** - * DOCUMENT ME! - * - * @return DOCUMENT ME! - */ - public Color getColor() - { - return color; - } + sequence = sequence.substring(0, i) + sequence.substring(i + 1); + } - public String getVamsasId() + /** + * DOCUMENT ME! + * + * @param i DOCUMENT ME! + * @param j DOCUMENT ME! + */ + public void deleteChars(int i, int j) + { + if (i >= sequence.length()) { - return vamsasId; + return; } - public void setVamsasId(String id) + if (j >= sequence.length()) { - vamsasId = id; + sequence = sequence.substring(0, i); } - - public void setDBRef(DBRefEntry [] dbref) + else { - dbrefs = dbref; + sequence = sequence.substring(0, i) + sequence.substring(j); } + } - public DBRefEntry [] getDBRef() - { - return dbrefs; - } + /** + * DOCUMENT ME! + * + * @param i DOCUMENT ME! + * @param c DOCUMENT ME! + * @param chop DOCUMENT ME! + */ + public void insertCharAt(int i, int length, char c) + { + StringBuffer tmp; - public void addDBRef(DBRefEntry entry) + if (i >= sequence.length()) { - if(dbrefs == null) - dbrefs = new DBRefEntry[0]; - - DBRefEntry [] temp = new DBRefEntry[dbrefs.length+1]; - System.arraycopy(dbrefs, 0, temp, 0, dbrefs.length); - - temp[temp.length-1] = entry; - - dbrefs = temp; + tmp = new StringBuffer(sequence); } + else + tmp = new StringBuffer(sequence.substring(0, i)); - public void setDatasetSequence(SequenceI seq) + while (length > 0) { - datasetSequence = seq; + tmp.append(c); + length--; } - public SequenceI getDatasetSequence() + if (i < sequence.length()) { - return datasetSequence; + tmp.append(sequence.substring(i)); } - public AlignmentAnnotation [] getAnnotation() - { - if(annotation==null) - return null; + sequence = tmp.toString(); + } - AlignmentAnnotation [] ret = new AlignmentAnnotation[annotation.size()]; - for(int r = 0; rsequence.length()) - end = sequence.length(); + int i, iSize = dbrefs.length; - if (start > 0) + for(i=0; isequence.length()) - end = sequence.length(); + public AlignmentAnnotation[] getAnnotation() + { + if (annotation == null) + return null; - if (start > 0) - { - newSeq.append(sequence.substring(0, start)); - } + AlignmentAnnotation[] ret = new AlignmentAnnotation[annotation.size()]; + for (int r = 0; r < ret.length; r++) + ret[r] = (AlignmentAnnotation) annotation.elementAt(r); - char nextChar; - for(int c=start; c