mgkit.workflow.hmmer2gff module

Script to convert HMMER results files (domain table) to a GFF file, the name of the profiles are expected to be now in the form GENEID_TAXONID_TAXON-NAME(-nr) by default, but any other profile name is accepted.

The profiles tested are those made from Kegg Orthologs, from the download_profiles script. If the –no-custom-profiles options is used, the script can be used with any profile name. The profile name will be used for gene_id, taxon_id and taxon_name in the GFF file.

It is possible to use seuqnces not translated using mgkit, no information on the frame is assumed, so this script can be used against a protein DB. For example Uniprot can be searched for profiles, in which case the –no-frame options must be used.

Note

for GENEID, old documentation points to KOID, it is the same

Warning

The compatibility with old data has been removed, meaning that old experiments must use the scripts from those versions. It is possible to use multiple environments, with virtualenv for this purpose. An examples is given in Installation.

Changes

Changed in version 0.1.15: adapted to new GFF module and specs

Changed in version 0.2.1: added options to customise output and filters and old restrictions

Changed in version 0.3.1: added –no-frame option for non mgkit-translated proteins, sequence headers are handled the same way as HMMER (truncated at the first space)

mgkit.workflow.hmmer2gff.get_aa_data(f_handle)

Load aminoacid seuqnces used by HMMER.

mgkit.workflow.hmmer2gff.main()

Main loop

mgkit.workflow.hmmer2gff.parse_domain_table_contigs(options)

Parse the HMMER result file

mgkit.workflow.hmmer2gff.set_parser()

Setup command line options