NAME
Bio::Tools::Run::Alignment::Blat - wrapper modules for the Blat program
SYNOPSIS
use Bio::Tools::Run::Alignment::Blat;my $factory = Bio::Tools::Run::Alignment::Blat->new(-db => $database);# $report is a Bio::SearchIO-compliant objectmy $report = $factory->run($seqobj);
DESCRIPTION
Wrapper module for Blat program. This newer version allows for all parameters to be set by passing them as an option to new().
Key bits not implemented yet (TODO):
- •
-
Implement all needed Bio::Tools::Run::WrapperBase methods
Missing are a few, including version().
- •
-
Re-implement using IPC::Run
Would like to get this running under something less reliant on OS-dependent changes within code.
- •
-
No .2bit or .nib conversions yet
These require callouts to faToNib or faTwoTwoBit, which may or may not be installed on a user's machine. We can possibly add functionality to check for faToTwoBit/faToNib and other UCSC tools in the future.
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 one of the Bioperl mailing lists. 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 the bugs and their resolution. Bug reports can be submitted via the web:
http://redmine.open-bio.org/projects/bioperl/
AUTHOR
Chris Fields - cjfields at bioperl dot orgOriginal author - Bala Email [email protected]
APPENDIX
The rest of the documentation details each of the object methods. Internal methods are usually preceded with a _
new
Title : newUsage : $blat->new( -db => '' )Function: Create a new Blat factoryReturns : A new Bio::Tools::Run::Alignment::Blat objectArgs : -db : Mandatory parameter. See db()-qsegment : see qsegment()-tsegment : see tsegment()Also, Blat parameters and flags are accepted: -t, -q, -minIdentity,-minScore, -out, -trimT, ...See Blat's manual for details.
program_name
Title : program_nameUsage : $factory->program_name()Function: Get the program nameReturns : stringArgs : None
program_dir
Title : program_dirUsage : $factory->program_dir(@params)Function: returns the program directory, obtained from ENV variable.Returns : stringArgs :
run
Title : run()Usage : $obj->run($query)Function: Run Blat.Returns : A Bio::SearchIO object that holds the resultsArgs : A Bio::PrimarySeqI object or a file of query sequences
align
Title : alignUsage : $obj->align($query)Function: Alias to run()
db
Title : dbUsage : $obj->db()Function: Get or set the file of database sequences (.fa , .nib or .2bit)Returns : Database filenameArgs : Database filename
qsegment
Title : qsegmentUsage : $obj->qsegment('sequence_a:0-1000')Function : pass in a B<UCSC-compliant> string for the query sequence(s)Returns : stringArgs : stringNote : Requires the sequence(s) in question be 2bit or nib formatReminder : UCSC segment/regions coordinates are 0-based half-open (sequencebegins at 0, but start isn't counted with length), whereas BioPerlcoordinates are 1-based closed (sequence begins with 1, both startand end are counted in the length of the segment). For example, asegment that is 'sequence_a:0-1000' will have BioPerl coordinates of'sequence_a:1-1000', both with the same length (1000).
tsegment
Title : tsegmentUsage : $obj->tsegment('sequence_a:0-1000')Function : pass in a B<UCSC-compliant> string for the target sequence(s)Returns : stringArgs : stringNote : Requires the sequence(s) in question be 2bit or nib formatReminder : UCSC segment/regions coordinates are 0-based half-open (sequencebegins at 0, but start isn't counted with length), whereas BioPerlcoordinates are 1-based closed (sequence begins with 1, both startand end are counted in the length of the segment). For example, asegment that is 'sequence_a:0-1000' will have BioPerl coordinates of'sequence_a:1-1000', both with the same length (1000).
outfile_name
Title : outfile_nameUsage : $obj->outfile_name('out.blat')Function : Get or set the name for the BLAT output fileReturns : stringArgs : string
searchio
Title : searchioUsage : $obj->searchio{-writer => $writer}Function : Pass in additional parameters to the returned Bio::SearchIO parserReturns : Hash reference with Bio::SearchIO parametersArgs : Hash referenceNote : Currently, this implementation overrides any passed -formatparameter based on whether the output is changed ('out'). Thismay change if requested, but we can't see the utility of doing so,as requesting mismatched output/parser combinations is just a recipefor disaster
Bio::ParameterBaseI-specific methods
These methods are part of the Bio::ParameterBaseI interface
set_parameters
Title : set_parametersUsage : $pobj->set_parameters(%params);Function: sets the parameters listed in the hash or arrayReturns : NoneArgs : [optional] hash or array of parameter/values. These can optionallybe hash or array referencesNote : This only sets parameters; to set methods use the method name
reset_parameters
Title : reset_parametersUsage : resets valuesFunction: resets parameters to either undef or value in passed hashReturns : noneArgs : [optional] hash of parameter-value pairs
validate_parameters
Title : validate_parametersUsage : $pobj->validate_parameters(1);Function: sets a flag indicating whether to validate parameters viaset_parameters() or reset_parameters()Returns : BoolArgs : [optional] value evaluating to True/FalseNote : NYI
parameters_changed
Title : parameters_changedUsage : if ($pobj->parameters_changed) {...}Function: Returns boolean true (1) if parameters have changedReturns : Boolean (0 or 1)Args : NoneNote : This module does not run state checks, so this always returns True
available_parameters
Title : available_parametersUsage : @params = $pobj->available_parameters()Function: Returns a list of the available parametersReturns : Array of parametersArgs : [optional] name of executable being used; defaults to returning allavailable parameters
get_parameters
Title : get_parametersUsage : %params = $pobj->get_parameters;Function: Returns list of set key-value pairs, parameter => valueReturns : List of key-value pairsArgs : none
to_* methods
All to_* methods are implementation-specific
to_exe_string
Title : to_exe_stringUsage : $string = $pobj->to_exe_string;Function: Returns string (command line string in this case)Returns : StringArgs :