REvolutionH-tl: Reconstruction of Evolutionary Histories tool
Project description
Bioinformatics tool for the reconstruction of evolutionary histories. Input: fasta files or sequence alignment hits, Output: orthology. event-labeled gene trees, and reconciliations.
Bioinformatics & complex networks lab
- José Antonio Ramírez-Rafael [jose.ramirezra@cinvestav.mx]
- Maribel Hernandez-Rosales [maribel.hr@cinvestav.mx ]
Steps
- Orthogroup & best hit selection. Input: alignment hits (generate this using
revolutionhtl.diamond
) . - Orthology and gene tree reconstruction. Input: best hits (generate this at step 1).
- Species tree reconstruction. Input: gene trees (generate this at step 2).
- Tree reconciliation. Input: gene and species trees (generate this at steps 2 and 3).
Install
pip install revolutionhtl
Requirements
If you want to run sequence alignments using revolutionhtl, then install Diamond.
Usage
python -m revolutionhtl <arguments>
See bellow the full list of arguments
Input data (Click to expand)
- -h --helpshow this help message and exit
-steps [int ...]
List of steps to run (default: 1 2 3 4).
-alignment_h str --alignment_hits str
Directory containing alignment hits, the input of step 1. (default: ./).
-best_h str --best_hits str
.tsv file containing best hits, the input of step 2. (default: use output of step 1).
-T str --gene_trees str
.tsv file containing gene trees, the input of steps 3 and 4. (default: use output of step 2).
-S str --species_tree str
.nhx file containing a species tree, an input of step 4. (default: use output of step 3).
File names (Click to expand)
-o str --output_prefix strPrefix used for output files (default "tl_project").
-og str --orthogroup_column str
Column in -best_h -T, and output files specifying orthogroups (default: OG).
-Nm int --N_max int
Indicates the maximum number of genes in a orthogroup, bigger orthogroups are splitted. If 0, no orthogroup is splitted. (default= 2000).
-k int --k_size_partition int
Integer indicatng how many best hit graphs will be processed in bunch:: first graphs with
Algorithm parameters (Click to expand)
-bh_heuristic string --besthit_heuristicIndicates how to normalize bit-score in step 1 (default: normal). Normal: no normalization, prt: use proteinortho auxiliary files, smallest: use length of the smallest sequence, target: use target sequence, query: use query sequence, directed: x->y hit, bidirectional: use x->y and y->x hits.
Options: normal, prt, smallest_bidirectional, smallest_directed, query_directed, target_directed, alignment_directed, query_bidirectional, target_bidirectional, alignment_bidirectional
-f float --f_value float
Real number between 0 and 1, a parameter of step 1. Defines the adaptative threshhold as: f\*max_bit_score (default: 0.95).
-bmg_h --bmg_heuristic
Comunity detection method, an heuristic of step 2. (default: Louvain).
Options: Mincut, BPMF, Karger, Greedy, Gradient_Walk, Louvain, Louvain_Obj
-bmgh_nb bool --bmgh_no_binary bool
Flag, specifies if force binary tree in step 2. (no flag: force binary, flag: do not force binary).
-stree_h --species_tree_heuristic
Comunity detection method, an heuristic of step 3. (default: louvain_weight).
Options: naive, louvain, mincut, louvain_weight
-streeh_repeats int --stree_heuristic_repeats int
integer, specifies how many times run the heuristic of step 3. (default: 3)
-streeh_b bool --streeh_binary bool
Flag, specifies if force binary tree in step 3. (no flag: do not force binary, flag: force binary).
-streeh_ndb bool --streeh_no_doble_build bool
Flag, specifies if run build algorithm twice to obtain less resolved tree in step 3. (no flag: double build, flag: single build).
Project details
Release history Release notifications | RSS feed
Download files
Download the file for your platform. If you're not sure which to choose, learn more about installing packages.
Source Distribution
revolutionhtl-1.0.12.tar.gz
(36.7 kB
view hashes)
Built Distribution
Close
Hashes for revolutionhtl-1.0.12-py3-none-any.whl
Algorithm | Hash digest | |
---|---|---|
SHA256 | a61394f898c91c965c17237b55aa6b874974efb7af7d7441e3f99bf4b6b0e87c |
|
MD5 | 8c4b125215c720f30c2ea5629de40eef |
|
BLAKE2b-256 | 30a50fc2999fb3ff1114193958a392aec05c2d037b2835baf25364f3a2e7a317 |