View source for AMPHORA

= AMPHORA2 =
'''AMPHORA has been updated. See [http://wolbachia.biology.virginia.edu/WuLab/Software.html] for AMPHORA2.'''

==AMPHORA==

AMPHORA is an Automated Phylogenomic Inference Pipeline for bacterial sequences. From a given a set of protein sequences, it automatically identifies 31 phylogenetic marker genes. It then generates high-quality multiple sequence alignments 
for these genes and make tree-based phylotype assignments.

==CITATION==

Please cite AMPHORA as: Martin Wu and Jonathan A Eisen. A simple, fast, and accurate method of phylogenomic inference Genome Biology 2008, 9:R151

==SYSTEM REQUIREMENT==

Linux OS (kernel version 2.6 or later)

The following software is required by the AMPHORA package. They need to be downloaded and installed separately from AMPHORA. 
# Perl 5.8.8 or later (http://www.perl.org) 
# Bioperl core package 1.5.2 or later  (http://www.bioperl.org) 
# HMMER (http://hmmer.janelia.org)
# WU BLAST (http://blast.wustl.edu)

The following software is included in the AMPHORA distribution. Their source codes have been slightly modified to suit the needs of AMPHORA. 
# seqboot (http://evolution.genetics.washington.edu/phylip) 
# quicktree (http://www.sanger.ac.uk/Software/analysis/quicktree) 
# raxml (http://icwww.epfl.ch/~stamatak/index-Dateien/Page443.htm)

==LICENSE==

AMPHORA Copyright 2008 by Martin Wu

AMPHORA is free software: you may redistribute it and/or modify its under the terms of the GNU General Public License as published by the Free Software Foundation; either version 2 of the License, or any later version.

AMPHORA is distributed in the hope that it will be useful, but WITHOUT ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the GNU General Public License for more details 
(http://www.gnu.org/licenses/).

==INSTALLATION==

# Download the AMPHORA package AMPHORA.tar.gz from http://bobcat.genomecenter.ucdavis.edu/AMPHORA/download.html
# Unpack the package <pre>tar -xzvf AMPHORA.tar.gz</pre>
# Install AMPHORA to a user specified directory. Make sure you have write permission to the directory.<pre>cd AMPHORA </pre><pre>perl INSTALL.pl -AMPHORA_home user-specified-directory -Bioperl_home path-to-bioperl</pre>

==PACKAGE CONTENTS==

After successful installation, there should be several folders in the home directory of AMPHORA
# '''Marker''' <br>It contains curated seed multiple sequence alignments of the phylogenetic markers in GDE format (with embedded masks) and associated Hidden Markov Models. Currently there are 31 protein marker genes. They are dnaG, frr, infC, nusA, pgk, pyrG, rplA, rplB, rplC, rplD, rplE, rplF, rplK, rplL, rplM, rplN, rplP, rplS, rplT, rpmA, rpoB, rpsB, rpsC, rpsE, rpsI, rpsJ, rpsK, rpsM, rpsS, smpB, tsf.
# '''Reference''' <br>It contains protein sequences of the marker genes from all complete bacterial genomes.  It also contains a bacterial genome tree that was made from the concatenated protein sequences of all the marker genes. Reference trees for each marker gene were derived from the genome tree by replacing the species names with their corresponding gene names and keeping the topology intact.
# '''Taxonomy''' <br>The NCBI taxonomy database (ftp://ftp.ncbi.nih.gov/pub/taxonomy/) with minor modifications. The changes are listed in the file change.note
# '''Scripts''' <br>Perl scripts for identifying markers, generating trimmed multiple sequence alignments, and assigning phylotypes based on phylogenetic inferences.
# '''bin''' <br>Helper programs

==USING AMPHORA==

===Phylotying bacterial sequences===
*Usage: 
  <code>AMPHORA_home/Scripts/Phylotyping.pl <options> protein-sequence-file output-file</code>
*Options:
  -Replicates: number of bootstrap replicates
  -BootstrapCutoff: normalized to 100 replicates (1-100)%  default 70 
*Output:
  foo.pep (identified maker sequences in fasta format, i.e. rpoB.pep)
  foo.aln (aligned and trimmed marker sequences, i.e., rpoB.aln)
  Trees/foo.tree (generated phylogenetic trees)
===Identify marker sequences===
*Usage:
  perl AMPHORA_home/Scripts/MarkerScanner.pl protein-sequence-file
*Output: 
  foo.pep (identified maker sequences in fasta format, i.e., rpoB.pep)
===Align and trim the marker sequences===
*Usage: 
  perl AMPHORA_home/Scripts/MarkerAlignTrim.pl
*Options: 
  -Partial:  query sequences contain partial genes
  -Trim:  trim the alignment using masks embedded within the marker database
  -Strict: use a conservative mask
  -Directory: the file directory where marker sequences are located. Default: current directory
  -Help:  print the help message
*Output:
   foo.aln (aligned and trimmed marker sequences, i.e., rpoB.aln)