NAME

       homGeneMapping - create summary of gene homology

SYNOPSIS

       homGeneMapping [options] --gtfs=gffilenames.tbl --halfile=aln.hal

DESCRIPTION

       homGeneMapping takes a set of gene predictions of different genomes and a hal alignment of
       the genomes and prints a summary for each gene, e.g.

       •   how many of its exons/introns are in agreement with genes of other genomes

       •   how many of its exons/introns are supported by extrinsic evidence from any of the
           genomes

       •   a list of geneids of homologous genes

OPTIONS

   Mandatory parameters
       --halfile=aln.hal
           input hal file

       --gtfs=gtffilenames.tbl
           a text file containing the locations of the input gene files and optionally the hints
           files (both in GTF format). The file is formatted as follows:

               name_of_genome_1  path/to/genefile/of/genome_1  path/to/hintsfile/of/genome_1
               name_of_genome_2  path/to/genefile/of/genome_2  path/to/hintsfile/of/genome_2
               ...
               name_of_genome_N  path/to/genefile/of/genome_N  path/to/hintsfile/of/genome_N

   Additional options
       --cpus=N
           N is the number of CPUs to use (default: 1)

       --noDupes
           do not map between duplications in hal graph (default: off)

       --details
           print detailed output (default: off)

       --halLiftover_exec_dir=DIR
           Directory that contains the executable halLiftover If not specified it must be in
           $PATH environment variable.

       --unmapped
           print a GTF attribute with a list of all genomes, that are not aligned to the
           corresponding gene feature, e.g. hgm_unmapped \"1,4,5\"; (default; off)

       --tmpdir=DIR
           a temporary file directory that stores lifted over files (default 'tmp/' in current
           directory)

       --outdir=DIR
           file directory that stores output gene files (default: current directory)

       --printHomologs=FILE
           prints disjunct sets of homologous transcripts to FILE, e.g.

               # 0     dana
               # 1     dere
               # 2     dgri
               # 3     dmel
               # 4     dmoj
               # 5     dper
               (0, jg4139.t1) (0, jg4140.t1) (1, jg7797.t1) (2, jg3247.t1) (4, jg6720.t1) (5, jg313.t1)
               (1, jg14269.t1) (3, jg89.t1) (5, jg290.t1)
               ...
               Two transcripts are in the same set, if all their exons/introns are homologs and their are
               no additional exons/introns.
               This option requires the Boost C++ Library

       --dbaccess=db
           retrieve hints from an SQLite database. In order to set up a database and populate it
           with hints a separate tool 'load2sqlitedb' is provided. For more information, see the
           documentation in docs/RUNNING-AUGUSTUS-IN-CGP-MODE.md (DATABASE ACCESS / SQLite) in
           the Augustus package. If both a database and hint files in 'gtffilenames.tbl' are
           specified, hints are retrieved from both sources.

EXAMPLE

       homGeneMapping --noDupes
       --halLiftover_exec_dir=~/tools/progressiveCactus/submodules/hal/bin
       --gtfs=gtffilenames.tbl --halfile=msca.hal homGeneMapping --gtfs=gtffilenames.tbl
       --halfile=aln.hal --outdir=outdir --cpus=4

AUTHORS

       AUGUSTUS was written by M. Stanke, O. Keller, S. König, L. Gerischer and L. Romoth.

ADDITIONAL DOCUMENTATION

       An exhaustive documentation can be found in the file /usr/share/doc/augustus/README.md.gz.

                                                                                HOMGENEMAPPING(1)