agentsclimarketplace

Macs3 peak calling starter

Skill ma-compbio-lab/SkillFoundry/skills/epigenomics-and-chromatin/macs3-peak-calling-starter

A framework for discovering, compiling, and validating reusable skills for scientific agents.

Install
npx -y skills add ma-compbio-lab/SkillFoundry --skill macs3-peak-calling-starter

Assembled from the repository path, not quoted from the project. Check it against their README if it does not work.

SKILL.md

1.2 KB, 321 tokens by cl100k_base, as published. Nobody here has run it

MACS3 Peak Calling Starter

Use this skill to run a deterministic toy peak-calling example with MACS3 callpeak.

What it does

  • Runs macs3 callpeak on a small BED file using --nomodel --extsize 75 so toy data can be verified reliably.
  • Writes the standard narrowPeak, summits.bed, and peaks.xls outputs to a work directory.
  • Summarizes the strongest peak into a compact JSON payload for tests and demos.

When to use it

  • You need a runnable starter for simple ChIP-seq or chromatin-style peak calling.
  • You want a deterministic local example before moving to real replicate-aware workflows.
  • You need a smoke-testable wrapper around MACS3.

Example

slurm/envs/genomics/bin/python skills/epigenomics-and-chromatin/macs3-peak-calling-starter/scripts/run_macs3_peak_calling.py \
  --treatment skills/epigenomics-and-chromatin/macs3-peak-calling-starter/examples/toy_treatment.bed \
  --summary-out scratch/epigenomics/macs3_peak_summary.json

Verification

  • Skill-local tests: python3 -m unittest discover -s skills/epigenomics-and-chromatin/macs3-peak-calling-starter/tests -p 'test_*.py'
  • Expected summary: peak_count == 1 and top_peak["name"] == "toy_peak_1"

Keep looking

Skills are one crate of 328,083. Ordering is by how many stacks a row turns up in, so the top of any crate is what has actually been picked rather than what has the most stars.