bamsignals

bamsignals extracts read count vectors and computes read and coverage profiles from indexed BAM files, including paired-end data, to quantify sequencing read distributions for genomic analyses such as gene expression, variant detection, and epigenetic studies.

Key Features:

• Count Vectors from Indexed BAM Files: Generates count vectors by counting reads within specified genomic ranges from indexed BAM files.
• Reads and Coverage Profiles: Computes per-region read and coverage profiles that describe sequencing depth and read distribution across the genome.
• Handling Paired-End Data: Supports processing of paired-end sequencing data and accounts for paired reads in count and coverage calculations.

Scientific Applications:

• Gene Expression Analysis: Quantifies gene- or region-level expression by counting reads overlapping specified genomic ranges from BAM files.
• Variant Calling and Genomic Studies: Provides coverage profiles that aid identification of regions with altered read depth to support variant detection and other genomic analyses.
• Epigenetic Research: Uses coverage and read distribution profiles to inform analyses of chromatin accessibility and methylation-associated patterns.

Methodology:

Operates within the Bioconductor framework and leverages R's statistical capabilities to process indexed BAM files and compute count vectors and coverage profiles.