Faster and More Accurate Sequence Alignment with SNAP

Abstract

We present the Scalable Nucleotide Alignment Program (SNAP), a new short and long read aligner that is both more accurate (i.e., aligns more reads with fewer errors) and 10-100× faster than state-of-the-art tools such as BWA. Unlike recent aligners based on the Burrows-Wheeler transform, SNAP uses a simple hash index of short seed sequences from the genome, similar to BLAST's. However, SNAP greatly reduces the number and cost of local alignment checks performed through several measures: it uses longer seeds to reduce the false positive locations considered, leverages larger memory capacities to speed index lookup, and excludes most candidate locations without fully computing their edit distance to the read. The result is an algorithm that scales well for reads from one hundred to thousands of bases long and provides a rich error model that can match classes of mutations (e.g., longer indels) that today's fast aligners ignore. We calculate that SNAP can align a dataset with 30× coverage of a human genome in less than an hour for a cost of $2 on Amazon EC2, with higher accuracy than BWA. Finally, we describe ongoing work to further improve SNAP.

Abattoir survey of bovine tuberculosis in tanta, centre of the Nile delta, with in silico analysis of gene mutations and protein–protein interactions of the involved mycobacteria

“…• Then, alignment process to the reference genome H37Rv was done using BWA and BCF tools labelled the small variants (Zaharia et al, 2011).…”

A Chemical Genetic Method for Monitoring Genome-Wide Dynamics of O-GlcNAc Turnover on Chromatin-Associated Proteins

“…After sequencing two independent replicates from each time point in the TC O-GlcNAc-seq experiment in both WT and OGA-null, we normalized the number of sequenced Drosophila reads using the spiked-in yeast controls by alignment to the yeast genome. 25,26 Genomic regions bearing O-GlcNAz modified proteins were determined by quantifying peaks, which represent a large number of overlapping sequencing reads observed at specific genomic loci, from the experimental samples at 0 h as compared to the vehicle fed Drosophila. Using this strategy, we identified a total of 1883 peaks present in both WT and OGA-null TC O-GlcNAc-seq Drosophila samples.…”

Increased Heat Resilience of Intraspecific Outbred Compared to Inbred Lineages in the Kelp Laminaria digitata: Physiology and Transcriptomics

“…The L. digitata transcriptome was assembled de novo from 33 libraries (n = 3) sampled on day 0 (10 • C, N = 12) and day 18 (10 • C, N = 12, and 20.5 • C, N = 9) of experiment 2 using rnaSPAdes (Bushmanova et al, 2019) with kmer length 37 on a Linux server (40 cores and 250Gb RAM). Quality assessment of the assembly was performed using Transrate v1.0.3 (Smith-Unna et al, 2016) with read-mapping metrics (SNAP; Zaharia et al, 2011;Salmon;Patro et al, 2015). A reference set of transcripts was built by predicting open reading frames (ORFs) from the high-quality contigs returned by Transrate using FragGeneScan (Rho et al, 2010).…”

“…As polyunsaturated fatty acids (PUFAs) and high temperature increase cell membrane fluidity, decreasing the expression of a lipid desaturase (linoleoyl-CoA desaturase) might maintain cell membrane rigidity in response to stress (Neidleman, 1987;Tatsuzawa et al, 1996;Maulucci et al, 2016). Modifications to sulfolipid (sulfoquinovosyltransferase) and glycolipid metabolism (1,3-propanediol dehydrogenase) indicate the maintenance of photosynthetic membranes in the chloroplasts (Boudière et al, 2014;Zhan et al, 2019).…”