[Bioperl-l] local blast
Bradley W. Langhorst
brad@langhorst.com
Sun, 5 Nov 2000 16:04:27 -0500
Okay
so im trying to do something that should be stupid easy.
I want to blast against a small local nucleotide database.
I've modified LocalBlast.pm with the following diffs
RCS file: /home/repository/bioperl/bioperl-
live/Bio/Tools/Blast/Run/LocalBlast.pm,v
retrieving revision 1.4
diff -r1.4 LocalBlast.pm
175c175
< @Blast_dbp_local = qw( YOUR LOCAL PEPTIDE BLAST
DATASETS HERE );
---
> @Blast_dbp_local = qw( ); #none
178c178
< @Blast_dbn_local = qw( YOUR LOCAL NUCLEOTIDE BLAST
DATASETS HERE );
---
> @Blast_dbn_local = qw( "~bwlang/tmp/unh/blast/etr2.blast.db");
181c181
< @Blast_matrix_local = qw( YOUR LOCAL SUBSTITUTION
SCORING MATRICES HERE );
---
> @Blast_matrix_local = qw( ); #na, nucleotide only
188c188
< my $program =''; # your default blast program
---
> my $program ='blastall'; # your default blast program
190c190
< my $gappedAlignmentFlag =''; # flag to indicate if gapping is
on or off
---
> my $gappedAlignmentFlag ='1'; # flag to indicate if gapping is
on or off
but when i run like this
%runParams = (-prog => 'blastn',
-method => 'local',
-database => "blast/etr2.blast.db",
-seqs => [$seqobj]);
#blast this file against the blastdb containing all the files
$blastObj = Bio::Tools::Blast->new( -run => \%runParams,
-parse => 1,
-signif => '1e-10',
-strict => 1,
);
I get
"blast_local" is not exported by the
Bio::Tools::Blast::Run::Webblast module at
./find_common_elements.pl line 49
so I go over there and make these changes
diff -r1.17 Blast.pm
1281,1282c1281,1282
< require Bio::Tools::Blast::Run::Webblast;
< Bio::Tools::Blast::Run::Webblast->import(qw(&blast_local));
---
> require Bio::Tools::Blast::Run::LocalBlast;
> Bio::Tools::Blast::Run::LocalBlast->import(qw(&blast_local));
But this fails very badly with inability to resolve _rearrange...
Is this local blasting stuff supposed to be working?
Did I make some error in calling the blast?
Or in the customizatoin of that LocalBlast.pm module?
thanks!
brad