Skip to content
alignment-level-QC logo

alignment-level-QC

Calculates technical mapping statistics for any aligned BAM file (ChIP or ATAC). It assesses the performance of the aligner itself by generating metrics on read depth, mapping quality, error rates, and read group data using samtools and Picard.Use this skill to check "how well the reads mapped" o...

SKILL.md

Full skill instructions

Alignment Quality Control for ChIP-seq/​ATAC-seq

Overview

Perform comprehensive preliminary alignment-level quality control for ChIP-seq and ATAC-seq BAM files using samtools, Picard, and MultiQC.

Main steps include:

  • Initialize the project directory.
  • 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.
  • Sort and add read groups if missing in the BAM file.
  • Run preliminary QC metrics
  • Generate MultiQC report

When to use this skill

  • Use skill when you want to perform alignment-level quality control for ChIP-seq or ATAC-seq BAM files.

Inputs & Outputs

Inputs

<sample>.bam

Outputs

alignment_qc/
  ${sample}.bam                        # Original input
  ${sample}.sorted.bam                 # (Optional) Created if sorting was needed
  ${sample}.RG.bam                     # (Optional) Created if RG was needed
  ${sample}.RG.bam.bai                 # Index file
  qc_results/
    ${sample}.flagstat.txt
    ${sample}.stats.txt
    ${sample}.insertsize_metrics.txt
    ${sample}.dup_metrics.txt
    alignment_qc_report.html      # Visual MultiQC report
    qc_summary.txt                # Pass/​Warn/​Fail table
  temp/
    ${sample}.markdup.bam            # Intermediate file (safe to delete later)
    ...

Decision Tree

Step 0: Initialize Project

Call:

  • mcp__project-init-tools__project_init

with:

  • sample: all
  • task: alignment_qc

The tool will:

  • Create ${sample}_alignment_qc directory.
  • Return the full path of the ${sample}_alignment_qc directory, which will be used as ${proj_dir}.

Step 1: Check and Fix BAMs

  • Ensure all BAM files are coordinate-sorted, have Read Groups, and are indexed. This tool will skip files that are already correct and only create temporary files when fixes are needed.

Call:

  • mcp__qc-tools__check_and_fix_bams

with:

  • bam_files: List of BAM files to process.
  • temp_dir: ${proj_dir}/​temp

Step 2: Run Alignment QC Metrics

Call:

  • mcp__qc-tools__run_bam_qc

with:

  • bam_files: List of BAM files to process.
  • qc_dir: ${proj_dir}/​qc_results
  • temp_dir: ${proj_dir}/​temp

Step 3: Generate Summary Report

Call:

  • mcp__qc-tools__generate_qc_report

with:

  • qc_dir: ${proj_dir}/​qc_results

Quality Assessment

Key QC Metrics

  • Total reads – overall sequencing depth
  • Mapped reads (%) – alignment efficiency
  • Properly paired (%) – valid pair fraction (paired-end)
  • Duplicate rate (%) – PCR duplication estimate
  • Mitochondrial reads (%) – mitochondrial contamination
  • Insert size distribution – fragment length profile

All metrics are derived from samtools/Picard and summarized by MultiQC.


Quality Thresholds

CategoryCriteriaInterpretation
PassAll metrics within recommended thresholdsSuitable for downstream analysis
WarnOne or more borderline metricsLikely acceptable; review recommended
FailCritical metrics outside acceptable rangesRe-sequencing or reprocessing suggested

Report Generation

After MultiQC completes, generate a sample-wise summary (PASS/​WARN/​FAIL) per thresholds in references/​qc_metrics.md and save it as:

qc_results/​qc_summary.txt

Resources

Use references/​qc_metrics.md for:

  • Metric definitions and recommended thresholds
  • Troubleshooting guidance
  • Readiness criteria for peak calling
  • Pointers to ENCODE/​nf-core QC standards

More skills from diegosouzapw

accessibility logo
diegosouzapw/awesome-omni-skill

accessibility

Accessibility guidelines for VS Code features — covers accessibility help dialogs, accessible views, verbosity settings, accessibility signals, ARIA alerts/status announcements, keyboard navigation, and ARIA labels/roles. Applies to both new interactive UI surfaces and updates to existing feature...

62 1
View
mermaid-diagrams logo
diegosouzapw/awesome-omni-skill

mermaid-diagrams

Mermaid diagram creation for flowcharts, sequences, ERDs, and more. Generate diagrams from text in markdown files. Use for documentation, architecture diagrams, and visual representations. Triggers on mermaid, flowchart, sequence diagram, ERD, entity relationship, gantt chart, pie chart, class di...

62 1
View
tauri2-react-rust logo
diegosouzapw/awesome-omni-skill

tauri2-react-rust

Guides development of cross-platform desktop apps with Tauri 2, TypeScript, React, and Rust. Use when building Tauri apps, implementing IPC, designing Rust backend or TypeScript/React frontend, when researching or cloning a website (open site, snapshot elements), when verifying local dev or built...

62 1
View
fp-ts-errors logo
diegosouzapw/awesome-omni-skill

fp-ts-errors

Handle errors as values using fp-ts Either and TaskEither for cleaner, more predictable TypeScript code. Use when implementing error handling patterns with fp-ts.

62 0
View
ring:pre-dev-research logo
diegosouzapw/awesome-omni-skill

ring:pre-dev-research

Gate 0 research phase for pre-dev workflow. Dispatches 3 parallel research agents to gather codebase patterns, external best practices, and framework documentation BEFORE creating PRD/TRD. Outputs research.md with file:line references.

62 0
View
create-assistant logo
diegosouzapw/awesome-omni-skill

create-assistant

Create and configure Vapi voice AI assistants with models, voices, transcribers, tools, hooks, and advanced settings. Use when building voice agents, phone bots, customer support assistants, or any conversational AI that handles phone or web calls.

62 0
View
codex-headless logo
diegosouzapw/awesome-omni-skill

codex-headless

Delegiere Aufgaben an OpenAI Codex CLI im Headless-Mode. Nutzt ChatGPT Subscription (KEIN API Key). Nur Workspace-Zugriff, KEIN System-Zugriff. Auto-Accept mit --full-auto. Web-Recherche mit --search.

62 0
View
social-media-scheduler logo
diegosouzapw/awesome-omni-skill

social-media-scheduler

Generate a full week of social media content for any topic. Outputs platform-optimized posts for Twitter/X, LinkedIn, and Instagram with hashtags and posting times.

62 0
View
Data Privacy Compliance logo
diegosouzapw/awesome-omni-skill

Data Privacy Compliance

Data privacy and regulatory compliance specialist for GDPR, CCPA, HIPAA, and international data protection laws. Use when implementing privacy controls, conducting data protection impact assessments, ensuring regulatory compliance, or managing data subject rights. Expert in consent management, da...

62 0
View
triple-review logo
diegosouzapw/awesome-omni-skill

triple-review

Reviews a task using the TRIPLE protocol, as the review agent.

62 0
View
Tech Stack & Architecture Decision logo
diegosouzapw/awesome-omni-skill

Tech Stack & Architecture Decision

Define and document the technology stack and architecture decisions for a project. Use when the user needs to choose a tech stack, make architecture decisions, define infrastructure choices, or document technology selections. Triggers on requests like "define the tech stack", "choose technologies...

62 0
View

Popular AI tools

Kaiber logo
Video

Kaiber

Generate, edit, and beat-sync AI video with leading models in one workspace.

Paid
View
Vimcal logo
Productivity

Vimcal

The world's fastest calendar for remote work

Free
View

Transform Your Design with AI Designer by ImgCreator.ai

Freemium
View
Akool AI logo
Content & writing

Akool AI

Revolutionizing Video Production with AI-Powered Creativity

Paid
View

Extend an image past the frame and let AI fill the new aspect ratio.

Freemium
View
StarByFace logo
Security

StarByFace

Discover your celebrity doppelgänger with StarByFace!

Free
View
C

ChainClarity explains 700+ crypto whitepapers in plain English, with layered summaries, comparisons, research tools, alerts, and a $4.99 Pro plan.

Freemium
View
Opus Clip logo
Coding & apps

Opus Clip

Opus.ai: Revolutionize Your Web Experience

Free
View