Abstract:Contigs assembled from the third-generation sequencing long reads are usually more complete than the second-generation short reads. However, the current algorithms still have difficulty in assembling the long reads into the ideal complete and accurate genome, or the theoretical best result [1]. To improve the long read contigs and with more and more fully sequenced genomes available, it could still be possible to use the similar genome-assisted reassembly method [2], which was initially proposed for the short reads making use of a closely related genome (similar genome) to the sequencing genome (target genome). The method aligns the contigs and reads to the similar genome, and then extends and refines the aligned contigs with the aligned reads. Here, we introduce AlignGraph2, a similar genome-assisted reassembly pipeline for the PacBio long reads. The AlignGraph2 pipeline is the second version of AlignGraph algorithm proposed by us but completely redesigned, can be inputted with either error-prone or HiFi long reads, and contains four novel algorithms: similarity-aware alignment algorithm and alignment filtration algorithm for alignment of the long reads and preassembled contigs to the similar genome, and reassembly algorithm and weight-adjusted consensus algorithm for extension and refinement of the preassembled contigs. In our performance tests on both error-prone and HiFi long reads, AlignGraph2 can align 5.7-27.2% more long reads and 7.3-56.0% more bases than some current alignment algorithm and is more efficient or comparable to the others. For contigs assembled with various de novo algorithms and aligned to similar genomes (aligned contigs), AlignGraph2 can extend 8.7-94.7% of them (extendable contigs), and obtain contigs of 7.0-249.6% larger N50 value and 5.2-87.7% smaller number of indels per 100 kbp (extended contigs). With genomes of decreased similarities, AlignGraph2 also has relatively stable performance. The AlignGraph2 software can be downloaded for free from this site: https://github.com/huangs001/AlignGraph2.

Champuru 2: Improved scoring of alignments and a user-friendly graphical interface

SeqKit2: A Swiss army knife for sequence and alignment processing

CHSalign: A Web Server That Builds upon Junction-Explorer and RNAJAG for Pairwise Alignment of RNA Secondary Structures with Coaxial Helical Stacking

Tracy: basecalling, alignment, assembly and deconvolution of sanger chromatogram trace files

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

Chrom-pro: A User-Friendly Toolkit for De-novo Chromosome Assembly and Genomic Analysis

AlignGraph2: similar genome-assisted reassembly pipeline for PacBio long reads

Mixed Sequence Reader: A Program for Analyzing DNA Sequences with Heterozygous Base Calling

Sequenceserver: A Modern Graphical User Interface for Custom BLAST Databases

CESAR 2.0 substantially improves speed and accuracy of comparative gene annotation

DiGAlign: Versatile and Interactive Visualization of Sequence Alignment for Comparative Genomics

AAscan, PCRdesign and MutantChecker: A Suite of Programs for Primer Design and Sequence Analysis for High-Throughput Scanning Mutagenesis

RealSeq2: a Software Integrated with UMI Identification, Error Correction, and Methylation Modifications Storing

STAR: ultrafast universal RNA-seq aligner

RNAprofiling 2.0: Enhanced cluster analysis of structural ensembles

New strategies to improve minimap2 alignment accuracy

GUIDANCE: a web server for assessing alignment confidence scores

UCHIME improves sensitivity and speed of chimera detection

Acceleration of short and long DNA read mapping without loss of accuracy using suffix array

CGAP-align: a high performance DNA short read alignment tool.

GSP: a web-based platform for designing genome-specific primers in polyploids.