X-Git-Url: http://source.jalview.org/gitweb/?a=blobdiff_plain;f=src%2Fjalview%2Fio%2FPIRFile.java;h=910d18929e48d60b65664bc33312104b556b1335;hb=ce886100c1eae01eac7d59c961e694f4493d4bcb;hp=d4b818b4d64c554c8433b687e811fc1168985c00;hpb=55e2e9b22b133db8b9ff0979b0338a33081fc8fd;p=jalview.git

diff --git a/src/jalview/io/PIRFile.java b/src/jalview/io/PIRFile.java
index d4b818b..910d189 100755
--- a/src/jalview/io/PIRFile.java
+++ b/src/jalview/io/PIRFile.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) 2007 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
@@ -21,23 +21,19 @@ package jalview.io;
 import java.io.*;
 import java.util.*;
 
-import jalview.analysis.*;
 import jalview.datamodel.*;
 
 public class PIRFile
     extends AlignFile
 {
+  public static boolean useModellerOutput = false;
+
   Vector words = new Vector(); //Stores the words in a line after splitting
 
   public PIRFile()
   {
   }
 
-  public PIRFile(String inStr)
-  {
-    super(inStr);
-  }
-
   public PIRFile(String inFile, String type)
       throws IOException
   {
@@ -45,54 +41,58 @@ public class PIRFile
   }
 
   public void parse()
+      throws IOException
   {
-    try
-    {
-      StringBuffer sequence;
-      String line = null;
+    StringBuffer sequence;
+    String line = null;
+    ModellerDescription md;
 
-      while ( (line = nextLine()) != null)
+    while ( (line = nextLine()) != null)
+    {
+      if (line.length() == 0)
       {
-        if(line.length()==0)
-        {
-          //System.out.println("blank line");
-          continue;
-        }
+        //System.out.println("blank line");
+        continue;
+      }
+      if (line.indexOf("C;") == 0 || line.indexOf("#") == 0)
+      {
+        continue;
+      }
+      Sequence newSeq = parseId(line.substring(line.indexOf(";") + 1));
 
-        Sequence newSeq = parseId(line.substring(line.indexOf(";") + 1));
+      sequence = new StringBuffer();
 
-        sequence = new StringBuffer();
+      newSeq.setDescription(nextLine()); // this is the title line
 
-        newSeq.setDescription(nextLine()); // this is the title line
+      boolean starFound = false;
 
-        boolean starFound = false;
+      while (!starFound)
+      {
+        line = nextLine();
+        sequence.append(line);
 
-        while(!starFound)
+        if (line == null)
         {
-          line = nextLine();
-          sequence.append(line);
-
-          if (line == null)
-            break;
-
-          if (line.indexOf("*") > -1)
-          {
-            starFound = true;
-          }
+          break;
         }
 
-        if(sequence.length()>0)
+        if (line.indexOf("*") > -1)
         {
-          sequence.setLength(sequence.length() - 1);
-
-          newSeq.setSequence(sequence.toString());
-          seqs.addElement(newSeq);
+          starFound = true;
         }
       }
-    }
-    catch (Exception ex)
-    {
-      ex.printStackTrace();
+
+      if (sequence.length() > 0)
+      {
+        sequence.setLength(sequence.length() - 1);
+        newSeq.setSequence(sequence.toString());
+
+        seqs.addElement(newSeq);
+
+        md = new ModellerDescription(newSeq.
+                                     getDescription());
+        md.updateSequenceI(newSeq);
+      }
     }
   }
 
@@ -107,20 +107,61 @@ public class PIRFile
     int len = 72;
     StringBuffer out = new StringBuffer();
     int i = 0;
+    ModellerDescription md;
 
     while ( (i < s.length) && (s[i] != null))
     {
-      String seq = s[i].getSequence();
+      String seq = s[i].getSequenceAsString();
       seq = seq + "*";
 
-      out.append(">P1;" + printId(s[i]) + "\n");
-
-      if(s[i].getDescription()!=null)
-        out.append(s[i].getDescription()+"\n");
+      if (is_NA)
+      {
+        // modeller doesn't really do nucleotides, so we don't do anything fancy
+        // Official tags area as follows, for now we'll use P1 and DL
+        // Protein (complete) P1
+        // Protein (fragment) F1
+        // DNA (linear) Dl
+        // DNA (circular) DC
+        // RNA (linear) RL
+        // RNA (circular) RC
+        // tRNA N3
+        // other functional RNA N1
+
+        out.append(">N1;" + s[i].getName() + "\n");
+        if (s[i].getDescription() == null)
+        {
+          out.append(s[i].getName() + " " +
+                     (s[i].getEnd() - s[i].getStart() + 1));
+          out.append(is_NA ? " bases\n" : " residues\n");
+        }
+        else
+        {
+          out.append(s[i].getDescription() + "\n");
+        }
+      }
       else
       {
-        out.append(s[i].getName()+" "+ (s[i].getEnd() - s[i].getStart() + 1));
-        out.append( is_NA ? " bases\n" : " residues\n");
+
+        if (useModellerOutput)
+        {
+          out.append(">P1;" + s[i].getName() + "\n");
+          md = new ModellerDescription(s[i]);
+          out.append(md.getDescriptionLine() + "\n");
+        }
+        else
+        {
+          out.append(">P1;" + printId(s[i]) + "\n");
+          if (s[i].getDescription() != null)
+          {
+            out.append(s[i].getDescription() + "\n");
+          }
+          else
+          {
+            out.append(s[i].getName() + " "
+                       + (s[i].getEnd() - s[i].getStart() + 1)
+                       + " residues\n");
+          }
+        }
       }
       int nochunks = (seq.length() / len) + 1;