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

hmmalign(1)							   HMMER Manual 						       hmmalign(1)

NAME
hmmalign - align sequences to a profile HMM SYNOPSIS
hmmalign [options] <hmmfile> <seqfile> DESCRIPTION
Perform a multiple sequence alignment of all the sequences in seqfile, by aligning them individually to the profile HMM in hmmfile. The new alignment is output to stdout in Stockholm format. The sequences in seqfile are aligned in unihit local alignment mode. Therefore they should already be known to contain a single domain; they should not contain more than one domain. They may be fragments. The optimal alignment may assign some residues as nonhomologous (N and C states), in which case these residues are still included in the resulting alignment, but shoved to the outer edges. To trim these nonhomologous residues from the result, see the --trim option. OPTIONS
-h Help; print a brief reminder of command line usage and all available options. -o <f> Direct the output alignment to file <f>, rather than to stdout. --allcol Include columns in the output alignment for every match (consensus) state in the hmmfile, even if it means having all-gap columns. This is useful in analysis pipelines that need to be able to maintain a predetermined profile HMM architecture (with an unchanging number of consensus columns) through an hmmalign step. --mapali <f> Merge the existing alignment in file <f> into the result, where <f> is exactly the same alignment that was used to build the model in hmmfile. This is done using a map of alignment columns to consensus profile positions that is stored in the hmmfile. The multi- ple alignment in <f> will be exactly reproduced in its consensus columns (as defined by the profile), but the displayed alignment in insert columns may be altered, because insertions relative to a profile are considered by convention to be unaligned data. --trim Trim nonhomologous residues (assigned to N and C states in the optimal alignments) from the resulting multiple alignment output. --amino Specify that all sequences in seqfile are proteins. By default, alphabet type is autodetected from looking at the residue composi- tion. --dna Specify that all sequences in seqfile are DNAs. --rna Specify that all sequences in seqfile are RNAs. --informat <s> Declare that the input seqfile is in format <s>. Accepted sequence file formats include FASTA, EMBL, Genbank, DDBJ, Uniprot, Stock- holm, and SELEX. Default is to autodetect the format of the file. --outformat <s> Specify that the msafile is in format <s>. Currently the accepted multiple alignment sequence file formats only include Stockholm and SELEX. Default is to autodetect the format of the file. SEE ALSO
See hmmer(1) for a master man page with a list of all the individual man pages for programs in the HMMER package. For complete documentation, see the user guide that came with your HMMER distribution (Userguide.pdf); or see the HMMER web page (@HMMER_URL@). COPYRIGHT
@HMMER_COPYRIGHT@ @HMMER_LICENSE@ For additional information on copyright and licensing, see the file called COPYRIGHT in your HMMER source distribution, or see the HMMER web page (@HMMER_URL@). AUTHOR
Eddy/Rivas Laboratory Janelia Farm Research Campus 19700 Helix Drive Ashburn VA 20147 USA http://eddylab.org HMMER
@HMMER_VERSION@ @HMMER_DATE@ hmmalign(1)
Man Page