NAME
Bio::SearchIO::infernal - SearchIO-based Infernal parser
SYNOPSIS
my $parser = Bio::SearchIO->new(-format => 'infernal',-file => 'purine.inf');while( my $result = $parser->next_result ) {# general result info, such as model used, Infernal versionwhile( my $hit = $result->next_hit ) {while( my $hsp = $hit->next_hsp ) {# ...}}}
DESCRIPTION
This is a SearchIO-based parser for Infernal output from the cmsearch program. It currently parses cmsearch output for Infernal versions 0.7-1.1; older versions may work but will not be supported.
The latest version of Infernal is 1.1. The output has changed substantially relative to version 1.0. Versions 1.x are stable releases (and output has stabilized) therefore it is highly recommended that users upgrade to using the latest Infernal release. Support for the older pre-v.1 developer releases will be dropped for future core 1.6 releases.
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.
[email protected] - General discussionhttp://bioperl.org/wiki/Mailing_lists - About the mailing lists
Support
Please direct usage questions or support issues to the mailing list:
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 - Chris Fields
Email cjfields-at-uiuc-dot-edu
CONTRIBUTORS
Jeffrey Barrick, Michigan State UniversityPaul Cantalupo, University of Pittsburgh
APPENDIX
The rest of the documentation details each of the object methods. Internal methods are usually preceded with a _
new
Title : newUsage : my $obj = Bio::SearchIO::infernal->new();Function: Builds a new Bio::SearchIO::infernal objectReturns : Bio::SearchIO::infernalArgs : -fh/-file => cmsearch (infernal) filename-format => 'infernal'-model => query model (Rfam ID) (default undef)-database => database name (default undef)-query_acc => query accession, eg. Rfam accession RF####-query_desc => query description, eg. Rfam description-hsp_minscore => minimum HSP score cutoff-convert_meta => boolean, set to convert meta string to simple WUSS format-symbols => hash ref of structure symbols to use(default symbols in %STRUCTURE_SYMBOLS hash)
next_result
Title : next_resultUsage : my $hit = $searchio->next_result;Function: Returns the next Result from a searchReturns : Bio::Search::Result::ResultI objectArgs : none
start_element
Title : start_elementUsage : $eventgenerator->start_elementFunction: Handles a start element eventReturns : noneArgs : hashref with at least 2 keys 'Data' and 'Name'
end_element
Title : start_elementUsage : $eventgenerator->end_elementFunction: Handles an end element eventReturns : noneArgs : hashref with at least 2 keys, 'Data' and 'Name'
element
Title : elementUsage : $eventhandler->element({'Name' => $name, 'Data' => $str});Function: Convenience method that calls start_element, characters, end_elementReturns : noneArgs : Hash ref with the keys 'Name' and 'Data'
element_hash
Title : elementUsage : $eventhandler->element_hash({'Hsp_hit-from' => $start,'Hsp_hit-to' => $end,'Hsp_score' => $lastscore});Function: Convenience method that takes multiple simple data elements andmaps to appropriate parametersReturns : noneArgs : Hash ref with the mapped key (in %MAPPING) and value
characters
Title : charactersUsage : $eventgenerator->characters($str)Function: Send a character eventsReturns : noneArgs : string
within_element
Title : within_elementUsage : if( $eventgenerator->within_element($element) ) {}Function: Test if we are within a particular elementThis is different than 'in' because within can be testedfor a whole block.Returns : booleanArgs : string element name
in_element
Title : in_elementUsage : if( $eventgenerator->in_element($element) ) {}Function: Test if we are in a particular elementThis is different than 'within' because 'in' onlytests its immediate parent.Returns : booleanArgs : string element name
start_document
Title : start_documentUsage : $eventgenerator->start_documentFunction: Handle a start document eventReturns : noneArgs : none
end_document
Title : end_documentUsage : $eventgenerator->end_documentFunction: Handles an end document eventReturns : Bio::Search::Result::ResultI objectArgs : none
result_count
Title : result_countUsage : my $count = $searchio->result_countFunction: Returns the number of results we have processedReturns : integerArgs : none
model
Title : modelUsage : my $model = $parser->model();Function: Get/Set model; Infernal currently does not outputthe model name (Rfam ID)Returns : String (name of model)Args : [optional] String (name of model)
database
Title : databaseUsage : my $database = $parser->database();Function: Get/Set database; pre-v.1 versions of Infernal do not outputthe database nameReturns : String (database name)Args : [optional] String (database name)
algorithm
Title : algorithmUsage : my $algorithm = $parser->algorithm();Function: Get/Set algorithm; pre-v.1 versions of Infernal do not outputthe algorithm nameReturns : String (algorithm name)Args : [optional] String (algorithm name)
query_accession
Title : query_accessionUsage : my $acc = $parser->query_accession();Function: Get/Set query (model) accession; pre-v1.1 Infernal does not outputthe accession number (Rfam accession #)Returns : String (accession)Args : [optional] String (accession)
query_description
Title : query_descriptionUsage : my $acc = $parser->query_description();Function: Get/Set query (model) description; pre-v1.1 Infernal does not outputthe Rfam descriptionReturns : String (description)Args : [optional] String (description)
hsp_minscore
Title : hsp_minscoreUsage : my $cutoff = $parser->hsp_minscore();Function: Get/Set min bit score cutoff (for generating Hits/HSPs)Returns : score (number)Args : [optional] score (number)
convert_meta
Title : convert_metaUsage : $parser->convert_meta(1);Function: Get/Set boolean flag for converting Infernal WUSS formatto a simple bracketed format (simple WUSS by default)Returns : boolean flag (TRUE or FALSE)Args : [optional] boolean (eval's to TRUE or FALSE)
version
Title : versionUsage : $parser->version();Function: Set the Infernal cmsearch versionReturns : versionArgs : [optional] version
structure_symbols
Title : structure_symbolsUsage : my $hashref = $parser->structure_symbols();Function: Get/Set RNA structure symbolsReturns : Hash ref of delimiters (5' stem, 3' stem, single-strand, etc): default = < (5-prime)> (3-prime): (single-strand)? (unknown). (gap)Args : Hash ref of substitute delimiters, using above keys.
simple_meta
Title : simple_metaUsage : my $string = $parser->simple_meta($str);Function: converts more complex WUSS meta format into simple bracket formatusing symbols defined in structure_symbols()Returns : converted stringArgs : [required] string to convertNote : This is a very simple conversion method to get simple bracketedformat from Infernal data. If the convert_meta() flag is set,this is the method used to convert the strings.