Mercurial > repos > devteam > blast_datatypes
annotate blast.py @ 8:de11e1a921c4 draft
Uploaded v0.0.18, tweak display_data for running tests
author | peterjc |
---|---|
date | Tue, 21 Jan 2014 13:33:20 -0500 |
parents | a04cf51612f1 |
children | 2bda64d39931 |
rev | line source |
---|---|
3 | 1 """ |
2 BlastXml class | |
3 """ | |
4 | |
5 from galaxy.datatypes.data import get_file_peek | |
4 | 6 from galaxy.datatypes.data import Text, Data |
3 | 7 from galaxy.datatypes.xml import GenericXml |
4 | 8 from galaxy.datatypes.metadata import MetadataElement |
3 | 9 |
8
de11e1a921c4
Uploaded v0.0.18, tweak display_data for running tests
peterjc
parents:
6
diff
changeset
|
10 from time import sleep |
de11e1a921c4
Uploaded v0.0.18, tweak display_data for running tests
peterjc
parents:
6
diff
changeset
|
11 import os |
de11e1a921c4
Uploaded v0.0.18, tweak display_data for running tests
peterjc
parents:
6
diff
changeset
|
12 import logging |
de11e1a921c4
Uploaded v0.0.18, tweak display_data for running tests
peterjc
parents:
6
diff
changeset
|
13 |
de11e1a921c4
Uploaded v0.0.18, tweak display_data for running tests
peterjc
parents:
6
diff
changeset
|
14 log = logging.getLogger(__name__) |
5
b3a3ba0c1d47
Uploaded v0.0.15 which updates the BLAST database definitions.
peterjc
parents:
4
diff
changeset
|
15 |
3 | 16 class BlastXml( GenericXml ): |
17 """NCBI Blast XML Output data""" | |
18 file_ext = "blastxml" | |
19 | |
20 def set_peek( self, dataset, is_multi_byte=False ): | |
21 """Set the peek and blurb text""" | |
22 if not dataset.dataset.purged: | |
23 dataset.peek = get_file_peek( dataset.file_name, is_multi_byte=is_multi_byte ) | |
24 dataset.blurb = 'NCBI Blast XML data' | |
25 else: | |
26 dataset.peek = 'file does not exist' | |
27 dataset.blurb = 'file purged from disk' | |
4 | 28 |
3 | 29 def sniff( self, filename ): |
30 """ | |
31 Determines whether the file is blastxml | |
32 | |
33 >>> fname = get_test_fname( 'megablast_xml_parser_test1.blastxml' ) | |
34 >>> BlastXml().sniff( fname ) | |
35 True | |
36 >>> fname = get_test_fname( 'tblastn_four_human_vs_rhodopsin.xml' ) | |
37 >>> BlastXml().sniff( fname ) | |
38 True | |
39 >>> fname = get_test_fname( 'interval.interval' ) | |
40 >>> BlastXml().sniff( fname ) | |
41 False | |
42 """ | |
43 #TODO - Use a context manager on Python 2.5+ to close handle | |
44 handle = open(filename) | |
45 line = handle.readline() | |
46 if line.strip() != '<?xml version="1.0"?>': | |
47 handle.close() | |
48 return False | |
49 line = handle.readline() | |
50 if line.strip() not in ['<!DOCTYPE BlastOutput PUBLIC "-//NCBI//NCBI BlastOutput/EN" "http://www.ncbi.nlm.nih.gov/dtd/NCBI_BlastOutput.dtd">', | |
51 '<!DOCTYPE BlastOutput PUBLIC "-//NCBI//NCBI BlastOutput/EN" "NCBI_BlastOutput.dtd">']: | |
52 handle.close() | |
53 return False | |
54 line = handle.readline() | |
55 if line.strip() != '<BlastOutput>': | |
56 handle.close() | |
57 return False | |
58 handle.close() | |
59 return True | |
60 | |
61 def merge(split_files, output_file): | |
62 """Merging multiple XML files is non-trivial and must be done in subclasses.""" | |
63 if len(split_files) == 1: | |
64 #For one file only, use base class method (move/copy) | |
65 return Text.merge(split_files, output_file) | |
4 | 66 if not split_files: |
67 raise ValueError("Given no BLAST XML files, %r, to merge into %s" \ | |
68 % (split_files, output_file)) | |
3 | 69 out = open(output_file, "w") |
70 h = None | |
71 for f in split_files: | |
8
de11e1a921c4
Uploaded v0.0.18, tweak display_data for running tests
peterjc
parents:
6
diff
changeset
|
72 if not os.path.isfile(f): |
de11e1a921c4
Uploaded v0.0.18, tweak display_data for running tests
peterjc
parents:
6
diff
changeset
|
73 log.warning("BLAST XML file %s missing, retry in 1s..." % f) |
de11e1a921c4
Uploaded v0.0.18, tweak display_data for running tests
peterjc
parents:
6
diff
changeset
|
74 sleep(1) |
de11e1a921c4
Uploaded v0.0.18, tweak display_data for running tests
peterjc
parents:
6
diff
changeset
|
75 if not os.path.isfile(f): |
de11e1a921c4
Uploaded v0.0.18, tweak display_data for running tests
peterjc
parents:
6
diff
changeset
|
76 log.error("BLAST XML file %s missing" % f) |
de11e1a921c4
Uploaded v0.0.18, tweak display_data for running tests
peterjc
parents:
6
diff
changeset
|
77 raise ValueError("BLAST XML file %s missing" % f) |
3 | 78 h = open(f) |
79 body = False | |
80 header = h.readline() | |
81 if not header: | |
82 out.close() | |
83 h.close() | |
8
de11e1a921c4
Uploaded v0.0.18, tweak display_data for running tests
peterjc
parents:
6
diff
changeset
|
84 #Retry, could be transient error with networked file system... |
de11e1a921c4
Uploaded v0.0.18, tweak display_data for running tests
peterjc
parents:
6
diff
changeset
|
85 log.warning("BLAST XML file %s empty, retry in 1s..." % f) |
de11e1a921c4
Uploaded v0.0.18, tweak display_data for running tests
peterjc
parents:
6
diff
changeset
|
86 sleep(1) |
de11e1a921c4
Uploaded v0.0.18, tweak display_data for running tests
peterjc
parents:
6
diff
changeset
|
87 h = open(f) |
de11e1a921c4
Uploaded v0.0.18, tweak display_data for running tests
peterjc
parents:
6
diff
changeset
|
88 header = h.readline() |
de11e1a921c4
Uploaded v0.0.18, tweak display_data for running tests
peterjc
parents:
6
diff
changeset
|
89 if not header: |
de11e1a921c4
Uploaded v0.0.18, tweak display_data for running tests
peterjc
parents:
6
diff
changeset
|
90 log.error("BLAST XML file %s was empty" % f) |
de11e1a921c4
Uploaded v0.0.18, tweak display_data for running tests
peterjc
parents:
6
diff
changeset
|
91 raise ValueError("BLAST XML file %s was empty" % f) |
3 | 92 if header.strip() != '<?xml version="1.0"?>': |
93 out.write(header) #for diagnosis | |
94 out.close() | |
95 h.close() | |
96 raise ValueError("%s is not an XML file!" % f) | |
97 line = h.readline() | |
98 header += line | |
99 if line.strip() not in ['<!DOCTYPE BlastOutput PUBLIC "-//NCBI//NCBI BlastOutput/EN" "http://www.ncbi.nlm.nih.gov/dtd/NCBI_BlastOutput.dtd">', | |
100 '<!DOCTYPE BlastOutput PUBLIC "-//NCBI//NCBI BlastOutput/EN" "NCBI_BlastOutput.dtd">']: | |
101 out.write(header) #for diagnosis | |
102 out.close() | |
103 h.close() | |
104 raise ValueError("%s is not a BLAST XML file!" % f) | |
105 while True: | |
106 line = h.readline() | |
107 if not line: | |
108 out.write(header) #for diagnosis | |
109 out.close() | |
110 h.close() | |
111 raise ValueError("BLAST XML file %s ended prematurely" % f) | |
112 header += line | |
113 if "<Iteration>" in line: | |
114 break | |
115 if len(header) > 10000: | |
116 #Something has gone wrong, don't load too much into memory! | |
117 #Write what we have to the merged file for diagnostics | |
118 out.write(header) | |
119 out.close() | |
120 h.close() | |
121 raise ValueError("BLAST XML file %s has too long a header!" % f) | |
122 if "<BlastOutput>" not in header: | |
123 out.close() | |
124 h.close() | |
125 raise ValueError("%s is not a BLAST XML file:\n%s\n..." % (f, header)) | |
126 if f == split_files[0]: | |
127 out.write(header) | |
128 old_header = header | |
129 elif old_header[:300] != header[:300]: | |
130 #Enough to check <BlastOutput_program> and <BlastOutput_version> match | |
131 out.close() | |
132 h.close() | |
133 raise ValueError("BLAST XML headers don't match for %s and %s - have:\n%s\n...\n\nAnd:\n%s\n...\n" \ | |
134 % (split_files[0], f, old_header[:300], header[:300])) | |
135 else: | |
136 out.write(" <Iteration>\n") | |
137 for line in h: | |
138 if "</BlastOutput_iterations>" in line: | |
139 break | |
140 #TODO - Increment <Iteration_iter-num> and if required automatic query names | |
141 #like <Iteration_query-ID>Query_3</Iteration_query-ID> to be increasing? | |
142 out.write(line) | |
143 h.close() | |
144 out.write(" </BlastOutput_iterations>\n") | |
145 out.write("</BlastOutput>\n") | |
146 out.close() | |
147 merge = staticmethod(merge) | |
148 | |
4 | 149 |
150 class _BlastDb(object): | |
151 """Base class for BLAST database datatype.""" | |
152 | |
153 def set_peek( self, dataset, is_multi_byte=False ): | |
154 """Set the peek and blurb text.""" | |
155 if not dataset.dataset.purged: | |
156 dataset.peek = "BLAST database (multiple files)" | |
157 dataset.blurb = "BLAST database (multiple files)" | |
158 else: | |
159 dataset.peek = 'file does not exist' | |
160 dataset.blurb = 'file purged from disk' | |
161 | |
162 def display_peek( self, dataset ): | |
163 """Create HTML content, used for displaying peek.""" | |
164 try: | |
165 return dataset.peek | |
166 except: | |
167 return "BLAST database (multiple files)" | |
168 | |
169 def display_data(self, trans, data, preview=False, filename=None, | |
170 to_ext=None, size=None, offset=None, **kwd): | |
8
de11e1a921c4
Uploaded v0.0.18, tweak display_data for running tests
peterjc
parents:
6
diff
changeset
|
171 """Documented as an old display method, but still gets called via tests etc |
4 | 172 |
173 This allows us to format the data shown in the central pane via the "eye" icon. | |
174 """ | |
8
de11e1a921c4
Uploaded v0.0.18, tweak display_data for running tests
peterjc
parents:
6
diff
changeset
|
175 if filename is not None and filename != "index": |
de11e1a921c4
Uploaded v0.0.18, tweak display_data for running tests
peterjc
parents:
6
diff
changeset
|
176 #Change nothing - important for the unit tests to access child files: |
de11e1a921c4
Uploaded v0.0.18, tweak display_data for running tests
peterjc
parents:
6
diff
changeset
|
177 return Data.display_data(self, trans, data, preview, filename, |
de11e1a921c4
Uploaded v0.0.18, tweak display_data for running tests
peterjc
parents:
6
diff
changeset
|
178 to_ext, size, offset, **kwd) |
de11e1a921c4
Uploaded v0.0.18, tweak display_data for running tests
peterjc
parents:
6
diff
changeset
|
179 if self.file_ext == "blastdbn": |
de11e1a921c4
Uploaded v0.0.18, tweak display_data for running tests
peterjc
parents:
6
diff
changeset
|
180 title = "This is a nucleotide BLAST database" |
de11e1a921c4
Uploaded v0.0.18, tweak display_data for running tests
peterjc
parents:
6
diff
changeset
|
181 elif self.file_ext =="blastdbp": |
de11e1a921c4
Uploaded v0.0.18, tweak display_data for running tests
peterjc
parents:
6
diff
changeset
|
182 title = "This is a protein BLAST database" |
de11e1a921c4
Uploaded v0.0.18, tweak display_data for running tests
peterjc
parents:
6
diff
changeset
|
183 else: |
de11e1a921c4
Uploaded v0.0.18, tweak display_data for running tests
peterjc
parents:
6
diff
changeset
|
184 #Error? |
de11e1a921c4
Uploaded v0.0.18, tweak display_data for running tests
peterjc
parents:
6
diff
changeset
|
185 title = "This is a BLAST database." |
de11e1a921c4
Uploaded v0.0.18, tweak display_data for running tests
peterjc
parents:
6
diff
changeset
|
186 msg = "" |
de11e1a921c4
Uploaded v0.0.18, tweak display_data for running tests
peterjc
parents:
6
diff
changeset
|
187 try: |
de11e1a921c4
Uploaded v0.0.18, tweak display_data for running tests
peterjc
parents:
6
diff
changeset
|
188 #Try to use any text recorded in the dummy index file: |
de11e1a921c4
Uploaded v0.0.18, tweak display_data for running tests
peterjc
parents:
6
diff
changeset
|
189 handle = open(data.file_name, "rU") |
de11e1a921c4
Uploaded v0.0.18, tweak display_data for running tests
peterjc
parents:
6
diff
changeset
|
190 msg = handle.read().strip() |
de11e1a921c4
Uploaded v0.0.18, tweak display_data for running tests
peterjc
parents:
6
diff
changeset
|
191 handle.close() |
de11e1a921c4
Uploaded v0.0.18, tweak display_data for running tests
peterjc
parents:
6
diff
changeset
|
192 except Exception, err: |
de11e1a921c4
Uploaded v0.0.18, tweak display_data for running tests
peterjc
parents:
6
diff
changeset
|
193 #msg = str(err) |
de11e1a921c4
Uploaded v0.0.18, tweak display_data for running tests
peterjc
parents:
6
diff
changeset
|
194 pass |
de11e1a921c4
Uploaded v0.0.18, tweak display_data for running tests
peterjc
parents:
6
diff
changeset
|
195 if not msg: |
de11e1a921c4
Uploaded v0.0.18, tweak display_data for running tests
peterjc
parents:
6
diff
changeset
|
196 msg = title |
de11e1a921c4
Uploaded v0.0.18, tweak display_data for running tests
peterjc
parents:
6
diff
changeset
|
197 #Galaxy assumes HTML for the display of composite datatypes, |
de11e1a921c4
Uploaded v0.0.18, tweak display_data for running tests
peterjc
parents:
6
diff
changeset
|
198 return "<html><head><title>%s</title></head><body><pre>%s</pre></body></html>" % (title, msg) |
4 | 199 |
200 def merge(split_files, output_file): | |
201 """Merge BLAST databases (not implemented for now).""" | |
202 raise NotImplementedError("Merging BLAST databases is non-trivial (do this via makeblastdb?)") | |
203 | |
204 def split( cls, input_datasets, subdir_generator_function, split_params): | |
205 """Split a BLAST database (not implemented for now).""" | |
206 if split_params is None: | |
207 return None | |
208 raise NotImplementedError("Can't split BLAST databases") | |
209 | |
210 | |
211 class BlastNucDb( _BlastDb, Data ): | |
212 """Class for nucleotide BLAST database files.""" | |
213 file_ext = 'blastdbn' | |
8
de11e1a921c4
Uploaded v0.0.18, tweak display_data for running tests
peterjc
parents:
6
diff
changeset
|
214 allow_datatype_change = False |
de11e1a921c4
Uploaded v0.0.18, tweak display_data for running tests
peterjc
parents:
6
diff
changeset
|
215 composite_type = 'basic' |
4 | 216 |
5
b3a3ba0c1d47
Uploaded v0.0.15 which updates the BLAST database definitions.
peterjc
parents:
4
diff
changeset
|
217 def __init__(self, **kwd): |
4 | 218 Data.__init__(self, **kwd) |
5
b3a3ba0c1d47
Uploaded v0.0.15 which updates the BLAST database definitions.
peterjc
parents:
4
diff
changeset
|
219 self.add_composite_file('blastdb.nhr', is_binary=True) # sequence headers |
b3a3ba0c1d47
Uploaded v0.0.15 which updates the BLAST database definitions.
peterjc
parents:
4
diff
changeset
|
220 self.add_composite_file('blastdb.nin', is_binary=True) # index file |
b3a3ba0c1d47
Uploaded v0.0.15 which updates the BLAST database definitions.
peterjc
parents:
4
diff
changeset
|
221 self.add_composite_file('blastdb.nsq', is_binary=True) # nucleotide sequences |
b3a3ba0c1d47
Uploaded v0.0.15 which updates the BLAST database definitions.
peterjc
parents:
4
diff
changeset
|
222 self.add_composite_file('blastdb.nal', is_binary=False, optional=True) # alias ( -gi_mask option of makeblastdb) |
b3a3ba0c1d47
Uploaded v0.0.15 which updates the BLAST database definitions.
peterjc
parents:
4
diff
changeset
|
223 self.add_composite_file('blastdb.nhd', is_binary=True, optional=True) # sorted sequence hash values ( -hash_index option of makeblastdb) |
b3a3ba0c1d47
Uploaded v0.0.15 which updates the BLAST database definitions.
peterjc
parents:
4
diff
changeset
|
224 self.add_composite_file('blastdb.nhi', is_binary=True, optional=True) # index of sequence hash values ( -hash_index option of makeblastdb) |
b3a3ba0c1d47
Uploaded v0.0.15 which updates the BLAST database definitions.
peterjc
parents:
4
diff
changeset
|
225 self.add_composite_file('blastdb.nnd', is_binary=True, optional=True) # sorted GI values ( -parse_seqids option of makeblastdb and gi present in the description lines) |
b3a3ba0c1d47
Uploaded v0.0.15 which updates the BLAST database definitions.
peterjc
parents:
4
diff
changeset
|
226 self.add_composite_file('blastdb.nni', is_binary=True, optional=True) # index of GI values ( -parse_seqids option of makeblastdb and gi present in the description lines) |
b3a3ba0c1d47
Uploaded v0.0.15 which updates the BLAST database definitions.
peterjc
parents:
4
diff
changeset
|
227 self.add_composite_file('blastdb.nog', is_binary=True, optional=True) # OID->GI lookup file ( -hash_index or -parse_seqids option of makeblastdb) |
b3a3ba0c1d47
Uploaded v0.0.15 which updates the BLAST database definitions.
peterjc
parents:
4
diff
changeset
|
228 self.add_composite_file('blastdb.nsd', is_binary=True, optional=True) # sorted sequence accession values ( -hash_index or -parse_seqids option of makeblastdb) |
b3a3ba0c1d47
Uploaded v0.0.15 which updates the BLAST database definitions.
peterjc
parents:
4
diff
changeset
|
229 self.add_composite_file('blastdb.nsi', is_binary=True, optional=True) # index of sequence accession values ( -hash_index or -parse_seqids option of makeblastdb) |
6
a04cf51612f1
Uploaded v0.0.16, MIT License, development moved to GitHub, nucleotide database definition aware of MegaBLAST index superheader
peterjc
parents:
5
diff
changeset
|
230 # self.add_composite_file('blastdb.00.idx', is_binary=True, optional=True) # first volume of the MegaBLAST index generated by makembindex |
a04cf51612f1
Uploaded v0.0.16, MIT License, development moved to GitHub, nucleotide database definition aware of MegaBLAST index superheader
peterjc
parents:
5
diff
changeset
|
231 # The previous line should be repeated for each index volume, with filename extensions like '.01.idx', '.02.idx', etc. |
a04cf51612f1
Uploaded v0.0.16, MIT License, development moved to GitHub, nucleotide database definition aware of MegaBLAST index superheader
peterjc
parents:
5
diff
changeset
|
232 self.add_composite_file('blastdb.shd', is_binary=True, optional=True) # MegaBLAST index superheader (-old_style_index false option of makembindex) |
5
b3a3ba0c1d47
Uploaded v0.0.15 which updates the BLAST database definitions.
peterjc
parents:
4
diff
changeset
|
233 # self.add_composite_file('blastdb.naa', is_binary=True, optional=True) # index of a WriteDB column for e.g. mask data |
b3a3ba0c1d47
Uploaded v0.0.15 which updates the BLAST database definitions.
peterjc
parents:
4
diff
changeset
|
234 # self.add_composite_file('blastdb.nab', is_binary=True, optional=True) # data of a WriteDB column |
b3a3ba0c1d47
Uploaded v0.0.15 which updates the BLAST database definitions.
peterjc
parents:
4
diff
changeset
|
235 # self.add_composite_file('blastdb.nac', is_binary=True, optional=True) # multiple byte order for a WriteDB column |
6
a04cf51612f1
Uploaded v0.0.16, MIT License, development moved to GitHub, nucleotide database definition aware of MegaBLAST index superheader
peterjc
parents:
5
diff
changeset
|
236 # The previous 3 lines should be repeated for each WriteDB column, with filename extensions like ('.nba', '.nbb', '.nbc'), ('.nca', '.ncb', '.ncc'), etc. |
4 | 237 |
5
b3a3ba0c1d47
Uploaded v0.0.15 which updates the BLAST database definitions.
peterjc
parents:
4
diff
changeset
|
238 |
4 | 239 class BlastProtDb( _BlastDb, Data ): |
240 """Class for protein BLAST database files.""" | |
241 file_ext = 'blastdbp' | |
8
de11e1a921c4
Uploaded v0.0.18, tweak display_data for running tests
peterjc
parents:
6
diff
changeset
|
242 allow_datatype_change = False |
de11e1a921c4
Uploaded v0.0.18, tweak display_data for running tests
peterjc
parents:
6
diff
changeset
|
243 composite_type = 'basic' |
4 | 244 |
5
b3a3ba0c1d47
Uploaded v0.0.15 which updates the BLAST database definitions.
peterjc
parents:
4
diff
changeset
|
245 def __init__(self, **kwd): |
4 | 246 Data.__init__(self, **kwd) |
5
b3a3ba0c1d47
Uploaded v0.0.15 which updates the BLAST database definitions.
peterjc
parents:
4
diff
changeset
|
247 # Component file comments are as in BlastNucDb except where noted |
b3a3ba0c1d47
Uploaded v0.0.15 which updates the BLAST database definitions.
peterjc
parents:
4
diff
changeset
|
248 self.add_composite_file('blastdb.phr', is_binary=True) |
b3a3ba0c1d47
Uploaded v0.0.15 which updates the BLAST database definitions.
peterjc
parents:
4
diff
changeset
|
249 self.add_composite_file('blastdb.pin', is_binary=True) |
b3a3ba0c1d47
Uploaded v0.0.15 which updates the BLAST database definitions.
peterjc
parents:
4
diff
changeset
|
250 self.add_composite_file('blastdb.psq', is_binary=True) # protein sequences |
b3a3ba0c1d47
Uploaded v0.0.15 which updates the BLAST database definitions.
peterjc
parents:
4
diff
changeset
|
251 self.add_composite_file('blastdb.phd', is_binary=True, optional=True) |
b3a3ba0c1d47
Uploaded v0.0.15 which updates the BLAST database definitions.
peterjc
parents:
4
diff
changeset
|
252 self.add_composite_file('blastdb.phi', is_binary=True, optional=True) |
b3a3ba0c1d47
Uploaded v0.0.15 which updates the BLAST database definitions.
peterjc
parents:
4
diff
changeset
|
253 self.add_composite_file('blastdb.pnd', is_binary=True, optional=True) |
b3a3ba0c1d47
Uploaded v0.0.15 which updates the BLAST database definitions.
peterjc
parents:
4
diff
changeset
|
254 self.add_composite_file('blastdb.pni', is_binary=True, optional=True) |
b3a3ba0c1d47
Uploaded v0.0.15 which updates the BLAST database definitions.
peterjc
parents:
4
diff
changeset
|
255 self.add_composite_file('blastdb.pog', is_binary=True, optional=True) |
b3a3ba0c1d47
Uploaded v0.0.15 which updates the BLAST database definitions.
peterjc
parents:
4
diff
changeset
|
256 self.add_composite_file('blastdb.psd', is_binary=True, optional=True) |
b3a3ba0c1d47
Uploaded v0.0.15 which updates the BLAST database definitions.
peterjc
parents:
4
diff
changeset
|
257 self.add_composite_file('blastdb.psi', is_binary=True, optional=True) |
b3a3ba0c1d47
Uploaded v0.0.15 which updates the BLAST database definitions.
peterjc
parents:
4
diff
changeset
|
258 # self.add_composite_file('blastdb.paa', is_binary=True, optional=True) |
b3a3ba0c1d47
Uploaded v0.0.15 which updates the BLAST database definitions.
peterjc
parents:
4
diff
changeset
|
259 # self.add_composite_file('blastdb.pab', is_binary=True, optional=True) |
b3a3ba0c1d47
Uploaded v0.0.15 which updates the BLAST database definitions.
peterjc
parents:
4
diff
changeset
|
260 # self.add_composite_file('blastdb.pac', is_binary=True, optional=True) |
b3a3ba0c1d47
Uploaded v0.0.15 which updates the BLAST database definitions.
peterjc
parents:
4
diff
changeset
|
261 # The last 3 lines should be repeated for each WriteDB column, with filename extensions like ('.pba', '.pbb', '.pbc'), ('.pca', '.pcb', '.pcc'), etc. |