Abstract:Accurate multiple sequence alignment (MSA) is imperative for the comprehensive analysis of biological sequences. However, a notable challenge arises as no single MSA tool consistently outperforms its counterparts across diverse datasets. Users often have to try multiple MSA tools to achieve optimal alignment results, which can be time-consuming and memory-intensive. While the overall accuracy of certain MSA results may be lower, there could be local regions with the highest alignment scores, prompting researchers to seek a tool capable of merging these locally optimal results from multiple initial alignments into a globally optimal alignment. In this study, we introduce Two Pointers Meta-Alignment (TPMA), a novel tool designed for the integration of nucleic acid sequence alignments. TPMA employs two pointers to partition the initial alignments into blocks containing identical sequence fragments. It selects blocks with the high sum of pairs (SP) scores to concatenate them into an alignment with an overall SP score superior to that of the initial alignments. Through tests on simulated and real datasets, the experimental results consistently demonstrate that TPMA outperforms M-Coffee in terms of aSP, Q, and total column (TC) scores across most datasets. Even in cases where TPMA's scores are comparable to M-Coffee, TPMA exhibits significantly lower running time and memory consumption. Furthermore, we comprehensively assessed all the MSA tools used in the experiments, considering accuracy, time, and memory consumption. We propose accurate and fast combination strategies for small and large datasets, which streamline the user tool selection process and facilitate large-scale dataset integration. The dataset and source code of TPMA are available on GitHub (https://github.com/malabz/TPMA). Accurate multiple sequence alignment (MSA) is vital for comprehensive biological sequence analysis. However, as no single MSA tool consistently outperforms others across diverse datasets, researchers must invest significant time exploring multiple tools to identify the most suitable one for their specific dataset. To address this, researchers seek tools that can merge locally optimal results from diverse initial alignments into a globally optimal alignment. Our novel approach, Two Pointers Meta-Alignment (TPMA), employs a two-pointer to partition initial alignments into blocks, selecting those with the higher sum of pairs (SP) scores for integration into a globally optimal alignment. TPMA consistently outperforms M-Coffee, demonstrating superior aSP, Q, and total column (TC) scores, coupled with faster running times and lower memory consumption. We present comprehensive assessments of various MSA tools, proposing efficient combination strategies for diverse datasets. Our tool, TPMA, and associated resources are publicly available on GitHub (https://github.com/malabz/TPMA), offering a valuable contribution to the field of evolutionary biology and streamlining the selection process for users dealing with large-scale datasets.

Minimap2: pairwise alignment for nucleotide sequences

New strategies to improve minimap2 alignment accuracy

Minimap and miniasm: fast mapping and de novo assembly for noisy long sequences

MinimapR: A Parallel Alignment Tool for the Analysis of Large-Scale Third-Generation Sequencing Data

GMAP and GSNAP for Genomic Sequence Alignment: Enhancements to Speed, Accuracy, and Functionality

mm2-gb: GPU Accelerated Minimap2 for Long Read DNA Mapping

An efficient Burrows-Wheeler transform-based aligner for short read mapping

Efficient end-to-end long-read sequence mapping using minimap2-fpga integrated with hardware accelerated chaining

SaAlign: Multiple DNA/RNA Sequence Alignment and Phylogenetic Tree Construction Tool for Ultra-Large Datasets and Ultra-Long Sequences Based on Suffix Array

Accelerating Minimap2 for Accurate Long Read Alignment on GPUs

GenMPI: Cluster Scalable Variant Calling for Short/Long Reads Sequencing Data

Umap: Use Unique Sequence for Alignment of Short Sequence Reads and SNP Detection

Short Read Alignment Based On Maximal Approximate Match Seeds

pathMap: a path-based mapping tool for long noisy reads with high sensitivity

Protein-to-genome alignment with miniprot

kngMap: Sensitive and Fast Mapping Algorithm for Noisy Long Reads Based on the K-Mer Neighborhood Graph

Sensitive Long-Indel-Aware Alignment of Sequencing Reads

ReadsMap: a new tool for high precision mapping of DNAseq and RNAseq read sequences

TPMA: A two pointers meta-alignment tool to ensemble different multiple nucleic acid sequence alignments

STAR: ultrafast universal RNA-seq aligner