identify

Identify marker genes in genome(s). The following heuristic is used to establish the translation table used by Prodigal: use table 11 unless the coding density using table 4 is 5% higher than when using table 11 and the coding density under table 4 is >70%. Distinguishing between tables 4 and 25 is challenging so GTDB-Tk does not attempt to distinguish between these two tables. If you know the correct translation table for your genomes this can be provided to GTDB-Tk in the –batchfile.

Arguments

## Files output

• [prefix].log

• [prefix].json

• [prefix].warnings.log

• identify

  • [prefix].[domain].markers_summary.tsv

  • [prefix].translation_table_summary.tsv

  • [prefix].failed_genomes.tsv

  • intermediate_results/marker_genes/[genome_id]/

  • [genome_id]_pfam_tophit.tsv

  • [genome_id]_pfam.tsv

  • [genome_id]_protein.faa

  • [genome_id]_protein.fna

  • [genome_id]_protein.gff

  • [genome_id]_tigrfam.out

  • [genome_id]_tigrfam_tophit.tsv

  • [genome_id]_tigrfam.tsv

  • prodigal_translation_table.tsv

Example

Input

gtdbtk identify --genome_dir genomes/ --out_dir identify_output --cpus 3

Output

[2022-04-11 11:48:59] INFO: GTDB-Tk v2.0.0
[2022-04-11 11:48:59] INFO: gtdbtk identify --genome_dir /tmp/gtdbtk/genomes --out_dir /tmp/gtdbtk/identify --extension gz --cpus 2
[2022-04-11 11:48:59] INFO: Using GTDB-Tk reference data version r207: /srv/db/gtdbtk/official/release207
[2022-04-11 11:48:59] INFO: Identifying markers in 2 genomes with 2 threads.
[2022-04-11 11:48:59] TASK: Running Prodigal V2.6.3 to identify genes.
[2022-04-11 11:49:10] INFO: Completed 2 genomes in 10.94 seconds (5.47 seconds/genome).
[2022-04-11 11:49:10] TASK: Identifying TIGRFAM protein families.
[2022-04-11 11:49:16] INFO: Completed 2 genomes in 5.78 seconds (2.89 seconds/genome).
[2022-04-11 11:49:16] TASK: Identifying Pfam protein families.
[2022-04-11 11:49:16] INFO: Completed 2 genomes in 0.42 seconds (4.81 genomes/second).
[2022-04-11 11:49:16] INFO: Annotations done using HMMER 3.1b2 (February 2015).
[2022-04-11 11:49:16] TASK: Summarising identified marker genes.
[2022-04-11 11:49:16] INFO: Completed 2 genomes in 0.05 seconds (40.91 genomes/second).
[2022-04-11 11:49:16] INFO: Done.