(md5sum: 6ea9b1768bd7085756798d9f96610430 )
Instructions for running seqMMLigner
Synopsis:seqmmligner [sequence-1 fasta file] [sequence-2 fasta] [OPTION]...
|optimal: (default)||finds optimal alignment (and infers associated params) under the information-theoretic measure (i.e., finds the alignment that maximizes the joint probability of the alignment and the two input sequences)|
|marginal:||finds the marginal probability (and infers associated params) that the two input sequences are related, under the information-theoretic measure|
|marginal-interactive:||same as above, but allows interactive exploration of the marginal probability landscape and probe competing sequence alignments|
Score any SEQUENCE_ALIGNMENT_FILE (in fasta format) using I-value. Changes default operation: an alignment is not computed. Instead, SEQUENCE_ALIGNMENT_FILE (FASTA format) is scored using the I-value measure.
Some example command line runs:
Generating an alignment using seqMMLigner:./seqmmligner seq1.fa seq2.fa
./seqmmligner seq1.fa seq2.fa --criterion marginal
./seqmmligner --ivalue alignment.afasta
Instructions for building seqMMLigner (v1.1)
Dependencies: GNUMake or equivalent. A modern C++ compiler. seqMMLigner is known to build with g++ (GCC) >= 4.1.2. If these dependencies are met, follow these instructions:
seqMMLigner is free software: you can redistribute it and/or modify it under the terms of the GNU General Public License as published by the Free Software Foundation, either version 3 of the License. seqMMLigner 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. You should have received a copy of the GNU General Public License along with seqMMLigner. If not, see http://www.gnu.org/licenses/.
- Arun Konagurthu <arun DOT konagurthu AT monash DOT edu>
- Dinithi Sumanaweera <dinithi DOT sumanaweera AT monash DOT edu>
Main and supplementary PDFs
main+sm (combined PDF)
SCOP domain sequence pairs used to infer Dirichlet priors: (click here)
Distribution of #alignments vs sequence-distance parameter (i.e. n of PAM-n): (click here)
Inferred Dirichlet (priors') parameters for sequence-distance parameter n in [1,1000]: (click here)
A selection of marginal probability landscapes: (click here)
Benchmark statistics across the programs seqMMLigner, ClustalW, CONTRAlign, KAlign, MAFFT, MUSCLE, ProbCons, T-Coffee: