Abstract:MOTIVATIONNext-generation sequencing (NGS) technologies generate large amounts of short read data for many different organisms. The fact that NGS reads are generally short makes it challenging to assemble the reads and reconstruct the original genome sequence. For clustering genomes using such NGS data, word-count based alignment-free sequence comparison is a promising approach, but for this approach, the underlying expected word counts are essential.A plausible model for this underlying distribution of word counts is given through modeling the DNA sequence as a Markov chain (MC). For single long sequences, efficient statistics are available to estimate the order of MCs and the transition probability matrix for the sequences. As NGS data do not provide a single long sequence, inference methods on Markovian properties of sequences based on single long sequences cannot be directly used for NGS short read data.RESULTSHere we derive a normal approximation for such word counts. We also show that the traditional Chi-square statistic has an approximate gamma distribution ,: using the Lander-Waterman model for physical mapping. We propose several methods to estimate the order of the MC based on NGS reads and evaluate those using simulations. We illustrate the applications of our results by clustering genomic sequences of several vertebrate and tree species based on NGS reads using alignment-free sequence dissimilarity measures. We find that the estimated order of the MC has a considerable effect on the clustering results ,: and that the clustering results that use a N: MC of the estimated order give a plausible clustering of the species.AVAILABILITY AND IMPLEMENTATIONOur implementation of the statistics developed here is available as R package 'NGS.MC' at http://www-rcf.usc.edu/∼fsun/Programs/NGS-MC/NGS-MC.htmlCONTACTfsun@usc.eduSUPPLEMENTARY INFORMATIONSupplementary data are available at Bioinformatics online.

Comparison of Next-Generation Sequencing Samples Using Compression-Based Distances and Its Application to Phylogenetic Reconstruction

Alignment-Free Sequence Analysis and Applications

New Developments Of Alignment-Free Sequence Comparison: Measures, Statistics And Next-Generation Sequencing

Alignment-Free Sequence Comparison Based on Next Generation Sequencing Reads: Extended Abstract.

Evaluation and Comparison of Multiple Aligners for Next-Generation Sequencing Data Analysis.

Alignment-free sequence comparison based on next-generation sequencing reads.

Comparison of the Two Up-to-date Sequencing Technologies for Genome Assembly: HiFi Reads of Pacific Biosciences Sequel II System and Ultralong Reads of Oxford Nanopore

Co-Phylog: an Assembly-Free Phylogenomic Approach for Closely Related Organisms

Novel Computational Technologies for Next-Generation Sequencing Data Analysis and Their Applications.

A Practical Comparison Of De Novo Genome Assembly Software Tools For Next-Generation Sequencing Technologies

Comparing DNA sequence collections by direct comparison of compressed text indexes

Research of NGS Alignment Algorthims

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

Comparison of the Experimental Methods in Haplotype Sequencing Via Next Generation Sequencing.

[Alignment-Free Biomolecular Sequence Comparison Method].

HUGO: Hierarchical Multi-Reference Genome Compression for Aligned Reads.

Inference of Markovian Properties of Molecular Sequences from NGS Data and Applications to Comparative Genomics

Comparison of high-throughput sequencing data compression tools

Distangsd: Fast and Accurate Inference of Genetic Distances for Next-Generation Sequencing Data

Recent advancement in Next Generation Sequencing techniques and its computational analysis

Comparison of the Two Major Classes of Assembly Algorithms: Overlap-Layout-consensus and De-Bruijn-graph