-
Notifications
You must be signed in to change notification settings - Fork 0
Workflow
This document describes the analytical workflow implemented in PoODLE.
The pipeline processes sequencing data using a series of modular steps implemented in Nextflow DSL2. Each step is executed independently and automatically parallelized when possible.
The workflow consists of several stages that transform sequencing data into interpretable results for genomic surveillance and linkage investigation.

-
Reference-based SNP calling (Snippy)
- Runs
Snippyper sample within each cluster - Automatically detects if results already exist
- Reuses previous results when the same reference was used → This is critical for surveillance workflows where clusters grow over time
- Runs
-
Core genome alignment & SNP distances
Recombination events can distort phylogenetic inference. When enabled, recombinant regions are identified and removed.
-
Gubbinsmasks recombinant regions - SNP-only alignment is rebuilt with
snp-sites - Phylogenetic tree and SNP distances are generated
Gene presence and absence across genomes is analyzed with Panaroo to characterize genomic diversity. Hamming distances based on gene presence are calculated as a metric to compare genome similarity including accessory genes.
A rapid whole-genome comparison can be performed using assembly-based distance estimation with MashTree.
Results from all analysis steps are summarized into an interactive report made with Quarto and several R packages including phylocanvas, plottly and tidyverse.
Outputs include:
- cluster-specific HTML reports
- phylogenetic trees
- SNP and gene distance matrices
- pangenome visualizations
- linkage classification tables
These reports are designed for epidemiological interpretation and outbreak investigation.