line |
stmt |
bran |
cond |
sub |
pod |
time |
code |
1
|
|
|
|
|
|
|
# $Id: FASTAwriter.pm,v 1.1 2004/04/28 15:03:43 aphillip Exp $ |
2
|
|
|
|
|
|
|
|
3
|
|
|
|
|
|
|
package TIGR::FASTA::Writer; |
4
|
|
|
|
|
|
|
{ |
5
|
|
|
|
|
|
|
|
6
|
|
|
|
|
|
|
=head1 NAME |
7
|
|
|
|
|
|
|
|
8
|
|
|
|
|
|
|
TIGR::FASTA::Writer - TIGR::FASTA::Writer class for writing TIGR::FASTA::Record |
9
|
|
|
|
|
|
|
objects to a file |
10
|
|
|
|
|
|
|
|
11
|
|
|
|
|
|
|
=head1 SYNOPSIS |
12
|
|
|
|
|
|
|
|
13
|
|
|
|
|
|
|
use TIGR::FASTA:Writer; |
14
|
|
|
|
|
|
|
my $obj_instance = new TIGR::FASTA::Writer ($tigr_foundation_obj, |
15
|
|
|
|
|
|
|
$output_file_name); |
16
|
|
|
|
|
|
|
|
17
|
|
|
|
|
|
|
=head1 DESCRIPTION |
18
|
|
|
|
|
|
|
|
19
|
|
|
|
|
|
|
This module provides an object definition for a TIGR::FASTA::Writer. |
20
|
|
|
|
|
|
|
The TIGR::FASTA::Writer object accepts TIGR::FASTA::Record objects for |
21
|
|
|
|
|
|
|
printing to an output file. |
22
|
|
|
|
|
|
|
|
23
|
|
|
|
|
|
|
=cut |
24
|
|
|
|
|
|
|
|
25
|
|
|
|
|
|
|
BEGIN { |
26
|
1
|
|
|
1
|
|
1508
|
require 5.006_00; |
27
|
|
|
|
|
|
|
} |
28
|
|
|
|
|
|
|
|
29
|
1
|
|
|
1
|
|
4
|
use strict; |
|
1
|
|
|
|
|
2
|
|
|
1
|
|
|
|
|
19
|
|
30
|
1
|
|
|
1
|
|
4
|
use IO::File; |
|
1
|
|
|
|
|
1
|
|
|
1
|
|
|
|
|
142
|
|
31
|
1
|
|
|
1
|
|
9
|
use TIGR::FASTA::Record; |
|
1
|
|
|
|
|
3
|
|
|
1
|
|
|
|
|
844
|
|
32
|
|
|
|
|
|
|
|
33
|
|
|
|
|
|
|
|
34
|
|
|
|
|
|
|
|
35
|
|
|
|
|
|
|
## internal variables and identifiers |
36
|
|
|
|
|
|
|
|
37
|
|
|
|
|
|
|
our $REVISION = (qw$Revision: 1.1 $)[-1]; |
38
|
|
|
|
|
|
|
our $VERSION = '1.0'; |
39
|
|
|
|
|
|
|
our $VERSION_STRING = "$VERSION (Build $REVISION)"; |
40
|
|
|
|
|
|
|
our @DEPEND = (); |
41
|
|
|
|
|
|
|
|
42
|
|
|
|
|
|
|
my $SYS_ERR = 0; # this flag specifies non-user related error |
43
|
|
|
|
|
|
|
my $USR_ERR = 1; # this flag specifies user related error |
44
|
|
|
|
|
|
|
|
45
|
|
|
|
|
|
|
#debugging scheme |
46
|
|
|
|
|
|
|
# |
47
|
|
|
|
|
|
|
# Debugging via the TIGR Foundation uses increasing log levels based on |
48
|
|
|
|
|
|
|
# nesting. 'MAIN' starts at level 1. Every nest increments the level by 1. |
49
|
|
|
|
|
|
|
# Subroutines always start nesting at level 2. As debugging levels |
50
|
|
|
|
|
|
|
# increase, logging is more verbose. This makes sense as you log at |
51
|
|
|
|
|
|
|
# greater depth (ie. deeper branching). |
52
|
|
|
|
|
|
|
# |
53
|
|
|
|
|
|
|
# The following definitions help emphasize the debugging in the program. |
54
|
|
|
|
|
|
|
# |
55
|
|
|
|
|
|
|
my $DEBUG_LEVEL_1 = 1; |
56
|
|
|
|
|
|
|
my $DEBUG_LEVEL_2 = 2; |
57
|
|
|
|
|
|
|
my $DEBUG_LEVEL_3 = 3; |
58
|
|
|
|
|
|
|
my $DEBUG_LEVEL_4 = 4; |
59
|
|
|
|
|
|
|
my $DEBUG_LEVEL_5 = 5; |
60
|
|
|
|
|
|
|
my $DEBUG_LEVEL_6 = 6; |
61
|
|
|
|
|
|
|
my $DEBUG_LEVEL_7 = 7; |
62
|
|
|
|
|
|
|
my $DEBUG_LEVEL_8 = 8; |
63
|
|
|
|
|
|
|
my $DEBUG_LEVEL_9 = 9; |
64
|
|
|
|
|
|
|
|
65
|
|
|
|
|
|
|
## prototypes |
66
|
|
|
|
|
|
|
|
67
|
|
|
|
|
|
|
sub new(;$$$); |
68
|
|
|
|
|
|
|
sub open($;$); |
69
|
|
|
|
|
|
|
sub close(); |
70
|
|
|
|
|
|
|
sub write($); |
71
|
|
|
|
|
|
|
sub _errorHandler($$$); |
72
|
|
|
|
|
|
|
|
73
|
|
|
|
|
|
|
## implementation |
74
|
|
|
|
|
|
|
|
75
|
|
|
|
|
|
|
=over |
76
|
|
|
|
|
|
|
|
77
|
|
|
|
|
|
|
=item $obj_instance = new TIGR::FASTA::Writer ($foundation_object, |
78
|
|
|
|
|
|
|
$error_array_ref, $output_file); |
79
|
|
|
|
|
|
|
|
80
|
|
|
|
|
|
|
This method returns a new instance of a TIGR::FASTA::Writer object. It takes |
81
|
|
|
|
|
|
|
three optional parameters: a TIGR::Foundation object (C<$foundation_object>), |
82
|
|
|
|
|
|
|
a reference to an array for logging user error messages (C<$error_array_ref>), |
83
|
|
|
|
|
|
|
and an output file name, C<$output_file>, as parameters. A new object instance |
84
|
|
|
|
|
|
|
is returned on success and successful opening of a specified output |
85
|
|
|
|
|
|
|
file. |
86
|
|
|
|
|
|
|
If the file supplied cannot be opened, this method returns undefined. |
87
|
|
|
|
|
|
|
This method also returns undefined if the parameters supplied are invalid. |
88
|
|
|
|
|
|
|
Writing errors are written to the array at C<$error_array_ref> and the |
89
|
|
|
|
|
|
|
log file. |
90
|
|
|
|
|
|
|
|
91
|
|
|
|
|
|
|
=cut |
92
|
|
|
|
|
|
|
|
93
|
|
|
|
|
|
|
|
94
|
|
|
|
|
|
|
sub new(;$$$) { |
95
|
|
|
|
|
|
|
|
96
|
0
|
|
|
0
|
1
|
|
my $pkg = shift; |
97
|
0
|
|
|
|
|
|
my @method_args = @_; |
98
|
|
|
|
|
|
|
|
99
|
0
|
|
|
|
|
|
my $error_condition = 0; |
100
|
0
|
|
|
|
|
|
my $self = {}; |
101
|
0
|
|
|
|
|
|
bless $self, $pkg; |
102
|
|
|
|
|
|
|
|
103
|
0
|
0
|
0
|
|
|
|
if ( ( scalar (@method_args) > 0 ) && |
104
|
|
|
|
|
|
|
( ( ref ($method_args[0]) ) =~ /foundation/i ) ) { |
105
|
0
|
|
|
|
|
|
$self->{foundation} = shift @method_args; |
106
|
0
|
|
|
|
|
|
$self->_errorHandler("Got TIGR::Foundation in new()", $DEBUG_LEVEL_3, |
107
|
|
|
|
|
|
|
$SYS_ERR); |
108
|
|
|
|
|
|
|
} |
109
|
|
|
|
|
|
|
else { |
110
|
0
|
|
|
|
|
|
$self->{foundation} = undef; |
111
|
0
|
|
|
|
|
|
$self->_errorHandler("No TIGR::Foundation in new()", $DEBUG_LEVEL_3, |
112
|
|
|
|
|
|
|
$SYS_ERR); |
113
|
|
|
|
|
|
|
} |
114
|
|
|
|
|
|
|
|
115
|
0
|
0
|
0
|
|
|
|
if ( ( scalar (@method_args) > 0 ) && |
116
|
|
|
|
|
|
|
( ( ref ($method_args[0]) ) =~ /array/i ) ) { |
117
|
0
|
|
|
|
|
|
$self->{error_ref} = shift @method_args; |
118
|
|
|
|
|
|
|
|
119
|
0
|
|
|
|
|
|
$self->_errorHandler("Got Error ARRAY in new()", $DEBUG_LEVEL_3, |
120
|
|
|
|
|
|
|
$SYS_ERR); |
121
|
|
|
|
|
|
|
} |
122
|
|
|
|
|
|
|
else { |
123
|
0
|
|
|
|
|
|
$self->{error_ref} = undef; |
124
|
0
|
|
|
|
|
|
$self->_errorHandler("No Error ARRAY in new()", $DEBUG_LEVEL_3, |
125
|
|
|
|
|
|
|
$SYS_ERR); |
126
|
|
|
|
|
|
|
} |
127
|
|
|
|
|
|
|
|
128
|
0
|
0
|
0
|
|
|
|
if ( ( scalar (@method_args) > 0 ) && |
129
|
|
|
|
|
|
|
( ! ref ($method_args[0]) ) ) { |
130
|
0
|
|
|
|
|
|
my $filename = shift @method_args; |
131
|
0
|
0
|
|
|
|
|
if(defined($filename)) { |
132
|
0
|
|
|
|
|
|
$self->{db_file_name} = $filename ; |
133
|
0
|
|
|
|
|
|
$self->_errorHandler("Got file name in new()", $DEBUG_LEVEL_4, |
134
|
|
|
|
|
|
|
$SYS_ERR); |
135
|
|
|
|
|
|
|
} |
136
|
|
|
|
|
|
|
else { |
137
|
0
|
|
|
|
|
|
$self->{db_file_name} = "" ; |
138
|
0
|
|
|
|
|
|
$self->_errorHandler("undef passed as filename", $DEBUG_LEVEL_4, |
139
|
|
|
|
|
|
|
$USR_ERR); |
140
|
|
|
|
|
|
|
} |
141
|
|
|
|
|
|
|
} |
142
|
|
|
|
|
|
|
else { |
143
|
0
|
|
|
|
|
|
$self->{db_file_name} = undef; |
144
|
0
|
|
|
|
|
|
$self->_errorHandler("No file name in new()", $DEBUG_LEVEL_3, |
145
|
|
|
|
|
|
|
$SYS_ERR); |
146
|
|
|
|
|
|
|
} |
147
|
|
|
|
|
|
|
|
148
|
|
|
|
|
|
|
# check for invocation errors |
149
|
0
|
0
|
0
|
|
|
|
if ( ( scalar (@method_args) > 0 ) ) { |
|
|
0
|
|
|
|
|
|
150
|
0
|
|
|
|
|
|
$error_condition = 1; |
151
|
0
|
|
|
|
|
|
$self->_errorHandler("Too many parameters passed to new() method", |
152
|
|
|
|
|
|
|
$DEBUG_LEVEL_3, $SYS_ERR); |
153
|
|
|
|
|
|
|
} |
154
|
|
|
|
|
|
|
elsif ((defined ( $self->{db_file_name} )) && |
155
|
|
|
|
|
|
|
(! defined ( $self->open($self->{db_file_name}, "w") ) )) { |
156
|
|
|
|
|
|
|
# the error message is logged via the open() routine |
157
|
0
|
|
|
|
|
|
$self = undef; |
158
|
|
|
|
|
|
|
} |
159
|
|
|
|
|
|
|
|
160
|
0
|
0
|
|
|
|
|
return ( $error_condition == 0 ) ? $self : undef; |
161
|
|
|
|
|
|
|
|
162
|
|
|
|
|
|
|
} |
163
|
|
|
|
|
|
|
|
164
|
|
|
|
|
|
|
|
165
|
|
|
|
|
|
|
=item $result = $obj_instance->open($file_name, $flag); |
166
|
|
|
|
|
|
|
|
167
|
|
|
|
|
|
|
This method opens a FASTA file for writing or appending. The file, |
168
|
|
|
|
|
|
|
F<$file_name>, is opened using the C flags specified by C<$flag>. |
169
|
|
|
|
|
|
|
Supported flags include 'w' and 'a'. On success, this method returns 1. |
170
|
|
|
|
|
|
|
The default C method is 'w', or truncated open. If the file cannot |
171
|
|
|
|
|
|
|
be opened, this method returns undefined. |
172
|
|
|
|
|
|
|
|
173
|
|
|
|
|
|
|
=cut |
174
|
|
|
|
|
|
|
|
175
|
|
|
|
|
|
|
|
176
|
|
|
|
|
|
|
sub open($;$) { |
177
|
0
|
|
|
0
|
1
|
|
my $self = shift; |
178
|
0
|
|
|
|
|
|
my $db_file_name = shift; |
179
|
0
|
|
|
|
|
|
my $open_flags = shift; |
180
|
|
|
|
|
|
|
|
181
|
0
|
|
|
|
|
|
my $error_condition = 0; |
182
|
|
|
|
|
|
|
|
183
|
0
|
0
|
0
|
|
|
|
if ( ( ! defined ($open_flags) ) || |
|
|
|
0
|
|
|
|
|
184
|
|
|
|
|
|
|
( ( $open_flags !~ /^w$/i ) && |
185
|
|
|
|
|
|
|
( $open_flags !~ /^a$/i ) ) ) { |
186
|
0
|
|
|
|
|
|
$open_flags = "w"; |
187
|
|
|
|
|
|
|
} |
188
|
|
|
|
|
|
|
|
189
|
|
|
|
|
|
|
# close a previously open file |
190
|
0
|
0
|
|
|
|
|
if ( defined ($self->{db_handle}) ) { |
191
|
0
|
|
|
|
|
|
$self->close(); |
192
|
|
|
|
|
|
|
} |
193
|
|
|
|
|
|
|
|
194
|
0
|
0
|
0
|
|
|
|
if (!( |
|
|
|
0
|
|
|
|
|
|
|
|
0
|
|
|
|
|
195
|
|
|
|
|
|
|
( defined ( $db_file_name ) ) && |
196
|
|
|
|
|
|
|
( $self->{db_file_name} = $db_file_name ) && |
197
|
|
|
|
|
|
|
( defined ( $self->{db_file_name} )) && |
198
|
|
|
|
|
|
|
( defined ( $self->{db_handle} = |
199
|
|
|
|
|
|
|
new IO::File $self->{db_file_name}, $open_flags )) |
200
|
|
|
|
|
|
|
) ) { |
201
|
0
|
|
|
|
|
|
$error_condition = 1; |
202
|
0
|
|
|
|
|
|
$self->_errorHandler( |
203
|
|
|
|
|
|
|
"Cannot open file \'$self->{db_file_name}\'", $DEBUG_LEVEL_3, |
204
|
|
|
|
|
|
|
$USR_ERR); |
205
|
|
|
|
|
|
|
} |
206
|
0
|
0
|
|
|
|
|
return ($error_condition == 1) ? undef : 1; |
207
|
|
|
|
|
|
|
} |
208
|
|
|
|
|
|
|
|
209
|
|
|
|
|
|
|
|
210
|
|
|
|
|
|
|
=item $result = $obj_instance->close(); |
211
|
|
|
|
|
|
|
|
212
|
|
|
|
|
|
|
This method closes the object file stream and resets all internal data |
213
|
|
|
|
|
|
|
structures. The result of the operation is returned. If the file stream |
214
|
|
|
|
|
|
|
is closed successfully, this object returns true (1), otherwise false |
215
|
|
|
|
|
|
|
(undefined). |
216
|
|
|
|
|
|
|
|
217
|
|
|
|
|
|
|
=cut |
218
|
|
|
|
|
|
|
|
219
|
|
|
|
|
|
|
|
220
|
|
|
|
|
|
|
sub close() { |
221
|
0
|
|
|
0
|
1
|
|
my $self = shift; |
222
|
0
|
|
|
|
|
|
my $return_val = undef; |
223
|
|
|
|
|
|
|
|
224
|
0
|
0
|
|
|
|
|
if ( defined ( $self->{db_handle} ) ) { |
225
|
0
|
|
|
|
|
|
$return_val = $self->{db_handle}->close(); |
226
|
0
|
0
|
|
|
|
|
if ( ! defined ($return_val) ) { |
227
|
0
|
|
|
|
|
|
$return_val = undef; |
228
|
0
|
|
|
|
|
|
$self->_errorHandler( |
229
|
|
|
|
|
|
|
"Error closing FASTA file: $self->{db_file_name}", |
230
|
|
|
|
|
|
|
$DEBUG_LEVEL_4, $USR_ERR); |
231
|
|
|
|
|
|
|
} |
232
|
|
|
|
|
|
|
} |
233
|
0
|
|
|
|
|
|
$self->{db_file_name} = undef; |
234
|
0
|
|
|
|
|
|
$self->{db_handle} = undef; |
235
|
0
|
|
|
|
|
|
return $return_val; |
236
|
|
|
|
|
|
|
} |
237
|
|
|
|
|
|
|
|
238
|
|
|
|
|
|
|
|
239
|
|
|
|
|
|
|
=item $result = $obj_instance->write($fasta_obj); |
240
|
|
|
|
|
|
|
|
241
|
|
|
|
|
|
|
This method takes a TIGR::FASTA::Record object, C<$fasta_obj>, and writes it |
242
|
|
|
|
|
|
|
to the file specified in C or C. On success, this method |
243
|
|
|
|
|
|
|
returns true (1). On error, this method returns false (undefined) and logs |
244
|
|
|
|
|
|
|
an error message. |
245
|
|
|
|
|
|
|
|
246
|
|
|
|
|
|
|
=cut |
247
|
|
|
|
|
|
|
|
248
|
|
|
|
|
|
|
|
249
|
|
|
|
|
|
|
sub write($) { |
250
|
0
|
|
|
0
|
1
|
|
my $self = shift; |
251
|
0
|
|
|
|
|
|
my $fasta_obj = shift; |
252
|
0
|
|
|
|
|
|
my $return_val = 1; |
253
|
|
|
|
|
|
|
|
254
|
0
|
0
|
0
|
|
|
|
if ( ( defined ($fasta_obj) ) && |
|
|
|
0
|
|
|
|
|
255
|
|
|
|
|
|
|
( ( ref($fasta_obj) ) =~ /fastarecord/i ) && |
256
|
|
|
|
|
|
|
( defined ($self->{db_handle}) ) ) { |
257
|
|
|
|
|
|
|
|
258
|
0
|
0
|
|
|
|
|
if ( $self->{db_handle}->print($fasta_obj->toString()) ) { |
259
|
0
|
|
|
|
|
|
$return_val = 1; |
260
|
|
|
|
|
|
|
} |
261
|
|
|
|
|
|
|
else { |
262
|
|
|
|
|
|
|
|
263
|
0
|
|
|
|
|
|
$return_val = undef; |
264
|
0
|
|
|
|
|
|
$self->_errorHandler( |
265
|
|
|
|
|
|
|
"Error printing to FASTA file: $self->{db_file_name}", |
266
|
|
|
|
|
|
|
$DEBUG_LEVEL_3, $USR_ERR); |
267
|
|
|
|
|
|
|
} |
268
|
|
|
|
|
|
|
} |
269
|
|
|
|
|
|
|
else { |
270
|
0
|
|
|
|
|
|
$return_val = undef; |
271
|
0
|
|
|
|
|
|
$self->_errorHandler( |
272
|
|
|
|
|
|
|
"Invalid method of initialization for " . |
273
|
|
|
|
|
|
|
"TIGR::FASTA::Writer", $DEBUG_LEVEL_3, $USR_ERR); |
274
|
|
|
|
|
|
|
} |
275
|
0
|
|
|
|
|
|
return $return_val; |
276
|
|
|
|
|
|
|
} |
277
|
|
|
|
|
|
|
|
278
|
|
|
|
|
|
|
# $message = $obj_instance->_errorHandler($message, $tf_level, |
279
|
|
|
|
|
|
|
# $internal_log_flag); |
280
|
|
|
|
|
|
|
|
281
|
|
|
|
|
|
|
#This method handles logging to the TIGR::Foundation module and |
282
|
|
|
|
|
|
|
#internal error record reference array. The C<$message> argument is logged |
283
|
|
|
|
|
|
|
#to the appropriate service. The C<$tf_level> parameter specifies the |
284
|
|
|
|
|
|
|
#logging level for TIGR::Foundation, while the C<$internal_log_flag> parameter |
285
|
|
|
|
|
|
|
#specifies if C<$message> should be written to the internal array reference |
286
|
|
|
|
|
|
|
#specified in C. If a TIGR::Foundation instance does not exist, |
287
|
|
|
|
|
|
|
#no logging to that facility occurs. This method returns C<$message>. |
288
|
|
|
|
|
|
|
|
289
|
|
|
|
|
|
|
|
290
|
|
|
|
|
|
|
|
291
|
|
|
|
|
|
|
sub _errorHandler($$$) { |
292
|
|
|
|
|
|
|
|
293
|
0
|
|
|
0
|
|
|
my $self = shift; |
294
|
|
|
|
|
|
|
|
295
|
0
|
|
|
|
|
|
my ( $message, $tf_level, $log_facility ) = @_; |
296
|
|
|
|
|
|
|
|
297
|
0
|
0
|
0
|
|
|
|
if ( defined ($message) && |
|
|
|
0
|
|
|
|
|
298
|
|
|
|
|
|
|
defined ($tf_level) && |
299
|
|
|
|
|
|
|
defined ($log_facility) ) { |
300
|
|
|
|
|
|
|
|
301
|
0
|
0
|
|
|
|
|
if ( defined ($self->{foundation}) ) { |
302
|
0
|
0
|
|
|
|
|
if ( $log_facility != $USR_ERR ) { # all user errors go to .error |
303
|
0
|
|
|
|
|
|
$self->{foundation}->logLocal($message, $tf_level); |
304
|
|
|
|
|
|
|
} |
305
|
|
|
|
|
|
|
else { |
306
|
0
|
|
|
|
|
|
$self->{foundation}->logError($message); |
307
|
|
|
|
|
|
|
} |
308
|
|
|
|
|
|
|
} |
309
|
|
|
|
|
|
|
|
310
|
0
|
0
|
0
|
|
|
|
if ( ( defined ($self->{error_ref}) ) && |
311
|
|
|
|
|
|
|
( $log_facility == $USR_ERR ) ) { |
312
|
0
|
|
|
|
|
|
push @{$self->{error_ref}}, $message; |
|
0
|
|
|
|
|
|
|
313
|
|
|
|
|
|
|
} |
314
|
|
|
|
|
|
|
} |
315
|
0
|
|
|
|
|
|
return $message; |
316
|
|
|
|
|
|
|
} |
317
|
|
|
|
|
|
|
|
318
|
|
|
|
|
|
|
|
319
|
|
|
|
|
|
|
=back |
320
|
|
|
|
|
|
|
|
321
|
|
|
|
|
|
|
=head1 USAGE |
322
|
|
|
|
|
|
|
|
323
|
|
|
|
|
|
|
To use this module, load the C and C |
324
|
|
|
|
|
|
|
modules with the C |
325
|
|
|
|
|
|
|
via the C method, as shown below. There are several invocations |
326
|
|
|
|
|
|
|
possible for this method since all parameters to C are optional. |
327
|
|
|
|
|
|
|
An example script using this module follows. The C |
328
|
|
|
|
|
|
|
module is included for completeness but does not have to be used. |
329
|
|
|
|
|
|
|
|
330
|
|
|
|
|
|
|
#!/usr/local/bin/perl -w |
331
|
|
|
|
|
|
|
|
332
|
|
|
|
|
|
|
# This example uses the TIGR::FASTA::Writer object to write |
333
|
|
|
|
|
|
|
# a simple TIGR::FASTA::Record object to a file specified with |
334
|
|
|
|
|
|
|
# the '-o' option to this script. |
335
|
|
|
|
|
|
|
# Writing errors are collected to the '@errors_list' array. |
336
|
|
|
|
|
|
|
|
337
|
|
|
|
|
|
|
use strict; |
338
|
|
|
|
|
|
|
use TIGR::Foundation; |
339
|
|
|
|
|
|
|
use TIGR::FASTA::Record; |
340
|
|
|
|
|
|
|
use TIGR::FASTA::Writer; |
341
|
|
|
|
|
|
|
|
342
|
|
|
|
|
|
|
MAIN: |
343
|
|
|
|
|
|
|
{ |
344
|
|
|
|
|
|
|
my $tf_object = new TIGR::Foundation; |
345
|
|
|
|
|
|
|
my @errors_list = (); |
346
|
|
|
|
|
|
|
my $output_file = undef; |
347
|
|
|
|
|
|
|
|
348
|
|
|
|
|
|
|
my $getopts_result = undef; |
349
|
|
|
|
|
|
|
|
350
|
|
|
|
|
|
|
$getopts_result = $tf_object->TIGR_GetOptions( "o=s" => \$output_file ); |
351
|
|
|
|
|
|
|
|
352
|
|
|
|
|
|
|
if ( $getopts_result != 1 ) { |
353
|
|
|
|
|
|
|
$tf_object->bail("Invalid command line option."); |
354
|
|
|
|
|
|
|
} |
355
|
|
|
|
|
|
|
|
356
|
|
|
|
|
|
|
if ( ! defined ( $output_file ) ) { |
357
|
|
|
|
|
|
|
$tf_object->bail("Must specify an output file with the '-o' option"); |
358
|
|
|
|
|
|
|
} |
359
|
|
|
|
|
|
|
|
360
|
|
|
|
|
|
|
my $header = "ORF00001"; |
361
|
|
|
|
|
|
|
my $data = "ATGC"; |
362
|
|
|
|
|
|
|
|
363
|
|
|
|
|
|
|
my $fasta_record = new TIGR::FASTA::Record $header, $data; |
364
|
|
|
|
|
|
|
if ( ! defined ( $fasta_record ) ) { |
365
|
|
|
|
|
|
|
$tf_object->bail("Cannot create TIGR::FASTA::Record object"); |
366
|
|
|
|
|
|
|
} |
367
|
|
|
|
|
|
|
|
368
|
|
|
|
|
|
|
# Create a TIGR::FASTA::Writer instance using TIGR::Foundation and |
369
|
|
|
|
|
|
|
# an error message list. |
370
|
|
|
|
|
|
|
|
371
|
|
|
|
|
|
|
my $fasta_writer = new TIGR::FASTA::Writer $tf_object, \@errors_list; |
372
|
|
|
|
|
|
|
|
373
|
|
|
|
|
|
|
$fasta_writer->open($output_file) or |
374
|
|
|
|
|
|
|
$tf_object->logLocal("Cannot open output file $output_file", |
375
|
|
|
|
|
|
|
$DEBUG_LEVEL_1); |
376
|
|
|
|
|
|
|
|
377
|
|
|
|
|
|
|
if ( scalar(@errors_list) > 0 ) { # are there parse errors? |
378
|
|
|
|
|
|
|
while ( @errors_list ) { # get the messages from the list |
379
|
|
|
|
|
|
|
my $message = shift @errors_list; |
380
|
|
|
|
|
|
|
print STDERR $message, "\n"; |
381
|
|
|
|
|
|
|
} |
382
|
|
|
|
|
|
|
} |
383
|
|
|
|
|
|
|
|
384
|
|
|
|
|
|
|
$fasta_writer->write($fasta_record ) or |
385
|
|
|
|
|
|
|
$tf_object->logLocal("Cannot write FASTA record to $output_file", |
386
|
|
|
|
|
|
|
$DEBUG_LEVEL_1); |
387
|
|
|
|
|
|
|
} |
388
|
|
|
|
|
|
|
|
389
|
|
|
|
|
|
|
=cut |
390
|
|
|
|
|
|
|
|
391
|
|
|
|
|
|
|
} |
392
|
|
|
|
|
|
|
|
393
|
|
|
|
|
|
|
1; |