The London Perl and Raku Workshop takes place on 26th Oct 2024. If your company depends on Perl, please consider sponsoring and/or attending.

NAME

Bio::Tools::PrositeScan - Parser for ps_scan result

SYNOPSIS

use Bio::Tools::PrositeScan;

my $factory = Bio::Tools::PrositeScan->new(
    -file   => 'out.PrositeScan',
    -format => 'fasta'
);

while(my $match = $factory->next_prediction){
    #  $match is a Bio::SeqFeature::FeaturePair

    # Sequence ID
    my $seq_id = $match->seq_id;

    # PROSITE accession number
    my $psac = $match->hseq_id;

    # Coordinates
    my @coords = ( $match->start, $match->end );

    # Subsequence
    my $seq = $match->feature1->seq;
}

DESCRIPTION

This is a parser of the output of the ps_scan program. It takes either a file handle or a file name, and returns a Bio::SeqFeature::FeaturePair object.

Note that the current implementation parses the entire file at once.

AUTHOR

Juguang Xiao, juguang@tll.org.sg

SEE ALSO

new

Title   : new
Usage   : Bio::Tools::PrositeScan->new(-file => 'out.PrositeScan');
          Bio::Tools::PrositeScan->new(-fh => \*FH);
Returns : L<Bio::Tools::PrositeScan>
Args    : -format => string representing the format type for the
                     ps_scan output, REQUIRED

The -format argument must currently be set to fasta since this is the only parser implemented. This corresponds with using the ps_scan arguments -o fasta.

next_prediction

Title   : new
Usage   :
    while($result = $factory->next_prediction){
        ;
    }

Returns : a Bio::SeqFeature::FeaturePair object where
          feature1 is the matched subsequence and
          feature2 is the PROSITE accession number.
          See <http://prosite.expasy.org/prosuser.html#conv_ac>.