GPTomics/bio-reporting-automated-qc-reports
reporting/automated-qc-reports/SKILL.md
Generates standardized quality control reports by aggregating metrics from FastQC, alignment, and other tools using MultiQC. Use when summarizing QC metrics across samples, creating shareable quality reports, or building automated QC pipelines.
$ install --global
skillsauthnpx skillsauth add GPTomics/bioSkills bio-reporting-automated-qc-reportsInstall this skill globally with one command. Works with Claude Code, Cursor, and Windsurf.
Security Scan Results
3 of 9 scanners reported clean
Some scanners were skipped, did not run, or reported a non-clean status. Review each row below.
3
Scanners Passed
9
Scanners in report
Clean
TrivyContainer and dependency vulnerability scanner
95%
Clean
SemgrepStatic code analysis for vulnerabilities
95%
Clean
mcp-scan (Snyk)Model Context Protocol security validation
95%
Skipped
Snyk (dep)Open source security scanning
50%
Skipped
Socket.devSupply chain security analysis
50%
Skipped
VirusTotalMulti-engine malware detection
50%
Skipped
CrowdStrikeAdvanced threat intelligence
50%
Skipped
OSV-ScannerOpen Source Vulnerability database check
50%
Skipped
OWASP Dep-Check
50%
Last scanned: Jun 5, 2026, 2:14 AM155.6s3 files scanned
SKILL.md
- name:
- bio-reporting-automated-qc-reports
- description:
- Generates standardized quality control reports by aggregating metrics from FastQC, alignment, and other tools using MultiQC. Use when summarizing QC metrics across samples, creating shareable quality reports, or building automated QC pipelines.
- tool_type:
- cli
- primary_tool:
- multiqc
Version Compatibility
Reference examples tested with: Cell Ranger 8.0+, FastQC 0.12+, GATK 4.5+, HISAT2 2.2.1+, MultiQC 1.21+, STAR 2.7.11+, Subread 2.0+, bcftools 1.19+, fastp 0.23+, kallisto 0.50+
Before using code patterns, verify installed versions match. If versions differ:
- Python:
pip show <package>thenhelp(module.function)to check signatures - CLI:
<tool> --versionthen<tool> --helpto confirm flags
If code throws ImportError, AttributeError, or TypeError, introspect the installed package and adapt the example to match the actual API rather than retrying.
Automated QC Reports with MultiQC
"Aggregate QC results into one report" -> Combine outputs from FastQC, samtools, Picard, and other tools into a single interactive HTML report.
- CLI:
multiqc .(scans current directory for recognized tool outputs)
Basic Usage
# Aggregate all QC outputs in directory
multiqc results/ -o qc_report/
# Specify output name
multiqc results/ -n my_project_qc
# Include specific tools only
multiqc results/ --module fastqc --module star
Supported Tools
MultiQC recognizes outputs from 100+ bioinformatics tools:
| Category | Tools | |----------|-------| | Read QC | FastQC, fastp, Cutadapt | | Alignment | STAR, HISAT2, BWA, Bowtie2 | | Quantification | featureCounts, Salmon, kallisto | | Variant Calling | bcftools, GATK | | Single-cell | CellRanger, STARsolo |
Configuration
Create multiqc_config.yaml:
title: "RNA-seq QC Report"
subtitle: "Project XYZ"
intro_text: "QC metrics for all samples"
# Custom sample name cleaning
extra_fn_clean_exts:
- '.sorted'
- '.dedup'
# Report sections to include
module_order:
- fastqc
- star
- featurecounts
# Highlight samples
table_cond_formatting_rules:
pct_mapped:
fail: [{lt: 50}]
warn: [{lt: 70}]
Custom Data
# Add custom data file
# File format: sample\tmetric1\tmetric2
multiqc results/ --data-format tsv --custom-data-file custom_metrics.tsv
Python API
from multiqc import run as multiqc_run
# Run programmatically
multiqc_run(analysis_dir='results/', outdir='qc_report/')
Related Skills
- read-qc/quality-reports - Generate input FastQC reports
- read-qc/fastp-workflow - Preprocessing QC
- workflows/rnaseq-to-de - Full workflow with QC
Related Skills
GPTomics/phasing-imputation/foundations
tools
VerifiedTrustedCommunity
--- name: bio-phasing-imputation-foundations description: Frames the phasing/imputation pipeline before any tool runs: phasing and imputation are one Li-Stephens copying HMM (recombination is the transition, mutation the emission, the genetic map and Ne set the rates), imputation's honest output is a dosage with a self-estimated quality (INFO/R2/DR2) not a hard genotype, and the stages are ordered and each fails silently (QC, align build and strand to the panel, phase, impute per chromosome, fil
894SKILL.mdUpdated Jun 13, 2026
GPTomics/bio-pathway-enrichment-foundations
tools
VerifiedTrustedCommunity
Chooses the enrichment generation before any tool runs, mapping the input shape to a method class - a pre-selected gene list plus a background to over-representation analysis (ORA, hypergeometric), a ranked statistic for all genes to gene set enrichment (GSEA), a signed signaling topology to pathway-topology (SPIA) - then making the null explicit (competitive vs self-contained, gene vs subject sampling) and running a trustworthiness checklist (testable-gene universe, FDR, redundancy collapse, leading-edge check, version reporting). Covers why every clusterProfiler GSEA is the inter-gene-correlation-uncorrected competitive null, why the background not the gene list decides ORA significance, and why no method is universally best. Use when deciding ORA vs GSEA vs topology, which gene-set DB, whether a result is trustworthy, or which null a tool computes. For ORA see go-enrichment, GSEA see gsea, databases kegg-pathways/reactome-pathways/wikipathways; the ranking comes from differential-expression/de-results.
894SKILL.mdUpdated Jun 13, 2026
GPTomics/bio-workflows-gwas-pipeline
testing
VerifiedTrustedCommunity
End-to-end GWAS workflow from VCF to association results. Covers PLINK QC, population structure correction, and association testing for case-control or quantitative traits. Use when running genome-wide association studies.
894SKILL.mdUpdated Jun 11, 2026
GPTomics/bio-workflows-expression-to-pathways
development
VerifiedTrustedCommunity
Orchestrates the full path from differential expression results to redundancy-collapsed functional enrichment: choose ORA vs GSEA, convert gene IDs per method, run enrichGO/enrichKEGG/enrichPathway/enrichWP or gseGO/gseKEGG (clusterProfiler, ReactomePA, rWikiPathways), and visualize. Routes the ORA-vs-GSEA generation fork and the null/universe/reproducibility theory to pathway-analysis/enrichment-foundations. Use when a DESeq2/edgeR/limma result must become enriched GO terms, KEGG/Reactome/WikiPathways pathways, or a GSEA leading edge; when deciding whether a ranking exists for all genes (GSEA, named decreasing vector) or only a pre-selected list (ORA plus a defensible background universe); or when assembling DE-to-pathway end to end. The DE list and ranking statistic come from differential-expression/de-results; per-method nuance lives in the pathway-analysis skills.
894SKILL.mdUpdated Jun 10, 2026