From: amwaterhouse Date: Tue, 18 Jul 2006 15:30:37 +0000 (+0000) Subject: getSelectionAsNewSequence added X-Git-Tag: Release_2_1~277 X-Git-Url: http://source.jalview.org/gitweb/?a=commitdiff_plain;h=79cc40807a44006b6510d6533e57be05a66b760c;p=jalview.git getSelectionAsNewSequence added --- diff --git a/src/jalview/datamodel/SequenceGroup.java b/src/jalview/datamodel/SequenceGroup.java index ecdf2bc..19e0da0 100755 --- a/src/jalview/datamodel/SequenceGroup.java +++ b/src/jalview/datamodel/SequenceGroup.java @@ -113,6 +113,62 @@ public class SequenceGroup endRes = end; } + public SequenceI [] getSelectionAsNewSequences(AlignmentI align) + { + int iSize = sequences.size(); + SequenceI [] seqs = new SequenceI[iSize]; + SequenceI [] inorder = getSequencesInOrder(align); + + char ch; + int sres, eres; + + for (int i = 0; i < iSize; i++) + { + SequenceI seq = inorder[i]; + + //FIND START RES + //Returns residue following index if gap + sres = seq.findPosition(startRes); + + //FIND END RES + //Need to find the residue preceeding index if gap + eres = 0; + + for (int j = 0; j < endRes + 1 && j < seq.getLength(); j++) + { + ch = seq.getCharAt(j); + if (!jalview.util.Comparison.isGap( (ch))) + { + eres++; + } + } + + if (eres > 0) + { + eres += seq.getStart() - 1; + } + + seqs[i] = new Sequence(seq.getName(), + seq.getSequence(startRes, endRes + 1), + sres, + eres); + seqs[i].setDescription(seq.getDescription()); + seqs[i].setDBRef(seq.getDBRef()); + seqs[i].setSequenceFeatures(seq.getSequenceFeatures()); + if (seqs[i].getDatasetSequence() != null) + seqs[i].setDatasetSequence(seq.getDatasetSequence()); + + if(seq.getAnnotation()!=null) + { + for(int a=0; a