DOKK / manpages / debian 12 / libbio-perl-perl / Bio::Tools::PrositeScan.3pm.en
Bio::Tools::PrositeScan(3pm) User Contributed Perl Documentation Bio::Tools::PrositeScan(3pm)

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

  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;
  }

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.

Juguang Xiao, juguang@tll.org.sg

  • ps_scan software <ftp://ftp.expasy.org/databases/prosite/ps_scan>
  • PROSITE User Manual <http://prosite.expasy.org/prosuser.html>

  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".

  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>.
2021-08-15 perl v5.32.1