annotate variant_effect_predictor/Bio/SeqIO/fasta.pm @ 0:1f6dce3d34e0

Uploaded
author mahtabm
date Thu, 11 Apr 2013 02:01:53 -0400
parents
children
Ignore whitespace changes - Everywhere: Within whitespace: At end of lines:
rev   line source
0
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
1 # $Id: fasta.pm,v 1.41.2.4 2003/09/18 02:43:16 jason Exp $
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
2 # BioPerl module for Bio::SeqIO::fasta
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
3 #
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
4 # Cared for by Ewan Birney <birney@ebi.ac.uk>
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
5 # and Lincoln Stein <lstein@cshl.org>
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
6 #
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
7 # Copyright Ewan Birney & Lincoln Stein
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
8 #
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
9 # You may distribute this module under the same terms as perl itself
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
10 # _history
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
11 # October 18, 1999 Largely rewritten by Lincoln Stein
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
12
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
13 # POD documentation - main docs before the code
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
14
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
15 =head1 NAME
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
16
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
17 Bio::SeqIO::fasta - fasta sequence input/output stream
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
18
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
19 =head1 SYNOPSIS
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
20
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
21 Do not use this module directly. Use it via the Bio::SeqIO class.
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
22
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
23 =head1 DESCRIPTION
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
24
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
25 This object can transform Bio::Seq objects to and from fasta flat
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
26 file databases.
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
27
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
28 A method L<preferred_id_type()> can be used to specify the type of ID
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
29 we would like to parse from the fasta line. By default 'display' is
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
30 used, which means it parses everything from the '>' to the first space
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
31 and makes that the 'display_id' for the sequence.
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
32
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
33 Can be one of:
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
34 - accession
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
35 - accession.version
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
36 - display
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
37 - primary
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
38
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
39 =head1 FEEDBACK
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
40
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
41 =head2 Mailing Lists
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
42
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
43 User feedback is an integral part of the evolution of this and other
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
44 Bioperl modules. Send your comments and suggestions preferably to one
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
45 of the Bioperl mailing lists. Your participation is much appreciated.
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
46
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
47 bioperl-l@bioperl.org - General discussion
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
48 http://bioperl.org/MailList.shtml - About the mailing lists
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
49
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
50 =head2 Reporting Bugs
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
51
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
52 Report bugs to the Bioperl bug tracking system to help us keep track
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
53 the bugs and their resolution. Bug reports can be submitted via the
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
54 web:
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
55
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
56 http://bugzilla.bioperl.org/
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
57
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
58 =head1 AUTHORS - Ewan Birney & Lincoln Stein
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
59
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
60 Email: birney@ebi.ac.uk
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
61 lstein@cshl.org
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
62
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
63 =head1 CONTRIBUTORS
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
64
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
65 Jason Stajich, jason-at-bioperl.org
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
66
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
67 =head1 APPENDIX
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
68
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
69 The rest of the documentation details each of the object
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
70 methods. Internal methods are usually preceded with a _
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
71
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
72 =cut
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
73
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
74 # Let the code begin...
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
75
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
76 package Bio::SeqIO::fasta;
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
77 use vars qw(@ISA $WIDTH @SEQ_ID_TYPES $DEFAULT_SEQ_ID_TYPE);
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
78 use strict;
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
79 # Object preamble - inherits from Bio::Root::Object
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
80
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
81 use Bio::SeqIO;
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
82 use Bio::Seq::SeqFactory;
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
83 use Bio::Seq::SeqFastaSpeedFactory;
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
84
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
85 @ISA = qw(Bio::SeqIO);
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
86
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
87 @SEQ_ID_TYPES = qw(accession accession.version display primary);
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
88 $DEFAULT_SEQ_ID_TYPE = 'display';
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
89
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
90 BEGIN { $WIDTH = 60}
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
91
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
92 sub _initialize {
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
93 my($self,@args) = @_;
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
94 $self->SUPER::_initialize(@args);
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
95 my ($width) = $self->_rearrange([qw(WIDTH)], @args);
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
96 $width && $self->width($width);
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
97 unless ( defined $self->sequence_factory ) {
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
98 $self->sequence_factory(Bio::Seq::SeqFastaSpeedFactory->new());
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
99 }
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
100 }
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
101
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
102 =head2 next_seq
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
103
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
104 Title : next_seq
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
105 Usage : $seq = $stream->next_seq()
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
106 Function: returns the next sequence in the stream
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
107 Returns : Bio::Seq object
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
108 Args : NONE
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
109
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
110 =cut
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
111
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
112 sub next_seq {
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
113 my( $self ) = @_;
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
114 my $seq;
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
115 my $alphabet;
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
116 local $/ = "\n>";
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
117 return unless my $entry = $self->_readline;
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
118
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
119 chomp($entry);
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
120 if ($entry =~ m/\A\s*\Z/s) { # very first one
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
121 return unless $entry = $self->_readline;
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
122 chomp($entry);
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
123 }
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
124 $entry =~ s/^>//;
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
125
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
126 my ($top,$sequence) = split(/\n/,$entry,2);
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
127 defined $sequence && $sequence =~ s/>//g;
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
128 # my ($top,$sequence) = $entry =~ /^>?(.+?)\n+([^>]*)/s
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
129 # or $self->throw("Can't parse fasta entry");
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
130
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
131 my ($id,$fulldesc);
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
132 if( $top =~ /^\s*(\S+)\s*(.*)/ ) {
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
133 ($id,$fulldesc) = ($1,$2);
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
134 }
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
135
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
136 if (defined $id && $id eq '') {$id=$fulldesc;} # FIX incase no space
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
137 # between > and name \AE
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
138 defined $sequence && $sequence =~ s/\s//g; # Remove whitespace
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
139
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
140 # for empty sequences we need to know the mol.type
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
141 $alphabet = $self->alphabet();
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
142 if(defined $sequence && length($sequence) == 0) {
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
143 if(! defined($alphabet)) {
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
144 # let's default to dna
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
145 $alphabet = "dna";
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
146 }
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
147 } else {
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
148 # we don't need it really, so disable
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
149 $alphabet = undef;
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
150 }
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
151
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
152 $seq = $self->sequence_factory->create(
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
153 -seq => $sequence,
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
154 -id => $id,
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
155 # Ewan's note - I don't think this healthy
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
156 # but obviously to taste.
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
157 #-primary_id => $id,
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
158 -desc => $fulldesc,
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
159 -alphabet => $alphabet,
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
160 -direct => 1,
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
161 );
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
162
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
163
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
164
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
165
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
166 # if there wasn't one before, set the guessed type
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
167 unless ( defined $alphabet ) {
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
168 $self->alphabet($seq->alphabet());
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
169 }
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
170 return $seq;
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
171
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
172 }
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
173
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
174 =head2 write_seq
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
175
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
176 Title : write_seq
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
177 Usage : $stream->write_seq(@seq)
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
178 Function: writes the $seq object into the stream
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
179 Returns : 1 for success and 0 for error
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
180 Args : array of 1 to n Bio::PrimarySeqI objects
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
181
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
182
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
183 =cut
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
184
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
185 sub write_seq {
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
186 my ($self,@seq) = @_;
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
187 my $width = $self->width;
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
188 foreach my $seq (@seq) {
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
189 $self->throw("Did not provide a valid Bio::PrimarySeqI object")
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
190 unless defined $seq && ref($seq) && $seq->isa('Bio::PrimarySeqI');
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
191
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
192 my $str = $seq->seq;
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
193 my $top;
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
194
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
195 # Allow for different ids
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
196 my $id_type = $self->preferred_id_type;
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
197 if( $id_type =~ /^acc/i ) {
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
198 $top = $seq->accession_number();
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
199 if( $id_type =~ /vers/i ) {
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
200 $top .= "." . $seq->version();
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
201 }
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
202 } elsif($id_type =~ /^displ/i ) {
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
203 $top = $seq->display_id();
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
204 } elsif($id_type =~ /^pri/i ) {
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
205 $top = $seq->primary_id();
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
206 }
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
207
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
208 if ($seq->can('desc') and my $desc = $seq->desc()) {
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
209 $desc =~ s/\n//g;
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
210 $top .= " $desc";
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
211 }
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
212 if(length($str) > 0) {
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
213 $str =~ s/(.{1,$width})/$1\n/g;
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
214 } else {
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
215 $str = "\n";
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
216 }
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
217 $self->_print (">",$top,"\n",$str) or return;
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
218 }
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
219
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
220 $self->flush if $self->_flush_on_write && defined $self->_fh;
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
221 return 1;
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
222 }
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
223
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
224 =head2 width
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
225
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
226 Title : width
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
227 Usage : $obj->width($newval)
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
228 Function: Get/Set the line width for FASTA output
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
229 Returns : value of width
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
230 Args : newvalue (optional)
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
231
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
232
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
233 =cut
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
234
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
235 sub width{
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
236 my ($self,$value) = @_;
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
237 if( defined $value) {
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
238 $self->{'width'} = $value;
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
239 }
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
240 return $self->{'width'} || $WIDTH;
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
241 }
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
242
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
243 =head2 preferred_id_type
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
244
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
245 Title : preferred_id_type
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
246 Usage : $obj->preferred_id_type('accession')
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
247 Function: Get/Set the preferred type of identifier to use in the ">ID" position
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
248 for FASTA output.
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
249 Returns : string, one of values defined in @Bio::SeqIO::fasta::SEQ_ID_TYPES.
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
250 Default = $Bio::SeqIO::fasta::DEFAULT_SEQ_ID_TYPE ('display').
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
251 Args : string when setting. This must be one of values defined in
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
252 @Bio::SeqIO::fasta::SEQ_ID_TYPES. Allowable values:
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
253 accession, accession.version, display, primary
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
254 Throws : fatal exception if the supplied id type is not in @SEQ_ID_TYPES.
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
255
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
256 =cut
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
257
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
258 sub preferred_id_type {
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
259 my ($self,$type) = @_;
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
260 if( defined $type ) {
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
261 if( ! grep lc($type) eq $_, @SEQ_ID_TYPES) {
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
262 $self->throw(-class=>'Bio::Root::BadParameter',
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
263 -text=>"Invalid ID type \"$type\". Must be one of: @SEQ_ID_TYPES");
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
264 }
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
265 $self->{'_seq_id_type'} = lc($type);
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
266 # print STDERR "Setting preferred_id_type=$type\n";
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
267 }
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
268 $self->{'_seq_id_type'} || $DEFAULT_SEQ_ID_TYPE;
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
269 }
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
270
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
271 1;