Fast and sensitive protein alignment using DIAMOND

Benjamin Buchfink,Chao Xie,Daniel H Huson

DOI: https://doi.org/10.1038/nmeth.3176

IF: 48

2014-11-17

Nature Methods

Abstract:The open-source DIAMOND software provides protein alignment that is 20,000 times faster on short reads than BLASTX at similar sensitivity, for rapid analysis of large metagenomics data sets on a desktop computer.

biochemical research methods

What problem does this paper attempt to address?

Sensitive protein alignments at tree-of-life scale using DIAMOND

Benjamin Buchfink,Klaus Reuter,Hajk-Georg Drost

DOI: https://doi.org/10.1038/s41592-021-01101-x

IF: 48

2021-04-01

Nature Methods

Abstract:Abstract We are at the beginning of a genomic revolution in which all known species are planned to be sequenced. Accessing such data for comparative analyses is crucial in this new age of data-driven biology. Here, we introduce an improved version of DIAMOND that greatly exceeds previous search performances and harnesses supercomputing to perform tree-of-life scale protein alignments in hours, while matching the sensitivity of the gold standard BLASTP.

biochemical research methods
ProteinInferencer: Confident protein identification and multiple experiment comparison for large scale proteomics projects

Yaoyang Zhang,Tao Xu,Bing Shan,Jonathan Hart,Aaron Aslanian,Xuemei Han,Nobel Zong,Haomin Li,Howard Choi,Dong Wang,Lipi Acharya,Lisa Du,Peter K. Vogt,Peipei Ping,John R. Yates

DOI: https://doi.org/10.1016/j.jprot.2015.07.006

IF: 3.855

2015-01-01

Journal of Proteomics

Abstract:Shotgun proteomics generates valuable information from large-scale and target protein characterizations, including protein expression, protein quantification, protein post-translational modifications (PTMs), protein localization, and protein–protein interactions. Typically, peptides derived from proteolytic digestion, rather than intact proteins, are analyzed by mass spectrometers because peptides are more readily separated, ionized and fragmented. The amino acid sequences of peptides can be interpreted by matching the observed tandem mass spectra to theoretical spectra derived from a protein sequence database. Identified peptides serve as surrogates for their proteins and are often used to establish what proteins were present in the original mixture and to quantify protein abundance. Two major issues exist for assigning peptides to their originating protein. The first issue is maintaining a desired false discovery rate (FDR) when comparing or combining multiple large datasets generated by shotgun analysis and the second issue is properly assigning peptides to proteins when homologous proteins are present in the database. Herein we demonstrate a new computational tool, ProteinInferencer, which can be used for protein inference with both small- or large-scale data sets to produce a well-controlled protein FDR. In addition, ProteinInferencer introduces confidence scoring for individual proteins, which makes protein identifications evaluable.
Fast and accurate protein structure search with Foldseek

Michel van Kempen,Stephanie S Kim,Charlotte Tumescheit,Milot Mirdita,Jeongjae Lee,Cameron L M Gilchrist,Johannes Söding,Martin Steinegger

DOI: https://doi.org/10.1038/s41587-023-01773-0

Abstract:As structure prediction methods are generating millions of publicly available protein structures, searching these databases is becoming a bottleneck. Foldseek aligns the structure of a query protein against a database by describing tertiary amino acid interactions within proteins as sequences over a structural alphabet. Foldseek decreases computation times by four to five orders of magnitude with 86%, 88% and 133% of the sensitivities of Dali, TM-align and CE, respectively.
Rapid and Sensitive Protein Complex Alignment with Foldseek-Multimer

Woosub Kim,Milot Mirdita,Eli Levy Karin,Cameron L.M. Gilchrist,Hugo Schweke,Johannes Soeding,Emmanuel D. Levy,Martin Steinegger

DOI: https://doi.org/10.1101/2024.04.14.589414

2024-10-28

Abstract:Advances in computational structure prediction will vastly augment the hundreds of thousands of currently-available protein complex structures. Translating these into discoveries requires aligning them, which is computationally prohibitive. Foldseek-Multimer computes complex alignments from compatible chain-to-chain alignments, identified by efficiently clustering their superposition vectors. Foldseek-Multimer is 3-4 orders of magnitudes faster than the gold standard, while producing comparable alignments; allowing it to compare billions of complex-pairs in 11 hours. Foldseek-Multimer is open-source software: https://github.com/steineggerlab/foldseek, webserver: https://search.foldseek.com and the BFMD database.

Bioinformatics
Fast-Part: Fast and Accurate Data Partitioning for Biological Sequence Analysis

Shafayat Ahmed,Muhit Islam Emon,Nazifa Ahmed Moumi,LIQING ZHANG

DOI: https://doi.org/10.1101/2024.11.13.623463

2024-11-15

Abstract:Developing effective machine learning models for classifications of biological sequences depends heavily on the quality of the training and test datasets split. Existing tools are either computationally expensive, unable to maintain the desired level of similarity between the training and test datasets, or unable to retain training: test ratio stratification. Here, we present Fast-Part, a fast and accurate sequence data partitioning tool that ensures strict homology separation between the training and test datasets and the best possible training: test stratification ratio, and at the same time is computationally fast. Evaluation of Fast-Part on multiple protein sequence datasets shows that it performs data partitioning with exceptional speed and maintains strict partitioning compared to the existing tools. Fast-Part can handle massive datasets like CD-HIT[1] and MMseq[2] and maintain strict homology partitioning like GraphPart[3].

Bioinformatics
Lapidary: Identifying and reporting amino acid sequences in metagenomes using sequence reads and Diamond

Samuel J Bloomfield,Aldert L Zomer,Alison E Mather

DOI: https://doi.org/10.1101/2024.03.25.586564

2024-03-28

Abstract:Genome and metagenome comparisons rely on identifying genetic elements that differ or are in common between samples. These genetic elements can be identified by assembling sequenced reads and identifying the genetic element in the assembly, or by aligning nucleotide sequences in the reads to the nucleotide sequences of a reference genetic element. The first relies on the complete assembly of the genetic element of interest, and the second relies on a reference sequence represented in nucleotides. This is particularly challenging with metagenome data, where the genetic elements, including genes, are often fragmented because sequences are shared between different species in the metagenomic data, resulting in contig breaks in or around genetic elements. This presents a difficulty when identifying genetic elements through the first approach. A common approach with metagenomes is to map reads against reference nucleotide sequences and extract the depth and coverage from those reference sequences. However, currently no software exists to identity and report genetic elements using DNA-protein alignments in metagenomes. We have developed the software Lapidary to identify the identity, coverage, depth, and most likely sequence of amino acid sequences from both genome and metagenome read files. We tested the effectiveness of the method against simulated, genomic and metagenomic read datasets. Lapidary is more sensitive than assembly methods for metagenomic data that often have fragmented assemblies but is less sensitive when assemblies are more complete, as is the case with genomic data.

Bioinformatics
Fast, sensitive detection of protein homologs using deep dense retrieval

Liang Hong,Zhihang Hu,Siqi Sun,Xiangru Tang,Jiuming Wang,Qingxiong Tan,Liangzhen Zheng,Sheng Wang,Sheng Xu,Irwin King,Mark Gerstein,Yu Li

DOI: https://doi.org/10.1038/s41587-024-02353-6

2024-08-09

Abstract:The identification of protein homologs in large databases using conventional methods, such as protein sequence comparison, often misses remote homologs. Here, we offer an ultrafast, highly sensitive method, dense homolog retriever (DHR), for detecting homologs on the basis of a protein language model and dense retrieval techniques. Its dual-encoder architecture generates different embeddings for the same protein sequence and easily locates homologs by comparing these representations. Its alignment-free nature improves speed and the protein language model incorporates rich evolutionary and structural information within DHR embeddings. DHR achieves a >10% increase in sensitivity compared to previous methods and a >56% increase in sensitivity at the superfamily level for samples that are challenging to identify using alignment-based approaches. It is up to 22 times faster than traditional methods such as PSI-BLAST and DIAMOND and up to 28,700 times faster than HMMER. The new remote homologs exclusively found by DHR are useful for revealing connections between well-characterized proteins and improving our knowledge of protein evolution, structure and function.
Sample-Align-D: A High Performance Multiple Sequence Alignment System using Phylogenetic Sampling and Domain Decomposition

Fahad Saeed,Ashfaq Khokhar

DOI: https://doi.org/10.48550/arXiv.0901.2742

2009-01-19

Abstract:Multiple Sequence Alignment (MSA) is one of the most computationally intensive tasks in Computational Biology. Existing best known solutions for multiple sequence alignment take several hours (in some cases days) of computation time to align, for example, 2000 homologous sequences of average length 300. Inspired by the Sample Sort approach in parallel processing, in this paper we propose a highly scalable multiprocessor solution for the MSA problem in phylogenetically diverse sequences. Our method employs an intelligent scheme to partition the set of sequences into smaller subsets using kmer count based similarity index, referred to as k-mer rank. Each subset is then independently aligned in parallel using any sequential approach. Further fine tuning of the local alignments is achieved using constraints derived from a global ancestor of the entire set. The proposed Sample-Align-D Algorithm has been implemented on a cluster of workstations using MPI message passing library. The accuracy of the proposed solution has been tested on standard benchmarks such as PREFAB. The accuracy of the alignment produced by our methods is comparable to that of well known sequential MSA techniques. We were able to align 2000 randomly selected sequences from the Methanosarcina acetivorans genome in less than 10 minutes using Sample-Align-D on a 16 node cluster, compared to over 23 hours on sequential MUSCLE system running on a single cluster node.

Distributed, Parallel, and Cluster Computing,Genomics,Quantitative Methods
Fast and accurate short read alignment with hybrid hash-tree data structure

Junichiro Makino,Toshikazu Ebisuzaki,Ryutaro Himeno,Yoshihide Hayashizaki

DOI: https://doi.org/10.1101/2024.02.20.581311

2024-02-23

Abstract:Rapidly increasing amount of short read data generated by NGSs (new-generation sequencers) calls for the development of fast and accurate read alignment programs. The programs based on hash table (BLAST) and Burrows-Wheeler transform (bwa-mem) are used, and the latter is known to give superior performance. We here present a new algorithm, a hybrid of hash table and suffix tree, which we designed to speed up the alignment of short reads against large reference sequences such as human genome. The total turnaround time for processing one human genome sample (read depth of 30) is just 31 minutes with our system while that was more than 25 hours with bwa-mem/gatk. The time for aligner alone is 28 minutes for our system but around 2 hours for bwa-mem. Our new algorithm is 4.4 times faster than bwa-mem while achieving similar accuracy. Variant calling and other downstream analyses after the alignment can be done with open-source tools such as SAMtools and Genome Analysis Toolkit (gatk) packages, as well as our own fast variant caller, which is well parallelized and much faster than gatk.

Bioinformatics
DIAMOND2GO: A rapid Gene Ontology assignment and enrichment tool for functional genomics

Christopher Golden,David John Studholme,Rhys A Farrer

DOI: https://doi.org/10.1101/2024.08.19.608700

2024-08-20

Abstract:DIAMOND2GO (D2GO) is a new toolset to rapidly assign Gene Ontology (GO) terms to genes or proteins based on sequence similarity searches. D2GO uses DIAMOND for alignment, which is 100 - 20,000 X faster than BLAST. D2GO leverages GO-terms already assigned to sequences in the NCBI non-redundant database to achieve rapid GO-term assignment on large sets of query sequences. In one test, 98% of the 130,184 predicted human proteins and splice variants were assigned GO-terms (>2 million in total) in < 13 minutes on a laptop computer. D2GO also features the ability to perform enrichment analysis between subsets of data, thereby allowing rapid assignment and detection of over-represented GO-terms in novel sets of sequences. D2GO is freely available under the MIT licence from https://github.com/rhysf/DIAMOND2GO

Bioinformatics
MUSCLE: a multiple sequence alignment method with reduced time and space complexity

Robert C Edgar

DOI: https://doi.org/10.1186/1471-2105-5-113

IF: 3.307

2004-01-01

BMC Bioinformatics

Abstract:BackgroundIn a previous paper, we introduced MUSCLE, a new program for creating multiple alignments of protein sequences, giving a brief summary of the algorithm and showing MUSCLE to achieve the highest scores reported to date on four alignment accuracy benchmarks. Here we present a more complete discussion of the algorithm, describing several previously unpublished techniques that improve biological accuracy and / or computational complexity. We introduce a new option, MUSCLE-fast, designed for high-throughput applications. We also describe a new protocol for evaluating objective functions that align two profiles.ResultsWe compare the speed and accuracy of MUSCLE with CLUSTALW, Progressive POA and the MAFFT script FFTNS1, the fastest previously published program known to the author. Accuracy is measured using four benchmarks: BAliBASE, PREFAB, SABmark and SMART. We test three variants that offer highest accuracy (MUSCLE with default settings), highest speed (MUSCLE-fast), and a carefully chosen compromise between the two (MUSCLE-prog). We find MUSCLE-fast to be the fastest algorithm on all test sets, achieving average alignment accuracy similar to CLUSTALW in times that are typically two to three orders of magnitude less. MUSCLE-fast is able to align 1,000 sequences of average length 282 in 21 seconds on a current desktop computer.ConclusionsMUSCLE offers a range of options that provide improved speed and / or alignment accuracy compared with currently available programs. MUSCLE is freely available at http://www.drive5.com/muscle.

biochemical research methods,biotechnology & applied microbiology,mathematical & computational biology
Fast Structural Alignment of Biomolecules Using a Hash Table, N-Grams and String Descriptors

Raphael André Bauer,Kristian Rother,Peter Moor,Knut Reinert,Thomas Steinke,Janusz Bujnicki,Robert Preissner

DOI: https://doi.org/10.3390/a2020692

2009-04-21

Algorithms

Abstract:This work presents a generalized approach for the fast structural alignment of thousands of macromolecular structures. The method uses string representations of a macromolecular structure and a hash table that stores n-grams of a certain size for searching. To this end, macromolecular structure-to-string translators were implemented for protein and RNA structures. A query against the index is performed in two hierarchical steps to unite speed and precision. In the first step the query structure is translated into n-grams, and all target structures containing these n-grams are retrieved from the hash table. In the second step all corresponding n-grams of the query and each target structure are subsequently aligned, and after each alignment a score is calculated based on the matching n-grams of query and target. The extendable framework enables the user to query and structurally align thousands of protein and RNA structures on a commodity machine and is available as open source from http://lajolla.sf.net.
Sequence alignment using large protein structure alphabets improves sensitivity to remote homologs

Robert C Edgar

DOI: https://doi.org/10.1101/2024.05.24.595840

2024-06-09

Abstract:Recent breakthroughs in protein fold prediction from amino acid sequences have unleashed a deluge of new structures, raising new opportunities for expanding insights into the universe of proteins and pursuing practical applications in bio-engineering and therapeutics while also presenting new challenges to protein search and analysis algorithms. Here, I describe Reseek, a protein alignment algorithm which improves sensitivity in protein homolog detection compared to state-of-the-art methods including DALI, TM-align and Foldseek, with improved speed over Foldseek, the fastest previous method. Reseek is based on alignment of sequences where each residue in the protein backbone is represented by a letter in a novel "mega-alphabet" of 85,899,345,920 (∼ 1011) distinct states. Code is vailable at https://github.com/rcedgar/reseek.

Bioinformatics
Protein-to-genome alignment with miniprot

Heng Li

DOI: https://doi.org/10.48550/arXiv.2210.08052

2022-12-29

Abstract:Motivation: Protein-to-genome alignment is critical to annotating genes in non-model organisms. While there are a few tools for this purpose, all of them were developed over ten years ago and did not incorporate the latest advances in alignment algorithms. They are inefficient and could not keep up with the rapid production of new genomes and quickly growing protein databases. Results: Here we describe miniprot, a new aligner for mapping protein sequences to a complete genome. Miniprot integrates recent techniques such as k-mer sketch and SIMD-based dynamic programming. It is tens of times faster than existing tools while achieving comparable accuracy on real data. Availability and implementation: <a class="link-external link-https" href="https://github.com/lh3/miniprot" rel="external noopener nofollow">this https URL</a>

Genomics
fastMSA: Accelerating Multiple Sequence Alignment with Dense Retrieval on Protein Language

Liang Hong,Siqi Sun,Liangzhen Zheng,Qingxiong Tan,Yu Li

DOI: https://doi.org/10.1101/2021.12.20.473431

2021-01-01

bioRxiv

Abstract:Evolutionarily related sequences provide information for the protein structure and function. Multiple sequence alignment, which includes homolog searching from large databases and sequence alignment, is efficient to dig out the information and assist protein structure and function prediction, whose efficiency has been proved by AlphaFold. Despite the existing tools for multiple sequence alignment, searching homologs from the entire UniProt is still time-consuming. Considering the success of AlphaFold, foreseeably, large-scale multiple sequence alignments against massive databases will be a trend in the field. It is very desirable to accelerate this step. Here, we propose a novel method, fastMSA, to improve the speed significantly. Our idea is orthogonal to all the previous accelerating methods. Taking advantage of the protein language model based on BERT, we propose a novel dual encoder architecture that can embed the protein sequences into a low-dimension space and filter the unrelated sequences efficiently before running BLAST. Extensive experimental results suggest that we can recall most of the homologs with a 34-fold speed-up. Moreover, our method is compatible with the downstream tasks, such as structure prediction using AlphaFold. Using multiple sequence alignments generated from our method, we have little performance compromise on the protein structure prediction with much less running time. fastMSA will effectively assist protein sequence, structure, and function analysis based on homologs and multiple sequence alignment. ### Competing Interest Statement The authors have declared no competing interest.
CRAFT: Compact genome Representation toward large-scale Alignment-Free daTabase

Yang Young Lu,Jiaxing Bai,Yiwen Wang,Ying Wang,Fengzhu Sun

DOI: https://doi.org/10.1093/bioinformatics/btaa699

IF: 5.8

2020-08-07

Bioinformatics

Abstract:Abstract Motivation Rapid developments in sequencing technologies have boosted generating high volumes of sequence data. To archive and analyze those data, one primary step is sequence comparison. Alignment-free sequence comparison based on k-mer frequencies offers a computationally efficient solution, yet in practice, the k-mer frequency vectors for large k of practical interest lead to excessive memory and storage consumption. Results We report CRAFT, a general genomic/metagenomic search engine to learn compact representations of sequences and perform fast comparison between DNA sequences. Specifically, given genome or high throughput sequencing data as input, CRAFT maps the data into a much smaller embedding space and locates the best matching genome in the archived massive sequence repositories. With 102−104-fold reduction of storage space, CRAFT performs fast query for gigabytes of data within seconds or minutes, achieving comparable performance as six state-of-the-art alignment-free measures. Availability and implementation CRAFT offers a user-friendly graphical user interface with one-click installation on Windows and Linux operating systems, freely available at https://github.com/jiaxingbai/CRAFT. Supplementary information Supplementary data are available at Bioinformatics online.

biochemical research methods,biotechnology & applied microbiology,mathematical & computational biology
Kraken: ultrafast metagenomic sequence classification using exact alignments

Derrick E Wood,Steven L Salzberg

DOI: https://doi.org/10.1186/gb-2014-15-3-r46

IF: 17.906

2014-01-01

Genome Biology

Abstract:Kraken is an ultrafast and highly accurate program for assigning taxonomic labels to metagenomic DNA sequences. Previous programs designed for this task have been relatively slow and computationally expensive, forcing researchers to use faster abundance estimation programs, which only classify small subsets of metagenomic data. Using exact alignment of k-mers, Kraken achieves classification accuracy comparable to the fastest BLAST program. In its fastest mode, Kraken classifies 100 base pair reads at a rate of over 4.1 million reads per minute, 909 times faster than Megablast and 11 times faster than the abundance estimation program MetaPhlAn. Kraken is available at http://ccb.jhu.edu/software/kraken/.

genetics & heredity,biotechnology & applied microbiology
MM-align: a quick algorithm for aligning multiple-chain protein complex structures using iterative dynamic programming

Srayanta Mukherjee,Yang Zhang

DOI: https://doi.org/10.1093/nar/gkp318

IF: 14.9

2009-05-14

Nucleic Acids Research

Abstract:Structural comparison of multiple-chain protein complexes is essential in many studies of protein-protein interactions. We develop a new algorithm, MM-align, for sequence-independent alignment of protein complex structures. The algorithm is built on a heuristic iteration of a modified Needleman-Wunsch dynamic programming (DP) algorithm, with the alignment score specified by the inter-complex residue distances. The multiple chains in each complex are first joined, in every possible order, and then simultaneously aligned with cross-chain alignments prevented. The alignments of interface residues are enhanced by an interface-specific weighting factor. MM-align is tested on a large-scale benchmark set of 205 x 3897 non-homologous multiple-chain complex pairs. Compared with a naïve extension of the monomer alignment program of TM-align, the alignment accuracy of MM-align is significantly higher as judged by the average TM-score of the physically-aligned residues. MM-align is about two times faster than TM-align because of omitting the cross-alignment zone of the DP matrix. It also shows that the enhanced alignment of the interfaces helps in identifying biologically relevant protein complex pairs.

biochemistry & molecular biology
IMPACT: Interval-based Multi-pass Proteomic Alignment with Constant Traceback

Sahand Kashani,Stuart Byma,James R. Larus

DOI: https://doi.org/10.48550/arXiv.1902.03238

2019-02-10

Abstract:Darwin is a genomics co-processor that achieved a 15000x acceleration on long read assembly through innovative hardware and algorithm co-design. Darwins algorithms and hardware implementation were specifically designed for DNA analysis pipelines. This paper analyzes the feasibility of applying Darwins algorithms to the problem of protein sequence alignment. In addition to a behavioral analysis of Darwin when aligning proteins, we propose an algorithmic improvement to Darwins alignment algorithm, GACT, in the form of a multi-pass variant that increases its accuracy on protein sequence alignment. Concretely, our proposed multi-pass variant of GACT achieves on average 14\% better alignment scores.

Quantitative Methods
FastMG: a simple, fast, and accurate maximum likelihood procedure to estimate amino acid replacement rate matrices from large data sets

Cuong Cao Dang,Vinh Sy Le,Olivier Gascuel,Bart Hazes,Quang Si Le

DOI: https://doi.org/10.1186/1471-2105-15-341

IF: 3.307

2014-10-24

BMC Bioinformatics

Abstract:Background: Amino acid replacement rate matrices are a crucial component of many protein analysis systems such as sequence similarity search, sequence alignment, and phylogenetic inference. Ideally, the rate matrix reflects the mutational behavior of the actual data under study; however, estimating amino acid replacement rate matrices requires large protein alignments and is computationally expensive and complex. As a compromise, sub-optimal pre-calculated generic matrices are typically used for protein-based phylogeny. Sequence availability has now grown to a point where problem-specific rate matrices can often be calculated if the computational cost can be controlled. Results: The most time consuming step in estimating rate matrices by maximum likelihood is building maximum likelihood phylogenetic trees from protein alignments. We propose a new procedure, called FastMG, to overcome this obstacle. The key innovation is the alignment-splitting algorithm that splits alignments with many sequences into non-overlapping sub-alignments prior to estimating amino acid replacement rates. Experiments with different large data sets showed that the FastMG procedure was an order of magnitude faster than without splitting. Importantly, there was no apparent loss in matrix quality if an appropriate splitting procedure is used. Conclusions: FastMG is a simple, fast and accurate procedure to estimate amino acid replacement rate matrices from large data sets. It enables researchers to study the evolutionary relationships for specific groups of proteins or taxa with optimized, data-specific amino acid replacement rate matrices. The programs, data sets, and the new mammalian mitochondrial protein rate matrix are available at http://fastmg.codeplex.com.

Fast and sensitive protein alignment using DIAMOND

Sensitive protein alignments at tree-of-life scale using DIAMOND

ProteinInferencer: Confident protein identification and multiple experiment comparison for large scale proteomics projects

Fast and accurate protein structure search with Foldseek

Rapid and Sensitive Protein Complex Alignment with Foldseek-Multimer

Fast-Part: Fast and Accurate Data Partitioning for Biological Sequence Analysis

Lapidary: Identifying and reporting amino acid sequences in metagenomes using sequence reads and Diamond

Fast, sensitive detection of protein homologs using deep dense retrieval

Sample-Align-D: A High Performance Multiple Sequence Alignment System using Phylogenetic Sampling and Domain Decomposition

Fast and accurate short read alignment with hybrid hash-tree data structure

DIAMOND2GO: A rapid Gene Ontology assignment and enrichment tool for functional genomics

MUSCLE: a multiple sequence alignment method with reduced time and space complexity

Fast Structural Alignment of Biomolecules Using a Hash Table, N-Grams and String Descriptors

Sequence alignment using large protein structure alphabets improves sensitivity to remote homologs

Protein-to-genome alignment with miniprot

fastMSA: Accelerating Multiple Sequence Alignment with Dense Retrieval on Protein Language

CRAFT: Compact genome Representation toward large-scale Alignment-Free daTabase

Kraken: ultrafast metagenomic sequence classification using exact alignments

MM-align: a quick algorithm for aligning multiple-chain protein complex structures using iterative dynamic programming

IMPACT: Interval-based Multi-pass Proteomic Alignment with Constant Traceback

FastMG: a simple, fast, and accurate maximum likelihood procedure to estimate amino acid replacement rate matrices from large data sets