Comparing DNA sequence collections by direct comparison of compressed text indexes
From MaRDI portal
Abstract: Popular sequence alignment tools such as BWA convert a reference genome to an indexing data structure based on the Burrows-Wheeler Transform (BWT), from which matches to individual query sequences can be rapidly determined. However the utility of also indexing the query sequences themselves remains relatively unexplored. Here we show that an all-against-all comparison of two sequence collections can be computed from the BWT of each collection with the BWTs held entirely in external memory, i.e. on disk and not in RAM. As an application of this technique, we show that BWTs of transcriptomic and genomic reads can be compared to obtain reference-free predictions of splice junctions that have high overlap with results from more standard reference-based methods. Code to construct and compare the BWT of large genomic data sets is available at http://beetl.github.com/BEETL/ as part of the BEETL library.
Recommendations
Cited in
(9)- Detecting mutations by eBWT
- The Burrows-Wheeler transform between data compression and combinatorics on words
- When a dollar makes a BWT
- An external-memory algorithm for string graph construction
- Comparison of BWA-MEM and NovoAlign using human whole-genome next-generation sequencing reads
- BEETL
- An Experimental Study of Compressed Indexing and Local Alignments of DNA
- A head-to-head comparison of de Bruijn indices and names
- Generalized substring compression
This page was built for publication: Comparing DNA sequence collections by direct comparison of compressed text indexes
Report a bug (only for logged in users!)Click here to report a bug for this page (MaRDI item Q4916910)