feat: add holodeck modules (simulate, mutate, methylate, eval)#12176
Merged
Conversation
atrigila
requested changes
Jul 2, 2026
Comment on lines
+46
to
+54
| "##fileformat=VCFv4.4", | ||
| "##holodeckVersion=0.3.0", | ||
| "##FORMAT=<ID=GT,Number=1,Type=String,Description=\"Genotype\">", | ||
| "##FORMAT=<ID=MT,Number=.,Type=String,Description=\"Methylation state, top strand. Per-haplotype pipe-separated bitstring: 1=methylated, 0=unmethylated, .=haplotype carries REF or no owned CpG.\">", | ||
| "##FORMAT=<ID=MB,Number=.,Type=String,Description=\"Methylation state, bottom strand. Per-haplotype pipe-separated bitstring: 1=methylated, 0=unmethylated, .=haplotype carries REF or no owned CpG.\">", | ||
| "##contig=<ID=MT192765.1,length=29829>", | ||
| "#CHROM\tPOS\tID\tREF\tALT\tQUAL\tFILTER\tINFO\tFORMAT\tMETHYLATE", | ||
| "MT192765.1\t37\t.\tC\t.\t.\t.\t.\tMT:MB\t0|0:0|0", | ||
| "MT192765.1\t64\t.\tC\t.\t.\t.\t.\tMT:MB\t0|0:0|0", |
Member
Author
There was a problem hiding this comment.
This block is the decompressed VCF body — one record per CpG site with the simulated methylation state (MT/MB FORMAT) — which is deterministic for a fixed input, so it snapshots stably. I've now also filtered the ##holodeckVersion header line (alongside ##holodeckCommand) so the snapshot won't churn on holodeck releases. Happy to adjust if you had a different concern in mind.
Add an nf-core module wrapping `holodeck simulate`. Simulate sequencing reads from a reference genome. Emits paired-end FastQ plus optional ground-truth golden BAM/VCF, with an optional VCF input to embed variants into the reads. holodeck is a modern NGS read simulator written in Rust. The module is tested against the sarscov2 reference and includes a stub run.
Add an nf-core module wrapping `holodeck mutate`. Generate a VCF of random mutations from a reference genome. holodeck is a modern NGS read simulator written in Rust. The module is tested against the sarscov2 reference and includes a stub run.
Add an nf-core module wrapping `holodeck methylate`. Generate a methylation-annotated VCF (MT/MB FORMAT fields) from a reference genome, optionally annotating an input VCF of phased variants. holodeck is a modern NGS read simulator written in Rust. The module is tested against the sarscov2 reference and includes a stub run.
Add an nf-core module wrapping `holodeck eval`. Evaluate alignment placement accuracy of holodeck-simulated reads, comparing mapped positions against the simulated truth (golden BAM or encoded read names). holodeck is a modern NGS read simulator written in Rust. The module is tested against the sarscov2 reference and includes a stub run.
Member
Author
|
@atrigila thanks for suggestions, any chance you can take another look? |
atrigila
approved these changes
Jul 4, 2026
atrigila
left a comment
Contributor
There was a problem hiding this comment.
Thanks for addressing the comments
This file contains hidden or bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
Learn more about bidirectional Unicode characters
Sign up for free
to join this conversation on GitHub.
Already have an account?
Sign in to comment
Add this suggestion to a batch that can be applied as a single commit.This suggestion is invalid because no changes were made to the code.Suggestions cannot be applied while the pull request is closed.Suggestions cannot be applied while viewing a subset of changes.Only one suggestion per line can be applied in a batch.Add this suggestion to a batch that can be applied as a single commit.Applying suggestions on deleted lines is not supported.You must change the existing code in this line in order to create a valid suggestion.Outdated suggestions cannot be applied.This suggestion has been applied or marked resolved.Suggestions cannot be applied from pending reviews.Suggestions cannot be applied on multi-line comments.Suggestions cannot be applied while the pull request is queued to merge.Suggestion cannot be applied right now. Please check back later.
Description
Adds four new modules for holodeck, a modern NGS read simulator written in Rust:
holodeck/simulate— simulate paired-end reads from a reference (optional variants VCF, optional ground-truth golden BAM/VCF).holodeck/mutate— generate a VCF of random mutations from a reference.holodeck/methylate— generate a methylation-annotated VCF (MT/MB FORMAT fields).holodeck/eval— evaluate alignment placement accuracy of simulated reads.All tested against the sarscov2 reference with stub runs.
PR checklist
label.bioconda::holodeck=0.3.0; Seqera Wave community containers).nf-core modules test holodeck/<sub> --profile dockerpasses for all four.topic: versions.