[Bioperl-l] Bio::Tools::Run::Cap3 - Parameters
Kevin Brown
Kevin.M.Brown at asu.edu
Wed May 6 15:23:30 UTC 2009
BEGIN {
@PARAMS = qw(a b c d e f g m n o p s u v x);
$PROGRAMDIR = '/usr/local/bin';
# Authorize attribute fields
foreach my $attr (@PARAMS) { $OK_FIELD{$attr}++;
}
That is the list of params that Cap3 will accept in the BioPerl module.
I'm guessing if you add the y to that list that it might work.
> -----Original Message-----
> From: bioperl-l-bounces at lists.open-bio.org
> [mailto:bioperl-l-bounces at lists.open-bio.org] On Behalf Of
> Michael Stubbington
> Sent: Wednesday, May 06, 2009 7:39 AM
> To: bioperl-l at lists.open-bio.org
> Subject: [Bioperl-l] Bio::Tools::Run::Cap3 - Parameters
>
> Dear all,
>
>
>
> I am using the Bio::Tools::Run::Cap3 wrapper to the Cap3 assembly
> program. I have some reads that will only assemble if cap3 is
> used with
> the '-y 150' option. This is fine from the command line but I
> can't work
> out how to pass this option to the Cap3 factory object in my script.
>
>
>
> If I do the following
>
>
>
> my $params = "y 150" ;
>
> my $cap3Factory = Bio::Tools::Run::Cap3->new($params);
>
> my $assembly = $cap3Factory->run($file);
>
>
>
> Then I get an exception as follows:
>
>
>
> ------------- EXCEPTION: Bio::Root::Exception -------------
>
> MSG: Unallowed parameter: y !
>
> STACK: Error::throw
>
> STACK: Bio::Root::Root::throw
> /Users/mike/lib/perl5/site_perl/5.8.9/Bio/Root/Root.pm:357
>
> STACK: Bio::Tools::Run::Cap3::AUTOLOAD
> /Users/mike/lib/perl5/site_perl/5.8.9/Bio/Tools/Run/Cap3.pm:116
>
> STACK: Bio::Tools::Run::Cap3::new
> /Users/mike/lib/perl5/site_perl/5.8.9/Bio/Tools/Run/Cap3.pm:101
>
> STACK: /Users/mike/perlScripts/QGenotype.pl:150
>
>
>
> If I don't try to pass any parameters to Cap3 it runs fine but just
> fails to assemble the reads that need the -y 150 flag.
>
>
>
> I'd very much appreciate any help with this. I'm pretty new
> to bioperl,
> hope I haven't missed anything obvious!
>
>
>
> Thanks in advance,
>
>
>
> Mike
>
>
>
> --------------------------------------------------------------
> ----------
> ----
>
> Mike Stubbington
>
> Novel and Dangerous Pathogens
>
> Health Protection Agency
>
> Centre for Emergency Preparedness and Response
>
> Porton Down
>
> Salisbury
>
> SP4 0JG
>
>
>
> Tel: +44 1980 619812
>
>
>
>
>
> -----------------------------------------
> **************************************************************
> ************
> The information contained in the EMail and any attachments is
> confidential and intended solely and for the attention and use of
> the named addressee(s). It may not be disclosed to any other person
> without the express authority of the HPA, or the intended
> recipient, or both. If you are not the intended recipient, you must
> not disclose, copy, distribute or retain this message or any part
> of it. This footnote also confirms that this EMail has been swept
> for computer viruses, but please re-sweep any attachments before
> opening or saving. HTTP://www.HPA.org.uk
> **************************************************************
> ************
> _______________________________________________
> Bioperl-l mailing list
> Bioperl-l at lists.open-bio.org
> http://lists.open-bio.org/mailman/listinfo/bioperl-l
>
More information about the Bioperl-l
mailing list