2 * Jalview - A Sequence Alignment Editor and Viewer ($$Version-Rel$$)
3 * Copyright (C) $$Year-Rel$$ The Jalview Authors
5 * This file is part of Jalview.
7 * Jalview is free software: you can redistribute it and/or
8 * modify it under the terms of the GNU General Public License
9 * as published by the Free Software Foundation, either version 3
10 * of the License, or (at your option) any later version.
12 * Jalview is distributed in the hope that it will be useful, but
13 * WITHOUT ANY WARRANTY; without even the implied warranty
14 * of MERCHANTABILITY or FITNESS FOR A PARTICULAR
15 * PURPOSE. See the GNU General Public License for more details.
17 * You should have received a copy of the GNU General Public License
18 * along with Jalview. If not, see <http://www.gnu.org/licenses/>.
19 * The Jalview Authors are detailed in the 'AUTHORS' file.
23 import java.util.Locale;
25 import jalview.datamodel.Sequence;
26 import jalview.datamodel.SequenceI;
27 import jalview.util.Comparison;
28 import jalview.util.Format;
30 import java.io.IOException;
31 import java.util.ArrayList;
32 import java.util.Hashtable;
33 import java.util.List;
34 import java.util.StringTokenizer;
42 public class MSFfile extends AlignFile
46 * Creates a new MSFfile object.
53 * Creates a new MSFfile object.
63 public MSFfile(String inFile, DataSourceType type) throws IOException
68 public MSFfile(FileParse source) throws IOException
74 * Read and parse MSF sequence data
77 public void parse() throws IOException
79 boolean seqFlag = false;
80 List<String> headers = new ArrayList<String>();
81 Hashtable<String, StringBuilder> seqhash = new Hashtable<String, StringBuilder>();
86 while ((line = nextLine()) != null)
88 StringTokenizer str = new StringTokenizer(line);
91 while (str.hasMoreTokens())
93 String inStr = str.nextToken();
95 // If line has header information add to the headers vector
96 if (inStr.indexOf("Name:") != -1)
98 key = str.nextToken();
102 // if line has // set SeqFlag so we know sequences are coming
103 if (inStr.indexOf("//") != -1)
108 // Process lines as sequence lines if seqFlag is set
109 if ((inStr.indexOf("//") == -1) && seqFlag)
111 // sequence id is the first field
114 StringBuilder tempseq;
116 // Get sequence from hash if it exists
117 if (seqhash.containsKey(key))
119 tempseq = seqhash.get(key);
123 tempseq = new StringBuilder(64);
124 seqhash.put(key, tempseq);
127 // loop through the rest of the words
128 while (str.hasMoreTokens())
130 // append the word to the sequence
131 String sequenceBlock = str.nextToken();
132 tempseq.append(sequenceBlock);
137 } catch (IOException e)
139 jalview.bin.Console.errPrintln("Exception parsing MSFFile " + e);
143 this.noSeqs = headers.size();
145 // Add sequences to the hash
146 for (int i = 0; i < headers.size(); i++)
148 if (seqhash.get(headers.get(i)) != null)
150 String head = headers.get(i);
151 String seq = seqhash.get(head).toString();
153 if (maxLength < head.length())
155 maxLength = head.length();
159 * replace ~ (leading/trailing positions) with the gap character;
160 * use '.' as this is the internal gap character required by MSF
162 seq = seq.replace('~', '.');
164 Sequence newSeq = parseId(head);
166 newSeq.setSequence(seq);
168 seqs.addElement(newSeq);
173 .errPrintln("MSFFile Parser: Can't find sequence for "
185 * @return DOCUMENT ME!
187 public int checkSum(String seq)
190 String sequence = seq.toUpperCase(Locale.ROOT);
192 for (int i = 0; i < sequence.length(); i++)
197 int value = sequence.charAt(i);
200 check += (i % 57 + 1) * value;
202 } catch (Exception e)
205 .errPrintln("Exception during MSF Checksum calculation");
210 return check % 10000;
221 * @return DOCUMENT ME!
224 public String print(SequenceI[] sqs, boolean jvSuffix)
227 boolean is_NA = Comparison.isNucleotide(sqs);
229 SequenceI[] s = new SequenceI[sqs.length];
231 StringBuilder out = new StringBuilder(256);
232 out.append("!!").append(is_NA ? "NA" : "AA")
233 .append("_MULTIPLE_ALIGNMENT 1.0");
234 // TODO: JBPNote : Jalview doesn't remember NA or AA yet.
241 while ((i < sqs.length) && (sqs[i] != null))
244 * modify to MSF format: uses '.' for internal gaps,
245 * and '~' for leading or trailing gaps
247 String seqString = sqs[i].getSequenceAsString().replace('-', '.');
249 StringBuilder sb = new StringBuilder(seqString);
251 for (int ii = 0; ii < sb.length(); ii++)
253 if (sb.charAt(ii) == '.')
255 sb.setCharAt(ii, '~');
263 for (int ii = sb.length() - 1; ii > 0; ii--)
265 if (sb.charAt(ii) == '.')
267 sb.setCharAt(ii, '~');
274 s[i] = new Sequence(sqs[i].getName(), sb.toString(),
275 sqs[i].getStart(), sqs[i].getEnd());
277 if (sb.length() > max)
285 Format maxLenpad = new Format(
286 "%" + (new String("" + max)).length() + "d");
287 Format maxChkpad = new Format(
288 "%" + (new String("1" + max)).length() + "d");
292 int[] checksums = new int[s.length];
295 checksums[i] = checkSum(s[i].getSequenceAsString());
296 bigChecksum += checksums[i];
301 out.append(" MSF: " + s[0].getLength() + " Type: "
302 + (is_NA ? "N" : "P") + " Check: " + (bigChecksum % 10000)
308 String[] nameBlock = new String[s.length];
309 String[] idBlock = new String[s.length];
312 while ((i < s.length) && (s[i] != null))
315 nameBlock[i] = new String(" Name: " + printId(s[i], jvSuffix) + " ");
317 idBlock[i] = new String("Len: " + maxLenpad.form(s[i].getLength())
318 + " Check: " + maxChkpad.form(checksums[i])
319 + " Weight: 1.00" + newline);
321 if (s[i].getName().length() > maxid)
323 maxid = s[i].getName().length();
326 if (nameBlock[i].length() > maxNB)
328 maxNB = nameBlock[i].length();
344 Format nbFormat = new Format("%-" + maxNB + "s");
346 for (i = 0; (i < s.length) && (s[i] != null); i++)
348 out.append(nbFormat.form(nameBlock[i]) + idBlock[i]);
359 int nochunks = (max / len) + (max % len > 0 ? 1 : 0);
361 for (i = 0; i < nochunks; i++)
365 while ((j < s.length) && (s[j] != null))
367 String name = printId(s[j], jvSuffix);
369 out.append(new Format("%-" + maxid + "s").form(name + " "));
371 for (int k = 0; k < 5; k++)
373 int start = (i * 50) + (k * 10);
374 int end = start + 10;
376 int length = s[j].getLength();
377 if ((end < length) && (start < length))
379 out.append(s[j].getSequence(start, end));
394 out.append(s[j].getSequenceAsString().substring(start));
413 return out.toString();