SAMBLASTER : fast duplicate marking and structural variant read extraction

Abstract

Motivation: Illumina DNA sequencing is now the predominant source of raw genomic data, and data volumes are growing rapidly. Bioinformatic analysis pipelines are having trouble keeping pace. A common bottleneck in such pipelines is the requirement to read, write, sort and compress large BAM files multiple times. Results: We present SAMBLASTER, a tool that reduces the number of times such costly operations are performed. SAMBLASTER is designed to mark duplicates in read-sorted SAM files as a piped post-pass on DNA aligner output before it is compressed to BAM. In addition, it can simultaneously output into separate files the discordant read-pairs and/or split-read mappings used for structural variant calling. As an alignment post-pass, its own runtime overhead is negligible, while dramatically reducing overall pipeline complexity and runtime. As a stand-alone duplicate marking tool, it performs significantly better than PICARD or SAMBAMBA in terms of both speed and memory usage, while achieving nearly identical results. Availability and implementation: SAMBLASTER is open-source C++ code and freely available for download from https://github.com/GregoryFaust/samblaster. Contact: imh4y@virginia.edu

Full text

Its licence does not allow us to show the full text. Read it at PubMed Central:

References 3

1. 1.

   Faust GG, Hall IM. YAHA: fast and flexible long-read alignment with optimal breakpoint detection. Bioinformatics. 2012;28:2417-2424.

   Not in our library10.1093/bioinformatics/bts456

2. 2.

   Li H, et al. The Sequence Alignment/Map format and SAMtools. Bioinformatics. 2009;25:2078-2079.

   Not in our library10.1093/bioinformatics/btp352

3. 3.

   Quinlan AR, et al. Genome-wide mapping and assembly of structural variant breakpoints in the mouse genome. Genome Res. 2010;20:623-635.

   Not in our library10.1101/gr.102970.109

Cited by papers we hold 26

• A chromosome-scale genome assembly of the false clownfish, <i>Amphiprion ocellaris</i> · 2022
  Not checked
• Nuclear-embedded mitochondrial DNA sequences in 66,083 human genomes · 2022
  No notice
• Detection of copy number variants in African goats using whole genome sequence data · 2021
  Not checked
• Multimodal analysis of cell-free DNA whole-genome sequencing for pediatric cancers with low mutational burden · 2021
  No notice
• Samplot: a platform for structural variant visual validation and automated filtering · 2021
  No notice
• Lung cancer organoids analyzed on microwell arrays predict drug responses of patients within a week · 2021
  No notice
• Acquired resistance to DZNep-mediated apoptosis is associated with copy number gains of AHCY in a B-cell lymphoma model · 2020
  Not checked
• Unmasking Retinitis Pigmentosa complex cases by a whole genome sequencing algorithm based on open-access tools: hidden recessive inheritance and potential oligogenic variants · 2020
  Not checked
• Best practices for variant calling in clinical sequencing · 2020
  No notice
• Histone H2AK119 Mono-Ubiquitination Is Essential for Polycomb-Mediated Transcriptional Repression · 2020
  No notice
• Developing a High-Quality Linkage Map for the Atlantic Killifish Fundulus heteroclitus · 2019
  Not checked
• Increased ultra-rare variant load in an isolated Scottish population impacts exonic and regulatory regions. · 2019
  Not checked
• Integrative and comparative genomic analyses identify clinically relevant pulmonary carcinoid groups and unveil the supra-carcinoids · 2019
  No notice
• Clinical validation of the tempus xT next-generation targeted oncology sequencing assay · 2019
  No notice
• Large, three-generation human families reveal post-zygotic mosaicism and variability in germline mutation accumulation · 2019
  No notice
• The GenomeAsia 100K Project enables genetic discoveries across Asia · 2019
  No notice
• SyRI: finding genomic rearrangements and local sequence differences from whole-genome assemblies · 2019
  No notice
• Functional equivalence of genome sequencing analysis pipelines enables harmonized variant calling across human genetics projects · 2018
  No notice
• Defining the diverse spectrum of inversions, complex structural variation, and chromothripsis in the morbid human genome · 2017
  No notice
• Population scale mapping of transposable element diversity reveals links to gene regulation and epigenomic variation · 2016
  No notice

and 6 more

Reads

A person counts once a day per paper. We keep no cookie and no address: a reader is a hash that changes every day. Any month under five reads shows no number.