proteinortho5(1) orthology detection tool


proteinortho5 [,OPTIONS/] ,FASTA1 FASTA2 /[,FASTA/...]


Proteinortho is a stand-alone tool that is geared towards large datasets and makes use of distributed computing techniques when run on multi-core hardware. It implements an extended version of the reciprocal best alignment heuristic. Proteinortho was applied to compute orthologous proteins in the complete set of all 717 eubacterial genomes available at NCBI at the beginning of 2009. Authors succeeded identifying thirty proteins present in 99% of all bacterial proteomes.


E-value for blast [default: 1e-05]
blast program {blastn|blastp|blastn+|blastp+} [default: blastp+]
prefix for all result file names [default: myproject]
activate PoFF extension to separate similar sequences by contextual adjacencies (requires .gff for each .fasta)
PoFF: number of reiterations for adjacencies heuristic, to determine duplicated regions (default: 0)
PoFF: Size of a maximum common substring (MCS) for adjacency matches (default: 3)
PoFF: weight of adjacencies vs. sequence similarity (default: 0.5)
write description files (for NCBI FASTA input only)
stores temporary blast results for reuse
forces recalculation of blast results in any case
number of processors to use [default: auto]
apply selfblast, detects paralogs without orthologs
report singleton genes without any hit
min. percent identity of best blast hits [default: 25]
min. coverage of best blast alignments in % [default: 50]
min. algebraic connectivity [default: 0.1]
min. similarity for additional hits (0..1) [default: 0.95]
1 -> generate indices 2 -> run blast (and ff-adj, if -synteny is set) 3 -> clustering 0 -> all (default)
path to your local blast (if not installed globally)
keeps you informed about the progress
remove all unnecessary files after processing
generate .graph files (pairwise orthology relations)
gives detailed information for bug tracking

More specific blast parameters can be defined by

-blastParameters='[parameters]' (e.g. -blastParameters='-seg no')

In case jobs should be distributed onto several machines, use

-startat= File number to start with (default: 0)
-stopat= File number to end with (default: -1)