forked from bioperl/bioperl-live
-
Notifications
You must be signed in to change notification settings - Fork 1
/
README
258 lines (193 loc) · 10.3 KB
/
README
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
145
146
147
148
149
150
151
152
153
154
155
156
157
158
159
160
161
162
163
164
165
166
167
168
169
170
171
172
173
174
175
176
177
178
179
180
181
182
183
184
185
186
187
188
189
190
191
192
193
194
195
196
197
198
199
200
201
202
203
204
205
206
207
208
209
210
211
212
213
214
215
216
217
218
219
220
221
222
223
224
225
226
227
228
229
230
231
232
233
234
235
236
237
238
239
240
241
242
243
244
245
246
247
248
249
250
251
252
253
254
255
256
257
258
This is the README file for the BioPerl central distribution.
o Getting Started
Please see the the INSTALL or INSTALL.WIN documents for installation
instructions.
o About BioPerl
BioPerl is a package of public domain Perl tools for computational
molecular biology.
Our website, http://bioperl.org/, provides an online resource of
modules, scripts, and web links for developers of Perl-based software
for life science research.
o Contact info
BioPerl mailing list: [email protected]
There's quite a variety of tools available in BioPerl, and more are
added all the time. If the tool you're looking for isn't described in
the documentation please write us, it could be undocumented or in
process.
Project website : http://bioperl.org/
Bug reports : https://redmine.open-bio.org/projects/bioperl/
Please send us bugs, in particular about documentation which you
think is unclear or problems in installation. We are also very
interested in functions which don't work the way you think they
do!
o The directory structure
The BioPerl directory structure is organized as follows:
- Bio/ - BioPerl modules
- doc/ - Documentation utilities
- examples/ - Scripts demonstrating the many uses of BioPerl
- ide/ - files for developing BioPerl using an IDE
- maintenance/ - BioPerl housekeeping scripts
- models/ - DIA drawing program generated OO UML for BioPerl classes
(these are quite out-of-date)
- scripts/ - Useful production-quality scripts with POD documentation
- t/ - Perl built-in tests, tests are divided into subdirectories
based on the specific classes being tested
- t/data/ - Data files used for the tests, provides good example data
o Documentation
For documentation on BioPerl see the HOWTO documents and tutorials
online at http://bioperl.org.
Useful documentation in the form of example code can also be found in
the examples/ and scripts/ directories. The current collection
includes scripts that run BLAST, index flat files, parse PDB
structure files, make primers, retrieve ESTs based on tissue, align
protein to nucleotide sequence, run GENSCAN on multiple sequences,
and much more! See bioscripts.pod for a complete listing.
Individual *.pm modules have their own embedded POD documentation as
well. A complete set of hyperlinked POD, or module, documentation is
available at http://www.bioperl.org/.
Remember that 'perldoc' is your friend. You can use it to read any
file containing POD formatted documentation without needing any type
of translator (e.g. 'perldoc Bio::SeqIO').
If you used the Build.PL installation, and depending on your
platform, you may have documentation installed as man pages, which
can be accessed in the usual way.
o Releases
BioPerl releases are always available from the website at
http://www.bioperl.org/DIST or in CPAN. The latest code can be found
at https://github.com/bioperl.
BioPerl formerly used a numbering scheme to indicate stable release
series vs. development release series. A release number is a three
digit number like 1.2.0. The first digit indicates the major release
- the idea being that all the API calls in a major release are
reasonably consistent. The second number is the release series. This
is probably the most important number.
From the 1.0 release until the 1.6 release, even numbers (1.0, 1.2
etc) indicated stable releases. Stable releases were well tested and
recommended for most uses. Odd numbers (1.1, 1.3 etc) were development
releases which one would only use if one were interested in the
latest and greatest features. The final number (e.g. 1.2.0, 1.2.1) is
the bug fix release. The higher the number the more bug fixes has
been incorporated. In theory you can upgrade from one bug fix release
to the next with no changes to your own code (for production cases,
obviously check things out carefully before you switch over).
The 1.6 release will be the last release series to utilize the
alternating 'stable'/'developer' convention. Starting immediately
after the 1.6 branch, we will start splitting BioPerl into several
smaller easier-to-manage distributions, including a developer
distribution for cutting-edge (in development) code, untested
modules, and alternative implementations.
o Caveats and warnings
When you run the tests ("./Build test") some tests may issue warnings
messages or even fail. Sometimes this is because we didn't have
anyone to test the test system on the combination of your operating
system, version of perl, and associated libraries and other modules.
Because BioPerl depends on several outside libraries we may not be
able to test every single combination so if there are warnings you
may find that the package is still perfectly useful.
If you install the bioperl-run system and run tests when you don't
have the program installed you'll get messages like 'program XXX not
found, skipping tests'. That's okay, BioPerl is doing what it is
supposed to do. If you wanted to run the program you'd need to
install it first.
Not all scripts in the examples/ directory are correct and up-to-date.
We need volunteers to help maintain these so if you find they do not
submit a bug report to https://redmine.open-bio.org/projects/bioperl/
and consider helping out in their maintenance.
If you are confused about what modules are appropriate when you try
and solve a particular issue in bioinformatics we urge you to look at
HOWTO documents first.
o A simple module summary
Here is a quick summary of many of the useful modules and how the
toolkit is laid out:
All modules are in the Bio/ namespace,
- Perl is for newbies and gives a functional interface to the main
parts of the package
- Seq is for Sequences (protein and DNA).
o Bio::PrimarySeq is a plain sequence (sequence data + identifiers)
o Bio::Seq is a PrimarySeq plus it has a Bio::Annotation::Collection
and Bio::SeqFeatureI objects attached
(via Bio::FeatureHolderI).
o Bio::Seq::RichSeq is all of the above plus it has slots for
extra information specific to GenBank/EMBL/SwissProt files.
o Bio::Seq::LargeSeq is for sequences which are too big for
fitting into memory.
- SeqIO is for reading and writing Sequences, it is a front end
module for separate driver modules supporting the different
sequence formats
- SeqFeature - start/stop/strand annotations of sequences
o Bio::SeqFeature::Generic is basic catchall
o Bio::SeqFeature::Similarity a similarity sequence feature
o Bio::SeqFeature::FeaturePair a sequence feature which is pairwise
such as query/hit pairs
- SearchIO is for reading and writing pairwise alignment reports like
BLAST or FASTA
- Search is where the alignment objects are defined
o Bio::Search::Result::GenericResult is the result object (a blast
query is a Result object)
o Bio::Search::Hit::GenericHit is the Hit object (a query will have
0 -> many hits in a database)
o Bio::Search::HSP::GenericHSP is the High-scoring Segment Pair
object defining the alignment(s) of the query and hit.
- SimpleAlign is for multiple sequence alignments
- AlignIO is for reading and writing multiple sequence alignment
formats
- Assembly provides the start of an infrastructure for assemblies and
Assembly::IO IO converters for them
- DB is the namespace for all the database query objects
o Bio::DB::GenBank/GenPept are two modules which query NCBI entrez
for sequences
o Bio::DB::SwissProt/EMBL query various EMBL and SwissProt
repositories for a sequences
o Bio::DB::GFF is Lincoln Stein's fast, lightweight feature and
sequence database which is the backend to his GBrowse system (see
www.gmod.org)
o Bio::DB::Flat is a fast implementation of the OBDA flat-file
indexing system (cross-language and cross-platform supported by
O|B|F projects see http://obda.open-bio.org).
o Bio::DB::BioFetch/DBFetch for OBDA, Web (HTTP) access to remote
databases.
o Bio::DB::InMemoryCache/FileCache (fast local caching of sequences
from remote dbs to speed up your access).
o Bio::DB::Registry interface to the OBDA specification for remote
data sources
o Bio::DB::Biblio for access to remote bibliographic databases.
o Bio::DB::EUtilities is the initial set of modules used for
generic queried using NCBI's eUtils.
- Annotation collection of annotation objects (comments, DBlinks,
References, and misc key/value pairs)
- Coordinate is a system for mapping between different coordinate
systems such as DNA to protein or between assemblies
- Index is for locally indexed flatfiles with BerkeleyDB
- Tools contains many miscellaneous parsers and function for
different bioinformatics needs
o Gene prediction parser (Genscan, MZEF, Grail, Genemark)
o Annotation format (GFF)
o Enumerate codon tables and valid sequences symbols (CodonTable,
IUPAC)
o Phylogenetic program parsing (PAML, Molphy, Phylip)
- Map genetic and physical map representations
- Structure - parse and represent protein structure data
- TreeIO is for reading and writing Tree formats
- Tree is the namespace for all the associated Tree objects
o Bio::Tree::Tree is the basic tree object
o Bio::Tree::Node are the nodes which make up the tree
o Bio::Tree::Statistics is for computing statistics for a tree
o Bio::Tree::TreeFunctionsI is where specific tree functions are
implemented (like is_monophyletic and lca)
- Bio::Biblio is where bibliographic data and database access objects
are kept
- Variation represent sequences with mutations and variations applied
so one can compare and represent wild-type and mutation versions of
a sequence.
- Root, basic objects for the internals of BioPerl
o Upgrading from an older version
If you have a previously installed version of BioPerl on your system
some of these notes may help you.
Some modules have been removed because they have been superceded by
new development efforts. They are documented in the DEPRECATED file
that is included in the release. In addition some methods, or the
Application Programming Interface (API), have changed or been
removed. You may find that scripts which worked with BioPerl 1.4 may
give you warnings or may not work at all (although we have tried very
hard to minimize this!). Send an email to the list and we'll be happy
to give you pointers.