Templates and tools

This page is generated from the app's own catalog (version 2026.17.1), so it always matches what ships: 14 templates and 85 tools.

Templates

A template builds a whole pipeline from your files. The settings that depend on your data (adapter, strandedness, genome size) are asked before the workflow is created.

TemplateFieldLevelReadsStepsTools
Genome assembly from long reads (Flye and QUAST)
Put long Nanopore or PacBio reads of a microbial genome together without a reference, then measure how good the assembly is.
Genome assembly Intermediate Long reads or signal 2 Flye, QUAST
Genome assembly from short reads (SPAdes and QUAST)
Put Illumina reads of a microbial genome together without a reference, then measure how good the assembly is.
Genome assembly Intermediate Single- or paired-end reads 2 SPAdes, QUAST
ATAC-seq open chromatin peaks with Genrich
Find open chromatin regions from paired-end ATAC-seq reads: align, call peaks in ATAC mode, merge them and make a signal track.
Epigenomics Advanced Paired-end reads 8 Bowtie2 build, Bowtie2, samtools sort, Genrich, bedtools merge, samtools index, deepTools bamCoverage
Basic read quality check
Check the quality of sequencing reads, trim them, check again and get one combined report.
Read quality Beginner Single-end reads 4 FastQC, fastp, MultiQC
ChIP-seq peak calling with MACS3
Find where a protein binds the genome: trim, align a ChIP sample and its input control, call peaks, draw the signal and say which genes the peaks lie near.
Epigenomics Advanced Single-end reads 14 fastp, Bowtie2 build, Bowtie2, samtools sort, samtools index, MACS3 callpeak, deepTools bamCoverage, deepTools plotFingerprint, HOMER annotatePeaks
Who is in my sample? Kraken2 and Bracken
Find which organisms a metagenomic sample contains and how abundant they are: trim, classify the reads with Kraken2, estimate abundances with Bracken and gather the results in one report.
Metagenomics Intermediate Single-end reads 4 fastp, Kraken2, Bracken, MultiQC
Nanopore reads from FASTQ to coverage
Filter Nanopore reads, check their quality, map them to a genome and measure coverage. Works on every computer.
Long reads Intermediate Long reads or signal 6 chopper, NanoPlot, minimap2 (long reads), samtools sort, samtools index, mosdepth
Nanopore from raw signal to coverage
Start from the raw POD5 signal files of a Nanopore run: call the bases with Dorado, filter, check the reads, map them to a genome and measure coverage.
Long reads Advanced Long reads or signal 10 pod5 merge, pod5 inspect, Dorado basecaller, samtools fastq, chopper, NanoPlot, minimap2 (long reads), samtools sort, samtools index, mosdepth
Ribo-seq with UMIs (riboWaltz)
Ribosome profiling with UMIs: trim, extract UMIs, remove rRNA/tRNA/ncRNA, map to transcripts with STAR, deduplicate, riboWaltz report and MultiQC.
RNA sequencing Advanced Single-end reads 20 FastQC, Cutadapt, UMI-tools extract, Bowtie2 build, Remove reads matching a database (bowtie2), STAR genome index, STAR (Ribo-seq, transcriptome BAM), samtools view, samtools sort, samtools index, UMI-tools dedup, riboWaltz report, MultiQC
RNA-seq alignment and gene counts
Align RNA-seq reads to a genome with HISAT2 and count the reads per gene, ready for differential expression.
RNA sequencing Intermediate Single-end reads 7 fastp, HISAT2 build, HISAT2, samtools sort, samtools index, featureCounts, MultiQC
RNA-seq quantification with Salmon
Estimate how much of each transcript is present in an RNA-seq sample, without aligning the reads to a genome.
RNA sequencing Intermediate Single-end reads 4 fastp, Salmon index, Salmon quant, MultiQC
Germline variant calling, paired-end reads (bcftools)
Find variants in one paired-end sample against a reference genome, with alignment statistics and a filtered VCF.
DNA sequencing Intermediate Paired-end reads 10 BWA index, BWA MEM, samtools sort, samtools markdup, samtools index, samtools flagstat, bcftools call, bcftools filter, bcftools stats, MultiQC
Germline variant calling with bcftools
Find single-base differences between one sample and a reference genome, from trimming to a filtered, summarised VCF.
DNA sequencing Intermediate Single-end reads 10 fastp, BWA index, BWA MEM, samtools sort, samtools markdup, samtools index, bcftools call, bcftools filter, bcftools stats, MultiQC
Germline variant calling with GATK
Find variants in one sample with GATK HaplotypeCaller, the standard caller for human and other diploid genomes.
DNA sequencing Advanced Single-end reads 11 fastp, BWA index, BWA MEM, samtools sort, Picard MarkDuplicates, samtools index, samtools faidx, GATK CreateSequenceDictionary, GATK HaplotypeCaller, bcftools stats, MultiQC

Catalog tools

Each entry names every file it reads and writes, installs the tool at a pinned version and was run against the real program. Jump to: Quality control | Trimming | Reference preparation | Alignment | BAM processing | Variant calling | Variant annotation | Quantification | Peak calling and epigenomics | Signal tracks and plots | Genome intervals | FASTA and FASTQ tools | Assembly | Metagenomics | Nanopore signal data

Quality control 9

ToolWhat it doesReadsMakesInstalled as
FastQC
Read quality
Quality report for raw sequencing reads. fastq html, zip fastqc 0.13.0 (bioconda)
MultiQC
Report summaries
Combines the reports of many tools into one HTML page. any html multiqc 1.35 (bioconda)
mosdepth
Coverage
Fast coverage summary: the average depth per chromosome and in windows along the genome. bam, bai txt, bed mosdepth 0.3.13 (bioconda)
Qualimap BamQC
Alignment QC
Quality report of an alignment: coverage, mapping quality, insert size and base content. bam html, txt qualimap 2.3 (bioconda)
RSeQC infer experiment
RNA-seq QC
Works out whether an RNA-seq library is stranded, and which way, from where reads fall on genes. bam, bed txt rseqc 5.0.5 (bioconda)
RSeQC read distribution
RNA-seq QC
Shows how the reads divide between exons, introns and the regions around genes. bam, bed txt rseqc 5.0.5 (bioconda)
deepTools plotFingerprint
Enrichment
Shows how well a ChIP-seq or ATAC-seq sample is enriched: a good ChIP puts most of its reads into a small part of the genome, an input control spreads them evenly. bam, bai image, tsv deeptools 3.5.6 (bioconda)
NanoPlot
Long-read QC
Summarises a Nanopore or PacBio run: how many reads, how long and how accurate, with plots. Reads FASTQ, a Dorado BAM or a sequencing summary. fastq, fasta, bam, tsv html, txt nanoplot 1.48.0 (bioconda)
seqkit stats
Sequence statistics
Counts the sequences in FASTA or FASTQ files and reports their total length, shortest, longest, N50 and average quality, one line per file. fastq, fasta tsv seqkit 2.14.0 (bioconda)

Trimming 6

ToolWhat it doesReadsMakesInstalled as
fastp
Trim and filter reads
Fast adapter and quality trimming for FASTQ reads. fastq fastq, html, json fastp 1.3.7 (bioconda)
Cutadapt
Trim and filter reads
Removes adapter sequences and low-quality ends from reads. fastq fastq, txt cutadapt 5.2 (bioconda)
Filtlong
Long-read filtering
Keeps the best long reads of a Nanopore or PacBio run: drops reads that are too short or too inaccurate, and can keep only as many bases as you need. Writes a compressed FASTQ. fastq fastq filtlong 0.3.1 (bioconda)
chopper
Long-read filtering
Filters and trims long reads by their average quality and length. Fast, and the usual first step for Nanopore reads. Writes a compressed FASTQ. fastq fastq chopper 0.14.1 (bioconda)
UMI-tools extract
UMI handling
Moves each read's UMI (unique molecular identifier) into its name and cuts it, and any spacer next to it, out of the sequence. Run it before aligning so duplicates can be removed later. fastq fastq, txt umi_tools 1.1.6 (bioconda)
Remove reads matching a database (bowtie2)
Remove contaminant reads
Aligns reads to sequences you do not want (rRNA, tRNA, small RNAs, PhiX, a contaminant genome) and keeps the reads that do not match. Matches are counted in the log, not saved. index, fastq fastq, txt bowtie2 2.5.5 (bioconda)

Reference preparation 10

ToolWhat it doesReadsMakesInstalled as
BWA index
Aligner indexes
Builds the index BWA MEM needs to align reads to a genome. fasta index, fasta bwa 0.7.19 (bioconda)
Bowtie2 build
Aligner indexes
Builds the genome index Bowtie2 needs to align reads. fasta index bowtie2 2.5.5 (bioconda)
samtools faidx
Reference indexes
Indexes a reference FASTA (makes the .fai file variant callers need). fasta fasta, fai samtools 1.24 (bioconda)
GATK CreateSequenceDictionary
Reference indexes
Makes the sequence dictionary (.dict) that GATK needs next to the reference. fasta fasta, fai, dict gatk4 4.6.2.0 (bioconda)
HISAT2 build
Aligner indexes
Builds the genome index HISAT2 needs to align RNA-seq reads. fasta index hisat2 2.2.3 (bioconda)
Salmon index
Aligner indexes
Builds the index Salmon quant needs from transcript sequences. fasta index salmon 2.8.0 (bioconda)
STAR genome index
Aligner indexes
Builds the genome index STAR needs to align RNA-seq reads, with the known splice junctions of an annotation. fasta, gtf index star 2.7.10b (bioconda)
kallisto index
Aligner indexes
Builds the index kallisto quant needs from transcript sequences. fasta index kallisto 0.52.0 (bioconda)
gffread transcripts
Annotation files
Cuts the transcript sequences (exons joined together) out of a genome using an annotation. fasta, gtf fasta gffread 0.12.9 (bioconda)
Annotation to BED12
Annotation files
Converts a gene annotation (GTF or GFF3) to the BED12 gene model that RSeQC reads. gtf bed gffread 0.12.9 (bioconda)

Alignment 7

ToolWhat it doesReadsMakesInstalled as
BWA MEM
DNA aligners
Aligns reads to a reference genome. Writes SAM. fasta, fastq sam bwa 0.7.19 (bioconda)
Bowtie2
DNA aligners
Aligns reads to an indexed genome, single-end or as read pairs. Writes SAM. index, fastq sam bowtie2 2.5.5 (bioconda)
minimap2 (short reads)
DNA aligners
Aligns short reads (the sr preset) to a reference genome. Writes SAM with a read group. fasta, fastq sam minimap2 2.30 (bioconda)
STAR
RNA-seq aligners
Splice-aware RNA-seq aligner, single-end or paired-end. Writes a coordinate-sorted BAM. index, fastq bam, txt star 2.7.10b (bioconda)
HISAT2
RNA-seq aligners
Splice-aware RNA-seq aligner, single-end or paired-end. Writes SAM. index, fastq sam, txt hisat2 2.2.3 (bioconda)
minimap2 (long reads)
Long-read aligners
Aligns long reads (Nanopore, PacBio HiFi or cDNA) to a reference genome. Pick the preset that matches how the reads were made. Writes SAM. fasta, fastq sam minimap2 2.30 (bioconda)
STAR (Ribo-seq, transcriptome BAM)
RNA-seq aligners
Aligns footprint reads to the genome and writes them in transcript coordinates, as riboWaltz needs. Defaults suit short footprints; reads that did not align stay in the BAM, marked unmapped. index, fastq bam, txt star 2.7.10b (bioconda)

BAM processing 10

ToolWhat it doesReadsMakesInstalled as
samtools sort
Alignment files
Sorts alignments by coordinate (or by read name). sam, bam bam samtools 1.24 (bioconda)
samtools index
Alignment files
Builds the .bai index of a coordinate-sorted BAM. bam bai samtools 1.24 (bioconda)
samtools view
Alignment files
Filters alignments and converts them to BAM. sam, bam bam samtools 1.24 (bioconda)
samtools flagstat
Alignment files
Counts reads by alignment flag. Writes a text summary. sam, bam txt samtools 1.24 (bioconda)
samtools fastq
Alignment files
Turns an unaligned BAM, such as the output of the Dorado basecaller, into a FASTQ file of reads. Writes compressed FASTQ. bam fastq samtools 1.24 (bioconda)
samtools merge
Alignment files
Combines several sorted BAM files into one. bam bam samtools 1.24 (bioconda)
samtools markdup
Duplicates
Finds PCR duplicates and marks (or removes) them. Runs fixmate and sort first, so plain aligner output works. sam, bam bam, txt samtools 1.24 (bioconda)
Picard MarkDuplicates
Duplicates
Finds PCR duplicates in a sorted BAM and marks (or removes) them. Writes a duplication report. bam bam, txt picard 3.4.0 (bioconda)
samtools depth
Coverage
Writes how many reads cover each position of the genome, one column per BAM file. bam tsv samtools 1.24 (bioconda)
UMI-tools dedup
Duplicates
Removes PCR duplicates: reads with the same start and UMI (allowing for UMI errors) count once. Needs a sorted BAM with its index, and read names that end in the UMI. bam, bai bam, txt umi_tools 1.1.6 (bioconda)

Variant calling 8

ToolWhat it doesReadsMakesInstalled as
bcftools call
SNPs and small indels
Calls SNPs and indels from a sorted BAM (mpileup then call). fasta, bam vcf bcftools 1.24 (bioconda)
FreeBayes
SNPs and small indels
Calls SNPs and small indels from one or more sorted, indexed BAM files. fasta, bam vcf freebayes 1.3.10 (bioconda)
GATK HaplotypeCaller
SNPs and small indels
Calls SNPs and small indels by re-assembling the reads around each variant. fasta, bam, bai vcf gatk4 4.6.2.0 (bioconda)
bcftools filter
VCF tools
Keeps the variants that are good enough (enough quality and depth), or marks the others. vcf vcf bcftools 1.24 (bioconda)
bcftools norm
VCF tools
Writes variants in one standard form: indels left-aligned, optionally one line per alternative allele. vcf, fasta vcf bcftools 1.24 (bioconda)
bcftools stats
VCF tools
Counts the variants in a VCF: SNPs, indels, transitions and transversions. MultiQC reads it. vcf txt bcftools 1.24 (bioconda)
bgzip
VCF tools
Compresses a VCF, BED or other table so it can be indexed with tabix. vcf, bed, gtf, tsv vcf, bed, gtf, tsv htslib 1.24 (bioconda)
tabix
VCF tools
Indexes a compressed, sorted VCF, BED or GFF so a region can be read without reading the whole file. vcf, bed, gtf, tsv vcf, bed, gtf, tsv, tbi htslib 1.24 (bioconda)

Variant annotation 2

ToolWhat it doesReadsMakesInstalled as
snpEff build database
Custom databases
Builds a snpEff database for your own genome from its FASTA and gene annotation (GTF). fasta, gtf database snpeff 5.4.0c (bioconda)
snpEff annotate
Variant effects
Says what each variant does to the genes: missense, synonymous, intron, intergenic and so on. Needs a snpEff database. vcf, database vcf, html, tsv snpeff 5.4.0c (bioconda)

Quantification 5

ToolWhat it doesReadsMakesInstalled as
Salmon quant
Transcript abundance
Quantifies transcript abundance from single-end reads. index, fastq tsv salmon 2.8.0 (bioconda)
featureCounts
Gene counts
Counts reads per gene from alignments and an annotation. bam, sam, gtf tsv, txt subread 2.1.1 (bioconda)
kallisto quant
Transcript abundance
Quantifies transcript abundance from single-end or paired-end reads, without aligning them. index, fastq tsv, txt kallisto 0.52.0 (bioconda)
StringTie assemble
Transcript assembly
Assembles transcripts from RNA-seq alignments, optionally guided by a gene annotation, and estimates their abundance. bam, gtf gtf, tsv stringtie 3.0.3 (bioconda)
StringTie quantify
Transcript assembly
Measures the abundance of the transcripts of an annotation (no new transcripts are made). bam, gtf gtf, tsv stringtie 3.0.3 (bioconda)

Peak calling and epigenomics 3

ToolWhat it doesReadsMakesInstalled as
MACS3 callpeak
Peak calling
Finds the places in the genome where a ChIP-seq (or ATAC-seq) sample has more reads than expected, with an optional input control. Makes sharp peaks, summits and a signal track. bam peaks, bed, tsv, bedgraph macs3 3.0.5 (bioconda)
Genrich
Peak calling
Finds open chromatin in ATAC-seq data (or enriched regions in ChIP-seq). It reads all replicates together and takes the Tn5 cut sites into account. bam, bed peaks, bed genrich 0.6.1 (bioconda)
HOMER annotatePeaks
Peak annotation
Tells you where each peak lies: in a promoter, an exon, an intron or between genes, and which gene is closest. Uses the genome FASTA and gene annotation you connect. Needs a Genome FASTA and gene annotation. bed, peaks, fasta, gtf tsv, txt homer 5.1 (bioconda)

Signal tracks and plots 4

ToolWhat it doesReadsMakesInstalled as
deepTools bamCoverage
Signal tracks
Turns alignments into a signal track (bigWig) that shows how many reads cover each stretch of the genome, scaled so that samples can be compared. bam, bai, bed bigwig deeptools 3.5.6 (bioconda)
deepTools computeMatrix
Signal around regions
Collects the signal around one point of every region (the centre of each peak, or the start of each gene) into a table. plotHeatmap draws it. bigwig, bed, peaks, gtf matrix, bed, tsv deeptools 3.5.6 (bioconda)
deepTools plotHeatmap
Plots
Draws a heatmap of the signal around your regions, one row per region, with an average profile on top. matrix image deeptools 3.5.6 (bioconda)
riboWaltz report
Ribosome profiling
Finds where the ribosome sits in each footprint (the P-site offset) and writes one HTML report: lengths, offsets, 3-nt periodicity, start and stop profiles, regions, codon usage, plus TSV tables. bam, gtf, fasta html, tsv ribowaltz 2.0 (bioconda)

Genome intervals 4

ToolWhat it doesReadsMakesInstalled as
bedtools intersect
Compare regions
Compares two sets of regions: keeps the ones that overlap, the ones that do not, or just the overlapping stretches. For example, which peaks of one sample are also peaks in another. bed, peaks, gtf, vcf bed, peaks bedtools 2.31.1 (bioconda)
bedtools merge
Combine regions
Joins regions that overlap (or lie close together) into single regions. Give it the peaks of several samples to get one combined set. bed, peaks bed bedtools 2.31.1 (bioconda)
bedtools genomecov
Coverage
Measures how deeply every stretch of the genome is covered by reads and writes it as a bedGraph track. bam bedgraph bedtools 2.31.1 (bioconda)
bedtools getfasta
Sequences
Cuts the DNA sequence of every region out of a genome FASTA, for example the sequences under your peaks for a motif search. fasta, bed, peaks, gtf fasta bedtools 2.31.1 (bioconda)

FASTA and FASTQ tools 2

ToolWhat it doesReadsMakesInstalled as
seqkit seq (filter or convert)
Filter and convert
Keeps the sequences within a length or quality range, and can reverse-complement or upper-case them. Writes a compressed file. fastq, fasta fastq, fasta seqkit 2.14.0 (bioconda)
seqkit grep (pick sequences)
Pick sequences
Picks sequences by name or by a stretch of sequence, from a list of names or a pattern you type. Can also leave those out instead. fastq, fasta, txt, tsv fastq, fasta seqkit 2.14.0 (bioconda)

Assembly 3

ToolWhat it doesReadsMakesInstalled as
Flye
Long-read assemblers
Assembles a genome from long reads (Nanopore or PacBio) without a reference. Writes the contigs, the assembly graph and a table of contig lengths and coverage. fastq, fasta fasta, gfa, tsv, txt flye 2.9.6 (bioconda)
SPAdes
Short-read assemblers
Assembles a genome from short Illumina reads without a reference. Writes the contigs, the scaffolds and the assembly graph. fastq fasta, gfa, txt spades 4.2.0 (bioconda)
QUAST
Assembly quality
Judges genome assemblies: how many contigs, how long, and (with a reference) how complete and correct. Several assemblies go into one comparison. fasta tsv, html, txt quast 5.3.0 (bioconda)

Metagenomics 4

ToolWhat it doesReadsMakesInstalled as
Kraken2 build database
Databases
Builds a Kraken2 database from your own genomes, for example the few species you expect in a sample. Needs a taxonomy that names every genome. fasta, taxonomy database kraken2 2.17.1 (bioconda)
Kraken2
Classification
Names the organism each read most likely comes from (its taxonomy) by matching it against a database of genomes. Writes a summary table by species and a line for every read. Needs a Kraken2 database. fastq, fasta, database tsv, txt kraken2 2.17.1 (bioconda)
Bracken prepare database
Databases
Gets a Kraken2 database ready for Bracken. Makes a new folder that links to your database and adds the files Bracken needs for one read length. Your database is not changed. Needs a Kraken2 database. database database bracken 2.6.1 (bioconda)
Bracken
Abundance
Turns a Kraken2 summary into better abundance estimates. Reads that Kraken2 could only place on a genus are shared out over the species of that genus. Needs a Bracken database. tsv, txt, database tsv, txt bracken 2.6.1 (bioconda)

Nanopore signal data 8

ToolWhat it doesReadsMakesInstalled as
pod5 view
POD5 files
Lists the reads in POD5 files as a table: read ID, channel, length in samples, start time and more. Use it to pick reads for pod5 filter or pod5 subset. pod5 tsv pod5 0.3.48 (bioconda)
pod5 inspect
POD5 files
Reports what a POD5 file holds: how many reads it has, how they are grouped into batches, and the version of the file format. pod5 txt pod5 0.3.48 (bioconda)
pod5 filter
POD5 files
Copies only the reads you list into a new POD5 file. Use it to basecall or share a small part of a run. pod5, txt, tsv pod5 pod5 0.3.48 (bioconda)
pod5 subset (split by column)
POD5 files
Splits POD5 files into several files by a value in a table, for example one file per channel. The table can come from pod5 view or from a sequencing summary. pod5, tsv pod5 pod5 0.3.48 (bioconda)
pod5 merge
POD5 files
Joins several POD5 files into one, for example the many files MinKNOW writes during a run. pod5 pod5 pod5 0.3.48 (bioconda)
pod5 convert from FAST5
POD5 files
Converts old FAST5 signal files to one POD5 file, the format Dorado reads. fast5 pod5 pod5 0.3.48 (bioconda)
Dorado basecaller
Basecalling
Turns the raw signal of a Nanopore run (POD5) into DNA sequence. Writes an unaligned BAM. Needs a graphics card or an Apple chip for anything but tests. Needs a Basecalling model. pod5 bam dorado 1.4.0 (downloaded from the vendor, Oxford Nanopore Technologies Public License)
Dorado summary
Basecalling
Makes a table with one line per read from a Dorado BAM: read ID, channel, length and quality. NanoPlot and pod5 subset read it. bam tsv dorado 1.4.0 (downloaded from the vendor, Oxford Nanopore Technologies Public License)