How do I install the Bio Microbiome Amplicon Processing skill?

Run npx skillmds@latest add diegosouzapw/bio-microbiome-amplicon-processing in your terminal (requires Node.js), paste this page's agent-chat prompt into Claude, Cursor, or any MCP-connected agent, or download the SKILL.md file and copy it into your agent's skills directory.

Is Bio Microbiome Amplicon Processing safe to use?

This skill has not completed SkillMD's automated safety review yet. Independent scanners report: SkillSpector: PASS, Skill Scanner: PASS. Capability flags: docs only. SkillMD never runs a skill's scripts for you; review the SKILL.md before installing.

Which AI agents work with Bio Microbiome Amplicon Processing?

This skill is tagged as working with Claude Code, Claude.ai, OpenAI Codex. SKILL.md is an open format, so most agents that read a skills directory can load it too.

Is Bio Microbiome Amplicon Processing free to use?

Yes. Installing skills from SkillMD is free, and the skill stays under its author's original license.

Who published Bio Microbiome Amplicon Processing?

diegosouzapw (@diegosouzapw) published this skill. Their other Agent Skills are listed on their SkillMD profile.

Bio Microbiome Amplicon Processing by diegosouzapw

Name: Bio Microbiome Amplicon Processing
Rating: 4.2 (5 reviews)
Author: diegosouzapw

Bio Microbiome Amplicon Processing

Amplicon sequence variant (ASV) inference from 16S rRNA or ITS amplicon sequencing using DADA2. Covers quality filtering, error learning, denoising, and chimera removal. Use when processing demultiplexed amplicon FASTQ files to generate an ASV table for downstream analysis.

diegosouzapw Updated Sep 8, 2026 54 repo stars

File contents

Amplicon Processing with DADA2

Complete DADA2 Workflow

library(dada2)

path <- 'raw_reads'
fnFs <- sort(list.files(path, pattern = '_R1_001.fastq.gz', full.names = TRUE))
fnRs <- sort(list.files(path, pattern = '_R2_001.fastq.gz', full.names = TRUE))
sample_names <- sapply(strsplit(basename(fnFs), '_'), `[`, 1)

# Quality profiles
plotQualityProfile(fnFs[1:2])
plotQualityProfile(fnRs[1:2])

Quality Filtering and Trimming

filtFs <- file.path('filtered', paste0(sample_names, '_F_filt.fastq.gz'))
filtRs <- file.path('filtered', paste0(sample_names, '_R_filt.fastq.gz'))
names(filtFs) <- sample_names
names(filtRs) <- sample_names

# Filter parameters depend on amplicon region and read length
out <- filterAndTrim(fnFs, filtFs, fnRs, filtRs,
                     truncLen = c(240, 160),      # Trim to quality scores
                     maxN = 0,                     # No ambiguous bases
                     maxEE = c(2, 2),              # Max expected errors
                     truncQ = 2,                   # Truncate at first Q <= 2
                     rm.phix = TRUE,               # Remove PhiX
                     compress = TRUE,
                     multithread = TRUE)

Error Rate Learning

errF <- learnErrors(filtFs, multithread = TRUE)
errR <- learnErrors(filtRs, multithread = TRUE)

# Visualize error rates
plotErrors(errF, nominalQ = TRUE)

Sample Inference (Denoising)

dadaFs <- dada(filtFs, err = errF, multithread = TRUE)
dadaRs <- dada(filtRs, err = errR, multithread = TRUE)

# Check results
dadaFs[[1]]

Merge Paired Reads

mergers <- mergePairs(dadaFs, filtFs, dadaRs, filtRs, verbose = TRUE)

# Check merge success
head(mergers[[1]])

Construct Sequence Table

seqtab <- makeSequenceTable(mergers)
dim(seqtab)

# Check length distribution
table(nchar(getSequences(seqtab)))

Remove Chimeras

seqtab_nochim <- removeBimeraDenovo(seqtab, method = 'consensus',
                                     multithread = TRUE, verbose = TRUE)

# Percentage retained
sum(seqtab_nochim) / sum(seqtab)

Track Reads Through Pipeline

getN <- function(x) sum(getUniques(x))
track <- cbind(out, sapply(dadaFs, getN), sapply(dadaRs, getN),
               sapply(mergers, getN), rowSums(seqtab_nochim))
colnames(track) <- c('input', 'filtered', 'denoisedF', 'denoisedR', 'merged', 'nonchim')
rownames(track) <- sample_names
track

ITS-Specific Processing

# For ITS, use cutadapt to remove primers first (variable length amplicons)
# Then skip truncLen (don't truncate ITS to fixed length)

out_its <- filterAndTrim(fnFs, filtFs, fnRs, filtRs,
                         maxN = 0, maxEE = c(2, 2), truncQ = 2,
                         minLen = 50,  # Minimum length
                         rm.phix = TRUE, compress = TRUE, multithread = TRUE)

Related Skills

taxonomy-assignment - Assign taxonomy to ASVs
read-qc/quality-reports - Pre-DADA2 quality assessment
diversity-analysis - Analyze ASV table

---name: bio-microbiome-amplicon-processingdescription: Amplicon sequence variant (ASV) inference from 16S rRNA or ITS amplicon sequencing using DADA2. Covers quality filtering, error learning, denoising, and chimera removal. Use when processing demultiplexed amplicon FASTQ files to generate an ASV table for downstream analysis.---# Amplicon Processing with DADA2## Complete DADA2 Workflow```rlibrary(dada2)path <- 'raw_reads'fnFs <- sort(list.files(path, pattern = '_R1_001.fastq.gz', full.names = TRUE))fnRs <- sort(list.files(path, pattern = '_R2_001.fastq.gz', full.names = TRUE))sample_names <- sapply(strsplit(basename(fnFs), '_'), `[`, 1)# Quality profilesplotQualityProfile(fnFs[1:2])plotQualityProfile(fnRs[1:2])```## Quality Filtering and Trimming```rfiltFs <- file.path('filtered', paste0(sample_names, '_F_filt.fastq.gz'))filtRs <- file.path('filtered', paste0(sample_names, '_R_filt.fastq.gz'))names(filtFs) <- sample_namesnames(filtRs) <- sample_names# Filter parameters depend on amplicon region and read lengthout <- filterAndTrim(fnFs, filtFs, fnRs, filtRs, truncLen = c(240, 160), # Trim to quality scores maxN = 0, # No ambiguous bases maxEE = c(2, 2), # Max expected errors truncQ = 2, # Truncate at first Q <= 2 rm.phix = TRUE, # Remove PhiX compress = TRUE, multithread = TRUE)```## Error Rate Learning```rerrF <- learnErrors(filtFs, multithread = TRUE)errR <- learnErrors(filtRs, multithread = TRUE)# Visualize error ratesplotErrors(errF, nominalQ = TRUE)```## Sample Inference (Denoising)```rdadaFs <- dada(filtFs, err = errF, multithread = TRUE)dadaRs <- dada(filtRs, err = errR, multithread = TRUE)# Check resultsdadaFs[[1]]```## Merge Paired Reads```rmergers <- mergePairs(dadaFs, filtFs, dadaRs, filtRs, verbose = TRUE)# Check merge successhead(mergers[[1]])```## Construct Sequence Table```rseqtab <- makeSequenceTable(mergers)dim(seqtab)# Check length distributiontable(nchar(getSequences(seqtab)))```## Remove Chimeras```rseqtab_nochim <- removeBimeraDenovo(seqtab, method = 'consensus', multithread = TRUE, verbose = TRUE)# Percentage retainedsum(seqtab_nochim) / sum(seqtab)```## Track Reads Through Pipeline```rgetN <- function(x) sum(getUniques(x))track <- cbind(out, sapply(dadaFs, getN), sapply(dadaRs, getN), sapply(mergers, getN), rowSums(seqtab_nochim))colnames(track) <- c('input', 'filtered', 'denoisedF', 'denoisedR', 'merged', 'nonchim')rownames(track) <- sample_namestrack```## ITS-Specific Processing```r# For ITS, use cutadapt to remove primers first (variable length amplicons)# Then skip truncLen (don't truncate ITS to fixed length)out_its <- filterAndTrim(fnFs, filtFs, fnRs, filtRs, maxN = 0, maxEE = c(2, 2), truncQ = 2, minLen = 50, # Minimum length rm.phix = TRUE, compress = TRUE, multithread = TRUE)```## Related Skills- taxonomy-assignment - Assign taxonomy to ASVs- read-qc/quality-reports - Pre-DADA2 quality assessment- diversity-analysis - Analyze ASV table

Frequently asked questions