SomatoSim

SomatoSim simulates somatic single nucleotide variants (SNVs) within SAM/BAM files to generate controlled datasets for benchmarking and evaluating somatic variant calling algorithms.

Key Features:

• Targeted somatic SNV simulation: Simulates somatic single nucleotide variants (SNVs) directly within sequence alignment map (SAM/BAM) files.
• Customizable simulation parameters: Allows specification of exact variant positions, number of variants, variant allele fractions, depth of coverage, read quality, base quality, and other simulation parameters.
• Three-stage process: Performs Variant Selection, Variant Simulation, and Variant Evaluation as discrete computational stages.
• Read-level mutation: Selects reads from the input BAM and mutates bases to introduce desired variants at specified locations.
• BAM input/output integration: Accepts an analysis-ready BAM as input and outputs a BAM file containing the simulated variants.

Scientific Applications:

• Benchmarking variant callers: Generates controlled datasets to evaluate and compare the performance of somatic variant calling tools and algorithms.
• Pipeline development: Enables development and refinement of bioinformatics pipelines for somatic variant detection.
• Sequencing parameter studies: Facilitates studies of the impact of variant allele fraction, sequencing depth, read quality, and base quality on variant calling accuracy.

Methodology:

Performs a three-stage computational workflow consisting of Variant Selection based on user-defined criteria, Variant Simulation by selecting and mutating reads to introduce SNVs at specified positions, and Variant Evaluation that summarizes the simulation results, operating on an analysis-ready BAM to produce an output BAM with simulated variants.