pysam

pysam provides programmatic access to read, manipulate, and analyze sequencing alignment and variant files (SAM, BAM, VCF, BCF) for genomic analyses.

Key Features:

• Support for Multiple Formats: Read and write SAM, BAM, VCF, and BCF files, including alignments for short and long reads up to 128 Mbp.
• Efficient Random Access: Provide efficient random access and retrieval from large alignment and variant files for high-throughput analyses.
• Integration with SAMtools: Expose SAMtools functionality via Python, including indexing, variant calling, and alignment viewing.

Scientific Applications:

• Genome Alignment Analysis: Read and manipulate alignment files to study how sequencing reads map to reference genomes.
• Variant Calling and Analysis: Identify genetic variants from sequence data by integrating SAMtools variant-calling functions.
• Data Visualization and Exploration: Support alignment viewing for visual inspection and exploration of sequencing data.
• Population-scale Variant Analysis (1000 Genomes Project): Process and analyze large population genomic datasets such as the 1000 Genomes Project to identify genetic variation across populations.

Methodology:

Provide a Python interface to read and manipulate SAM, BAM, VCF, and BCF files, enabling scripted automation of genomic data processing.