ATACseq-QC

Performs ATAC-specific biological validation. It calculates metrics unique to chromatin accessibility assays, such as TSS enrichment scores and fragment size distributions (nucleosome banding patterns). Use this skill when you have filtered BAM file and have called peak for the file. Do NOT use this skill for ChIP-seq data or general alignment statistics.

16 stars

Best use case

ATACseq-QC is best used when you need a repeatable AI agent workflow instead of a one-off prompt.

Performs ATAC-specific biological validation. It calculates metrics unique to chromatin accessibility assays, such as TSS enrichment scores and fragment size distributions (nucleosome banding patterns). Use this skill when you have filtered BAM file and have called peak for the file. Do NOT use this skill for ChIP-seq data or general alignment statistics.

Teams using ATACseq-QC should expect a more consistent output, faster repeated execution, less prompt rewriting.

When to use this skill

  • You want a reusable workflow that can be run more than once with consistent structure.

When not to use this skill

  • You only need a quick one-off answer and do not need a reusable workflow.
  • You cannot install or maintain the underlying files, dependencies, or repository context.

Installation

Claude Code / Cursor / Codex

$curl -o ~/.claude/skills/atacseq-qc/SKILL.md --create-dirs "https://raw.githubusercontent.com/diegosouzapw/awesome-omni-skill/main/skills/data-ai/atacseq-qc/SKILL.md"

Manual Installation

  1. Download SKILL.md from GitHub
  2. Place it in .claude/skills/atacseq-qc/SKILL.md inside your project
  3. Restart your AI agent — it will auto-discover the skill

How ATACseq-QC Compares

Feature / AgentATACseq-QCStandard Approach
Platform SupportNot specifiedLimited / Varies
Context Awareness High Baseline
Installation ComplexityUnknownN/A

Frequently Asked Questions

What does this skill do?

Performs ATAC-specific biological validation. It calculates metrics unique to chromatin accessibility assays, such as TSS enrichment scores and fragment size distributions (nucleosome banding patterns). Use this skill when you have filtered BAM file and have called peak for the file. Do NOT use this skill for ChIP-seq data or general alignment statistics.

Where can I find the source code?

You can find the source code on GitHub using the link provided at the top of the page.

SKILL.md Source

# ATAC-seq Quality Control

## Overview

This skill performs complete ATAC-seq data quality control from BAM and peak files.

Main steps include:
- Refer to the **Inputs & Outputs** section to check inputs and build the output architecture. All the output file should located in `${proj_dir}` in Step 0.
- **Always prompt user** for genome assembly used. Never decide by yourself. 
- Generate TSS files according to genome assembly.
- Compute TSS enrichment, fragment distribution and FRiP. 

---

## Inputs & Outputs

### Inputs

```bash
${sample}.bam # filtered bam files
${sample}.narrowPeak
```

### Outputs

```bash
all_atac_qc/
    ${sample}_qc_results/
        ataqv_metrics.json
        ataqv_report.html/
    temp/
```

---

## Decision Tree

### Step 0: Initialize Project

Call:

- `mcp__project-init-tools__project_init`

with:

- `sample`: all
- `task`: atac_qc
- `genome`: provided by user

The tool will:

- Create`all_atac_qc` directory.
- Return the full path of the `all_atac_qc` directory, which will be used as `${proj_dir}`.

### Step 1: Detect the name logic of the chromosomes in BAM file (have "chr" as prefix or not)

`samtools view <sample>.bam | head -n 10 | cut -f 3`

### Step 2: Generate reference files

Call:
- mcp__qc-tools__generate_reference

with:
- `genome`: Genome name (e.g., hg38), provided by user
- `temp_dir`: ${proj_dir}/temp
- `bam_uses_chr`: True if BAM uses 'chr' prefix (chr1), False if not (1).

### Step 3: Peform quality control for the ATAC-seq data

Call:
- mcp__qc-tools__run_ataqv_qc

- `bam_file`: Path to filtered BAM file
- `peak_file`: Path to peak file (narrowPeak) corresponding to the BAM file
- `tss_file`: ${proj_dir}/temp/${genome}.tss
- `species`: Species used, choose from (fly, human, mouse, rat, worm, yeast)
- `bam_uses_chr`: True if BAM uses 'chr' prefix (chr1), False if not (1).
- `output_dir`: ${proj_dir}/${sample}_qc_results
- `autosomal_ref_path`: Provided if `bam_uses_chr` is False, ${proj_dir}/temp/${genome}.autosomal.ref

Related Skills

bio-workflows-atacseq-pipeline

16
from diegosouzapw/awesome-omni-skill

End-to-end ATAC-seq workflow from FASTQ files to differential accessibility and TF footprinting. Covers alignment, peak calling with MACS3, QC metrics, and optional TOBIAS footprinting. Use when running end-to-end ATAC-seq analysis from FASTQ to differential accessibility.

bgo

10
from diegosouzapw/awesome-omni-skill

Automates the complete Blender build-go workflow, from building and packaging your extension/add-on to removing old versions, installing, enabling, and launching Blender for quick testing and iteration.

Coding & Development

large-data-with-dask

16
from diegosouzapw/awesome-omni-skill

Specific optimization strategies for Python scripts working with larger-than-memory datasets via Dask.

langsmith-fetch

16
from diegosouzapw/awesome-omni-skill

Debug LangChain and LangGraph agents by fetching execution traces from LangSmith Studio. Use when debugging agent behavior, investigating errors, analyzing tool calls, checking memory operations, or examining agent performance. Automatically fetches recent traces and analyzes execution patterns. Requires langsmith-fetch CLI installed.

langchain-tool-calling

16
from diegosouzapw/awesome-omni-skill

How chat models call tools - includes bind_tools, tool choice strategies, parallel tool calling, and tool message handling

langchain-notes

16
from diegosouzapw/awesome-omni-skill

LangChain 框架学习笔记 - 快速查找概念、代码示例和最佳实践。包含 Core components、Middleware、Advanced usage、Multi-agent patterns、RAG retrieval、Long-term memory 等主题。当用户询问 LangChain、Agent、RAG、向量存储、工具使用、记忆系统时使用此 Skill。

langchain-js

16
from diegosouzapw/awesome-omni-skill

Builds LLM-powered applications with LangChain.js for chat, agents, and RAG. Use when creating AI applications with chains, memory, tools, and retrieval-augmented generation in JavaScript.

langchain-agents

16
from diegosouzapw/awesome-omni-skill

Expert guidance for building LangChain agents with proper tool binding, memory, and configuration. Use when creating agents, configuring models, or setting up tool integrations in LangConfig.

lang-python

16
from diegosouzapw/awesome-omni-skill

Python 3.13+ development specialist covering FastAPI, Django, async patterns, data science, testing with pytest, and modern Python features. Use when developing Python APIs, web applications, data pipelines, or writing tests.

kramme:agents-md

16
from diegosouzapw/awesome-omni-skill

This skill should be used when the user asks to "update AGENTS.md", "add to AGENTS.md", "maintain agent docs", or needs to add guidelines to agent instructions. Guides discovery of local skills and enforces structured, keyword-based documentation style.

kontent-ai-automation

16
from diegosouzapw/awesome-omni-skill

Automate Kontent AI tasks via Rube MCP (Composio). Always search tools first for current schemas.

kitt-create-slash-commands

16
from diegosouzapw/awesome-omni-skill

Expert guidance for creating slash commands. Use when working with slash commands, creating custom commands, understanding command structure, or learning YAML configuration.