[Bioperl-l] Issues on Bio::Tools::Primer3->next_primer

John M.C. Ma manchunjohn-ma at uiowa.edu
Mon Jul 28 20:45:57 UTC 2008


Hi,

I am trying a generic script to parse outputs for Primer3 results from
Bio::Tools::Run::Primer3:

use Bio::Tools::Run::Primer3;
use Bio::Tools::Primer3;
use Bio::Seq::PrimedSeq;
[...]
my $primer3_handle=Bio::Tools::Run::Primer3->new(-path=>'foobar');
[...]
my $primer3_results=$primer3_handle->run();
my $primer_test=$primer3_results->next_primer;
[...]

Which is similar to what was on the docs, but when Perl gives the following
error at running next_primer:

------------- EXCEPTION: Bio::Root::Exception -------------
MSG: No target sequence
STACK: Error::throw
STACK: Bio::Root::Root::throw
/usr/lib/perl5/site_perl/5.10.0/Bio/Root/Root.pm:359
STACK: Bio::Tools::Primer3::next_primer
/usr/lib/perl5/site_perl/5.10.0/Bio/Tools/Primer3.pm:306
STACK: /home/johnma/workspace/NCBI-test/NCBI-test.pl:208
-----------------------------------------------------------

Whereas line 208 is the line that contained the next_primer.

I tried to access the results manually by primer_results() and hard-calling
fields, but the output of primer_results() is a hash and not very convenient
for further processing. Is there a way that I can go back to next_primer...?

Cheers,

John



More information about the Bioperl-l mailing list