omics.co.in logo omics.co.in
← Back to Directory Hub (All 18 Disciplines)
← Prev: Stage 10 | Next: Stage 12 →
Discipline 11 of 18 36 Tools Stage 11: Single-Molecule Real-Time Sequencing

⚡ 11. Long-Read Sequencing (PacBio HiFi & Oxford Nanopore)

Direct neural network basecalling, modified base detection (5mC/5hmC), full-length isoform sequencing, and long-range structural variant resolution across ONT and PacBio platforms.

Essential Software: Dorado Guppy pbmm2 Clair3 modkit IsoSeq
Comparative Analysis Production Benchmarks

Algorithm Architecture & Decision Matrices

Detailed comparison cards highlighting algorithmic trade-offs, strengths, weaknesses, and recommended production setups.

Basecalling, Platform Utilities & QC
Dorado (ONT Standard) Guppy (Legacy) SMRT Link (PacBio) pbtk NanoPlot NanoStat LongQC Filtlong modkit Nanopolish Remora pb-CpG-tools
Direct pod5 GPU basecalling, duplex read pairing, raw Nanopore signal QC, and 5mC/5hmC direct detection.
Alignment, Phasing & Assembly
minimap2 Winnowmap2 pbmm2 NGMLR Clair3 DeepVariant Longshot WhatsHap HapCUT2 hifiasm Flye Polypolish NextPolish Inspector
Error-tolerant spliced long-read alignment, tandem repeat mapping, haplotype phasing, and long-read neural variant calling.
Isoforms & Structural Variation
IsoSeq FLAIR StringTie2 Bambu IsoQuant SQANTI3 Sniffles2 cuteSV SVIM pbsv
Full-length cDNA/direct RNA sequencing, novel isoform validation, TSS/poly(A) site mapping, and structural variation resolution.
Hardware & Execution Profiles

Primary Software & Algorithms in Long-Read Sequencing Technologies

Direct command-line recipes, Conda environment specifications, memory footprints, and license metadata.

Genomic Alignment
🚀 32GB+ RAM C++ • AVX-512 / SSE4.1

BWA-MEM2

Docs ↗ By Heng Li & Vasimuddin Md (Intel Labs)
(License: MIT)

SIMD hardware-accelerated rewrite of BWA-MEM. Employs 8-way and 16-way vectorization to align short DNA reads 2x faster to human reference genomes (GRCh38/T2T).

Formats: Input: FASTQ, FASTA • Output: SAM, BAM (via samtools)
💻 8-32 CPU Cores • 32-45 GB RAM (Index resident)
conda install -c bioconda bwa-mem2
Long-Read & Splice Aligner
💻 16-32GB RAM C • Approximate Seed Matches

Minimap2

Docs ↗ By Heng Li (Dana-Farber / Harvard)
(License: MIT)

Fast pairwise aligner for PacBio HiFi, Oxford Nanopore, spliced direct RNA, and full-genome assembly synteny. Industry standard for third-generation sequencing.

Formats: Input: FASTQ, FASTA • Output: PAF, SAM, BAM
💻 8-16 CPU Cores • 16-32 GB RAM
conda install -c bioconda minimap2
Spliced RNA Aligner
🚀 32GB+ RAM C++ • Suffix Array / MMP

STAR

Docs ↗ By Alexander Dobin (CSHL)
(License: GPL-3.0)

Ultrafast universal spliced RNA-seq aligner based on sequential Maximum Mappable Prefix matching. High sensitivity for non-canonical splice junctions and chimeric gene fusions.

Formats: Input: FASTQ • Output: Sorted BAM, Gene Counts, Chimeric SAM
💻 8-16 CPU Cores • 32-40 GB RAM (Human genome index)
conda install -c bioconda star
Fast Read Alignment
💻 < 8GB RAM C++ • FM-Index

Bowtie 2

Docs ↗ By Ben Langmead & Steven Salzberg (JHU)
(License: GPL-3.0)

Memory-efficient Burrows-Wheeler transform tool for aligning sequencing reads to long reference genomes. Ideal for ChIP-seq, ATAC-seq, and bacterial assemblies.

Formats: Input: FASTQ, FASTA • Output: SAM, BAM
💻 4-8 CPU Cores • 4-8 GB RAM
conda install -c bioconda bowtie2
Graph-based Spliced Aligner
💻 16-32GB RAM C++ • Hierarchical Graph FM

HISAT2

Docs ↗ By Daehwan Kim & Steven Salzberg (JHU)
(License: GPL-3.0)

Hierarchical Graph FM index aligner capable of mapping RNA-seq reads across population-scale variant graphs, reducing reference allele alignment bias.

Formats: Input: FASTQ • Output: SAM, BAM
💻 8 CPU Cores • 8-16 GB RAM
conda install -c bioconda hisat2
Quasi-Mapping Quantifier
💻 16-32GB RAM C++ • Dual-phase Optimization

Salmon

Docs ↗ By Rob Patro Lab (UMD)
(License: GPL-3.0)

Ultra-fast transcript abundance quantification from RNA-seq reads using quasi-mapping with GC-bias, sequence-specific bias, and fragment-length correction.

Formats: Input: FASTQ • Output: quant.sf, TPM, Bootstrap Estimates
💻 4-8 CPU Cores • 8-16 GB RAM
conda install -c bioconda salmon
Pseudoalignment
💻 < 8GB RAM C++ • de Bruijn Graph

kallisto

Docs ↗ By Nicolas Bray & Lior Pachter (Caltech)
(License: GPL-3.0)

Near-instantaneous pseudoalignment algorithm for quantifying transcript abundances from RNA-Seq reads by matching k-mers to colored de Bruijn graph paths without base-to-base alignment.

Formats: Input: FASTQ • Output: abundance.h5, abundance.tsv
💻 4 CPU Cores • 4-8 GB RAM
conda install -c bioconda kallisto
Pangenome Graph Aligner
🚀 32GB+ RAM C++ • Variation Graphs / GBWT

vg / Giraffe

Docs ↗ By Human Pangenome Reference Consortium
(License: Apache-2.0)

Haplotype-aware pangenome read aligner. Maps sequencing reads to variation graphs as fast as linear aligners, eliminating reference bias and uncovering structural alleles.

Formats: Input: FASTQ, GFA, GBZ • Output: GAM, GAF, BAM
💻 16+ CPU Cores • 32-64 GB RAM
conda install -c bioconda vg
Germline & Somatic Caller
💻 16-32GB RAM Java • Apache Spark

GATK4

Docs ↗ By Broad Institute
(License: BSD-3-Clause)

Gold-standard genome analysis toolkit for germline (HaplotypeCaller) and somatic (Mutect2) variant discovery, Base Quality Score Recalibration (BQSR), and VQSR.

Formats: Input: BAM, CRAM, FASTA • Output: VCF, gVCF
💻 8-16 CPU Cores • 16-32 GB RAM
conda install -c bioconda gatk4
Deep Learning CNN Variant Caller
⚡ GPU Required Python / C++ • TensorFlow / GPU

DeepVariant

Docs ↗ By Google Health & Genomics
(License: BSD-3-Clause)

Deep convolutional neural network that calls genetic variants by transforming aligned BAM reads into multi-channel pileup images. Achieves precision scores >99.9%.

Formats: Input: BAM, CRAM, FASTA • Output: VCF, gVCF
💻 16 CPU Cores • 32 GB RAM • NVIDIA GPU (Optional)
docker pull google/deepvariant:latest
VCF/BCF Manipulation
💻 < 8GB RAM C • HTSlib

bcftools

Docs ↗ By Petr Danecek & Heng Li
(License: MIT / GPL-3)

Blazing-fast utility suite for calling, filtering, sorting, merging, and querying VCF/BCF variant call format files. Includes mpileup genotype caller and consensus builder.

Formats: Input: VCF, BCF, BAM • Output: Indexed VCF, BCF, FASTA
💻 2-4 CPU Cores • 2-4 GB RAM
conda install -c bioconda bcftools
Haplotype-based Caller
💻 < 8GB RAM C++ • Bayesian Haplotype Models

FreeBayes

Docs ↗ By Erik Garrison & Gabor Marth
(License: MIT)

Bayesian genetic variant detector designed to find SNPs, indels, MNPs, and complex haplotypes smaller than sequencing read length directly from short read alignments.

Formats: Input: BAM • Output: VCF
💻 4-8 CPU Cores • 8-16 GB RAM
conda install -c bioconda freebayes
Long-Read SV Caller
💻 16-32GB RAM Python / C++ • Clustering Algorithms

Sniffles2

Docs ↗ By Smok et al. (Fritz Sedlazeck Lab)
(License: MIT)

Fast and accurate structural variant caller for long reads (PacBio HiFi and Oxford Nanopore). Calls deletions, duplications, inversions, and translocations with multi-sample genotype calling.

Formats: Input: BAM, CRAM • Output: VCF, .snf files
💻 8 CPU Cores • 16 GB RAM
conda install -c bioconda sniffles
Short-Read SV Discovery
💻 16-32GB RAM C++ / Python • Graph Assembly

Manta

Docs ↗ By Illumina
(License: GPL-3.0)

Somatic and germline structural variant and indel caller for paired-end short reads. Combines discordant paired-end alignment evidence with split-read graph assembly.

Formats: Input: BAM, CRAM • Output: VCF
💻 8-16 CPU Cores • 16-32 GB RAM
conda install -c bioconda manta
Variant Effect Predictor
💻 16-32GB RAM Perl / C • Cache Indexes

Ensembl VEP

Docs ↗ By William McLaren et al. (EMBL-EBI)
(License: Apache-2.0)

Determines the clinical consequence of genomic variants in coding and non-coding regions. Annotates SIFT, PolyPhen, gnomAD allele frequencies, ClinVar, and splice sites.

Formats: Input: VCF • Output: Annotated VCF, HTML, JSON
💻 4-8 CPU Cores • 16-32 GB RAM (Local cache)
conda install -c bioconda ensembl-vep
Variant Annotation & Filtering
💻 < 8GB RAM Java • Pre-built Genome DBs

SnpEff & SnpSift

Docs ↗ By Pablo Cingolani
(License: LGPL-3.0)

Rapid genomic variant annotation and functional effect prediction toolbox. Calculates consequence categories (missense, nonsense, frameshift) and filters multi-sample VCFs.

Formats: Input: VCF • Output: Annotated VCF, HTML summary
💻 4 CPU Cores • 8-16 GB RAM
conda install -c bioconda snpeff
Structural Variant Caller
💻 < 8GB RAM C++ • Paired-End & Split-Read

Delly

Docs ↗ By Tobias Rausch et al. (EMBL)
(License: BSD-3-Clause)

Integrated structural variant discovery tool utilizing paired-end discordant mapping and split-read clipping to resolve somatic and germline SV breakpoints at base-pair precision.

Formats: Input: BAM, CRAM • Output: BCF, VCF
💻 4-8 CPU Cores • 8-16 GB RAM
conda install -c bioconda delly
BAM/SAM/CRAM Utilities
💻 < 8GB RAM C • High-Throughput HTSlib

SAMtools / HTSlib

Docs ↗ By Heng Li, James Bonfield & Petr Danecek
(License: MIT / BSD)

The essential foundational software library and command-line utilities for reading, writing, sorting, indexing, viewing, and calculating depth coverage across SAM, BAM, and CRAM alignment files.

Formats: Input: SAM, BAM, CRAM • Output: Indexed BAM, CRAM, Stats
💻 2-8 CPU Cores • 4-8 GB RAM
conda install -c bioconda samtools htslib
Genomic Interval Arithmetic
💻 < 8GB RAM C++ • Interval Trees

bedtools

Docs ↗ By Aaron Quinlan Lab (Utah)
(License: GPL-2.0)

The Swiss-army knife for genomic interval arithmetic. Enables rapid intersection, merging, coverage calculation, subtracting, and windowing across BED, GFF, VCF, and BAM genome annotations.

Formats: Input: BED, VCF, GFF, BAM • Output: Filtered BED, Intersected intervals
💻 2-4 CPU Cores • 4-8 GB RAM
conda install -c bioconda bedtools
Genomic Browser & Alignment Inspection
💻 < 8GB RAM Java • Desktop & igv.js

Integrative Genomics Viewer (IGV)

Docs ↗ By James Robinson & Broad Institute
(License: MIT)

High-performance interactive visualization tool for exploring diverse large-scale genomic datasets (BAM pileups, VCF mutations, BigWig ChIP peaks, gene annotations) in smooth real time.

Formats: Input: BAM, CRAM, VCF, BigWig, BED, GTF • Output: High-res Snapshots
💻 2-4 CPU Cores • 4-8 GB RAM
conda install -c bioconda igv
Recommended Workflow Progression

Proceed to Stage 12: Targeted, Panel & Clinical Diagnostics

Sequencing pipelines proceed sequentially. Move forward to the next analytical domain to explore downstream algorithms and execution profiles.

Advance to Stage 12 →
← Stage 10: Population Genetics & Phylogenetics Directory Hub Stage 12: Clinical Genomics & Diagnostics →