[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