Description added
[jalview.git] / src / jalview / datamodel / SeqCigar.java
index 74a7812..7363e9e 100644 (file)
@@ -59,7 +59,7 @@ public class SeqCigar
   public String getSequenceString(char GapChar)
   {
     return (length == 0) ? "" :
-        (String) getSequenceAndDeletions(refseq.getSequence().substring(start, end), GapChar)[0];
+        (String) getSequenceAndDeletions(refseq.getSequenceAsString(start, end), GapChar)[0];
   }
 
   /**
@@ -73,7 +73,7 @@ public class SeqCigar
     {
       return null;
     }
-    Object[] edit_result = getSequenceAndDeletions(refseq.getSequence().substring(start,end),
+    Object[] edit_result = getSequenceAndDeletions(refseq.getSequenceAsString(start,end),
         GapChar);
     if (edit_result == null)
     {
@@ -86,6 +86,7 @@ public class SeqCigar
                        refseq.getStart() + start+((bounds[2]==0) ? -1 : bounds[2]));
     // seq.checkValidRange(); probably not needed
     seq.setDatasetSequence(refseq);
+    seq.setDescription(refseq.getDescription());
     return seq;
   }
 
@@ -113,7 +114,7 @@ public class SeqCigar
     }
     if (_s<0)
       throw new Error("Implementation Error: _s="+_s);
-    String seq_string = seq.getSequence();
+    String seq_string = seq.getSequenceAsString();
     if (_e==0 || _e<_s || _e>seq_string.length())
       _e=seq_string.length();
     // resolve start and end positions relative to ungapped reference sequence
@@ -375,7 +376,7 @@ public class SeqCigar
     for (int i = 0; i < alseqs.length; i++)
     {
       alseqs_string[i]=alseqs[i].getRefSeq().
-          getSequence().substring(alseqs[i].start,alseqs[i].end);
+          getSequenceAsString(alseqs[i].start,alseqs[i].end);
       gs_regions[i] = alseqs[i].getSequenceAndDeletions(alseqs_string[i], gapCharacter); // gapped sequence, {start, start col, end. endcol}, hidden regions {{start, end, col}})
       if (gs_regions[i] == null)
       {
@@ -478,8 +479,8 @@ public class SeqCigar
     if (!gen_sgapped_s.getSequence().equals(s_gapped.getSequence()))
     {
       System.err.println("Couldn't reconstruct sequence.\n" +
-                         gen_sgapped_s.getSequence() + "\n" +
-                         s_gapped.getSequence());
+                         gen_sgapped_s.getSequenceAsString() + "\n" +
+                         s_gapped.getSequenceAsString());
       return false;
     }
     return true;
@@ -551,7 +552,7 @@ public class SeqCigar
                        + "\nCigar String:" + sub_se_gp.getCigarstring() + "\n"
         );
     SequenceI ssgp = sub_se_gp.getSeq('-');
-    System.out.println("\t " + ssgp.getSequence());
+    System.out.println("\t " + ssgp.getSequenceAsString());
     for (int r = 0; r < 10; r++)
     {
       sub_se_gp = new SeqCigar(s_subsequence_gapped, 8, 48);
@@ -561,7 +562,7 @@ public class SeqCigar
       {
         int e = st + rs;
         sub_se_gp.deleteRange(st, e);
-        String ssgapedseq = sub_se_gp.getSeq('-').getSequence();
+        String ssgapedseq = sub_se_gp.getSeq('-').getSequenceAsString();
         System.out.println(st + "," + e + "\t:" + ssgapedseq);
         st -=3;
       }
@@ -577,7 +578,7 @@ public class SeqCigar
         System.out.println("" + al.getSequenceAt(i).getName() + "\t" +
                            al.getSequenceAt(i).getStart() + "\t" +
                            al.getSequenceAt(i).getEnd() + "\t" +
-                           al.getSequenceAt(i).getSequence());
+                           al.getSequenceAt(i).getSequenceAsString());
       }
     }
     {
@@ -593,7 +594,7 @@ public class SeqCigar
         System.out.println("" + al.getSequenceAt(i).getName() + "\t" +
                            al.getSequenceAt(i).getStart() + "\t" +
                            al.getSequenceAt(i).getEnd() + "\t" +
-                           al.getSequenceAt(i).getSequence());
+                           al.getSequenceAt(i).getSequenceAsString());
       }
     }
 //    if (!ssgapedseq.equals("ryas---dtqqwa----slchvh"))