0
|
1 # $Id: RichSeqI.pm,v 1.8 2002/10/22 07:38:40 lapp Exp $
|
|
2 #
|
|
3 # BioPerl module for Bio::Seq::RichSeqI
|
|
4 #
|
|
5 # Cared for by Ewan Birney <birney@ebi.ac.uk>
|
|
6 #
|
|
7 # Copyright Ewan Birney
|
|
8 #
|
|
9 # You may distribute this module under the same terms as perl itself
|
|
10
|
|
11 # POD documentation - main docs before the code
|
|
12
|
|
13 =head1 NAME
|
|
14
|
|
15 Bio::Seq::RichSeqI - RichSeq interface, mainly for database orientated sequences
|
|
16
|
|
17 =head1 SYNOPSIS
|
|
18
|
|
19 @secondary = $richseq->get_secondary_accessions;
|
|
20 $division = $richseq->division;
|
|
21 $mol = $richseq->molecule;
|
|
22 @dates = $richseq->get_dates;
|
|
23 $seq_version = $richseq->seq_version;
|
|
24 $pid = $richseq->pid;
|
|
25 $keywords = $richseq->keywords;
|
|
26
|
|
27 =head1 DESCRIPTION
|
|
28
|
|
29 This interface extends the Bio::SeqI interface to give additional functionality
|
|
30 to sequences with richer data sources, in particular from database sequences
|
|
31 (EMBL, GenBank and Swissprot).
|
|
32
|
|
33 =head1 FEEDBACK
|
|
34
|
|
35 =head2 Mailing Lists
|
|
36
|
|
37 User feedback is an integral part of the evolution of this
|
|
38 and other Bioperl modules. Send your comments and suggestions preferably
|
|
39 to one of the Bioperl mailing lists.
|
|
40 Your participation is much appreciated.
|
|
41
|
|
42 bioperl-l@bioperl.org - General discussion
|
|
43 http://bio.perl.org/MailList.html - About the mailing lists
|
|
44
|
|
45 =head2 Reporting Bugs
|
|
46
|
|
47 Report bugs to the Bioperl bug tracking system to help us keep track
|
|
48 the bugs and their resolution.
|
|
49 Bug reports can be submitted via email or the web:
|
|
50
|
|
51 bioperl-bugs@bio.perl.org
|
|
52 http://bugzilla.bioperl.org/
|
|
53
|
|
54 =head1 AUTHOR - Ewan Birney
|
|
55
|
|
56 Email birney@ebi.ac.uk
|
|
57
|
|
58 Describe contact details here
|
|
59
|
|
60 =head1 APPENDIX
|
|
61
|
|
62 The rest of the documentation details each of the object methods. Internal methods are usually preceded with a _
|
|
63
|
|
64 =cut
|
|
65
|
|
66
|
|
67 # Let the code begin...
|
|
68
|
|
69
|
|
70 package Bio::Seq::RichSeqI;
|
|
71 use vars qw(@ISA);
|
|
72 use strict;
|
|
73 use Bio::SeqI;
|
|
74
|
|
75 @ISA = ('Bio::SeqI');
|
|
76
|
|
77
|
|
78 =head2 get_secondary_accessions
|
|
79
|
|
80 Title : get_secondary_accessions
|
|
81 Usage :
|
|
82 Function: Get the secondary accessions for a sequence.
|
|
83 Example :
|
|
84 Returns : an array of strings
|
|
85 Args : none
|
|
86
|
|
87
|
|
88 =cut
|
|
89
|
|
90 sub get_secondary_accessions{
|
|
91 my ($self,@args) = @_;
|
|
92
|
|
93 $self->throw("hit get_secondary_accessions in interface definition - error");
|
|
94
|
|
95 }
|
|
96
|
|
97
|
|
98 =head2 division
|
|
99
|
|
100 Title : division
|
|
101 Usage :
|
|
102 Function: Get (and set, depending on the implementation) the divison for
|
|
103 a sequence.
|
|
104
|
|
105 Examples from GenBank are PLN (plants), PRI (primates), etc.
|
|
106 Example :
|
|
107 Returns : a string
|
|
108 Args :
|
|
109
|
|
110
|
|
111 =cut
|
|
112
|
|
113 sub division{
|
|
114 my ($self,@args) = @_;
|
|
115
|
|
116 $self->throw("hit division in interface definition - error");
|
|
117
|
|
118 }
|
|
119
|
|
120
|
|
121 =head2 molecule
|
|
122
|
|
123 Title : molecule
|
|
124 Usage :
|
|
125 Function: Get (and set, depending on the implementation) the molecule
|
|
126 type for the sequence.
|
|
127
|
|
128 This is not necessarily the same as Bio::PrimarySeqI::alphabet(),
|
|
129 because it is databank-specific.
|
|
130 Example :
|
|
131 Returns : a string
|
|
132 Args :
|
|
133
|
|
134
|
|
135 =cut
|
|
136
|
|
137 sub molecule{
|
|
138 my ($self,@args) = @_;
|
|
139
|
|
140 $self->throw("hit molecule in interface definition - error");
|
|
141 }
|
|
142
|
|
143 =head2 pid
|
|
144
|
|
145 Title : pid
|
|
146 Usage :
|
|
147 Function: Get (and set, depending on the implementation) the PID property
|
|
148 for the sequence.
|
|
149 Example :
|
|
150 Returns : a string
|
|
151 Args :
|
|
152
|
|
153
|
|
154 =cut
|
|
155
|
|
156 sub pid {
|
|
157 my ($self,@args) = @_;
|
|
158
|
|
159 $self->throw("hit pid in interface definition - error");
|
|
160 }
|
|
161
|
|
162 =head2 get_dates
|
|
163
|
|
164 Title : get_dates
|
|
165 Usage :
|
|
166 Function: Get (and set, depending on the implementation) the dates the
|
|
167 databank entry specified for the sequence
|
|
168 Example :
|
|
169 Returns : an array of strings
|
|
170 Args :
|
|
171
|
|
172
|
|
173 =cut
|
|
174
|
|
175 sub get_dates{
|
|
176 my ($self,@args) = @_;
|
|
177
|
|
178 $self->throw("hit get_dates in interface definition - error");
|
|
179
|
|
180 }
|
|
181
|
|
182
|
|
183 =head2 seq_version
|
|
184
|
|
185 Title : seq_version
|
|
186 Usage :
|
|
187 Function: Get (and set, depending on the implementation) the version string
|
|
188 of the sequence.
|
|
189 Example :
|
|
190 Returns : a string
|
|
191 Args :
|
|
192
|
|
193
|
|
194 =cut
|
|
195
|
|
196 sub seq_version{
|
|
197 my ($self,@args) = @_;
|
|
198
|
|
199 $self->throw("hit seq_version in interface definition - error");
|
|
200
|
|
201 }
|
|
202
|
|
203 =head2 keywords
|
|
204
|
|
205 Title : keywords
|
|
206 Usage : $obj->keywords($newval)
|
|
207 Function:
|
|
208 Returns : value of keywords (a string)
|
|
209 Args : newvalue (optional) (a string)
|
|
210
|
|
211
|
|
212 =cut
|
|
213
|
|
214 sub keywords {
|
|
215 my ($self) = @_;
|
|
216 $self->throw("hit keywords in interface definition - error");
|
|
217 }
|
|
218
|
|
219 1;
|