tooluniverse-fastq-qc

Run FastQC quality control on FASTQ reads and decide adapter or quality trimming with fastp.

1.7k|254|Updated Mar 3, 2025
One-click install
npx skills add https://github.com/mims-harvard/ToolUniverse --skill tooluniverse-fastq-qc
Or copy as Structured Prompt for Agent
Please help me install this Agent Skill.
Skill: tooluniverse-fastq-qc
Source: https://github.com/mims-harvard/ToolUniverse/tree/main/plugins/tooluniverse/skills/tooluniverse-fastq-qc
Command: npx skills add https://github.com/mims-harvard/ToolUniverse --skill tooluniverse-fastq-qc

SYSTEM DOCUMENTATION & REQUIREMENTS

💡 This Skill includes scripts (resource) and references (resource) components.

What problem does it solve?

Raw sequencing reads often contain adapters, low-quality bases, and artifacts that silently corrupt downstream analysis. This Skill runs real local QC tools (FastQC, MultiQC, fastp, Cutadapt, seqkit) on FASTQ files, interprets every FastQC module flag, and produces an evidence-based trim-or-proceed decision without ever fabricating QC numbers or overwriting raw reads.

Core Features & Use Cases

  • Automated QC orchestration: A bundled script preflights tool availability, runs FastQC and seqkit stats into an isolated workdir, and emits a conda/mamba install plan if required binaries are missing.
  • FastQC report interpretation: Maps each module (per-base quality, adapter content, overrepresented sequences, duplication, GC content, N content) to its meaning and a concrete action for your library type.
  • Opt-in trimming: Trimming with fastp or Cutadapt only happens on explicit request, writes new files to a separate workdir, and is verified by re-running FastQC on trimmed output.
  • Use Case: You receive paired-end RNA-seq reads and ask whether they need trimming. The Skill runs FastQC, finds a 3' adapter ramp and quality tail decay, recommends fastp trimming, executes it into a new directory, and confirms the fix with a post-trim report.

Quick Start

Run QC on my paired-end reads sample_R1.fastq.gz and sample_R2.fastq.gz, interpret the FastQC report, and tell me whether I need to trim adapters before alignment.

Frequently Asked Questions about tooluniverse-fastq-qc

High-intent search queries and answers about installing and using this skill.

FAQPage Schema
How do I run FastQC on paired-end FASTQ files?

Run the bundled script with both read files: python scripts/run_fastq_qc.py --fastq sample_R1.fastq.gz sample_R2.fastq.gz --workdir /tmp/fastq_qc_run. It preflights FastQC, writes reports into the workdir, and leaves raw reads untouched.

fastp vs Cutadapt for adapter trimming, which should I use?

Use fastp for general Illumina QC and trimming since it auto-detects adapters, especially for paired-end via read overlap. Use Cutadapt when you know the exact adapter or primer sequence, such as amplicon panels or custom library prep.

What does a FastQC FAIL on per base sequence quality mean?

It usually means quality drops at the 3' end of reads, which is common especially on R2. The recommended action is gentle 3' quality trimming with fastp; a few low bases at the very end are often safe to leave since aligners soft-clip them.

Does this skill work without FastQC or fastp installed?

No. It requires real local binaries installed via bioconda. If a required tool is missing, the script prints a mamba or conda install plan and stops rather than fabricating QC results.

Should I trim reads before or after UMI extraction?

Extract UMIs first with umi_tools before any trimming. Naive adapter or quality trimming corrupts the UMI sequence and breaks downstream deduplication.

When should I not use FASTQ-level QC and trimming?

Do not use it for differential expression, read alignment, variant calling, or single-cell per-cell QC, which need dedicated tools. Also avoid trimming already-trimmed reads, and never deduplicate raw FASTQs since PCR duplicate removal belongs after alignment.