A Nextflow pipeline for tumor–normal long-read whole-genome sequencing (WGS) analysis using PacBio HiFi data.
The workflow performs mapping, germline variant calling, haplotype phasing, somatic structural variant (SV) detection with multiple callers, and SV merging.
-
Mapping & QC
minimap2, samtools, pandepth -
Germline SV Calling
pbsv -
Germline SNV Calling
DeepVariant (PACBIO model) -
Phasing & Haplotagging
HiPhase -
Somatic SV Calling
- nanomonsv
- Severus
- SAVANA
- SVision-Pro
-
SV Merging
SURVIVOR (SVs supported by ≥2 callers)
- Tumor HiFi BAM
- Normal HiFi BAM
- Reference genome (T2T-CHM13)
- minimap2 index
- Tandem repeat annotation (for pbsv / Severus)
- Phased tumor and normal BAM files
- Germline SNV and SV VCFs
- Somatic SV VCFs from individual callers
- Merged somatic SV VCF (
merged.2caller.vcf)
nextflow run main.nf \
--normalinput normal.bam \
--tumorinput tumor.bam \
--reference chm13v2.0.fa \
--index minimap2_index