2 * Jalview - A Sequence Alignment Editor and Viewer (Version 2.5)
3 * Copyright (C) 2010 J Procter, AM Waterhouse, G Barton, M Clamp, S Searle
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 of the License, or (at your option) any later version.
11 * Jalview is distributed in the hope that it will be useful, but
12 * WITHOUT ANY WARRANTY; without even the implied warranty
13 * of MERCHANTABILITY or FITNESS FOR A PARTICULAR
14 * PURPOSE. See the GNU General Public License for more details.
16 * You should have received a copy of the GNU General Public License along with Jalview. If not, see <http://www.gnu.org/licenses/>.
29 * Read and write PileUp style MSF Files. This used to be the MSFFile class, and
30 * was written according to the EBI's idea of a subset of the MSF alignment
31 * format. But, that was updated to reflect current GCG style IO fashion, as
32 * found in Emboss (thanks David Martin!)
37 import jalview.datamodel.*;
38 import jalview.util.*;
40 public class PileUpfile extends MSFfile
44 * Creates a new MSFfile object.
51 * Creates a new MSFfile object.
61 public PileUpfile(String inFile, String type) throws IOException
66 public PileUpfile(FileParse source) throws IOException
74 * @return DOCUMENT ME!
78 return print(getSeqsAsArray());
81 public String print(SequenceI[] s)
83 StringBuffer out = new StringBuffer("PileUp\n\n");
90 int[] checksums = new int[s.length];
93 checksums[i] = checkSum(s[i].getSequenceAsString());
94 bigChecksum += checksums[i];
98 out.append(" MSF: " + s[0].getSequence().length
99 + " Type: P Check: " + bigChecksum % 10000
103 while ((i < s.length) && (s[i] != null))
105 String seq = s[i].getSequenceAsString();
106 out.append(" Name: " + printId(s[i]) + " oo Len: " + seq.length()
107 + " Check: " + checksums[i] + " Weight: 1.00\n");
109 if (seq.length() > max)
114 if (s[i].getName().length() > maxid)
116 maxid = s[i].getName().length();
128 out.append("\n\n//\n\n");
132 int nochunks = (max / len) + 1;
134 if ((max % len) == 0)
139 for (i = 0; i < nochunks; i++)
143 while ((j < s.length) && (s[j] != null))
145 String name = printId(s[j]);
147 out.append(new Format("%-" + maxid + "s").form(name + " "));
149 for (int k = 0; k < 5; k++)
151 int start = (i * 50) + (k * 10);
152 int end = start + 10;
154 if ((end < s[j].getSequence().length)
155 && (start < s[j].getSequence().length))
157 out.append(s[j].getSequence(start, end));
170 if (start < s[j].getSequence().length)
172 out.append(s[j].getSequenceAsString().substring(start));
191 return out.toString();