bio-reporting-automated-qc-reports
$
npx mdskill add GPTomics/bioSkills/bio-reporting-automated-qc-reportsGenerates standardized QC reports by aggregating bioinformatics tool outputs with MultiQC
- Summarizes quality control metrics across multiple sequencing samples
- Leverages FastQC, STAR, HISAT2, and other tools via MultiQC integration
- Automatically detects and parses recognized tool output formats
- Produces interactive HTML reports for easy visualization and sharing
SKILL.md
.github/skills/bio-reporting-automated-qc-reportsView on GitHub ↗
---
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>` then `help(module.function)` to check signatures
- CLI: `<tool> --version` then `<tool> --help` to 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
```bash
# 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`:
```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
```bash
# Add custom data file
# File format: sample\tmetric1\tmetric2
multiqc results/ --data-format tsv --custom-data-file custom_metrics.tsv
```
## Python API
```python
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
More from GPTomics/bioSkills
- bio-admet-predictionPredicts ADMET properties using ADMETlab 3.0 API or DeepChem models. Estimates bioavailability, CYP inhibition, hERG liability, and 119 toxicity endpoints with uncertainty quantification. Filters for PAINS and other structural alerts. Use when filtering compounds for drug-likeness or prioritizing leads by predicted safety.
- bio-alignment-amplicon-clippingTrim PCR primers from aligned reads in amplicon-panel BAMs using samtools ampliconclip. Use when processing SARS-CoV-2 ARTIC, hereditary cancer panels, ctDNA hot-spot panels, or any amplicon assay where primer-derived bases would falsely confirm reference at primer footprints.
- bio-alignment-filteringFilter alignments by flags, mapping quality, and regions using samtools view and pysam. Use when extracting specific reads, removing low-quality alignments, or subsetting to target regions.
- bio-alignment-indexingCreate and use BAI/CSI indices for BAM/CRAM files using samtools and pysam. Use when enabling random access to alignment files or fetching specific genomic regions.
- bio-alignment-ioRead, write, and convert multiple sequence alignment files using Biopython Bio.AlignIO. Supports Clustal, PHYLIP, Stockholm, FASTA, Nexus, and other alignment formats for phylogenetics and conservation analysis. Use when reading, writing, or converting alignment file formats.
- bio-alignment-msa-parsingParse and analyze multiple sequence alignments using Biopython. Extract sequences, identify conserved regions, analyze gaps, work with annotations, and manipulate alignment data for downstream analysis. Use when parsing or manipulating multiple sequence alignments.
- bio-alignment-msa-statisticsCalculate alignment statistics including sequence identity, conservation scores, substitution matrices, and similarity metrics. Use when comparing alignment quality, measuring sequence divergence, and analyzing evolutionary patterns.
- bio-alignment-multiplePerform multiple sequence alignment using MAFFT, MUSCLE5, ClustalOmega, or T-Coffee. Guides tool and algorithm selection based on dataset size, sequence divergence, and downstream application. Use when aligning three or more homologous sequences for phylogenetics, conservation analysis, or evolutionary studies.
- bio-alignment-pairwisePerform pairwise sequence alignment using Biopython Bio.Align.PairwiseAligner. Use when comparing two sequences, finding optimal alignments, scoring similarity, and identifying local or global matches between DNA, RNA, or protein sequences.
- bio-alignment-sortingSort alignment files by coordinate or read name using samtools and pysam. Use when preparing BAM files for indexing, variant calling, or paired-end analysis.