Bio-samtools

Bio-samtools provides a Ruby interface to SAMtools for programmatic access, parsing, and analysis of high-throughput sequence alignments in SAM and BAM formats, including MPileup and VCF-based variant analyses and alignment visualization.

Key Features:

• SAM/BAM integration: Programmatic access to information stored in SAM and BAM files via the SAMtools library.
• Extensibility: Allows incorporation of additional SAMtools methods as new SAMtools releases become available.
• MPileup and VCF handling: Provides Ruby classes for processing MPileup outputs and Variant Call Format (VCF) data.
• Variant analysis: Implements alternative allele calculation and allele frequency calling for single nucleotide polymorphisms (SNPs).
• Full SAMtools support: Exposes the full set of functionalities provided by the SAMtools library.
• Alignment visualization: Includes methods to generate direct visual representations of alignment data.

Scientific Applications:

• Sequence alignment management: Processing and manipulation of high-throughput alignment data in SAM and BAM formats for downstream analyses.
• Variant discovery and genotyping: MPileup and VCF handling with alternative allele and allele frequency calculations for SNP detection and genotyping.
• Allele frequency estimation: Calculation of allele frequencies from MPileup and VCF data for population- and sample-level analyses.
• Alignment visualization and interpretation: Visual representation of alignments to support inspection and validation of mapping and variant calls.
• Genome-wide association studies (GWAS): Generation and preprocessing of variant calls and allele frequency data for GWAS workflows.
• Personalized medicine: Variant calling and allele frequency information to support genotype-based analyses relevant to personalized medicine.
• Evolutionary biology: Comparative variant and allele frequency analyses to support evolutionary inference.

Methodology:

Provides Ruby bindings to SAMtools to access and parse SAM/BAM files; defines Ruby classes for MPileup and VCF processing; implements methods for alternative allele calculation and allele frequency calling for SNPs; and includes methods for direct visualization of alignment data.