Linux and UNIX Man Pages

Linux & Unix Commands - Search Man Pages

lastex(1) [debian man page]

LASTEX(1)							   User Commands							 LASTEX(1)

NAME
lastex - Calculate expected numbers of alignments for random sequences. DESCRIPTION
Usage: lastex [options] reference-counts-file query-counts-file [alignments-file] Calculate expected numbers of alignments for random sequences. Options (default settings): -h: show all options and their default settings -s: strands (2 for DNA, 1 for protein) -r: match score (DNA: 1, protein: blosum62) -q: mismatch cost (DNA: 1, protein: blosum62) -p: file for residue pair scores -a: gap existence cost (DNA: 7, pro- tein: 11) -b: gap extension cost (DNA: 1, protein: 2) -g: do calculations for gapless alignments -y: find the expected number of align- ments with score >= this -E: maximum expected number of alignments -z: calculate the expected number of alignments per: 0 = reference counts file / query counts file 1 = reference counts file / each query sequence 2 = each reference sequence / query counts file 3 = each reference sequence / each query sequence (0) REPORTING BUGS
Report bugs to: last (ATmark) cbrc (dot) jp LAST home page: http://last.cbrc.jp/ lastex 199 May 2012 LASTEX(1)

Check Out this Related Man Page

clustalo(1)							   USER COMMANDS						       clustalo(1)

NAME
clustalo - General purpose multiple sequence alignment program for proteins SYNOPSIS
clustalo [-h] DESCRIPTION
Clustal-Omega is a general purpose multiple sequence alignment (MSA) program for proteins. It produces high quality MSAs and is capable of handling data-sets of hundreds of thousands of sequences in reasonable time. In default mode, users give a file of sequences to be aligned and these are clustered to produce a guide tree and this is used to guide a "progressive alignment" of the sequences. There are also facilities for aligning existing alignments to each other, aligning a sequence to an alignment and for using a hidden Markov model (HMM) to help guide an alignment of new sequences that are homologous to the sequences used to make the HMM. This latter procedure is referred to as "external profile alignment" or EPA. Clustal-Omega uses HMMs for the alignment engine, based on the HHalign package from Johannes Soeding [1]. Guide trees are made using an enhanced version of mBed [2] which can cluster very large numbers of sequences in O(N*log(N)) time. Multiple alignment then proceeds by aligning larger and larger alignments using HHalign, following the clustering given by the guide tree. In its current form Clustal-Omega can only align protein sequences but not DNA/RNA sequences. It is envisioned that DNA/RNA will become available in a future version. USAGE
Tool usage is available in /usr/share/doc/clustalo/README. DEVELOPMENT
Headers and libraries are available in libclustalo-dev package. CITING
Sievers F, Wilm A, Dineen DG, Gibson TJ, Karplus K, Li W, Lopez R, McWilliam H, Remmert M, Soding J, Thompson JD, Higgins DG (2011). Fast, scalable generation of high-quality protein multiple sequence alignments using Clustal Omega. Mol Syst Biol 7. AUTHOR
Olivier Sallou (olivier.sallou (at) irisa.fr) - Man page and packaging Conway Institute UCD Dublin (clustalw (at) ucd.ie) - clustalo version 1.0.3 December 14, 2011 clustalo(1)
Man Page

3 More Discussions You Might Find Interesting

1. Shell Programming and Scripting

is not an identifier error

Hi gurus, I am trying to execute a shell script which connects to the oracle db and get the values from a table and then it exports the values on UNIX. For this, I am creating a temp file which stores the values returned from the select query and then executes this file. Below is the code which... (13 Replies)
Discussion started by: svajhala
13 Replies

2. Shell Programming and Scripting

Finding a String and Replacing it with an incremental value

Hi Friends, I have a text file which has about 200,000 records in it. we have a string which repeats in each and every record. So we have to write a script in ksh which finds that string on each line and replaces it with a new string(incremental value) for a set every four records. To be... (12 Replies)
Discussion started by: vpv0002
12 Replies

3. UNIX for Advanced & Expert Users

Some mails not scanned by SpamAssassin 3.4.1 w sendmail/procmail on Fedora 21

We have a user who has about a 50% missed rate on spam detection. I'm wondering if his user prefs or something is preventing scanning of all messages? SpamAssassin version 3.4.1, running on Perl version 5.20.3, sendmail Version 8.15.2 The contents of the user_prefs file: # How many points... (2 Replies)
Discussion started by: RobbieTheK
2 Replies