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 ).

Pybedgraph: a Python Package for Fast Operations on 1D Genomic Signal Tracks

Analysis and Visualization of ChIP-Seq and RNA-Seq Sequence Alignments Using Ngs.plot.

Megadepth: efficient coverage quantification for BigWigs and BAMs

ChIPbinner: An R package for analyzing broad histone marks binned in uniform windows from ChIP-Seq or CUT&RUN/TAG data

ChIPseeker: an R/Bioconductor package for ChIP peak annotation, comparison and visualization

PanDepth, an ultrafast and efficient genomic tool for coverage calculation

W-ChIPeaks: a Comprehensive Web Application Tool for Processing ChIP-chip and ChIP-seq Data

Pytometry: Flow and Mass Cytometry Analytics in Python

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

Computation for ChIP-seq and RNA-seq studies

BEDTools: The Swiss‐Army Tool for Genome Feature Analysis

Genotype Representation Graphs: Enabling Efficient Analysis of Biobank-Scale Data

An Integrated Software System for Analyzing ChIP-chip and ChIP-seq Data

ATACgraph: Profiling Genome-Wide Chromatin Accessibility From ATAC-seq

Advances in computational ChIA-PET data analysis

genomepy: genes and genomes at your fingertips

GraphBio: A shiny web app to easily perform popular visualization analysis for omics data

ChIAPoP: a new tool for ChIA-PET data analysis

QChIPat: a quantitative method to identify distinct binding patterns for two biological ChIP-seq samples in different experimental conditions

PCBS: an R package for fast and accurate analysis of bisulfite sequencing data

Pse-Analysis: a python package for DNA/RNA and protein/peptide sequence analysis based on pseudo components and kernel methods