bio tools run tribemcl 1.4
Sponsored Links
Sponsored Links
Secleted [ 0 ] software to compare
Results 1 - 15 of about 5515
Bio::Tools::Run::TribeMCL 1.4
Bio::Tools::Run::TribeMCL is a method for clustering proteins into related groups, which are termed protein families. more>>
Bio::Tools::Run::TribeMCL is a method for clustering proteins into related groups, which are termed protein families.
SYNOPSIS
use Bio::Tools::Run::TribeMCL;
use Bio::SearchIO;
# 3 methods to input the blast results
# straight forward raw blast output (NCBI or WU-BLAST)
my @params = (inputtype=>blastfile);
# OR
# markov program format
# protein_id1 protein_id2 evalue_magnitude evalue_factor
# for example:
# proteins ENSP00000257547 and ENSP00000261659
# with a blast score evalue of 1e-50
# and proteins O42187 and ENSP00000257547
# with a blast score evalue of 1e-119
# entry would be
my @array = [[qw(ENSP00000257547 ENSP00000261659 1 50)],
[qw(O42187 ENSP00000257547 1 119)]];
my @params = (pairs=>@array,I=>2.0);
# OR
# pass in a searchio object
# slowest of the 3 methods as it does more rigourous parsing
# than required for us here
my $sio = Bio::SearchIO->new(-format=>blast,
-file=>blast.out);
my @params=(inputtype=>searchio,I=>2.0);
# you can specify the path to the executable manually in the following way
my @params=(inputtype=>blastfile,I=>2.0,
mcl=>/home/shawn/software/mcl-02-150/src/shmcl/mcl,
matrix=>/home/shawn/software/mcl-02-150/src/contrib/tribe/tribe-matrix);
my $fact = Bio::Tools::Run::TribeMCL->new(@params);
# OR
$fact->matrix_executable(/home/shawn/software/mcl-02-150/src/contrib/tribe/tribe-matrix);
$fact->mcl_executable(/home/shawn/software/mcl-02-150/src/shmcl/mcl);
# to run
my $fact = Bio::Tools::Run::TribeMCL->new(@params);
# Run the program
# returns an array reference to clusters where members are the ids
# for example :2 clusters with 3 members per cluster:
# $fam = [ [mem1 mem2 mem3],[mem1 mem2 mem3]]
# pass in either the blastfile path/searchio obj/the array ref to scores
my $fam = $fact->run($sio);
# print out your clusters
for (my $i = 0; $i
print "Cluster $i t ".scalar(@{$fam->[$i]})." membersn";
foreach my $member (@{$fam->[$i]}){
print "t$membern";
}
}
This clustering is achieved by analysing similarity patterns between proteins in a given dataset, and using these patterns to assign proteins into related groups. In many cases, proteins in the same protein family will have similar functional properties.
Enhancements:
- Perl
<<lessSYNOPSIS
use Bio::Tools::Run::TribeMCL;
use Bio::SearchIO;
# 3 methods to input the blast results
# straight forward raw blast output (NCBI or WU-BLAST)
my @params = (inputtype=>blastfile);
# OR
# markov program format
# protein_id1 protein_id2 evalue_magnitude evalue_factor
# for example:
# proteins ENSP00000257547 and ENSP00000261659
# with a blast score evalue of 1e-50
# and proteins O42187 and ENSP00000257547
# with a blast score evalue of 1e-119
# entry would be
my @array = [[qw(ENSP00000257547 ENSP00000261659 1 50)],
[qw(O42187 ENSP00000257547 1 119)]];
my @params = (pairs=>@array,I=>2.0);
# OR
# pass in a searchio object
# slowest of the 3 methods as it does more rigourous parsing
# than required for us here
my $sio = Bio::SearchIO->new(-format=>blast,
-file=>blast.out);
my @params=(inputtype=>searchio,I=>2.0);
# you can specify the path to the executable manually in the following way
my @params=(inputtype=>blastfile,I=>2.0,
mcl=>/home/shawn/software/mcl-02-150/src/shmcl/mcl,
matrix=>/home/shawn/software/mcl-02-150/src/contrib/tribe/tribe-matrix);
my $fact = Bio::Tools::Run::TribeMCL->new(@params);
# OR
$fact->matrix_executable(/home/shawn/software/mcl-02-150/src/contrib/tribe/tribe-matrix);
$fact->mcl_executable(/home/shawn/software/mcl-02-150/src/shmcl/mcl);
# to run
my $fact = Bio::Tools::Run::TribeMCL->new(@params);
# Run the program
# returns an array reference to clusters where members are the ids
# for example :2 clusters with 3 members per cluster:
# $fam = [ [mem1 mem2 mem3],[mem1 mem2 mem3]]
# pass in either the blastfile path/searchio obj/the array ref to scores
my $fam = $fact->run($sio);
# print out your clusters
for (my $i = 0; $i
print "Cluster $i t ".scalar(@{$fam->[$i]})." membersn";
foreach my $member (@{$fam->[$i]}){
print "t$membern";
}
}
This clustering is achieved by analysing similarity patterns between proteins in a given dataset, and using these patterns to assign proteins into related groups. In many cases, proteins in the same protein family will have similar functional properties.
Enhancements:
- Perl
Download (0.81MB)
Added: 2007-02-22 License: Perl Artistic License Price:
975 downloads
Bio::Tools::Run::PiseWorkflow 1.4
Bio::Tools::Run::PiseWorkflow is a class to create a Pise workflow using Pise application objects as methods. more>>
Bio::Tools::Run::PiseWorkflow is a class to create a Pise workflow using Pise application objects as methods. A workflow is defined by a set of methods which all instanciate the class PiseApplication.
SYNOPSIS
# First, create a Bio::Tools::Run::AnalysisFactory::Pise object:
my $factory = new Bio::Tools::Run::AnalysisFactory::Pise();
# Then create the application objects (Pise::Run::Tools::PiseApplication):
my $clustalw = $factory->program(clustalw);
$clustalw->infile($my_alignment_file);
my $protpars = $factory->program(protpars);
# You can specify different servers for different applications :
my $protdist = $factory->program(protpars
-remote => http://kun.homelinux.com/cgi-bin/Pise/5.a//protpars.pl,
-email => your_email);
# Create a new workflow object :
my $workflow = Bio::Tools::Run::PiseWorkflow->new();
# Define the workflows methods using the application objects:
# the application method $protpars will receive the output of
# type readseq_ok_alig from the application method $clustalw.
$workflow->addpipe(-method => $clustalw,
-tomethod => $protpars,
-pipetype => readseq_ok_alig);
# The application method $clustalw will be piped to a second
# application method ($protdist) using the output of type readseq_ok_alig.
$workflow->addpipe(-method => $clustalw,
-tomethod => $protdist,
-pipetype => readseq_ok_alig);
# The application method $protpars will be piped to the application
# method $consense using the output of type phylip_tree.
my $consense = $factory->program(consense);
$workflow->addpipe(-method => $protpars,
-tomethod => $consense,
-pipetype => phylip_tree);
# Run the workflow.
$workflow->run();
<<lessSYNOPSIS
# First, create a Bio::Tools::Run::AnalysisFactory::Pise object:
my $factory = new Bio::Tools::Run::AnalysisFactory::Pise();
# Then create the application objects (Pise::Run::Tools::PiseApplication):
my $clustalw = $factory->program(clustalw);
$clustalw->infile($my_alignment_file);
my $protpars = $factory->program(protpars);
# You can specify different servers for different applications :
my $protdist = $factory->program(protpars
-remote => http://kun.homelinux.com/cgi-bin/Pise/5.a//protpars.pl,
-email => your_email);
# Create a new workflow object :
my $workflow = Bio::Tools::Run::PiseWorkflow->new();
# Define the workflows methods using the application objects:
# the application method $protpars will receive the output of
# type readseq_ok_alig from the application method $clustalw.
$workflow->addpipe(-method => $clustalw,
-tomethod => $protpars,
-pipetype => readseq_ok_alig);
# The application method $clustalw will be piped to a second
# application method ($protdist) using the output of type readseq_ok_alig.
$workflow->addpipe(-method => $clustalw,
-tomethod => $protdist,
-pipetype => readseq_ok_alig);
# The application method $protpars will be piped to the application
# method $consense using the output of type phylip_tree.
my $consense = $factory->program(consense);
$workflow->addpipe(-method => $protpars,
-tomethod => $consense,
-pipetype => phylip_tree);
# Run the workflow.
$workflow->run();
Download (0.81MB)
Added: 2007-02-22 License: Perl Artistic License Price:
974 downloads
Bio::Root::Object 1.4
Bio::Root::Object is a core Perl 5 object. more>>
Bio::Root::Object is a core Perl 5 object.
SYNOPSIS
# Use this module as the root of your inheritance tree.
Object Creation
require Bio::Root::Object;
$dad = new Bio::Root::Object();
$son = new Bio::Root::Object(-name => Junior,
-parent => $dad,
-make => full);
See the new() method for a complete description of parameters. See also the USAGE section .
Bio::Root::Object attempts to encapsulate the "core" Perl5 object: What are the key data and behaviors ALL (or at least most) Perl5 objects should have?
Rationale
Use of Bio::Root::Object.pm within the Bioperl framework facilitates operational consistency across the different modules defined within the Bio:: namespace. Not all objects need to derive from Bio::Root::Object.pm. However, when generating lots of different types of potentially complex objects which should all conform to a set of basic expectations, this module may be handy.
At the very least, this module saves you from re-writing the new() method for each module you develop. It also permits consistent and robust handling of -tag => value method arguments via the Bio::Root::RootI::_rearrange() method and provides a object-oriented way handle exceptions and warnings via the Bio::Root::Root::throw() and Bio::Root::Root::warn() methods.
See the APPENDIX section for some other handy methods.
Fault-Tolerant Objects
A major motivation for this module was to promote the creation of robust, fault-tolerant Perl5 objects. The Bio::Root::Root::throw() method relies on Perls built-in eval{}/die exception mechanism to generate fatal exceptions. The data comprising an exception is managed by the Bio::Root::Err.pm module, which essentially allows the data thrown by a die() event to be wrapped into an object that can be easily examined and possibly re-thrown.
The intent here is three-fold:
1 Detailed error reporting.
Allow objects to report detailed information about the error condition (who, what, where, why, how).
2 Handle complex errors in objects.
The goal is to make it relatively painless to detect and handle the wide variety of errors possible with a complex Perl object. Perls error handling mechanism is a might clunky when it comes to handling complex errors within complex objects, but it is improving.
3 Efficient & easy exception handling.
To enable robust exception handling without incurring a significant performance penalty in the resulting code. Ideally, exception handling code should be transparent to the cpu until and unless an exception arises.
These goals may at times be at odds and we are not claiming to have achieved the perfect balance. Ultimately, we want self- sufficient object-oriented systems able to deal with their own errors. This area should improve as the module, and Perl, evolve. One possible modification might be to utilize Graham Barrs Error.pm module or Torsten Ekedahls Experimental::Exception.pm module (see "Other Exception Modules").
Technologies such as these may eventually be incorporated into future releases of Perl. The exception handling used by Bio::Root::Object.pm can be expected to change as Perls exception handling mechanism evolves.
TERMINOLOGY NOTE: In this discussion and elsewhere in this module, the terms "Exception" and "Error" are used interchangeably to mean "something unexpected occurred" either as a result of incorrect user input or faulty internal processing.
<<lessSYNOPSIS
# Use this module as the root of your inheritance tree.
Object Creation
require Bio::Root::Object;
$dad = new Bio::Root::Object();
$son = new Bio::Root::Object(-name => Junior,
-parent => $dad,
-make => full);
See the new() method for a complete description of parameters. See also the USAGE section .
Bio::Root::Object attempts to encapsulate the "core" Perl5 object: What are the key data and behaviors ALL (or at least most) Perl5 objects should have?
Rationale
Use of Bio::Root::Object.pm within the Bioperl framework facilitates operational consistency across the different modules defined within the Bio:: namespace. Not all objects need to derive from Bio::Root::Object.pm. However, when generating lots of different types of potentially complex objects which should all conform to a set of basic expectations, this module may be handy.
At the very least, this module saves you from re-writing the new() method for each module you develop. It also permits consistent and robust handling of -tag => value method arguments via the Bio::Root::RootI::_rearrange() method and provides a object-oriented way handle exceptions and warnings via the Bio::Root::Root::throw() and Bio::Root::Root::warn() methods.
See the APPENDIX section for some other handy methods.
Fault-Tolerant Objects
A major motivation for this module was to promote the creation of robust, fault-tolerant Perl5 objects. The Bio::Root::Root::throw() method relies on Perls built-in eval{}/die exception mechanism to generate fatal exceptions. The data comprising an exception is managed by the Bio::Root::Err.pm module, which essentially allows the data thrown by a die() event to be wrapped into an object that can be easily examined and possibly re-thrown.
The intent here is three-fold:
1 Detailed error reporting.
Allow objects to report detailed information about the error condition (who, what, where, why, how).
2 Handle complex errors in objects.
The goal is to make it relatively painless to detect and handle the wide variety of errors possible with a complex Perl object. Perls error handling mechanism is a might clunky when it comes to handling complex errors within complex objects, but it is improving.
3 Efficient & easy exception handling.
To enable robust exception handling without incurring a significant performance penalty in the resulting code. Ideally, exception handling code should be transparent to the cpu until and unless an exception arises.
These goals may at times be at odds and we are not claiming to have achieved the perfect balance. Ultimately, we want self- sufficient object-oriented systems able to deal with their own errors. This area should improve as the module, and Perl, evolve. One possible modification might be to utilize Graham Barrs Error.pm module or Torsten Ekedahls Experimental::Exception.pm module (see "Other Exception Modules").
Technologies such as these may eventually be incorporated into future releases of Perl. The exception handling used by Bio::Root::Object.pm can be expected to change as Perls exception handling mechanism evolves.
TERMINOLOGY NOTE: In this discussion and elsewhere in this module, the terms "Exception" and "Error" are used interchangeably to mean "something unexpected occurred" either as a result of incorrect user input or faulty internal processing.
Download (4.7MB)
Added: 2006-06-16 License: Perl Artistic License Price:
1225 downloads
Bio::Tools::CodonTable 1.4
Bio::Tools::CodonTable is a bioperl codon table object. more>>
Bio::Tools::CodonTable is a bioperl codon table object.
SYNOPSIS
# This is a read-only class for all known codon tables. The IDs are
# the ones used by nucleotide sequence databases. All common IUPAC
# ambiguity codes for DNA, RNA and animo acids are recognized.
# to use
use Bio::Tools::CodonTable;
# defaults to ID 1 "Standard"
$myCodonTable = Bio::Tools::CodonTable->new();
$myCodonTable2 = Bio::Tools::CodonTable -> new ( -id => 3 );
# change codon table
$myCodonTable->id(5);
# examine codon table
print join ( , "The name of the codon table no.", $myCodonTable->id(4),
"is:", $myCodonTable->name(), "n");
# translate a codon
$aa = $myCodonTable->translate(ACU);
$aa = $myCodonTable->translate(act);
$aa = $myCodonTable->translate(ytr);
# reverse translate an amino acid
@codons = $myCodonTable->revtranslate(A);
@codons = $myCodonTable->revtranslate(Ser);
@codons = $myCodonTable->revtranslate(Glx);
@codons = $myCodonTable->revtranslate(cYS, rna);
#boolean tests
print "Is a startn" if $myCodonTable->is_start_codon(ATG);
print "Is a termianatorn" if $myCodonTable->is_ter_codon(tar);
print "Is a unknownn" if $myCodonTable->is_unknown_codon(JTG);
Codon tables are also called translation tables or genetics codes since that is what they try to represent. A bit more complete picture of the full complexity of codon usage in various taxonomic groups presented at the NCBI Genetic Codes Home page.
CodonTable is a BioPerl class that knows all current translation tables that are used by primary nucleotide sequence databases (GenBank, EMBL and DDBJ). It provides methods to output information about tables and relationships between codons and amino acids.
This class and its methods recognized all common IUPAC ambiguity codes for DNA, RNA and animo acids. The translation method follows the conventions in EMBL and TREMBL databases.
It is a nuisance to separate RNA and cDNA representations of nucleic acid transcripts. The CodonTable object accepts codons of both type as input and allows the user to set the mode for output when reverse translating. Its default for output is DNA.
Note: This class deals primarily with individual codons and amino acids. However in the interest of speed you can translate longer sequence, too. The full complexity of protein translation is tackled by Bio::PrimarySeqI::translate.
<<lessSYNOPSIS
# This is a read-only class for all known codon tables. The IDs are
# the ones used by nucleotide sequence databases. All common IUPAC
# ambiguity codes for DNA, RNA and animo acids are recognized.
# to use
use Bio::Tools::CodonTable;
# defaults to ID 1 "Standard"
$myCodonTable = Bio::Tools::CodonTable->new();
$myCodonTable2 = Bio::Tools::CodonTable -> new ( -id => 3 );
# change codon table
$myCodonTable->id(5);
# examine codon table
print join ( , "The name of the codon table no.", $myCodonTable->id(4),
"is:", $myCodonTable->name(), "n");
# translate a codon
$aa = $myCodonTable->translate(ACU);
$aa = $myCodonTable->translate(act);
$aa = $myCodonTable->translate(ytr);
# reverse translate an amino acid
@codons = $myCodonTable->revtranslate(A);
@codons = $myCodonTable->revtranslate(Ser);
@codons = $myCodonTable->revtranslate(Glx);
@codons = $myCodonTable->revtranslate(cYS, rna);
#boolean tests
print "Is a startn" if $myCodonTable->is_start_codon(ATG);
print "Is a termianatorn" if $myCodonTable->is_ter_codon(tar);
print "Is a unknownn" if $myCodonTable->is_unknown_codon(JTG);
Codon tables are also called translation tables or genetics codes since that is what they try to represent. A bit more complete picture of the full complexity of codon usage in various taxonomic groups presented at the NCBI Genetic Codes Home page.
CodonTable is a BioPerl class that knows all current translation tables that are used by primary nucleotide sequence databases (GenBank, EMBL and DDBJ). It provides methods to output information about tables and relationships between codons and amino acids.
This class and its methods recognized all common IUPAC ambiguity codes for DNA, RNA and animo acids. The translation method follows the conventions in EMBL and TREMBL databases.
It is a nuisance to separate RNA and cDNA representations of nucleic acid transcripts. The CodonTable object accepts codons of both type as input and allows the user to set the mode for output when reverse translating. Its default for output is DNA.
Note: This class deals primarily with individual codons and amino acids. However in the interest of speed you can translate longer sequence, too. The full complexity of protein translation is tackled by Bio::PrimarySeqI::translate.
Download (4.7MB)
Added: 2007-08-15 License: Perl Artistic License Price:
801 downloads
Bio::Tools::Run::PiseApplication::freak 1.4
Bio::Tools::Run::PiseApplication::freak is a Perl module. more>>
Bio::Tools::Run::PiseApplication::freak is a Perl module.
Bio::Tools::Run::PiseApplication::freak
Bioperl class for:
FREAK Residue/base frequency table or plot (EMBOSS)
Parameters:
(see also:
http://bioweb.pasteur.fr/seqanal/interfaces/freak.html
for available values):
freak (String)
init (String)
seqall (Sequence)
seqall -- any [sequences] (-seqall)
pipe: seqsfile
letters (String)
Residue letters (-letters)
step (Integer)
Stepping value (-step)
window (Integer)
Averaging window (-window)
plot (Switch)
Produce graphic (-plot)
graph (Excl)
graph (-graph)
outfile (OutFile)
outfile (-outfile)
auto (String)
psouput (String)
<<lessBio::Tools::Run::PiseApplication::freak
Bioperl class for:
FREAK Residue/base frequency table or plot (EMBOSS)
Parameters:
(see also:
http://bioweb.pasteur.fr/seqanal/interfaces/freak.html
for available values):
freak (String)
init (String)
seqall (Sequence)
seqall -- any [sequences] (-seqall)
pipe: seqsfile
letters (String)
Residue letters (-letters)
step (Integer)
Stepping value (-step)
window (Integer)
Averaging window (-window)
plot (Switch)
Produce graphic (-plot)
graph (Excl)
graph (-graph)
outfile (OutFile)
outfile (-outfile)
auto (String)
psouput (String)
Download (0.81MB)
Added: 2006-08-24 License: GPL (GNU General Public License) Price:
1156 downloads
Bio::Tools::Run::PiseApplication::fasta 1.4
Bio::Tools::Run::PiseApplication::fasta is a Bioperl class for sequence database search. more>>
Bio::Tools::Run::PiseApplication::fasta is a Bioperl class for sequence database search.
Parameters:
fasta (Excl)
Fasta program
query (Sequence)
Query sequence File
pipe: seqfile
seqtype (Excl)
Is it a DNA or protein sequence (-n)
protein_db (Excl)
Protein Database
nucleotid_db (Excl)
Nucleotid Database
break_long (Integer)
Break long library sequences into blocks (-N)
ktup (Integer)
ktup : sensitivity and speed of the search (protein:2, DNA:6)
optcut (Integer)
OPTCUT : the threshold for optimization. (-c)
gapinit (Integer)
Penalty for gap initiation (-12 by default for fasta with proteins, -16 for DNA) (-f)
gapext (Integer)
Penalty for gap extention (-2 by default for fasta with proteins, -4 for DNA) (-g)
high_expect (Float)
Maximal expectation value threshold for displaying scores and alignments (-E)
low_expect (Float)
Minimal expectation value threshold for displaying scores and alignments (-F)
nucleotid_match (Integer)
Reward for a nucleotid match (-r)
nucleotid_mismatch (Integer)
Penalty for a nucleotid mismatch (-r)
matrix (Excl)
Scoring matrix file (-s)
X_penalty (Integer)
Penalty for a match to X (independently of the PAM matrix) (-x)
frameshift (Integer)
Penalty for frameshift between codon (fast[xy]/tfast[xy])(-h)
frameshift_within (Integer)
Penalty for frameshift within a codon (fasty/tfasty)(-j)
threeframe (Switch)
Search only the three forward frames (tfasta) (-3)
invert (Switch)
Reverse complement the query sequence (all tfasta) (-i)
genetic_code (Excl)
Use genetic code for translation (tfasta/tfast[xy]/fast[xy]) (-t)
band (Integer)
band-width used for optimization (-y)
swalig (Switch)
unlimited Smith-Waterman alignment for DNA (-A)
noopt (Switch)
no limited optimization (-o)
stat (Excl)
Specify statistical calculation. (-z)
random (Switch)
Estimate stat parameters from shuffled copies of each library sequence (-z)
histogram (Switch)
No histogram (-H)
scores (Integer)
number of similarity scores to be shown (-b)
alns (Integer)
number of alignments to be shown (-d)
html_output (Switch)
HTML output (-m)
markx (Excl)
Alternate display of matches and mismatches in alignments
init1 (Switch)
sequences ranked by the z-score based on the init1 score (-1)
z_score_out (Excl)
Show normalize score as (-B)
showall (Switch)
both sequences are shown in their entirety in alignments (fasta only) (-a)
linlen (Integer)
output line length for sequence alignments (max. 200) (-w)
offsets (String)
Start numbering the aligned sequences at position x1 x2 (2 numbers) (-X)
info (Switch)
Display more information about the library sequence in the alignment (-L)
statfile (OutFile)
Write out the sequence identifier, superfamily number, and similarity scores to this file (-R)
filter (Switch)
Lower case filtering (-S)
outfile (OutFile)
pipe: mview_input
html_outfile (OutFile)
<<lessParameters:
fasta (Excl)
Fasta program
query (Sequence)
Query sequence File
pipe: seqfile
seqtype (Excl)
Is it a DNA or protein sequence (-n)
protein_db (Excl)
Protein Database
nucleotid_db (Excl)
Nucleotid Database
break_long (Integer)
Break long library sequences into blocks (-N)
ktup (Integer)
ktup : sensitivity and speed of the search (protein:2, DNA:6)
optcut (Integer)
OPTCUT : the threshold for optimization. (-c)
gapinit (Integer)
Penalty for gap initiation (-12 by default for fasta with proteins, -16 for DNA) (-f)
gapext (Integer)
Penalty for gap extention (-2 by default for fasta with proteins, -4 for DNA) (-g)
high_expect (Float)
Maximal expectation value threshold for displaying scores and alignments (-E)
low_expect (Float)
Minimal expectation value threshold for displaying scores and alignments (-F)
nucleotid_match (Integer)
Reward for a nucleotid match (-r)
nucleotid_mismatch (Integer)
Penalty for a nucleotid mismatch (-r)
matrix (Excl)
Scoring matrix file (-s)
X_penalty (Integer)
Penalty for a match to X (independently of the PAM matrix) (-x)
frameshift (Integer)
Penalty for frameshift between codon (fast[xy]/tfast[xy])(-h)
frameshift_within (Integer)
Penalty for frameshift within a codon (fasty/tfasty)(-j)
threeframe (Switch)
Search only the three forward frames (tfasta) (-3)
invert (Switch)
Reverse complement the query sequence (all tfasta) (-i)
genetic_code (Excl)
Use genetic code for translation (tfasta/tfast[xy]/fast[xy]) (-t)
band (Integer)
band-width used for optimization (-y)
swalig (Switch)
unlimited Smith-Waterman alignment for DNA (-A)
noopt (Switch)
no limited optimization (-o)
stat (Excl)
Specify statistical calculation. (-z)
random (Switch)
Estimate stat parameters from shuffled copies of each library sequence (-z)
histogram (Switch)
No histogram (-H)
scores (Integer)
number of similarity scores to be shown (-b)
alns (Integer)
number of alignments to be shown (-d)
html_output (Switch)
HTML output (-m)
markx (Excl)
Alternate display of matches and mismatches in alignments
init1 (Switch)
sequences ranked by the z-score based on the init1 score (-1)
z_score_out (Excl)
Show normalize score as (-B)
showall (Switch)
both sequences are shown in their entirety in alignments (fasta only) (-a)
linlen (Integer)
output line length for sequence alignments (max. 200) (-w)
offsets (String)
Start numbering the aligned sequences at position x1 x2 (2 numbers) (-X)
info (Switch)
Display more information about the library sequence in the alignment (-L)
statfile (OutFile)
Write out the sequence identifier, superfamily number, and similarity scores to this file (-R)
filter (Switch)
Lower case filtering (-S)
outfile (OutFile)
pipe: mview_input
html_outfile (OutFile)
Download (0.081MB)
Added: 2007-06-06 License: Perl Artistic License Price:
872 downloads
Bio::Tools::Run::PiseApplication::charge 1.4
Bio::Tools::Run::PiseApplication::charge is a Perl module. more>>
Bio::Tools::Run::PiseApplication::charge is a Perl module.
Bio::Tools::Run::PiseApplication::charge
Bioperl class for:
CHARGE Protein charge plot (EMBOSS)
Parameters:
(see also:
http://bioweb.pasteur.fr/seqanal/interfaces/charge.html
for available values):
charge (String)
init (String)
seqall (Sequence)
seqall -- protein [sequences] (-seqall)
pipe: seqsfile
window (Integer)
Window (-window)
aadata (String)
Amino acid property data file name (-aadata)
plot (Switch)
Produce graphic (-plot)
graph (Excl)
graph (-graph)
outfile (OutFile)
outfile (-outfile)
auto (String)
psouput (String)
<<lessBio::Tools::Run::PiseApplication::charge
Bioperl class for:
CHARGE Protein charge plot (EMBOSS)
Parameters:
(see also:
http://bioweb.pasteur.fr/seqanal/interfaces/charge.html
for available values):
charge (String)
init (String)
seqall (Sequence)
seqall -- protein [sequences] (-seqall)
pipe: seqsfile
window (Integer)
Window (-window)
aadata (String)
Amino acid property data file name (-aadata)
plot (Switch)
Produce graphic (-plot)
graph (Excl)
graph (-graph)
outfile (OutFile)
outfile (-outfile)
auto (String)
psouput (String)
Download (0.81MB)
Added: 2006-08-24 License: Perl Artistic License Price:
1156 downloads
Bio::Tools::Run::JavaRunner 1.5.2_100
Bio::Tools::Run::JavaRunner is a Perl module that can run java programs. more>>
Bio::Tools::Run::JavaRunner is a Perl module that can run java programs.
SYNOPSIS
my $runner = Bio::Tools::Run::JavaRunner->new(-jar => $jar);
$runner->run();
This module is probably incomplete. It is intended to be a wrapper for running java programs.
<<lessSYNOPSIS
my $runner = Bio::Tools::Run::JavaRunner->new(-jar => $jar);
$runner->run();
This module is probably incomplete. It is intended to be a wrapper for running java programs.
Download (0.92MB)
Added: 2007-06-05 License: Perl Artistic License Price:
872 downloads
Bio::Tools::Run::PiseApplication::align2model 1.4
Bio::Tools::Run::PiseApplication::align2model is a Bioperl class for align2model - create a multiple alignment of sequences... more>>
Bio::Tools::Run::PiseApplication::align2model is a Bioperl class for align2model - create a multiple alignment of sequences to an existing model.
Parameters:
align2model (String)
run (String)
Run name
db (Sequence)
Sequences to align (-db)
model_file (InFile)
Model (-i)
pipe: sam_model
id (String)
Sequence identifier(s) selection (separated by commas) (-id)
nscoreseq (Integer)
Maximum number of sequences to be read (-nscoreseq)
adpstyle (Excl)
Dynamic programming style (-adpstyle
SW (Excl)
Sequence scoring (-SW)
auto_fim (Switch)
Add FIMs automatically (-auto_fim)
jump_in_prob (Float)
Probability cost of jumping into the center of the model (-jump_in_prob)
jump_out_prob (Float)
Probability cost of jumping out the center of the model (-jump_out_prob)
a2mdots (Switch)
Print dots to fill space need for other sequences insertions (-a2mdots)
dump_parameters (Excl)
(-dump_parameters)
<<lessParameters:
align2model (String)
run (String)
Run name
db (Sequence)
Sequences to align (-db)
model_file (InFile)
Model (-i)
pipe: sam_model
id (String)
Sequence identifier(s) selection (separated by commas) (-id)
nscoreseq (Integer)
Maximum number of sequences to be read (-nscoreseq)
adpstyle (Excl)
Dynamic programming style (-adpstyle
SW (Excl)
Sequence scoring (-SW)
auto_fim (Switch)
Add FIMs automatically (-auto_fim)
jump_in_prob (Float)
Probability cost of jumping into the center of the model (-jump_in_prob)
jump_out_prob (Float)
Probability cost of jumping out the center of the model (-jump_out_prob)
a2mdots (Switch)
Print dots to fill space need for other sequences insertions (-a2mdots)
dump_parameters (Excl)
(-dump_parameters)
Download (0.81MB)
Added: 2007-06-06 License: Perl Artistic License Price:
870 downloads
Bio::Ontology::GOterm 1.4
Bio::Ontology::GOterm is a representation of GO terms. more>>
Bio::Ontology::GOterm is a representation of GO terms.
SYNOPSIS
$term = Bio::Ontology::GOterm->new
( -go_id => "GO:0016847",
-name => "1-aminocyclopropane-1-carboxylate synthase",
-definition => "Catalysis of ...",
-is_obsolete => 0,
-comment => "" );
$term->add_definition_references( @refs );
$term->add_secondary_GO_ids( @ids );
$term->add_aliases( @aliases );
foreach my $dr ( $term->each_definition_reference() ) {
print $dr, "n";
}
# etc.
This is "dumb" class for GO terms (it provides no functionality related to graphs). Implements Bio::Ontology::TermI.
<<lessSYNOPSIS
$term = Bio::Ontology::GOterm->new
( -go_id => "GO:0016847",
-name => "1-aminocyclopropane-1-carboxylate synthase",
-definition => "Catalysis of ...",
-is_obsolete => 0,
-comment => "" );
$term->add_definition_references( @refs );
$term->add_secondary_GO_ids( @ids );
$term->add_aliases( @aliases );
foreach my $dr ( $term->each_definition_reference() ) {
print $dr, "n";
}
# etc.
This is "dumb" class for GO terms (it provides no functionality related to graphs). Implements Bio::Ontology::TermI.
Download (4.7MB)
Added: 2006-10-10 License: Perl Artistic License Price:
1112 downloads
Bio::Tools::Run::PiseApplication::consensus 1.4
Bio::Tools::Run::PiseApplication::consensus is a Perl module. more>>
Bio::Tools::Run::PiseApplication::consensus is a Perl module.
Bioperl class for:
CONSENSUS Identification of consensus patterns in unaligned DNA and protein sequences (Hertz, Stormo)
References:
G.Z. Hertz and G.D. Stormo. Identification of consensus patterns in unaligned DNA and protein sequences: a large-deviation statistical basis for penalizing gaps. In: Proceedings of the Third International Conference on Bioinformatics and Genome Research (H.A. Lim, and C.R. Cantor, editors). World Scientific Publishing Co., Ltd. Singapore, 1995. pages 201--216.
Parameters:
(see also:
http://bioweb.pasteur.fr/seqanal/interfaces/consensus.html
for available values):
consensus (Excl)
Program to run
sequence (Sequence)
Sequences file (-f)
pipe: seqsfile
width (Integer)
Width of pattern (consensus only) (-L)
out (String)
consensus_matrix (String)
complement (Excl)
Complement of nucleic acid sequences (-c)
ascii_alphabet (InFile)
Alphabet and normalization information (if not DNA) (-a)
prior (Switch)
Use the designated prior probabilities of the letters to override the observed frequencies (-d)
dna (Switch)
Alphabet and normalization information for DNA
protein (Switch)
Alphabet and normalization information for protein
queue (Integer)
Maximum number of matrices to save between cycles of the program -- ie: queue size (-q)
standard_deviation (Float)
Number of standard deviations to lower the information content at each position before identifying information peaks (mandatory for wconsensus) (-s)
progeny (Excl)
Save the top progeny matrices (-pr1)
linearly (Switch)
Seed with the first sequence and proceed linearly through the list (-l)
max_cycle_nb (Integer)
Maximum repeat of the matrix building cycle (-n or -N)
max_cycle (Excl)
How many words per matrix for each sequence to contribute (-n or -N)
distance (Integer)
Minimum distance between the starting points of words within the same matrix pattern (-m)
terminate (Integer)
Terminate the program this number of cycles after the current most significant alignment is identified (-t)
terminal_gap (Excl)
Permit terminal gaps (-pg) (wconsensus only)
top_matrices (Integer)
Number of top matrices to print (-pt)
final_matrices (Integer)
Number of final matrices to print (-pf)
<<lessBioperl class for:
CONSENSUS Identification of consensus patterns in unaligned DNA and protein sequences (Hertz, Stormo)
References:
G.Z. Hertz and G.D. Stormo. Identification of consensus patterns in unaligned DNA and protein sequences: a large-deviation statistical basis for penalizing gaps. In: Proceedings of the Third International Conference on Bioinformatics and Genome Research (H.A. Lim, and C.R. Cantor, editors). World Scientific Publishing Co., Ltd. Singapore, 1995. pages 201--216.
Parameters:
(see also:
http://bioweb.pasteur.fr/seqanal/interfaces/consensus.html
for available values):
consensus (Excl)
Program to run
sequence (Sequence)
Sequences file (-f)
pipe: seqsfile
width (Integer)
Width of pattern (consensus only) (-L)
out (String)
consensus_matrix (String)
complement (Excl)
Complement of nucleic acid sequences (-c)
ascii_alphabet (InFile)
Alphabet and normalization information (if not DNA) (-a)
prior (Switch)
Use the designated prior probabilities of the letters to override the observed frequencies (-d)
dna (Switch)
Alphabet and normalization information for DNA
protein (Switch)
Alphabet and normalization information for protein
queue (Integer)
Maximum number of matrices to save between cycles of the program -- ie: queue size (-q)
standard_deviation (Float)
Number of standard deviations to lower the information content at each position before identifying information peaks (mandatory for wconsensus) (-s)
progeny (Excl)
Save the top progeny matrices (-pr1)
linearly (Switch)
Seed with the first sequence and proceed linearly through the list (-l)
max_cycle_nb (Integer)
Maximum repeat of the matrix building cycle (-n or -N)
max_cycle (Excl)
How many words per matrix for each sequence to contribute (-n or -N)
distance (Integer)
Minimum distance between the starting points of words within the same matrix pattern (-m)
terminate (Integer)
Terminate the program this number of cycles after the current most significant alignment is identified (-t)
terminal_gap (Excl)
Permit terminal gaps (-pg) (wconsensus only)
top_matrices (Integer)
Number of top matrices to print (-pt)
final_matrices (Integer)
Number of final matrices to print (-pf)
Download (0.81MB)
Added: 2006-10-10 License: Perl Artistic License Price:
1109 downloads
Bio::Biblio 1.4
Bio::Biblio is a Bibliographic Query Service module. more>>
Bio::Biblio is a Bibliographic Query Service module.
SYNOPSIS
use Bio::Biblio;
my $biblio = new Bio::Biblio;
print $biblio->find (perl)->get_count . "n";
my $collection = $biblio->find (brazma, authors);
while ( $collection->has_next ) {
print $collection->get_next;
}
#The new() method can get parameters, for example:
$biblio = Bio::Biblio
(-access => soap,
-location => http://industry.ebi.ac.uk/soap/openBQS,
-destroy_on_exit => 0);
# See below for some one-liners
This is a class whose instances can access bibliographic repositories. It allows to query a bibliographic database (such as MEDLINE) and then to retrieve resulting citations from it. The citations are returned in an XML format which is native to the repository but there are also supporting modules for converting them into Perl objects.
The detailed descriptions of all query and retrieval methods are in Bio::DB::BiblioI (an interface). All those methods should be called on instances of this (Bio::Biblio) module.
The module complies (with some simplifications) with the specification described in the OpenBQS project. Its home page is at http://industry.ebi.ac.uk/openBQS. There are also links to available servers providing access to the bibliographic repositories (namely to MEDLINE).
The module also gives an access to a set of controlled vocabularies and their values. It allows to introspect bibliographic repositories and to find what citation resource types (such as journal and book articles, patents or technical reports) are provided, and what attributes they have, eventually what attribute values are allowed.
Here are some one-liners:
perl -MBio::Biblio -e print new Bio::Biblio->get_by_id ("94033980")
perl -MBio::Biblio
-e print join ("n", @{ new Bio::Biblio->find ("brazma")->get_all_ids })
perl -MBio::Biblio
-e print new Bio::Biblio->find ("Java")->find ("perl")->get_count
<<lessSYNOPSIS
use Bio::Biblio;
my $biblio = new Bio::Biblio;
print $biblio->find (perl)->get_count . "n";
my $collection = $biblio->find (brazma, authors);
while ( $collection->has_next ) {
print $collection->get_next;
}
#The new() method can get parameters, for example:
$biblio = Bio::Biblio
(-access => soap,
-location => http://industry.ebi.ac.uk/soap/openBQS,
-destroy_on_exit => 0);
# See below for some one-liners
This is a class whose instances can access bibliographic repositories. It allows to query a bibliographic database (such as MEDLINE) and then to retrieve resulting citations from it. The citations are returned in an XML format which is native to the repository but there are also supporting modules for converting them into Perl objects.
The detailed descriptions of all query and retrieval methods are in Bio::DB::BiblioI (an interface). All those methods should be called on instances of this (Bio::Biblio) module.
The module complies (with some simplifications) with the specification described in the OpenBQS project. Its home page is at http://industry.ebi.ac.uk/openBQS. There are also links to available servers providing access to the bibliographic repositories (namely to MEDLINE).
The module also gives an access to a set of controlled vocabularies and their values. It allows to introspect bibliographic repositories and to find what citation resource types (such as journal and book articles, patents or technical reports) are provided, and what attributes they have, eventually what attribute values are allowed.
Here are some one-liners:
perl -MBio::Biblio -e print new Bio::Biblio->get_by_id ("94033980")
perl -MBio::Biblio
-e print join ("n", @{ new Bio::Biblio->find ("brazma")->get_all_ids })
perl -MBio::Biblio
-e print new Bio::Biblio->find ("Java")->find ("perl")->get_count
Download (4.7MB)
Added: 2007-06-11 License: Perl Artistic License Price:
866 downloads
Bio::Tools::AlignFactory 1.4
Bio::Tools::AlignFactory is a base object for alignment factories. more>>
Bio::Tools::AlignFactory is a base object for alignment factories.
SYNOPSIS
You wont be using this as an object, but using a dervied class like Bio::Tools::pSW
Holds common Alignment Factory attributes in place
<<lessSYNOPSIS
You wont be using this as an object, but using a dervied class like Bio::Tools::pSW
Holds common Alignment Factory attributes in place
Download (4.7MB)
Added: 2006-10-10 License: Perl Artistic License Price:
1111 downloads
Bio::Seq 1.4
Bio::Seq is a sequence object, with features. more>>
Bio::Seq is a sequence object, with features.
SYNOPSIS
# This is the main sequence object in Bioperl
# gets a sequence from a file
$seqio = Bio::SeqIO->new( -format => embl , -file => myfile.dat);
$seqobj = $seqio->next_seq();
# SeqIO can both read and write sequences; see Bio::SeqIO
# for more information and examples
# get from database
$db = Bio::DB::GenBank->new();
$seqobj = $db->get_Seq_by_acc(X78121);
# make from strings in script
$seqobj = Bio::Seq->new( -display_id => my_id,
-seq => $sequence_as_string);
# gets sequence as a string from sequence object
$seqstr = $seqobj->seq(); # actual sequence as a string
$seqstr = $seqobj->subseq(10,50); # slice in biological coordinates
# retrieves information from the sequence
# features must implement Bio::SeqFeatureI interface
@features = $seqobj->get_SeqFeatures(); # just top level
foreach my $feat ( @features ) {
print "Feature ",$feat->primary_tag," starts ",$feat->start," ends ",
$feat->end," strand ",$feat->strand,"n";
# features retain link to underlying sequence object
print "Feature sequence is ",$feat->seq->seq(),"n"
}
# sequences may have a species
if( defined $seq->species ) {
print "Sequence is from ",$species->binomial_name," [",$species->common_name,"]n";
}
# annotation objects are Bio::AnnotationCollectionIs
$ann = $seqobj->annotation(); # annotation object
# references is one type of annotations to get. Also get
# comment and dblink. Look at Bio::AnnotationCollection for
# more information
foreach my $ref ( $ann->get_Annotations(reference) ) {
print "Reference ",$ref->title,"n";
}
# you can get truncations, translations and reverse complements, these
# all give back Bio::Seq objects themselves, though currently with no
# features transfered
my $trunc = $seqobj->trunc(100,200);
my $rev = $seqobj->revcom();
# there are many options to translate - check out the docs
my $trans = $seqobj->translate();
# these functions can be chained together
my $trans_trunc_rev = $seqobj->trunc(100,200)->revcom->translate();
A Seq object is a sequence with sequence features placed on it. The Seq object contains a PrimarySeq object for the actual sequence and also implements its interface.
In Bioperl we have 3 main players that people are going to use frequently
Bio::PrimarySeq - just the sequence and its names, nothing else.
Bio::SeqFeatureI - a location on a sequence, potentially with a sequence
and annotation.
Bio::Seq - A sequence and a collection of sequence features
(an aggregate) with its own annotation.
Although Bioperl is not tied heavily to file formats these distinctions do map to file formats sensibly and for some bioinformaticians this might help
Bio::PrimarySeq - Fasta file of a sequence
Bio::SeqFeatureI - A single entry in an EMBL/GenBank/DDBJ feature table
Bio::Seq - A single EMBL/GenBank/DDBJ entry
By having this split we avoid a lot of nasty circular references (sequence features can hold a reference to a sequence without the sequence holding a reference to the sequence feature). See Bio::PrimarySeq and Bio::SeqFeatureI for more information.
Ian Korf really helped in the design of the Seq and SeqFeature system.
<<lessSYNOPSIS
# This is the main sequence object in Bioperl
# gets a sequence from a file
$seqio = Bio::SeqIO->new( -format => embl , -file => myfile.dat);
$seqobj = $seqio->next_seq();
# SeqIO can both read and write sequences; see Bio::SeqIO
# for more information and examples
# get from database
$db = Bio::DB::GenBank->new();
$seqobj = $db->get_Seq_by_acc(X78121);
# make from strings in script
$seqobj = Bio::Seq->new( -display_id => my_id,
-seq => $sequence_as_string);
# gets sequence as a string from sequence object
$seqstr = $seqobj->seq(); # actual sequence as a string
$seqstr = $seqobj->subseq(10,50); # slice in biological coordinates
# retrieves information from the sequence
# features must implement Bio::SeqFeatureI interface
@features = $seqobj->get_SeqFeatures(); # just top level
foreach my $feat ( @features ) {
print "Feature ",$feat->primary_tag," starts ",$feat->start," ends ",
$feat->end," strand ",$feat->strand,"n";
# features retain link to underlying sequence object
print "Feature sequence is ",$feat->seq->seq(),"n"
}
# sequences may have a species
if( defined $seq->species ) {
print "Sequence is from ",$species->binomial_name," [",$species->common_name,"]n";
}
# annotation objects are Bio::AnnotationCollectionIs
$ann = $seqobj->annotation(); # annotation object
# references is one type of annotations to get. Also get
# comment and dblink. Look at Bio::AnnotationCollection for
# more information
foreach my $ref ( $ann->get_Annotations(reference) ) {
print "Reference ",$ref->title,"n";
}
# you can get truncations, translations and reverse complements, these
# all give back Bio::Seq objects themselves, though currently with no
# features transfered
my $trunc = $seqobj->trunc(100,200);
my $rev = $seqobj->revcom();
# there are many options to translate - check out the docs
my $trans = $seqobj->translate();
# these functions can be chained together
my $trans_trunc_rev = $seqobj->trunc(100,200)->revcom->translate();
A Seq object is a sequence with sequence features placed on it. The Seq object contains a PrimarySeq object for the actual sequence and also implements its interface.
In Bioperl we have 3 main players that people are going to use frequently
Bio::PrimarySeq - just the sequence and its names, nothing else.
Bio::SeqFeatureI - a location on a sequence, potentially with a sequence
and annotation.
Bio::Seq - A sequence and a collection of sequence features
(an aggregate) with its own annotation.
Although Bioperl is not tied heavily to file formats these distinctions do map to file formats sensibly and for some bioinformaticians this might help
Bio::PrimarySeq - Fasta file of a sequence
Bio::SeqFeatureI - A single entry in an EMBL/GenBank/DDBJ feature table
Bio::Seq - A single EMBL/GenBank/DDBJ entry
By having this split we avoid a lot of nasty circular references (sequence features can hold a reference to a sequence without the sequence holding a reference to the sequence feature). See Bio::PrimarySeq and Bio::SeqFeatureI for more information.
Ian Korf really helped in the design of the Seq and SeqFeature system.
Download (4.7MB)
Added: 2006-10-10 License: Perl Artistic License Price:
1111 downloads
Bio::Ontology::Relationship 1.4
Bio::Ontology::Relationship is a relationship for an ontology. more>>
Bio::Ontology::Relationship is a relationship for an ontology.
SYNOPSIS
$rel = Bio::Ontology::Relationship->new( -identifier => "16847",
-subject_term => $subj,
-object_term => $obj,
-predicate_term => $pred );
This is a basic implementation of Bio::Ontology::RelationshipI.
The terminology we use here is the one commonly used for ontologies, namely the triple of (subject, predicate, object), which in addition is scoped in a namespace (ontology). It is called triple because it is a tuple of three ontology terms.
There are other terminologies in use for expressing relationships. For those who it helps to better understand the concept, the triple of (child, relationship type, parent) would be equivalent to the terminology chosen here, disregarding the question whether the notion of parent and child is sensible in the context of the relationship type or not. Especially in the case of ontologies with a wide variety of predicates the parent/child terminology and similar ones can quickly become ambiguous (e.g., A synthesises B), meaningless (e.g., A binds B), or even conflicting (e.g., A is-parent-of B), and are therefore strongly discouraged.
<<lessSYNOPSIS
$rel = Bio::Ontology::Relationship->new( -identifier => "16847",
-subject_term => $subj,
-object_term => $obj,
-predicate_term => $pred );
This is a basic implementation of Bio::Ontology::RelationshipI.
The terminology we use here is the one commonly used for ontologies, namely the triple of (subject, predicate, object), which in addition is scoped in a namespace (ontology). It is called triple because it is a tuple of three ontology terms.
There are other terminologies in use for expressing relationships. For those who it helps to better understand the concept, the triple of (child, relationship type, parent) would be equivalent to the terminology chosen here, disregarding the question whether the notion of parent and child is sensible in the context of the relationship type or not. Especially in the case of ontologies with a wide variety of predicates the parent/child terminology and similar ones can quickly become ambiguous (e.g., A synthesises B), meaningless (e.g., A binds B), or even conflicting (e.g., A is-parent-of B), and are therefore strongly discouraged.
Download (4.7MB)
Added: 2007-08-09 License: Perl Artistic License Price:
806 downloads
Secleted [ 0 ] software to compare
Copyright Notice:
Software piracy is theft, Using crack, password, serial numbers, registration codes, key generators is illegal and prevent future software development. The above bio tools run tribemcl 1.4 search only lists software in full, demo and trial versions for free download. Download links are directly from our mirror sites or publisher sites, torrent files or links from rapidshare.com, yousendit.com or megaupload.com are not allowed