Can I run samtools commands through this skill?

Yes, Pysam provides direct access to samtools and bcftools commands via Python wrappers, allowing you to sort, index, and view genomic files programmatically.

What coordinate system does Pysam use?

Pysam follows Python conventions by using a 0-based, half-open coordinate system for internal operations, though it supports 1-based strings for genomic region fetching.

Does Pysam require indexed files for data access?

While Pysam can read files sequentially, random access to specific genomic regions requires index files such as .bai for BAM, .tbi for VCF, and .fai for FASTA.

How does Pysam handle large genomic datasets?

Pysam is built on htslib and is optimized for performance, using iterators and index-based lookups to minimize memory overhead when processing massive sequencing files.

Pysam Genomic Toolkit

Name: Pysam Genomic Toolkit
Author: plurigrid

byplurigrid

•

데이터 과학 및 ML

Manipulates and analyzes genomic datasets including SAM, BAM, VCF, and FASTA files using a Pythonic interface to htslib.

The Pysam skill equips Claude with specialized knowledge for processing high-throughput sequencing (NGS) data. It provides a robust framework for reading, writing, and querying standard genomic formats like BAM/CRAM alignments, VCF/BCF variants, and FASTA/FASTQ sequences. By implementing this skill, developers can efficiently perform region-based data extraction, calculate genomic coverage via pileup analysis, and automate complex bioinformatics workflows. It bridges the gap between low-level htslib operations and high-level Python analysis, making it indispensable for variant calling, quality control, and sequence manipulation tasks.

주요 기능

018 GitHub stars

02Integrated wrappers for samtools and bcftools command-line utilities.

03Efficient region-based fetching and random access using genomic index files.

04Advanced pileup analysis for base-by-base coverage and read depth calculations.

05Support for 0-based coordinate systems and Pythonic handling of genomic intervals.

06Comprehensive read/write support for SAM, BAM, CRAM, VCF, BCF, FASTA, and FASTQ formats.

사용 사례

01Developing automated bioinformatics pipelines for NGS data processing and variant discovery.

02Performing quality control and coverage analysis on large-scale sequencing datasets.

03Extracting specific genomic regions and reference sequences for downstream statistical analysis.

What are Skills?·How to Install

Install with 🐟 Skill.Fish

npx skillfish add plurigrid/asi pysam

For use in Claude.ai and ChatGPT

Download Skill