[Bioperl-guts-l] [16535] bioperl-dev/branches/blastplus-dev-branch/lib/Bio/Tools/Run/ StandAloneBlastPlus/BlastMethods.pm: POD for BlastMethods

Mark Allen Jensen maj at dev.open-bio.org
Fri Dec 18 22:23:39 EST 2009


Revision: 16535
Author:   maj
Date:     2009-12-18 22:23:39 -0500 (Fri, 18 Dec 2009)
Log Message:
-----------
POD for BlastMethods

Modified Paths:
--------------
    bioperl-dev/branches/blastplus-dev-branch/lib/Bio/Tools/Run/StandAloneBlastPlus/BlastMethods.pm

Modified: bioperl-dev/branches/blastplus-dev-branch/lib/Bio/Tools/Run/StandAloneBlastPlus/BlastMethods.pm
===================================================================
--- bioperl-dev/branches/blastplus-dev-branch/lib/Bio/Tools/Run/StandAloneBlastPlus/BlastMethods.pm	2009-12-19 03:06:27 UTC (rev 16534)
+++ bioperl-dev/branches/blastplus-dev-branch/lib/Bio/Tools/Run/StandAloneBlastPlus/BlastMethods.pm	2009-12-19 03:23:39 UTC (rev 16535)
@@ -18,8 +18,23 @@
 
 =head1 SYNOPSIS
 
-Give standard usage here
+ # create a factory:
+ $fac = Bio::Tools::Run::StandAloneBlastPlus->new(
+    -db_name => 'testdb'
+ );
+ # get your results
+ $result = $fac->blastn( -query => 'query_seqs.fas',
+                         -outfile => 'query.bls',
+                         -method_args => [ '-num_alignments' => 10 ,
 
+ $result = $fac->tblastx( -query => $an_alignment_object,
+                          -outfile => 'query.bls',
+                          -outformat => 7 );
+ # do a bl2seq
+ $fac->bl2seq( -method => 'blastp',
+               -query => $seq_object_1,
+               -subject => $seq_object_2 );
+
 =head1 DESCRIPTION
 
 This module provides the BLAST methods (blastn, blastp, psiblast,
@@ -31,7 +46,7 @@
 L<Bio::Tools::Run::StandAloneBlastPlus> factory object. The object
 itself has extensive facilities for creating, formatting, and masking
 BLAST databases; please refer to
-L<Bio::Tools::Run::StandAloneBlastPlus> POD for these details
+L<Bio::Tools::Run::StandAloneBlastPlus> POD for these details.
 
 Given a C<StandAloneBlastPlus> factory, such as
 
@@ -47,18 +62,55 @@
 
 Here, C<$result> is a L<Bio::Search::Result::BlastResult> object.
 
-The blast output file can be named explicitly:
+Other details:
 
+=over
+
+=item * The blast output file can be named explicitly:
+
  $result = $fac->blastn( -query => 'query_seqs.fas',
                          -outfile => 'query.bls' );
 
+=item * The output format can be specified:
 
+ $result = $fac->blastn( -query => 'query_seqs.fas',
+                         -outfile => 'query.bls',
+                         -outformat => 7 ); #tabular
 
- 
+=item * Additional arguments to the method can be specified:
+
+ $result = $fac->blastn( -query => 'query_seqs.fas',
+                         -outfile => 'query.bls',
+                         -method_args => [ '-num_alignments' => 10 ,
+                                           '-evalue' => 100 ]);
+
+=item * To get the name of the blast output file, do 
+
+ $file = $fac->blast_out;
+
+=item * To clean up the temp files (you must do this explicitly):
+
+ $fac->cleanup;
+
+=back
+
+=head2 bl2seq()
+
+Running C<bl2seq> is similar, but both C<-query> and C<-subject> are
+required, and the attached database is ignored. The blast method must
+be specified explicitly with the C<-method> parameter:
+
+ $fac->bl2seq( -method => 'blastp',
+               -query => $seq_object_1,
+               -subject => $seq_object_2 );
+
+Other parameters ( C<-method_args>, C<-outfile>, and C<-outformat> ) are valid. 
+
 =head2 Return values
 
+The return value is always a L<Bio::Search::Result::BlastResult>
+object on success, undef on failure.
 
-
 =head1 SEE ALSO
 
 L<Bio::Tools::Run::StandAloneBlastPlus>, L<Bio::Tools::Run::BlastPlus>
@@ -293,6 +345,28 @@
 
 }
 
+=head2 blast_out()
+
+ Title   : blast_out
+ Usage   : $file = $fac->blast_out
+ Function: get the filename of the blast report file
+ Returns : scalar string
+ Args    : none
+
+=cut
+
+sub blast_output { shift->{_blastout} }
+
+# =head2 _demodernize()
+
+#  Title   : _demodernize
+#  Usage   : 
+#  Function: Ha! Wouldn't you like to know!
+#  Returns : 
+#  Args    : 
+
+# =cut
+
 sub _demodernize {
     my $file = shift;
     my $tf = File::Temp->new();



More information about the Bioperl-guts-l mailing list