Evaluation of sequencing reads at scale using rdeval
Giulio FormentiiD, Bonhwang KooiD, Marco SollittoiD, Jennifer BalaccoiD, Nadolina BrajukaiD, Richard C. Burhans, Erick DuarteiD, Alice Maria GianiiD, Kirsty McCaffreyiD, Jack A MedicoiD, Eugene Wimberly MyersiD, Patrik Smeds, Anton NekrutenkoiD, Erich David JarvisiD
Abstract
MotivationLarge sequencing data sets are produced and deposited into public archives at unprecedented rates. The availability of tools that can reliably and efficiently generate and store sequencing read summary statistics has become critical. ResultsAs part of the effort by the Vertebrate Genomes Project (VGP) to generate high-quality reference genomes at scale, we sought to address the community need for efficient sequencing data evaluation by developing rdeval, a standalone tool to quickly compute and dynamically display sequencing read metrics. Rdeval can either run on the fly or store key sequence data metrics in read sketches, with dramatic compression gains. Statistics can then be efficiently recalled from sketches for additional processing. Rdeval can convert fa*[.gz] files to and from other popular formats including BAM and CRAM for better compression. Overall, while CRAM achieves the best compression, the gain is marginal, and BAM achieves the best compromise between data compression and accessing speed. Rdeval also generates a detailed visual report with multiple data analytics that can be exported in various formats. We showcase rdevals functionalities using human and VGP read data from different sequencing platforms and species. For PacBio long-read sequencing, our analysis shows dramatic improvements both in read length and quality over time, and a benefit of additional coverage for genome assembly. Availability and implementationRdeval is implemented in C++ for data processivity and in R for data visualization. Precompiled releases (Linux, MacOS, Windows) and commented source code for rdeval are available under MIT license at https://github.com/vgl-hub/rdeval. Documentation is available using ReadTheDocs (https://rdeval-documentation.readthedocs.io). Rdeval is also available in Bioconda and in Galaxy (https://usegalaxy.org). An automated test workflow ensures the consistency of software updates. Supplementary informationSupplementary data are available at Bioinformatics online.
02 /Specification & Classification
Click any tag to browse matching papers- Citations
- 4
- FWCI
- 1.21
- Field percentile
- top 23%
- Organism
- Assay
- Research Focus
- Approach
Reproducibility ledger
Nobody has reported trying to reproduce this paper yet. Did you replicate an experiment or run the code? Your report saves the next lab a week.
Sign in to file a reproducibility report.
Related papers
- Engineering the Compression of Sequencing Reads
Tomasz Kowalski, Szymon Grabowski·bioRxiv·May 1, 2020
- Efficient storage and analysis of quantitative genomics data with the Dense Depth Data Dump (D4) format and d4tools.
Hao Hou, Brent S. Pedersen, Aaron R. Quinlan·Nature Computational SciencePeer-reviewed·Jun 21, 2021
- Quark enables semi-reference-based compression of RNA-seq data
Hirak Sarkar, Rob Patro·BioinformaticsPeer-reviewed·Jun 29, 2017
- PgRC: Pseudogenome based Read Compressor
Tomasz Kowalski, Szymon Grabowski·BioinformaticsPeer-reviewed·Dec 5, 2019
- NanoPack2: Population scale evaluation of long-read sequencing data
Wouter De Coster, Rosa Rademakers·BioinformaticsPeer-reviewed·May 1, 2023
- Quality Assessment of High-throughput DNA Sequencing Data via Range analysis
Ali Fotouhi, Mina Majidi, Muhammed Oğuzhan Külekci·Lecture notes in computer sciencePeer-reviewed·Jan 1, 2018