infer¶
Infer tree from multiple sequence alignment.
Arguments¶
usage: gtdbtk infer --msa_file MSA_FILE --out_dir OUT_DIR
[--prot_model {JTT,WAG,LG}] [--no_support] [--gamma]
[--prefix PREFIX] [--cpus CPUS] [--tmpdir TMPDIR]
[--debug] [-h]
required named arguments¶
- --msa_file
multiple sequence alignment in FASTA format
- --out_dir
directory to output files
Named Arguments¶
- --prot_model
Possible choices: JTT, WAG, LG
protein substitution model for tree inference
Default: “WAG”
- --no_support
do not compute local support values using the Shimodaira-Hasegawa test
- --gamma
rescale branch lengths to optimize the Gamma20 likelihood
- --prefix
prefix for all output files
Default: “gtdbtk”
- --cpus
number of CPUs to use
Default: 1
- --tmpdir
specify alternative directory for temporary files
Default: “/tmp”
- --debug
create intermediate files for debugging purposes
## Files output
- infer/intermediate_results/
Example¶
Input¶
gtdbtk infer --msa_file msa.faa --out_dir infer_out
Output¶
[2020-04-14 09:37:55] INFO: GTDB-Tk v1.1.0
[2020-04-14 09:37:55] INFO: gtdbtk infer --msa_file msa.faa --out_dir infer_out
[2020-04-14 09:37:55] INFO: Using GTDB-Tk reference data version r89: /release89
[2020-04-14 09:37:55] INFO: Inferring FastTree (WAG, +gamma, support) using a maximum of 1 CPUs.
[2020-04-14 09:37:55] INFO: FastTree version: 2.1.10
[2020-04-14 09:37:55] INFO: FastTree version: 2.1.10
[2020-04-14 09:37:55] INFO: Done.