Communitygithub.com

pacsomatic

Operator-Werkzeugkasten für nf-core/pacsomatic Matched-Tumor-Normal-Workflows ab BAM-Eingabedateien. Dieses Skill wird verwendet, wenn der Benutzer die Validierung von Run-Eingaben, die Erstellung pacsomatic-konformer Samplesheets, die Vorbereitung reproduzierbarer Nextflow-Launch-Artefakte, die lokale Ausführung oder die Einreichung an Scheduler (LSF/Slurm/PBS/SGE) sowie die Fehlerbehandlung bei Ausführungsfehlern benötigt. Es wird bei Anfragen zum Ausführen von pacsomatic, zur Erstellung von Launch-Befehlen/Skripten, zur Durchführung von Dry-Run-Prüfungen oder zur Fehlerbehebung bei Pipeline-Start- und Scheduler-Einreichungsfehlern ausgelöst.

Was ist pacsomatic?

Der Skill 'pacsomatic' ist ein Werkzeugkasten für die reproduzierbare Ausführung des nf-core/pacsomatic-Workflows für Matched-Tumor-Normal-Analysen, ausgehend von BAM-Dateien. Er zentralisiert die Verarbeitung über ein einziges Einstiegsskript (`scripts/run_pacsomatic.py`), das die Validierung der Eingabedaten, die Erstellung von Samplesheets und Parameter-YAML-Dateien sowie die Generierung von Startskripten übernimmt. Der Skill unterstützt sowohl die lokale Ausführung als auch die Einreichung an gängige Scheduler (LSF, Slurm, PBS, SGE) und bietet einen `--dry-run`-Modus für die vorherige Validierung. Er ist ideal für DevOps-orientierte Bioinformatiker und Agenten, die konsistente, rückverfolgbare Läufe ohne manuelle Nextflow-Konfiguration benötigen. Der Skill stellt sicher, dass alle Validierungsschritte durchlaufen werden, bevor die tatsächliche Ausführung startet, und gibt bei Fehlern konkrete nächste Schritte zur Fehlerbehebung an.

Funktioniert mit~Claude Code~Codex CLI~Cursor
npx skills add https://github.com/K-Dense-AI/scientific-agent-skills/tree/main/skills/pacsomatic

In Ihrer bevorzugten KI fragen

Öffnet einen neuen Chat, in dem dieser Agent-Skill bereits geladen ist.

Dokumentation

Was macht pacsomatic?

Übersicht

Dieser Skill bietet einen reproduzierbaren Ausführungsworkflow für nf-core/pacsomatic, der auf einem einzigen Helfer-Einstiegspunkt basiert, der die Validierung, Artefaktgenerierung und optionale Ausführung übernimmt.

Primärer Einstiegspunkt:

  • scripts/run_pacsomatic.py

Das Hilfsskript:

  • validiert erforderliche IDs, Dateien, Referenzmodus und Laufzeitvoraussetzungen
  • erstellt ein pacsomatic-kompatibles Samplesheet (patient,sample,status,bam,pbi)
  • generiert eine Parameter-YAML und ein Startskript für reproduzierbare Wiederholungsläufe
  • unterstützt Dry-Run-Validierung und Run/Submit-Ausführungspfade

Verwenden Sie diesen Skill als Standardpfad für pacsomatic-Operationen. Umgehen Sie ihn nicht mit manuell zusammengestellten nextflow run nf-core/pacsomatic-Befehlen, es sei denn, der Benutzer fordert explizit die manuelle Befehlserstellung an.

Wann dieser Skill verwendet werden sollte

Rufen Sie diesen Skill auf, wenn der Benutzer Folgendes verlangt:

  • Durchführung einer Matched-Tumor-Normal-Analyse ab BAM-Dateien
  • Erstellung oder Korrektur eines pacsomatic-Samplesheets und der Start-Artefakte
  • lokale Ausführung oder Einreichung an Scheduler (LSF/Slurm/PBS/SGE)
  • Durchführung einer Dry-Run-Validierung vor der Ausführung
  • Fehlerbehebung bei Startfehlern oder Zusammenfassung von Run-Ergebnissen

Verwenden Sie diesen Skill nicht für:

  • tiefgehende biologische Interpretation jenseits von Run-Level-Plausibilitätsprüfungen
  • Bearbeitung von Pipeline-Interna, es sei denn, dies wird explizit angefordert

Typische Auslösephrasen:

  • "führe nf-core/pacsomatic für dieses Tumor-Normal-Paar aus"
  • "erstelle pacsomatic-Samplesheet und Startskript"
  • "mach zuerst einen Dry Run und sag mir, was fehlt"
  • "reiche pacsomatic an slurm/lsf ein und gib die Job-ID zurück"
  • "warum ist die pacsomatic-Einreichung fehlgeschlagen"

Routing- und Ausführungsregeln

  1. Sammeln Sie immer zuerst die erforderlichen Run-Eingaben.
  2. Leiten Sie immer über scripts/run_pacsomatic.py für Validierung und Artefaktgenerierung.
  3. Verwenden Sie standardmäßig --dry-run, wenn der Benutzer nur Prüfungen/Validierung wünscht.
  4. Verwenden Sie --run nur, wenn der Benutzer die Ausführung/Einreichung wünscht.
  5. Fügen Sie für Scheduler-Modi ausführungsspezifische Ressourcenargumente hinzu und geben Sie die erkannte Job-ID zurück, wenn verfügbar.
  6. Melden Sie bei Ausführungsfehlern den ersten Fehlerpunkt und das nächste Ziel für die Fehlerbehebung (.nextflow.log, pipeline_info, fehlgeschlagene Task-Logs).

Erforderliche Eingaben

Erforderlich:

  • Tumor-BAM-Pfad
  • Normal-BAM-Pfad
  • Patienten-ID
  • Tumor-Proben-ID
  • Normal-Proben-ID
  • Ausgabeverzeichnis
  • genau ein Referenzmodus: --fasta oder --genome

Optional:

  • Profil, Ressourcen, Scheduler-Konto/Warteschlange
  • Pipeline-Version (-r)
  • Parameter-Datei, Resume/Report/DAG-Flags
  • --dry-run und/oder --run

Arbeitsablauf

  1. Identität und Eingabebeschränkungen validieren.
  2. Erforderliche lokale Pfade validieren (BAM, optionales PBI, optionales FASTA).
  3. Laufzeit- und Abhängigkeitsprüfungen durchführen.
  4. Samplesheet und generierte Parameter-YAML erstellen.
  5. Startskript für den ausgewählten Executor generieren.
  6. Wenn --dry-run und nicht --run, nach der Artefaktgenerierung anhalten.
  7. Wenn --run, lokal ausführen oder an Scheduler übergeben.
  8. Befehl/Skript-Pfad, Validierungsstatus und Job-ID (falls erkannt) zurückgeben.

Agenten-Antwortvertrag

Jede Antwort nach dem Aufruf sollte Folgendes enthalten:

  • den exakten verwendeten Befehl oder den Pfad des generierten Skripts
  • Bestätigung, dass die Validierungsprüfungen durchgeführt wurden
  • Run-Typ (dry-run vs run)
  • Scheduler-Job-ID, wenn verfügbar
  • einen konkreten nächsten Schritt zur Validierung/Fehlerbehebung

Kurzanleitung

Dry Run:

python scripts/run_pacsomatic.py \
  --tumor-bam /path/to/tumor.bam \
  --normal-bam /path/to/normal.bam \
  --patient-id P001 \
  --tumor-sample-id P001_T \
  --normal-sample-id P001_N \
  --outdir /path/to/output \
  --genome GRCh38 \
  --profile singularity,sanger \
  --dry-run

Beispiel für Scheduler-Ausführung (Slurm):

python scripts/run_pacsomatic.py \
  --tumor-bam /path/to/tumor.bam \
  --normal-bam /path/to/normal.bam \
  --patient-id P001 \
  --tumor-sample-id P001_T \
  --normal-sample-id P001_N \
  --outdir /path/to/output \
  --genome GRCh38 \
  --profile singularity,sanger \
  --executor slurm \
  --queue compute \
  --project my_account \
  --cpus 16 \
  --memory-gb 64 \
  --walltime 48:00 \
  --run

Konfiguration

Verwenden Sie config.yaml als Grundlage für Profil-/Executor-/Laufzeit-Standardeinstellungen. Überschreiben Sie diese beim Aufruf, wenn die Benutzeranforderungen abweichen.

Tests

Führen Sie Unit-Tests aus dem Skill-Stammverzeichnis aus:

python -m unittest discover -s tests/pacsomatic -v

Referenzen

  • references/agent-playbook.md
  • references/config-and-output.md
  • references/pacsomatic_guide.md
  • scripts/run_pacsomatic.py

Individual skills in this repo

This repo contains 20 individual skills — each has its own dedicated page.

adaptyv

How to use the Adaptyv Bio Foundry API and Python SDK for protein experiment design, submission, and results retrieval. Use this skill whenever the user mentions Adaptyv, Foundry API, protein binding assays, protein screening experiments, BLI/SPR assays, thermostability assays, or wants to submit protein sequences for experimental characterization. Also trigger when code imports `adaptyv`, `adaptyv_sdk`, or `FoundryClient`, or references `foundry-api-public.adaptyvbio.com`.

aeon

This skill should be used for time series machine learning tasks including classification, regression, clustering, forecasting, anomaly detection, segmentation, and similarity search. Use when working with temporal data, sequential patterns, or time-indexed observations requiring specialized algorithms beyond standard ML approaches. Particularly suited for univariate and multivariate time series analysis with scikit-learn compatible APIs.

alphagenome

Look up precomputed AlphaGenome Atlas effects for any GRCh38 single-nucleotide variant (AVI score with Phred and 18 SHAP feature attributions, plus raw and quantile scores for RNA-seq, DNase, ATAC, ChIP-TF, ChIP-histone, CAGE, PRO-cap, splicing, polyadenylation and contact-map tracks), score variants or scan windows on demand with the AlphaGenome model for human and mouse (variant scoring, in silico mutagenesis, REF-versus-ALT track prediction), and build Atlas website deep links. Use when the user mentions AlphaGenome, AlphaGenome Atlas, AVI or AlphaGenome Variant Impact, DeepMind variant effect prediction, or wants to prioritise or mechanistically interpret non-coding, regulatory, splicing, enhancer, promoter, or chromatin-accessibility effects of SNVs from a VCF, credible set, or region. Research use only; not a clinical tool.

analytical-method-validation

Plan, execute, and document validation, verification, and transfer of analytical procedures under the governing framework - ICH Q2(R2) and Q14, USP <1220>/<1225>/<1226>, ICH M10 bioanalytical, CLSI EP, or ISO/IEC 17025. Use for HPLC, LC-MS/MS, GC, CE, ICP-MS, dissolution, qNMR, qPCR, NIR, and ligand binding or cell-based assays whenever the question is whether a procedure is fit for its intended purpose. Triggers include

anndata

Data structure for annotated matrices in single-cell analysis. Use when working with .h5ad files or integrating with the scverse ecosystem. This is the data format skill—for analysis workflows use scanpy; for probabilistic models use scvi-tools; for population-scale queries use cellxgene-census.

arbor

Autonomously improve a real artifact (code, training recipe, agent harness, data pipeline, prompt) against an objective and an evaluator, using Hypothesis Tree Refinement (HTR) from the Arbor paper. Use this whenever someone wants to iteratively optimize something over many experiments without overfitting — e.g.

arboreto

Infer gene regulatory networks (GRNs) from gene expression data using scalable algorithms (GRNBoost2, GENIE3). Use when analyzing transcriptomics data (bulk RNA-seq, single-cell RNA-seq) to identify transcription factor-target gene relationships and regulatory interactions. Supports distributed computation for large-scale datasets.

astropy

Core Python library for astronomy and astrophysics workflows that need Astropy APIs, including units/quantities, coordinates, FITS I/O, tables, time systems, WCS, and cosmology. Use when implementing or debugging astronomical data analysis code with Astropy.

autoskill

Observe the user

benchling-integration

Benchling Python SDK and REST API integration for registry entities, inventory, ELN entries, workflows, Benchling Apps, and Data Warehouse queries. Use when automating lab data with benchling-sdk or the v2 API.

bgpt-paper-search

Search scientific papers and retrieve structured experimental data extracted from full-text studies via the BGPT MCP server. Returns 25+ fields per paper including methods, results, sample sizes, quality scores, and conclusions. Use for literature reviews, evidence synthesis, and finding experimental details not available in abstracts alone.

bids

>

biopython

Comprehensive molecular biology toolkit. Use for sequence manipulation, file parsing (FASTA/GenBank/PDB), phylogenetics, and programmatic NCBI/PubMed access (Bio.Entrez). Best for batch processing, custom bioinformatics pipelines, BLAST automation. For quick lookups use gget; for multi-service integration use bioservices.

bioservices

Unified Python interface to 40+ bioinformatics services. Use when querying multiple databases (UniProt, KEGG, ChEMBL, Reactome) in a single workflow with consistent API. Best for cross-database analysis, ID mapping across services. For quick single-database lookups use gget; for sequence/file manipulation use biopython.

bulk-rnaseq

End-to-end bulk RNA-seq orchestrator — takes raw FASTQ reads through QC and trimming (FastQC, fastp/Trim Galore), alignment and quantification (STAR, Salmon, featureCounts), assembles a gene-level counts matrix, then hands off to differential expression (pydeseq2), pathway/GSEA enrichment (pathway-enrichment), and publication figures (scientific-visualization). Use whenever the user has bulk RNA-seq reads or quant output and wants a complete, reproducible differential-expression workflow — e.g.

cellxgene-census

Query the CZ CELLxGENE Census programmatically for versioned public single-cell and spatial transcriptomics data. Use when you need population-scale cell metadata, gene expression slices, Census summary counts, source H5AD URIs/downloads, embeddings, spatial Census data, or reference atlas comparisons across organisms, tissues, diseases, assays, and cell types. For analyzing your own local single-cell data use scanpy, anndata, or scvi-tools.

cirq

Google quantum computing framework. Use when targeting Google Quantum AI hardware, designing noise-aware circuits, or running quantum characterization experiments. Best for Google hardware, noise modeling, and low-level circuit design. For IBM hardware use qiskit; for quantum ML with autodiff use pennylane; for physics simulations use qutip.

citation-management

Comprehensive citation management for academic research. Search OpenAlex, PubMed, and Google Scholar for papers, extract accurate metadata, validate citations, and generate properly formatted BibTeX entries. This skill should be used when you need to find papers, verify citation information, convert DOIs to BibTeX, or ensure reference accuracy in scientific writing.

clinical-decision-support

Prepare and validate research-only clinical decision-support evaluation, evidence-profile, cohort, survival, biomarker/model, privacy, and governance artifacts. Use for aggregate or synthetic research documentation and traceability—not patient care or live clinical operation.

clinical-reports

Create safety-bounded draft structures and run local deterministic checks for clinical case, diagnostic, trial, safety, and aggregate research reports. Use only with synthetic, de-identified, or aggregate inputs and verified source-fact manifests; every output requires qualified review.

Verwandte Skills