mpileup

mpileup processes BAM files to filter positions by coverage and detect single nucleotide polymorphisms (SNPs) to support variant calling and population genetics analyses.

Key Features:

• BAM parsing: Parses BAM-formatted alignments to extract per-base and per-position information.
• Coverage filtering: Filters reads or genomic positions using user-defined coverage/read-depth thresholds to retain regions with sufficient depth.
• SNP identification: Analyzes nucleotide variation across aligned reads to identify single nucleotide polymorphisms (SNPs).
• Filtered output: Produces filtered per-position data suitable for downstream variant-calling and comparative analyses.

Scientific Applications:

• Variant calling: Provides filtered alignments and SNP information as input for variant-calling workflows.
• Population genetics: Enables analysis of allele frequencies and polymorphism patterns across samples for population-genetics studies.
• Metagenomics: Applies coverage filtering and SNP detection principles to characterize genetic variation in metagenomic datasets.

Methodology:

Parses BAM files to extract genomic information, applies user-defined criteria for coverage and SNP detection, and outputs filtered data for downstream analysis.