Does Pysam handle coordinate systems automatically?

Pysam uses 0-based, half-open coordinates (Python style) for internal operations, but can interpret 1-based region strings common in samtools convention.

What genomic file formats does Pysam support?

Pysam supports all industry-standard formats including SAM/BAM/CRAM for alignments, VCF/BCF for variants, and FASTA/FASTQ for sequences.

Can I run samtools commands directly with this skill?

Yes, the skill provides a direct interface to execute samtools and bcftools commands from within your Python scripts.

Do I need index files to use Pysam effectively?

For random access to specific genomic regions, index files (like .bai, .fai, or .tbi) are required. Without them, files must be processed sequentially.

Is Pysam suitable for large-scale NGS pipelines?

Yes, it is highly optimized for performance as it is built on htslib and supports parallel processing of independent genomic regions.

Pysam Genomic Data Processing

Name: Pysam Genomic Data Processing
Author: ovachiever

byovachiever

•

데이터 과학 및 ML

Processes genomic datasets including alignment, variant, and sequence files using a Pythonic interface to htslib.

Pysam is a specialized skill designed to empower developers and bioinformaticians working with high-throughput sequencing data. It provides a robust interface for reading, writing, and manipulating standard genomic formats like SAM, BAM, CRAM, VCF, and FASTA by wrapping the powerful htslib C-API. This skill is essential for building bioinformatics pipelines, performing variant annotation, and conducting quality control on sequencing results, offering features like region-based queries, pileup analysis for coverage, and seamless integration with samtools and bcftools commands directly within Python.

주요 기능

01Perform base-by-base pileup analysis for coverage calculations

02Manipulate VCF/BCF variant files and genotype data

03Read and write SAM/BAM/CRAM alignment files

04Extract genomic sequences from FASTA/FASTQ files

0515 GitHub stars

06Execute integrated samtools and bcftools commands

사용 사례

01Calculating read depth and coverage statistics across specific genomic regions

02Developing end-to-end bioinformatics pipelines for NGS data analysis

03Validating genetic variants against reference sequences and aligned reads

What are Skills?·How to Install

Install with 🐟 Skill.Fish

npx skillfish add ovachiever/droid-tings pysam

For use in Claude.ai and ChatGPT

Download Skill