Abstract:The development of second-generation sequencing (SGS) technology has provided scientists with a myriad of opportunities as well as new challenges. SGS machines are capable of sequencing billions of short reads at a fraction of the cost and time in comparison to older technology. Often, the study of sequence data begins with the alignment of billions of short dna reads to the 3 billion base pair human reference genome, a daunting computational task, especially if the error-rate between the reads and reference is high. For this reason, PerM was developed to use periodic spaced seeds to efficiently and accurately provide highly sensitive ungapped alignment for Illumina and SOLiD reads. Inexact alignments are often the most interesting biologically, because mismatches between the read and reference are often the result of genetic variation. To accurately detect and discern variation from machine errors, we developed ComB, which iteratively applies Bayesian statistics to color or base alignment to accurately determine mutation probability. This allowed us to study a host of biological phenomena which result in rare nucleotide differences, including single nucleotide polymorphisms (SNPs), RNA-editing, and allele-specific expression. DNA-methylation of cytosine residues also produces single-base mismatches when dna is treated with sodium-bisulte which changes all unmethylated cytosine residues to thymine. To accurately estimate methylation rates from sodium bisulte treated dna we developed FadE, an algorithm which uses Newton-Raphson optimization to estimate the methylation rate at every cytosine residue in the genome. Finally, we have applied all our statistical tools to study human mRNA editing, and have shown that RNA editing in human brain tissue occurs at a much lower rate than previously thought.

FANSe2: a robust and cost-efficient alignment tool for quantitative next-generation sequencing applications.

FANSe: an accurate algorithm for quantitative mapping of large scale sequencing reads.

Application of the Hyper-Accurate Mapping Algorithm FANSe for Next-Generation Sequencing in Non-Model Organisms

A Bidirectional Fuzzy Index and Approximate Search Algorithm for Next Generation Sequencing

BaseFormer: Transformer Based Base-Caller for Fast and Accurate Next Generation Sequencing

Salmon-Easy: an Alignment-Free RNA-seq Quantitative Strategy Based on Partial Genome Indexing.

Low-cost, Low-bias and Low-input RNA-seq with High Experimental Verifiability Based on Semiconductor Sequencing

Ultrafast functional profiling of RNA-seq data for nonmodel organisms

STAR: ultrafast universal RNA-seq aligner

A fast read alignment method based on seed-and-vote for next generation sequencing

RNASEQR—a Streamlined and Accurate RNA-seq Sequence Analysis Program

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

Fast and Accurate Read Alignment for Resequencing.

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

ResSeq: Enhancing Short-Read Sequencing Alignment by Rescuing Error-Containing Reads

An Efficient Filtration Method Based on Variable-Length Seeds for Sequence Alignment.

Efficient Seeding for Error-Prone Sequences with SubseqHash2

Too many needles in this haystack: algorithms for the analysis of next generation sequence data

A Novel Multi-Alignment Pipeline for High-Throughput Sequencing Data.

RNA-SeQC: RNA-seq metrics for quality control and process optimization

RabbitQC: high-speed scalable quality control for sequencing data