Abstract:Abstract Background Next-generation sequencing allows genome-wide analysis of changes in chromatin states and gene expression. Data analysis of these increasingly used methods either requires multiple analysis steps, or extensive computational time. We sought to develop a tool for rapid quantification of sequencing peaks from diverse experimental sources and an efficient method to produce coverage tracks for accurate visualization that can be intuitively displayed and interpreted by experimentalists with minimal bioinformatics background. We demonstrate its strength and usability by integrating data from several types of sequencing approaches. Results We have developed BAMscale , a one-step tool that processes a wide set of sequencing datasets. To demonstrate the usefulness of BAMscale , we analyzed multiple sequencing datasets from chromatin immunoprecipitation sequencing data (ChIP-seq), chromatin state change data (assay for transposase-accessible chromatin using sequencing: ATAC-seq, DNA double-strand break mapping sequencing: END-seq), DNA replication data (Okazaki fragments sequencing: OK-seq, nascent-strand sequencing: NS-seq, single-cell replication timing sequencing: scRepli-seq) and RNA-seq data. The outputs consist of raw and normalized peak scores (multiple normalizations) in text format and scaled bigWig coverage tracks that are directly accessible to data visualization programs. BAMScale also includes a visualization module facilitating direct, on-demand quantitative peak comparisons that can be used by experimentalists. Our tool can effectively analyze large sequencing datasets (~ 100 Gb size) in minutes, outperforming currently available tools. Conclusions BAMscale accurately quantifies and normalizes identified peaks directly from BAM files, and creates coverage tracks for visualization in genome browsers. BAMScale can be implemented for a wide set of methods for calculating coverage tracks, including ChIP-seq and ATAC-seq, as well as methods that currently require specialized, separate tools for analyses, such as splice-aware RNA-seq, END-seq and OK-seq for which no dedicated software is available. BAMscale is freely available on github ( https://github.com/ncbi/BAMscale ).

Twelve years of SAMtools and BCFtools

HTSlib: C library for reading/writing high-throughput sequencing data

TBtools: An Integrative Toolkit Developed for Interactive Analyses of Big Biological Data

The Sequence Alignment/Map format and SAMtools

Over 1000 Tools Reveal Trends in the Single-Cell RNA-seq Analysis Landscape

Sam2bam: High-Performance Framework for NGS Data Preprocessing Tools

Biopython: freely available Python tools for computational molecular biology and bioinformatics

A spectrum of free software tools for processing the VCF variant call format: vcflib, bio-vcf, cyvcf2, hts-nim and slivar

Vcflib and tools for processing the VCF variant call format

BCFtools/liftover: an accurate and comprehensive tool to convert genetic variants across genome assemblies

TBtools-II: A "One for All, All for One" Bioinformatics Platform for Biological Big-data Mining

vcfpp: a C++ API for rapid processing of the Variant Call Format

GenoTools: An Open-Source Python Package for Efficient Genotype Data Quality Control and Analysis

A data-supported history of bioinformatics tools

Alfred: interactive multi-sample BAM alignment statistics, feature counting and feature annotation for long- and short-read sequencing

Maftools: efficient and comprehensive analysis of somatic variants in cancer

Bioconductor: open software development for computational biology and bioinformatics

BAMscale: quantification of next-generation sequencing peaks and generation of scaled coverage tracks

BEDTools: The Swiss‐Army Tool for Genome Feature Analysis

GMAP and GSNAP for Genomic Sequence Alignment: Enhancements to Speed, Accuracy, and Functionality

Cbrowse: A Sam/Bam-Based Contig Browser For Transcriptome Assembly Visualization And Analysis