Abstract:Background: Several problems exist with current methods used to align DNA sequences for comparative sequence analysis. Most dynamic programming algorithms assume that conserved sequence elements are collinear. This assumption appears valid when comparing orthologous protein coding sequences. Functional constraints on proteins provide strong selective pressure against sequence inversions, and minimize sequence duplications and feature shuffling. For non-coding sequences this collinearity assumption is often invalid. For example, enhancers contain clusters of transcription factor binding sites that change in number, orientation, and spacing during evolution yet the enhancer retains its activity. Dot plot analysis is often used to estimate non-coding sequence relatedness. Yet dot plots do not actually align sequences and thus cannot account well for base insertions or deletions. Moreover, they lack an adequate statistical framework for comparing sequence relatedness and are limited to pairwise comparisons. Lastly, dot plots and dynamic programming text outputs fail to provide an intuitive means for visualizing DNA alignments. Results: To address some of these issues, we created a stand alone, platform independent, graphic alignment tool for comparative sequence analysis (GATA http://gata.sourceforge.net/). GATA uses the NCBI-BLASTN program and extensive post-processing to identify all small sub-alignments above a low cut-off score. These are graphed as two shaded boxes, one for each sequence, connected by a line using the coordinate system of their parent sequence. Shading and colour are used to indicate score and orientation. A variety of options exist for querying, modifying and retrieving conserved sequence elements. Extensive gene annotation can be added to both sequences using a standardized General Feature Format (GFF) file. Conclusions: GATA uses the NCBI-BLASTN program in conjunction with post-processing to exhaustively align two DNA sequences. It provides researchers with a fine-grained alignment and visualization tool aptly suited for non-coding, 0-200 kb, pairwise, sequence analysis. It functions independent of sequence feature ordering or orientation, and readily visualizes both large and small sequence inversions, duplications, and segment shuffling. Since the alignment is visual and does not contain gaps, gene annotation can be added to both sequences to create a thoroughly descriptive picture of DNA conservation that is well suited for comparative sequence analysis.

BAD2matrix: Phylogenomic matrix concatenation, indel coding, and more

SequenceMatrix: concatenation software for the fast assembly of multi‐gene datasets with character set and codon information

Agalma: an automated phylogenomics workflow

AMAS: a fast tool for alignment manipulation and computing of summary statistics

Biological Data Analysis Program (BDAP): a Multitasking Biological Sequence Analysis Program

DAMBE5: A Comprehensive Software Package for Data Analysis in Molecular Biology and Evolution

ProSeq4: A user‐friendly multiplatform program for preparation and analysis of large‐scale DNA polymorphism datasets

AIR: A batch-oriented web program package for construction of supermatrices ready for phylogenomic analyses

GATA: a graphic alignment tool for comparative sequence analysis

A general substitution matrix for structural phylogenetics.

Phylobook: a tool for display, clade annotation and extraction of sequences from molecular phylogenies

DAMBE7: New and Improved Tools for Data Analysis in Molecular Biology and Evolution

Dadasnake, a Snakemake implementation of DADA2 to process amplicon sequencing data for microbial ecology

PhyloSuite: An integrated and scalable desktop platform for streamlined molecular sequence data management and evolutionary phylogenetics studies

Minimap2: pairwise alignment for nucleotide sequences

Please Mind the Gap: Indel-Aware Parsimony for Fast and Accurate Ancestral Sequence Reconstruction and Multiple Sequence Alignment including Long Indels

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

diverse-seq: an application for alignment-free selecting and clustering biological sequences

MAVID: Constrained Ancestral Alignment of Multiple Sequences

PyCoM: a python library for large-scale analysis of residue-residue coevolution data