pysam - Python module for reading and writing genomic datasets
Reads, writes, and manipulates SAM/BAM/CRAM alignments, VCF/BCF variants, and FASTA/FASTQ genomic sequences using a Pythonic interface to htslib.
Tags
Updated: 2026-09-21Capabilities
Typical Inputs
Typical Outputs
What this skill does
- Read and write alignment files
- Perform pileup coverage analysis
- Read and write variant files
- Query variants in genomic regions
- Extract reference sequences by coordinates
- Execute samtools and bcftools commands
Inputs
- Alignment files
- Variant files
- Sequence files
- Index files
- Genomic coordinates
Outputs
- Alignment files
- Variant files
- Extracted genomic sequences
- Coverage statistics
- Console output
Requirements
- Python environment
- pysam package installation
- Index files for random access
