You'll have to read-up on the SearchIO system for it to make any sense.
http://bioperl.org/HOWTOs/SearchIO/index.html
The HSPs are the "features" which are written back out with the
Bio::Tools::GFF module.
Most of the work is already done in the search2gff.PLS script for you
-- there is a lot of code in there to handle asking for query or hit
strand (you can only output one or the other) and filtering.
The code is in bioperl scripts/utilities directory or you can pull it
down here:
http://bioperl.org/SRC/bioperl-live/scripts/utilities/search2gff.PLS
So you want to run it like this (argument order doesn't matter)
perl search2gff -o myout.gff -i myinput.psl -f psl
-jason
--
Jason Stajich
jason.stajich at duke.edu
http://www.duke.edu/~jes12/
On Mar 4, 2005, at 1:07 PM, Thiago Motta Venancio wrote:
> Dear Jason.
> Thanks for replying me.
> I saw the documentation of your package before writing to the list,
> but i did not understand how to use it.
> Sorry about my low knowledge in Bioperl.
> Here is your code:
>> use Bio::SearchIO;
> my $parser = new Bio::SearchIO(-file => 'file.psl',
> -format => 'psl');
> while( my $result = $parser->next_result ) {
> }
>>> The question is where to specify the GFF format...
> Regards.
> Thiago
>> --
> Thiago Motta Venancio - PhD student in Bioinformatics
>>-------------- next part --------------
A non-text attachment was scrubbed...
Name: PGP.sig
Type: application/pgp-signature
Size: 186 bytes
Desc: This is a digitally signed message part
Url : http://portal.open-bio.org/pipermail/bioperl-l/attachments/20050304/107cffd0/PGP.bin