The script contains commands used to access functionality related to
taxonomy, without the need to write ad-hoc code for functionality that
can be part of a workflow. One example is access to the the last common
ancestor function contained in the
Last Common Ancestor (lca and lca_line)¶
These commands expose the functionality of
last_common_ancestor_multiple(), making it accessible via the command
line. They differ in the input file format and the choice of output files.
the lca command can be used to define the last common ancestor of contigs from the annotation in a GFF file. The command uses the taxon_ids from all annotations belonging to a contig/sequence, if they have a bitscore higher or equal to the one passed (50 by default). The default output of the command is a tab separated file where the first column is the contig/sequence name, the taxon_id of the last common ancestor, its scientific/common name and its lineage.
contig_21 172788 uncultured phototrophic eukaryote cellular organisms,environmental samples
If the -r is used, by passing the fasta file containing the nucleotide sequences the output file is a GFF where for each an annotation for the full contig length contains the same information of the tab separated file format.
The lca_line command accept as input a file where each line consist of a list of taxon_ids. The separator for the list can be changed and it defaults to TAB. The last common ancestor for all taxa on a line is searched. The ouput of this command is the same as the tab separated file of the lca command, with the difference that instead of the first column, which in this command becames a list of all taxon_ids that were used to find the last common ancestor for that line. The list of taxon_ids is separated by semicolon “;”.
Both also accept the -n option, to report the config/line and the taxon_ids that had no common ancestors. These are treated as errors and do not appear in the output file.
New in version 0.3.0.
The lca command supports the writing of a file compatible with Krona. The output file can be used with the ktImportText/ImportText.pl script included with KronaTools. Specifically, the output from taxon_utils will be a file with all the lineages found (tab separated), that can be used with:
$ ktImportText -q taxon_utils_ouput
Note the use of -q to make the script count the lineages. Sequences with no LCA found will be marked as No LCA in the graph, the -n is not required.
Please note that the output won’t include any sequence that didn’t have a hit with the software used. If that’s important, the -kt option can be used to add a number of Unknown lines at the end, to read the total supplied.
Filter by Taxon¶
The filter command of this script allows to filter a GFF file using the taxon_id attribute to include only some annotations, or exclude some. The filter is based on the mgkit.taxon.is_ancestor function, and the mgkit.filter.taxon.filter_taxon_by_id_list. It can also filter a table (tab separated values) when the first element is an ID and the second is a taxon_id. An example of a table of this sort is the output of the download-ncbi-taxa.sh and download-uniprot-taxa.sh, where each accession of a database is associated to a taxon_id.
Multiple taxon_id can be passed, either for inclusion or exclusion. If both exclusion and inclusion is used, the first check is on the inclusion and then on the exclusion. In alternative to passing taxon_id, taxon_names can be passed, with values such as ‘cellular organisms’ that needs to be quoted. Example:
$ taxon-utils filter -i 2 -in archaea -en prevotella -t taxonomy.pickle in.gff out.gff
Which will keep only line that are from Bacteria (taxon_id=2) and exclude those from the genus Prevotella. It will be also include Archaea.
Annotations with no taxon_id are not included in the output of both filters
Convert Taxa Tables to HDF5¶
This command is used to convert the taxa tables download from Uniprot and NCBI, using the scripts mentioned in download-data, download-uniprot-taxa.sh and download-ncbi-taxa into a HDF5 file that can be used with the addtaxa command in add-gff-info - Add informations to GFF annotations.
The advantage is a faster lookup of the IDs. The other is a smaller memory footprint when a great number of annotations are kept in memory.
Changed in version 0.3.4: changed interface and behaviour for filter, also now can filter tables; lca has changed the interface and allows the output of a 2 column table
Changed in version 0.3.1: added to_hdf command
Changed in version 0.3.1: added -j option to lca, which outputs a JSON file with the LCA results
Changed in version 0.3.0: added -k and -kt options for Krona output, lineage now includes the LCA also added -a option to select between lineages with only ranked taxa. Now it defaults to all components.
Changed in version 0.2.6: added feat-type option to lca command, added phylum output to nolca
New in version 0.2.5.
get_taxon_info(taxonomy, taxon_id, only_ranked)¶
write_json(lca_dict, seq_id, taxonomy, taxon_id, only_ranked)¶
write_krona(file_handle, taxonomy, taxon_id, only_ranked)¶
write_lca_gff(file_handle, seq_id, seq, taxon_id, taxon_name, lineage, feat_type)¶
write_lca_tab(file_handle, seq_id, taxon_id, taxon_name, rank, lineage)¶
write_lca_tab_simple(file_handle, seq_id, taxon_id)¶
write_no_lca(file_handle, seq_id, taxon_ids, extra=None)¶