Bio::S...er_pull(3pm)
π NAME π
Bio::SearchIO::hmmer_pull β bestie, it's a parser for HMMER output β¨
πΈ SYNOPSIS πΈ
# do not use this class directly it is available through Bio::SearchIO
use Bio::SearchIO;
my $in = Bio::SearchIO->new(-format => 'hmmer_pull',
-file => 't/data/hmmpfam.bigout');
while (my $result = $in->next_result) {
# this is a Bio::Search::Result::HmmpfamResult object
print $result->query_name(), " for HMM ", $result->hmm_name(), "\n";
while (my $hit = $result->next_hit) {
print $hit->name(), "\n";
while (my $hsp = $hit->next_hsp) {
print "length is ", $hsp->length(), "\n";
}
}
}
π«Ά okay so listen: do not use this class directly, it is available through Bio::SearchIO. That's the vibe, that's the rule.
β¨ DESCRIPTION β¨
This object implements a pull-parser for HMMER output π§ It's fast since it only does work on request (hence 'pull') β like, it's literally not gonna lift a finger until you ask it to, and honestly? iconic behaviour π
π FEEDBACK π
π Mailing Lists π
User feedback is an integral part of the evolution of this and other Bioperl modules. Send your comments and suggestions preferably to the Bioperl mailing list. Your participation is much appreciated. π₯Ή
bioperl-l@bioperl.org - General discussion
http://bioperl.org/wiki/Mailing_lists - About the mailing lists
π Support π
Please direct usage questions or support issues to the mailing list:
bioperl-l@bioperl.org
rather than to the module maintainer directly. Many experienced and reponsive experts will be able look at the problem and quickly address it. Please include a thorough description of the problem with code and data examples if at all possible. π«Ά
π Reporting Bugs π
Report bugs to the Bioperl bug tracking system to help us keep track of the bugs and their resolution. Bug reports can be submitted via the web:
https://github.com/bioperl/bioperl-live/issues
π AUTHOR - Sendu Bala π
Email bix@sendu.me.uk π
π APPENDIX π
The rest of the documentation details each of the object methods. Internal methods are usually preceded with a _ π
π new π
Title : new
Usage : my $obj = Bio::SearchIO::hmmer_pull->new();
Function: Builds a new Bio::SearchIO::hmmer_pull object
Returns : Bio::SearchIO::hmmer_pull
Args : -fh/-file => HMMER output filename
-format => 'hmmer_pull'
-evalue => float or scientific notation number to be used
as an evalue cutoff for hits
-score => integer or scientific notation number to be used
as a score value cutoff for hits
-hsps => integer minimum number of hsps (domains) a hit must have
-piped_behaviour => 'temp_file'|'memory'|'sequential_read'
-piped_behaviour defines what the parser should do if the input is
an unseekable filehandle (eg. piped input), see
Bio::PullParserI::chunk for details. Default is 'sequential_read'.
So like, you hand it a -fh/-file (that's your HMMER output filename, bestie), you tell it -format => 'hmmer_pull', and then there's the whole lil squad of cutoffs: -evalue (float or scientific notation number to be used as an evalue cutoff for hits), -score (integer or scientific notation number to be used as a score value cutoff for hits), and -hsps (integer minimum number of hsps (domains) a hit must have) π
And -piped_behaviour is 'temp_file'|'memory'|'sequential_read' β it defines what the parser should do if the input is an unseekable filehandle (eg. piped input), see Bio::PullParserI::chunk for details. Default is 'sequential_read'. π«Ά
π¦ next_result π¦
Title : next_result
Usage : my $hit = $searchio->next_result;
Function: Returns the next Result from a search
Returns : Bio::Search::Result::ResultI object
Args : none
β result_count β
Title : result_count
Usage : my $count = $searchio->result_count
Function: Returns the number of results we have processed.
Returns : integer
Args : none
π rewind π
Title : rewind
Usage : $searchio->rewind;
Function: Allow one to reset the Result iterator to the beginning, so that
next_result() will subsequently return the first result and so on.
NB: result objects are not cached, so you will get new result objects
each time you rewind. Also, note that result_count() counts the
number of times you have called next_result(), so will not be able
tell you how many results there were in the file if you use rewind().
Returns : n/a
Args : none
π¨ okay bestie, read the fine print: result objects are not cached, so you will get new result objects each time you rewind. Also, note thatresult_count()counts the number of times you have callednext_result(), so will not be able tell you how many results there were in the file if you userewind(). fr fr, it's giving "the count is a lie" π
and that's the whole page, babe π this was perl v5.30.0, dated 2019-10-29, from the Bio::S...er_pull(3pm) era β mwah, slay, go parse some HMMER β¨
π¦ where this page comes from
- debian bookworm: libbio-searchio-hmmer-perl 1.7.3-2 Β· Artistic-1.0-Perl OR GPL-1.0-or-later
- debian trixie: libbio-searchio-hmmer-perl 1.7.3-2 Β· Artistic-1.0-Perl OR GPL-1.0-or-later
- ubuntu jammy: libbio-searchio-hmmer-perl 1.7.3-2 Β· Artistic-1.0-Perl OR GPL-1.0-or-later
- ubuntu noble: libbio-searchio-hmmer-perl 1.7.3-2 Β· Artistic-1.0-Perl OR GPL-1.0-or-later
- π girlified by deepseek-ai/DeepSeek-V4.1-Flash on : a language model's rewrite of the original, under its terms