Linux and UNIX Man Pages

Linux & Unix Commands - Search Man Pages

glam2format(1) [debian man page]

GLAM2FORMAT(1)							   glam2 Manual 						    GLAM2FORMAT(1)

NAME
glam2format - converts GLAM2 motifs to FASTA or MSF format SYNOPSIS
glam2format [options] my_format my_motif.glam2 Formats: fasta, msf. DESCRIPTION
glam2format reads in a motif found by glam2, and writes it in a standard alignment format (FASTA-with-gaps or MSF). This enables the alignment to be passed to third-party software, including graphical visualization tools such as Kalignvu, Boxshade, and WebLogo. On the other hand, not all the motif information is preserved: in particular, the key positions are lost. Only the top motif in glam2 output is converted. OPTIONS (DEFAULT SETTINGS) -o Output file (stdout). -c Make a compact alignment. By default, residues that are inserted between key positions are written as unaligned with each other. This best reflects glam2's intention, but it can make the alignment large and full of gaps. With -c, inserted residues are written as arbitrarily aligned with each other, just as they appear in the glam2 output. -f Sequence file to make a "global" alignment by adding flanking sequences from the original FASTA-format sequence file. The flanking sequences will be written as either unaligned with each other or arbitrarily aligned, depending on the -c option. The sequences should have unique names and their order should be unchanged. SEE ALSO
boxshade(1), glam2(1), glam2mask(1), glam2-purge(1), glam2scan(1) The full Hypertext documentation of GLAM2 is available online at http://bioinformatics.org.au/glam2/ or on this computer in /usr/share/doc/glam2/. REFERENCE
If you use GLAM2, please cite: MC Frith, NFW Saunders, B Kobe, TL Bailey (2008) Discovering sequence motifs with arbitrary insertions and deletions, PLoS Computational Biology (in press). AUTHORS
Martin Frith Author of GLAM2. Timothy Bailey Author of GLAM2. Charles Plessy <plessy@debian.org> Formatted this manpage in DocBook XML for the Debian distribution. COPYRIGHT
The source code and the documentation of GLAM2 are released in the public domain. GLAM2 1056 05/19/2008 GLAM2FORMAT(1)

Check Out this Related Man Page

HHCONSENSUS(1)							   User Commands						    HHCONSENSUS(1)

NAME
hhconsensus - calculate the consensus sequence for an A3M/FASTA input file SYNOPSIS
hhconsensus -i <file> [options] DESCRIPTION
HHconsensus version 2.0.15 (June 2012) Calculate the consensus sequence for an A3M/FASTA input file. (C) Johannes Soeding, Michael Rem- mert, Andreas Biegert, Andreas Hauser Remmert M, Biegert A, Hauser A, and Soding J. HHblits: Lightning-fast iterative protein sequence searching by HMM-HMM alignment. Nat. Methods 9:173-175 (2011). -i <file> query alignment (A2M, A3M, or FASTA), or query HMM Output options: -s <file> append consensus sequence in FASTA (default=<infile.seq>) -o <file> write alignment with consensus sequence in A3M -oa3m <file> same -oa2m <file> write alignment with consensus sequence in A2M -ofas <file> write alignment with consensus sequence in FASTA -v <int> verbose mode: 0:no screen output 1:only warings 2: verbose Filter input alignment (options can be combined): -id [0,100] maximum pairwise sequence identity (%) (def=100) -diff [0,inf[ filter most diverse set of sequences, keeping at least this many sequences in each block of >50 columns (def=0) -cov [0,100] minimum coverage with query (%) (def=0) -qid [0,100] minimum sequence identity with query (%) (def=0) -qsc [0,100] minimum score per column with query (def=-20.0) Input alignment format: -M a2m use A2M/A3M (default): upper case = Match; lower case = Insert; '-' = Delete; '.' = gaps aligned to inserts (may be omitted) -M first use FASTA: columns with residue in 1st sequence are match states -M [0,100] use FASTA: columns with fewer than X% gaps are match states Other options: -addss add predicted secondary structure information from PSIPRED Example: hhconsensus -i stdin -s stdout hhconsensus 2.0.15 June 2012 HHCONSENSUS(1)
Man Page