ExactMapping is in pdbchain
[jalview.git] / src / jalview / analysis / AlignSeq.java
index 60ad85c..2fd1773 100755 (executable)
@@ -69,8 +69,8 @@ public class AlignSeq
     int maxj;\r
     int[] aseq1;\r
     int[] aseq2;\r
-    String astr1 = "";\r
-    String astr2 = "";\r
+    public String astr1 = "";\r
+    public String astr2 = "";\r
 \r
     /** DOCUMENT ME!! */\r
     public int seq1start;\r
@@ -94,7 +94,7 @@ public class AlignSeq
     int[][] lookup = ResidueProperties.getBLOSUM62();\r
     String[] intToStr = pep;\r
     int defInt = 23;\r
-    String output = "";\r
+    StringBuffer output = new StringBuffer();\r
     String type;\r
     Runtime rt;\r
 \r
@@ -169,7 +169,7 @@ public class AlignSeq
      */\r
     public String getOutput()\r
     {\r
-        return output;\r
+        return output.toString();\r
     }\r
 \r
     /**\r
@@ -263,7 +263,7 @@ public class AlignSeq
         }\r
         else\r
         {\r
-            output = output + ("Wrong type = dna or pep only");\r
+            output.append("Wrong type = dna or pep only");\r
             System.exit(0);\r
         }\r
 \r
@@ -406,7 +406,7 @@ public class AlignSeq
     /**\r
      * DOCUMENT ME!\r
      */\r
-    public void printAlignment()\r
+    public void printAlignment(java.io.PrintStream os)\r
     {\r
         // Find the biggest id length for formatting purposes\r
         int maxid = s1.getName().length();\r
@@ -420,39 +420,33 @@ public class AlignSeq
         int nochunks = ((aseq1.length - count) / len) + 1;\r
         pid = 0;\r
 \r
-        output = output + ("Score = " + score[maxi][maxj] + "\n");\r
-        output = output +\r
-            ("Length of alignment = " + (aseq1.length - count) + "\n");\r
-        output = output + ("Sequence ");\r
-        output = output + (new Format("%" + maxid + "s").form(s1.getName()));\r
-        output = output +\r
-            (" :  " + seq1start + " - " + seq1end + " (Sequence length = " +\r
+        output.append("Score = " + score[maxi][maxj] + "\n");\r
+        output.append("Length of alignment = " + (aseq1.length - count) + "\n");\r
+        output.append("Sequence ");\r
+        output.append(new Format("%" + maxid + "s").form(s1.getName()));\r
+        output.append(" :  " + seq1start + " - " + seq1end + " (Sequence length = " +\r
             s1str.length() + ")\n");\r
-        output = output + ("Sequence ");\r
-        output = output + (new Format("%" + maxid + "s").form(s2.getName()));\r
-        output = output +\r
-            (" :  " + seq2start + " - " + seq2end + " (Sequence length = " +\r
+        output .append("Sequence ");\r
+        output.append(new Format("%" + maxid + "s").form(s2.getName()));\r
+        output .append(" :  " + seq2start + " - " + seq2end + " (Sequence length = " +\r
             s2str.length() + ")\n\n");\r
 \r
         for (int j = 0; j < nochunks; j++)\r
         {\r
             // Print the first aligned sequence\r
-            output = output +\r
-                (new Format("%" + (maxid) + "s").form(s1.getName()) + " ");\r
+            output.append(new Format("%" + (maxid) + "s").form(s1.getName()) + " ");\r
 \r
             for (int i = 0; i < len; i++)\r
             {\r
                 if ((count + i + (j * len)) < aseq1.length)\r
                 {\r
-                    output = output +\r
-                        (new Format("%s").form(intToStr[aseq1[count + i +\r
+                    output.append(new Format("%s").form(intToStr[aseq1[count + i +\r
                             (j * len)]]));\r
                 }\r
             }\r
 \r
-            output = output + ("\n");\r
-            output = output +\r
-                (new Format("%" + (maxid) + "s").form(" ") + " ");\r
+            output.append("\n");\r
+            output.append(new Format("%" + (maxid) + "s").form(" ") + " ");\r
 \r
             // Print out the matching chars\r
             for (int i = 0; i < len; i++)\r
@@ -464,7 +458,7 @@ public class AlignSeq
                             !intToStr[aseq1[count + i + (j * len)]].equals("-"))\r
                     {\r
                         pid++;\r
-                        output = output + ("|");\r
+                        output.append("|");\r
                     }\r
                     else if (type.equals("pep"))\r
                     {\r
@@ -472,40 +466,42 @@ public class AlignSeq
                                     intToStr[aseq1[count + i + (j * len)]],\r
                                     intToStr[aseq2[count + i + (j * len)]]) > 0)\r
                         {\r
-                            output = output + (".");\r
+                            output.append(".");\r
                         }\r
                         else\r
                         {\r
-                            output = output + (" ");\r
+                            output.append(" ");\r
                         }\r
                     }\r
                     else\r
                     {\r
-                        output = output + (" ");\r
+                        output.append(" ");\r
                     }\r
                 }\r
             }\r
 \r
             // Now print the second aligned sequence\r
-            output = output + ("\n");\r
-            output = output +\r
-                (new Format("%" + (maxid) + "s").form(s2.getName()) + " ");\r
+            output = output.append("\n");\r
+            output = output.append(new Format("%" + (maxid) + "s").form(s2.getName()) + " ");\r
 \r
             for (int i = 0; i < len; i++)\r
             {\r
                 if ((count + i + (j * len)) < aseq1.length)\r
                 {\r
-                    output = output +\r
-                        (new Format("%s").form(intToStr[aseq2[count + i +\r
+                    output .append(new Format("%s").form(intToStr[aseq2[count + i +\r
                             (j * len)]]));\r
                 }\r
             }\r
 \r
-            output = output + ("\n\n");\r
+            output = output .append("\n\n");\r
         }\r
 \r
         pid = pid / (float) (aseq1.length - count) * 100;\r
-        output = output + (new Format("Percentage ID = %2.2f\n\n").form(pid));\r
+        output = output.append(new Format("Percentage ID = %2.2f\n\n").form(pid));\r
+\r
+        try{\r
+          os.println(output.toString());\r
+        }catch(Exception ex){}\r
     }\r
 \r
     /**\r
@@ -659,14 +655,14 @@ public class AlignSeq
     public static String extractGaps(String gapChar, String seq)\r
     {\r
         StringTokenizer str = new StringTokenizer(seq, gapChar);\r
-        String newString = "";\r
+        StringBuffer newString = new StringBuffer();\r
 \r
         while (str.hasMoreTokens())\r
         {\r
-            newString = newString + str.nextToken();\r
+            newString.append( str.nextToken() );\r
         }\r
 \r
-        return newString;\r
+        return newString.toString();\r
     }\r
 \r
     /**\r