vamsasDemo new branch
[jalview.git] / src / jalview / io / BLCFile.java
1 /*\r
2 * Jalview - A Sequence Alignment Editor and Viewer\r
3 * Copyright (C) 2005 AM Waterhouse, J Procter, G Barton, M Clamp, S Searle\r
4 *\r
5 * This program is free software; you can redistribute it and/or\r
6 * modify it under the terms of the GNU General Public License\r
7 * as published by the Free Software Foundation; either version 2\r
8 * of the License, or (at your option) any later version.\r
9 *\r
10 * This program is distributed in the hope that it will be useful,\r
11 * but WITHOUT ANY WARRANTY; without even the implied warranty of\r
12 * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the\r
13 * GNU General Public License for more details.\r
14 *\r
15 * You should have received a copy of the GNU General Public License\r
16 * along with this program; if not, write to the Free Software\r
17 * Foundation, Inc., 51 Franklin Street, Fifth Floor, Boston, MA  02110-1301, USA\r
18 */\r
19 package jalview.io;\r
20 \r
21 import jalview.datamodel.*;\r
22 \r
23 import java.io.*;\r
24 \r
25 import java.util.*;\r
26 \r
27 \r
28 /**\r
29  * DOCUMENT ME!\r
30  *\r
31  * @author $author$\r
32  * @version $Revision$\r
33  */\r
34 public class BLCFile extends AlignFile\r
35 {\r
36     Vector titles;\r
37 \r
38     /**\r
39      * Creates a new BLCFile object.\r
40      */\r
41     public BLCFile()\r
42     {\r
43     }\r
44 \r
45     /**\r
46      * Creates a new BLCFile object.\r
47      *\r
48      * @param inStr DOCUMENT ME!\r
49      */\r
50     public BLCFile(String inStr)\r
51     {\r
52         super(inStr);\r
53     }\r
54 \r
55     /**\r
56      * Creates a new BLCFile object.\r
57      *\r
58      * @param inFile DOCUMENT ME!\r
59      * @param type DOCUMENT ME!\r
60      *\r
61      * @throws IOException DOCUMENT ME!\r
62      */\r
63     public BLCFile(String inFile, String type) throws IOException\r
64     {\r
65         super(inFile, type);\r
66     }\r
67 \r
68     /**\r
69      * DOCUMENT ME!\r
70      */\r
71     public void initData()\r
72     {\r
73         super.initData();\r
74         titles = new Vector();\r
75     }\r
76 \r
77     /**\r
78      * DOCUMENT ME!\r
79      */\r
80     public void parse()\r
81     {\r
82         boolean idsFound = false;\r
83         Vector ids = new Vector();\r
84         StringBuffer[] seqstrings;\r
85         Vector starts = new Vector();\r
86         Vector ends = new Vector();\r
87 \r
88         String line = null;\r
89 \r
90         try\r
91         {\r
92             do\r
93             {\r
94                 line = nextLine();\r
95 \r
96                 // seek end of ids\r
97                 if (line.indexOf("*") > -1)\r
98                 {\r
99                     idsFound = true;\r
100 \r
101                     break;\r
102                 }\r
103 \r
104                 int abracket = line.indexOf(">");\r
105 \r
106                 if (abracket > -1)\r
107                 {\r
108                   if (line.indexOf(" ") > -1) //\r
109                   {\r
110                     line = line.substring(abracket + 1,\r
111                                           line.indexOf(" ", abracket + 1));\r
112                   }\r
113                   else\r
114                     line = line.substring(abracket+1);\r
115 \r
116                   Sequence seq = parseId(line);\r
117                   ids.addElement(seq.getName());\r
118                   starts.addElement(seq.getStart() + "");\r
119                   ends.addElement(seq.getEnd() + "");\r
120                 }\r
121             }\r
122             while (!idsFound);\r
123 \r
124             int starCol = line.indexOf("*");\r
125             seqstrings = new StringBuffer[ids.size()];\r
126 \r
127             for (int i = 0; i < ids.size(); i++)\r
128             {\r
129                 if (seqstrings[i] == null)\r
130                 {\r
131                     seqstrings[i] = new StringBuffer();\r
132                 }\r
133             }\r
134 \r
135             while ((line = nextLine()).indexOf("*") == -1)\r
136             {\r
137                 for (int i = 0; i < ids.size(); i++)\r
138                 {\r
139                     if (line.length() > (i + starCol))\r
140                     {\r
141                         seqstrings[i].append(line.charAt(i + starCol));\r
142                     }\r
143                 }\r
144             }\r
145 \r
146             for (int i = 0; i < ids.size(); i++)\r
147             {\r
148                 Sequence newSeq = new Sequence(ids.elementAt(i).toString(),\r
149                         seqstrings[i].toString(),\r
150                         Integer.parseInt(starts.elementAt(i).toString()),\r
151                         Integer.parseInt(ends.elementAt(i).toString()));\r
152                 seqs.addElement(newSeq);\r
153             }\r
154         }\r
155         catch (Exception ex)\r
156         {\r
157             ex.printStackTrace();\r
158         }\r
159     }\r
160 \r
161     /**\r
162      * DOCUMENT ME!\r
163      *\r
164      * @return DOCUMENT ME!\r
165      */\r
166     public String print()\r
167     {\r
168         return print(getSeqsAsArray());\r
169     }\r
170 \r
171     /**\r
172      * DOCUMENT ME!\r
173      *\r
174      * @param s DOCUMENT ME!\r
175      *\r
176      * @return DOCUMENT ME!\r
177      */\r
178     public String print(SequenceI[] s)\r
179     {\r
180         StringBuffer out = new StringBuffer();\r
181         /**\r
182          * A general parser for ids. Will look for dbrefs in\r
183          * Uniprot format source|id\r
184          * And also Jalview /start-end\r
185          *\r
186          * @String id Id to be parsed\r
187      */\r
188         int i = 0;\r
189         int max = -1;\r
190 \r
191         while ((i < s.length) && (s[i] != null))\r
192         {\r
193             out.append(">" + printId(s[i]) +"\n");\r
194 \r
195             if (s[i].getSequence().length() > max)\r
196             {\r
197                 max = s[i].getSequence().length();\r
198             }\r
199 \r
200             i++;\r
201         }\r
202 \r
203         out.append("* iteration 1\n");\r
204 \r
205         for (int j = 0; j < max; j++)\r
206         {\r
207             i = 0;\r
208 \r
209             while ((i < s.length) && (s[i] != null))\r
210             {\r
211                 if (s[i].getSequence().length() > j)\r
212                 {\r
213                     out.append(s[i].getSequence().substring(j, j + 1));\r
214                 }\r
215                 else\r
216                 {\r
217                     out.append("-");\r
218                 }\r
219 \r
220                 i++;\r
221             }\r
222 \r
223             out.append("\n");\r
224         }\r
225 \r
226         out.append("*\n");\r
227 \r
228         return out.toString();\r
229     }\r
230 }\r