lalign

Name lalign
Description

lalign is part of the Fasta3 package. This version replaces that from the Fasta2 package.

While programs such as fasta and ssearch report only the best alignment between the query sequence and the library sequence, lalign reports a number of non-overlapping alignments between sequences.

Graphical versions of lalign output can be generated using the programs lav2ps or lav2svg. To do this, use the -m 11 flag to generate lav format output, and pipe the lalign results directly to the appropriate lav2X program. E.g. to generate postscript output, try the command:

lalign -m11 inputSeq1.fasta inputSeq2.fasta | lav2ps > lalignResult.ps

lalign the Huang and Miller SIM algorithm (X. Huang and W. Miller (1991) Adv. Appl. Math. 12:337-357), This is a linear-space version of an algorithm described by M. S. Waterman and M. Eggert (J. Mol. Biol. 197:723-728).


In the Bio-Linux package, the threaded versions of the fasta programs are the default.


The programs available in the Fasta3 package are:

  • fasta - scan a protein or DNA sequence library for similar sequences.
  • fastx - compare a DNA sequence to a protein sequence database, comparing the translated DNA sequence in forward and reverse frames.
  • tfastx - compare a protein sequence to a DNA sequence database, calculating similarities with frameshifts to the forward and reverse orientations.
  • fasty - compare a DNA sequence to a protein sequence database, comparing the translated DNA sequence in forward and reverse frames.
  • tfasty - compare a protein sequence to a DNA sequence database, calculating similarities with frameshifts to the forward and reverse orientations.
  • fasts - compare unordered peptides to a protein sequence database
  • tfasts - compare unordered peptides to a translated DNA sequence database
  • fastm - compare ordered peptides (or short DNA sequences) to a protein (DNA) sequence database
  • fastm - compare ordered peptides (or short DNA sequences) to a translated DNA sequence database
  • fastf - compare mixed peptides to a protein sequence database
  • tfastf - compare mixed peptides to a translated DNA sequence database
  • ssearch - compare a protein or DNA sequence to a sequence database using the Smith-Waterman algorithm.
  • ggsearch - compare a protein or DNA sequence to a sequence database using a global alignment (Needleman-Wunsch)
  • lalign - produce multiple non-overlapping alignments for protein and DNA sequences using the Huang and Miller SIM algorithm for the Waterman-Eggert algorithm. This version of lalign replaces that from the Fasta2 package.
  • prss - (discontinued, replaced in the fasta35 release by new versions of ssearch and fastx) estimate statistical significance of an alignment by comparing the score to the distribution of similarity scores generated by shuffling the second sequence. prss35 uses Smith-Waterman. prfx35 uses the fastx algorithm.

References:
Pearson, W.R. Flexible sequence similarity searching with the FASTA3 program package. Methods Mol Biol. 2000;132:185-219 [Entrez]

Pearson, W.R. Empirical statistical estimates for sequence similarity searches. J Mol Biol. 1998 Feb 13;276(1):71-84 [Entrez]

Pearson WR, Wood T, Zhang Z, Miller W. Comparison of DNA sequences with protein sequences. Genomics. 1997 Nov 15;46(1):24-36. [Entrez]


Homepage http://www.people.virginia.edu/~wrp/pearson.html  
Remote Documentation http://www.people.virginia.edu/~wrp/papers/ismb2000.pdf