<!doctype refentry PUBLIC "-//OASIS//DTD DocBook V4.1//EN" [
<!-- Process this file with docbook-to-man to generate an nroff manual
page: `docbook-to-man manpage.sgml > manpage.1'. You may view
the manual page with: `docbook-to-man manpage.sgml | nroff -man |
less'. A typical entry in a Makefile or Makefile.am is:
manpage.1: manpage.sgml
docbook-to-man $< > $@
The docbook-to-man binary is found in the docbook-to-man package.
Please remember that if you create the nroff version in one of the
debian/rules file targets (such as build), you will need to include
docbook-to-man in your Build-Depends control field.
-->
<!-- Fill in your name for FIRSTNAME and SURNAME. -->
<!ENTITY dhfirstname "<firstname>Steffen</firstname>">
<!ENTITY dhsurname "<surname>Möller</surname>">
<!-- Please adjust the date whenever revising the manpage. -->
<!ENTITY dhdate "<date>November 10, 2004</date>">
<!ENTITY dhsection "<manvolnum>1</manvolnum>">
<!ENTITY dhemail "<email>moeller@pzr.uni-rostock.de</email>">
<!ENTITY dhusername "Steffen Moeller">
<!ENTITY dhucpackage "<refentrytitle>TIGR-GLIMMER</refentrytitle>">
<!ENTITY dhpackage "tigr-glimmer">
<!ENTITY debian "<productname>Debian</productname>">
<!ENTITY gnu "<acronym>GNU</acronym>">
<!ENTITY gpl "&gnu; <acronym>GPL</acronym>">
]>
<refentry>
<refentryinfo>
<address>
&dhemail;
</address>
<author>
&dhfirstname;
&dhsurname;
</author>
<copyright>
<year>2003</year>
<holder>&dhusername;</holder>
</copyright>
&dhdate;
</refentryinfo>
<refmeta>
&dhucpackage;
&dhsection;
</refmeta>
<refnamediv>
<refname>&dhpackage;</refname>
<refpurpose>
Find/Score potential genes in genome-file using the probability model in icm-file
</refpurpose>
</refnamediv>
<refsynopsisdiv>
<cmdsynopsis>
<command>tigr-glimmer3</command>
<arg><option><replaceable>genome-file</replaceable></option></arg>
<arg><option><replaceable>icm-file</replaceable></option></arg>
<arg><option><replaceable>[options]</replaceable></option></arg>
</cmdsynopsis>
</refsynopsisdiv>
<refsect1>
<title>DESCRIPTION</title>
<para>
<command>&dhpackage;</command> is a system for finding genes in microbial DNA, especially the genomes of bacteria and archaea. <command>&dhpackage;</command> (Gene Locator and Interpolated Markov Modeler) uses interpolated Markov models (IMMs) to identify the coding regions and distinguish them from noncoding DNA. The IMM approach, described in our Nucleic Acids Research paper on <command>&dhpackage;</command> 1.0 and in our subsequent paper on <command>&dhpackage;</command> 2.0, uses a combination of Markov models from 1st through 8th-order, weighting each model according to its predictive power. <command>&dhpackage;</command> 1.0 and 2.0 use 3-periodic nonhomogenous Markov models in their IMMs.
</para><para>
<command>&dhpackage;</command> is the primary microbial gene finder at TIGR, and has been used to annotate the complete genomes of B. burgdorferi (Fraser et al., Nature, Dec. 1997), T. pallidum (Fraser et al., Science, July 1998), T. maritima, D. radiodurans, M. tuberculosis, and non-TIGR projects including C. trachomatis, C. pneumoniae, and others. Its analyses of some of these genomes and others is available at the TIGR microbial database site.
</para><para>
A special version of <command>&dhpackage;</command> designed for small eukaryotes, GlimmerM, was used to find the genes in chromosome 2 of the malaria parasite, P. falciparum.. GlimmerM is described in S.L. Salzberg, M. Pertea, A.L. Delcher, M.J. Gardner, and H. Tettelin, "Interpolated Markov models for eukaryotic gene finding," Genomics 59 (1999), 24-31. Click here (http://www.tigr.org/software/glimmerm/) to visit the GlimmerM site, which includes information on how to download the GlimmerM system.
</para><para>
The <command>&dhpackage;</command> system consists of two main programs. The first of these is the training program, build-imm. This program takes an input set of sequences and builds and outputs the IMM for them. These sequences can be complete genes or just partial orfs. For a new genome, this training data can consist of those genes with strong database hits as well as very long open reading frames that are statistically almost certain to be genes. The second program is glimmer, which uses this IMM to identify putative genes in an entire genome. <command>&dhpackage;</command> automatically resolves conflicts between most overlapping genes by choosing one of them. It also identifies genes that are suspected to truly overlap, and flags these for closer inspection by the user. These ``suspect'' gene candidates have been a very small percentage of the total for all the genomes analyzed thus far.
<command>&dhpackage;</command> is a program that...</para>
</refsect1>
<refsect1>
<title>OPTIONS</title>
<variablelist>
<varlistentry>
<term><option>-C <replaceable>n</replaceable></option></term>
<listitem>
<para>Use n as GC percentage of independent model</para>
<para>Note: n should be a percentage, e.g., -C 45.2</para>
</listitem>
</varlistentry>
<varlistentry>
<term>-f</term><listitem><para>Use ribosome-binding energy to choose start codon</para></listitem>
</varlistentry>
<varlistentry>
<term><option>+f</option></term><listitem><para>Use first codon in orf as start codon</para></listitem>
</varlistentry>
<varlistentry>
<term><option>-g <replaceable>n</replaceable></option></term><listitem><para>Set minimum gene length to n</para></listitem>
</varlistentry>
<varlistentry>
<term><option>-i <replaceable>filename</replaceable></option></term>
<listitem>
<para>Use <option><replaceable>filename</replaceable></option>
to select regions of bases that are off
limits, so that no bases within that area will be examined
</para>
</listitem>
</varlistentry>
<varlistentry>
<term><option>-l</option></term>
<listitem><para>Assume linear rather than circular genome, i.e., no wraparound</para></listitem>
</varlistentry>
<varlistentry>
<term><option>-L <replaceable>filename</replaceable></option></term>
<listitem><para>Use filename to specify a list of orfs that should
be scored separately, with no overlap rules
</para></listitem>
</varlistentry>
<varlistentry>
<term><option>-M</option></term>
<listitem><para>Input is a multifasta file of separate genes to be scored
separately, with no overlap rules
</para>
</listitem>
</varlistentry>
<varlistentry>
<term><option>-o <replaceable>n</replaceable></option></term>
<listitem>
<para>Set minimum overlap length to n. Overlaps shorter than this
are ignored.
</para></listitem>
</varlistentry>
<varlistentry>
<term><option>-p <replaceable>n</replaceable></option></term>
<listitem>
<para>
Set minimum overlap percentage to n%. Overlaps shorter than this percentage of *both* strings are ignored.
</para>
</listitem>
</varlistentry>
<varlistentry>
<term><option>-q <replaceable>n</replaceable></option></term>
<listitem>
<para>Set the maximum length orf that can be rejected because of
the independent probability score column to (n - 1)
</para>
</listitem>
</varlistentry>
<varlistentry>
<term><option>-r</option></term>
<listitem>
<para>
Don't use independent probability score column
</para>
</listitem>
</varlistentry>
<varlistentry>
<term><option>+r</option></term>
<listitem><para>
Use independent probability score column
</para>
</listitem>
</varlistentry>
<varlistentry>
<term><option>-r</option></term>
<listitem>
<para>
Don't use independent probability score column
</para> </listitem> </varlistentry> <varlistentry>
<term><option>-s <replaceable>s</replaceable></option></term>
<listitem><para> Use string s as the ribosome binding pattern to find start codons.</para>
</listitem>
</varlistentry>
<varlistentry>
<term><option>+S</option></term>
<listitem>
<para>
Do use stricter independent intergenic model that doesn't
give probabilities to in-frame stop codons. (Option is obsolete
since this is now the only behaviour
</para> </listitem>
</varlistentry>
<varlistentry>
<term><option>-t <replaceable>n</replaceable></option></term>
<listitem><para>
Set threshold score for calling as gene to n. If the in-frame
score >= n, then the region is given a number and considered
a potential gene.
</para> </listitem>
</varlistentry>
<varlistentry>
<term><option>-w <replaceable>n</replaceable> </option></term>
<listitem><para>
Use "weak" scores on tentative genes n or longer. Weak
scores ignore the independent probability score.
</para></listitem>
</varlistentry>
</variablelist>
</refsect1>
<refsect1>
<title>SEE ALSO</title>
<para>
tigr-adjust (1),
tigr-anomaly (1),
tigr-build-icm (1),
tigr-check (1),
tigr-codon-usage (1),
tigr-compare-lists (1),
tigr-extract (1),
tigr-generate (1),
tigr-get-len (1),
tigr-get-putative (1),
tigr-glimmer3 (1),
tigr-long-orfs (1)
</para>
<para>
http://www.tigr.org/software/glimmer/
</para>
<para>Please see the readme in /usr/share/doc/glimmer for a description on how to use Glimmer.</para>
</refsect1>
<refsect1>
<title>AUTHOR</title>
<para>This manual page was quickly copied from the glimmer web site by &dhusername; &dhemail; for
the &debian; system.
</para>
</refsect1>
</refentry>
<!-- Keep this comment at the end of the file
Local variables:
mode: sgml
sgml-omittag:t
sgml-shorttag:t
sgml-minimize-attributes:nil
sgml-always-quote-attributes:t
sgml-indent-step:2
sgml-indent-data:t
sgml-parent-document:nil
sgml-default-dtd-file:nil
sgml-exposed-tags:nil
sgml-local-catalogs:nil
sgml-local-ecat-files:nil
End:
-->