Sign in

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.

PDFDOICode

02 /Specification & Classification

Click any tag to browse matching papers

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.

Similarity Map
SelectedSame fieldCross-field
THISKowalski '20Hou '21Sarkar '17Kowalski '19Coster '23Fotouhi '18
61% semantic matchBioinformatics
Tomasz Kowalski, Szymon Grabowski·bioRxiv·May 1, 2020· 3 citations