make new dataset seq if residues deleted
authoramwaterhouse <Andrew Waterhouse>
Wed, 31 Jan 2007 14:47:35 +0000 (14:47 +0000)
committeramwaterhouse <Andrew Waterhouse>
Wed, 31 Jan 2007 14:47:35 +0000 (14:47 +0000)
src/jalview/datamodel/Sequence.java

index 02fc543..a0a858a 100755 (executable)
@@ -20,6 +20,7 @@ package jalview.datamodel;
 
 
 import java.util.*;
+import jalview.analysis.AlignSeq;
 
 
 /**
@@ -568,9 +569,31 @@ public class Sequence implements SequenceI
       System.arraycopy(sequence,j,tmp,i,sequence.length-j);
     }
 
+    if (this.datasetSequence != null)
+    {
+      for (int s = i; s < j; s++)
+      {
+        if (jalview.schemes.ResidueProperties.aaIndex[sequence[s]] != 23)
+        {
+
+          Sequence ds = new Sequence(name,
+                                     AlignSeq.extractGaps(
+                                         jalview.util.Comparison.GapChars,
+                                         this.getSequenceAsString()
+                                     ),
+                                     start,
+                                     end);
+          ds.setDescription(description);
+        }
+        break;
+      }
+    }
+
     sequence = tmp;
+
   }
 
+
   /**
    * DOCUMENT ME!
    *