[Bioperl-l] Bio::Tools::Genewise.pm problems

Jason Stajich jason.stajich at duke.edu
Thu Nov 18 11:01:23 EST 2004


You need -genesf.  If you want the protein and DNA sequence names 
parsed out add -para as well.

So -genesf -para

Bio::SearchIO::wise may be more forgiving, I don't remember - we 
probably shouldn't be throwing errors with the parser as you saw, but 
we basically get all the data from the -genesf option so without it I 
don't think the parser can do much.  We don't make any attempt to parse 
the wise alignment block.

You can always look at the tests in t/ to get a sense of how an object 
is used (there will be a reference to a datafile in the test script 
usually)

-jason
On Nov 18, 2004, at 10:50 AM, michael watson (IAH-C) wrote:

> Hi
>
> I'm on Bioperl 1.4 on linux.
>
> Are there any particular output options I need to specify when using
> genewise and the BioPerl parser?
>
> I have a genewise output, using default options.  When I use the sample
> code in Bio::Tools::Genewise, nothing is printed out.  When I use the
> sample code with an example output from the wise2 package
> (basic_genomic.out) I get a error trace:
>
> EXCEPTION
> MSG: Need a start
> STACK Bio::Tools::Genewise::_get_strand Genewise.pm::116
>
> Is anyone using this parser without any errors?  If so, any tips on how
> to get it to work?
>
> Thanks
>
> Mick
>
> _______________________________________________
> Bioperl-l mailing list
> Bioperl-l at portal.open-bio.org
> http://portal.open-bio.org/mailman/listinfo/bioperl-l
>
>
--
Jason Stajich
jason.stajich at duke.edu
http://www.duke.edu/~jes12/



More information about the Bioperl-l mailing list