epi2me-labs / wf-somatic-variation

Other
15 stars 8 forks source link

Erroneous detection of basecaller model #34

Closed sabackman closed 2 months ago

sabackman commented 3 months ago

Operating System

Centos-7

Other Linux

No response

Workflow Version

v1.3.1

Workflow Execution

CLI Cluster

Other workflow execution

No response

EPI2ME Version

No response

CLI command run

nextflow run /proj/sens2023615/wf-somatic-variation-1.3.1 \ --bam_normal /path/to/normal.bam \ --bam_tumor /path/to/tumor.bam \ --ref /path/to/genome.fa \ --sv --snv --mod --sample_name SAMPLENAME \ --override_basecaller_cfg dna_r10.4.1_e8.2_400bps_sup@v4.2.0 \ --qv 20 \ --ubam_map_threads 12 \ --normal_min_coverage 0 \ --tumor_min_coverage 0 \ -profile samuel -resume

Workflow Execution - CLI Execution Profile

None

What happened?

I'm able to successfully run the workflow on some of my samples. For other samples, multiple basecaller models are detected, distinguished only by a suffix. An error is called, and the execution halted.

Relevant log output

Note that NXF_HOME is set to $HOME/.nextflow
Please change NXF_HOME to a place in your project directory (export NXF_HOME=yourprojectfolder)

Please use   export NXF_VER=nextflow_version    to select different version
https://www.nextflow.io/docs/latest/config.html#environment-variables

curl: (6) Could not resolve host: www.nextflow.io; Unknown error

 N E X T F L O W   ~  version 24.04.4

WARN: It appears you have never run this project before -- Option `-resume` is ignored
Launching `/proj/sens2023615/wf-somatic-variation-1.3.1/main.nf` [lonely_fermi] DSL2 - revision: adf816389f

WARN: Found unexpected parameters:
* --qv: 20
- Ignore this warning: params.schema_ignore_params = "qv" 

||||||||||   _____ ____ ___ ____  __  __ _____      _       _
||||||||||  | ____|  _ \_ _|___ \|  \/  | ____|    | | __ _| |__  ___
|||||       |  _| | |_) | |  __) | |\/| |  _| _____| |/ _` | '_ \/ __|
|||||       | |___|  __/| | / __/| |  | | |__|_____| | (_| | |_) \__ \
||||||||||  |_____|_|  |___|_____|_|  |_|_____|    |_|\__,_|_.__/|___/
||||||||||  wf-somatic-variation v1.3.1
--------------------------------------------------------------------------------
Core Nextflow options
  runName                : lonely_fermi
  containerEngine        : singularity
  container              : [withLabel:wf_somatic_snv:file:///proj/sens2023615/wf1.3.1_images/wf-somatic-snp.img, withLabel:wf_somatic_snv_to:file:///proj/sens2023615/wf1.3.1_images/clairs-to.img, withLabel:wf_somatic_mod:file:///proj/sens2023615/wf1.3.1_images/modkit.img, withLabel:dss:file:///proj/sens2023615/wf1.3.1_images/dss.img, withLabel:wf_somatic_sv:file:///proj/sens2023615/wf1.3.1_images/wf-somatic-sv.img, withLabel:wf_common:file:///proj/sens2023615/wf1.3.1_images/wf-common.img, withLabel:snpeff_annotation:file:///proj/sens2023615/wf1.3.1_images/snpeff.img, default:file:///proj/sens2023615/wf1.3.1_images/wf-somatic-variation.img]
  launchDir              : /proj/sens2023615/RUN/v1.3.1_GRCh38/SINET-7_forbugreport
  workDir                : /proj/sens2023615/RUN/v1.3.1_GRCh38/SINET-7_forbugreport/work
  projectDir             : /proj/sens2023615/wf-somatic-variation-1.3.1
  userName               : samuelba
  profile                : samuel
  configFiles            : /proj/sens2023615/wf-somatic-variation-1.3.1/nextflow.config

Workflow Options
  snv                    : true
  sv                     : true
  mod                    : true

Main options
  sample_name            : SINET_7
  bam_normal             : /proj/nobackup/sens2023615/wharf/fredra/fredra-sens2023615/op_089_007.bam
  bam_tumor              : /proj/nobackup/sens2023615/wharf/fredra/fredra-sens2023615/op_077_008.bam
  ref                    : /sw/data/igenomes/Homo_sapiens/UCSC/hg38/Sequence/WholeGenomeFasta/genome.fa

Quality Control Options
  tumor_min_coverage     : 0
  normal_min_coverage    : 0

Multiprocessing Options
  ubam_map_threads       : 12

Advanced Options
  override_basecaller_cfg: dna_r10.4.1_e8.2_400bps_sup@v4.2.0

!! Only displaying parameters that differ from the pipeline defaults !!
--------------------------------------------------------------------------------
If you use epi2me-labs/wf-somatic-variation for your analysis please cite:

* The nf-core framework
  https://doi.org/10.1038/s41587-020-0439-x

--------------------------------------------------------------------------------
This is epi2me-labs/wf-somatic-variation v1.3.1.
--------------------------------------------------------------------------------
[-        ] cram_cache         -
[-        ] getVersions        -
[-        ] getVersions_somvar -
[-        ] getParams          -

[-        ] cram_cache                     -
[-        ] getVersions                    -
[-        ] getVersions_somvar             -
[-        ] getParams                      -
[-        ] ingress_normal:checkBamHeaders -
Searching input for [.bam, .ubam] files.

[-        ] cram_cache                     -
[-        ] getVersions                    -
[-        ] getVersions_somvar             -
[-        ] getParams                      -
[-        ] ingress_normal:checkBamHeaders -
[-        ] ingress_normal:sortBam         -
[-        ] ingress_normal:mergeBams       -
[-        ] ingress_normal:catSortBams     -
[-        ] ingress_normal:validateIndex   -
[-        ] ingress_normal:samtools_index  -
Searching input for [.bam, .ubam] files.

[-        ] cram_cache                     -
[-        ] getVersions                    -
[-        ] getVersions_somvar             -
[-        ] getParams                      -
[-        ] ingress_normal:checkBamHeaders -
[-        ] ingress_normal:sortBam         -
[-        ] ingress_normal:mergeBams       -
[-        ] ingress_normal:catSortBams     -
[-        ] ingress_normal:validateIndex   -
[-        ] ingress_normal:samtools_index  -
[-        ] ing…normal:check_for_alignment -
[-        ] ing…_normal:minimap2_alignment -
Searching input for [.bam, .ubam] files.
Searching input for [.bam, .ubam] files.

[-        ] cram_cache                     -
[-        ] getVersions                    -
[-        ] getVersions_somvar             -
[-        ] getParams                      -
[-        ] ingress_normal:checkBamHeaders -
[-        ] ingress_normal:sortBam         -
[-        ] ingress_normal:mergeBams       -
[-        ] ingress_normal:catSortBams     -
[-        ] ingress_normal:validateIndex   -
[-        ] ingress_normal:samtools_index  -
[-        ] ing…normal:check_for_alignment -
[-        ] ing…_normal:minimap2_alignment -
[-        ] ingress_tumor:checkBamHeaders  -
[-        ] ingress_tumor:sortBam          -
[-        ] ingress_tumor:mergeBams        -
[-        ] ingress_tumor:catSortBams      -
[-        ] ingress_tumor:validateIndex    -
[-        ] ingress_tumor:samtools_index   -
[-        ] ing…_tumor:check_for_alignment -
Searching input for [.bam, .ubam] files.
Searching input for [.bam, .ubam] files.

[....]

Detected basecaller_model: dna_r10.4.1_e8.2_5khz_400bps_sup@v4.2.0
Detected basecaller_model: dna_r10.4.1_e8.2_5khz_400bps_sup@v4.2.0-102E7C0C
Detected basecaller_model: dna_r10.4.1_e8.2_5khz_400bps_sup@v4.2.0-10487036
Detected basecaller_model: dna_r10.4.1_e8.2_5khz_400bps_sup@v4.2.0-10738655
Detected basecaller_model: dna_r10.4.1_e8.2_5khz_400bps_sup@v4.2.0-114FBD25
Detected basecaller_model: dna_r10.4.1_e8.2_5khz_400bps_sup@v4.2.0-11ED0433
Detected basecaller_model: dna_r10.4.1_e8.2_5khz_400bps_sup@v4.2.0-124F6F05
Detected basecaller_model: dna_r10.4.1_e8.2_5khz_400bps_sup@v4.2.0-12673C44
Detected basecaller_model: dna_r10.4.1_e8.2_5khz_400bps_sup@v4.2.0-13CE7C29
Detected basecaller_model: dna_r10.4.1_e8.2_5khz_400bps_sup@v4.2.0-1479D126
Detected basecaller_model: dna_r10.4.1_e8.2_5khz_400bps_sup@v4.2.0-14918A8B
Detected basecaller_model: dna_r10.4.1_e8.2_5khz_400bps_sup@v4.2.0-153922DD
Detected basecaller_model: dna_r10.4.1_e8.2_5khz_400bps_sup@v4.2.0-155EC44B
Detected basecaller_model: dna_r10.4.1_e8.2_5khz_400bps_sup@v4.2.0-164B763
Detected basecaller_model: dna_r10.4.1_e8.2_5khz_400bps_sup@v4.2.0-16D2427F
Detected basecaller_model: dna_r10.4.1_e8.2_5khz_400bps_sup@v4.2.0-170119CD
Detected basecaller_model: dna_r10.4.1_e8.2_5khz_400bps_sup@v4.2.0-17BC2B6D
Detected basecaller_model: dna_r10.4.1_e8.2_5khz_400bps_sup@v4.2.0-17EBBAD5
Detected basecaller_model: dna_r10.4.1_e8.2_5khz_400bps_sup@v4.2.0-180A11F4
Detected basecaller_model: dna_r10.4.1_e8.2_5khz_400bps_sup@v4.2.0-1AB5602F
Detected basecaller_model: dna_r10.4.1_e8.2_5khz_400bps_sup@v4.2.0-1C03EB88
Detected basecaller_model: dna_r10.4.1_e8.2_5khz_400bps_sup@v4.2.0-1DC39BD
Detected basecaller_model: dna_r10.4.1_e8.2_5khz_400bps_sup@v4.2.0-1F41D19D
Detected basecaller_model: dna_r10.4.1_e8.2_5khz_400bps_sup@v4.2.0-1F972023
Detected basecaller_model: dna_r10.4.1_e8.2_5khz_400bps_sup@v4.2.0-1FDF9C0
Detected basecaller_model: dna_r10.4.1_e8.2_5khz_400bps_sup@v4.2.0-1FDFE777
Detected basecaller_model: dna_r10.4.1_e8.2_5khz_400bps_sup@v4.2.0-2198C901
Detected basecaller_model: dna_r10.4.1_e8.2_5khz_400bps_sup@v4.2.0-2329653D
Detected basecaller_model: dna_r10.4.1_e8.2_5khz_400bps_sup@v4.2.0-247F6304
Detected basecaller_model: dna_r10.4.1_e8.2_5khz_400bps_sup@v4.2.0-257B99EA
Detected basecaller_model: dna_r10.4.1_e8.2_5khz_400bps_sup@v4.2.0-26D6F255
Detected basecaller_model: dna_r10.4.1_e8.2_5khz_400bps_sup@v4.2.0-27AE60D2
Detected basecaller_model: dna_r10.4.1_e8.2_5khz_400bps_sup@v4.2.0-27BE9E92
Detected basecaller_model: dna_r10.4.1_e8.2_5khz_400bps_sup@v4.2.0-2980FD57
Detected basecaller_model: dna_r10.4.1_e8.2_5khz_400bps_sup@v4.2.0-2A9455DD
Detected basecaller_model: dna_r10.4.1_e8.2_5khz_400bps_sup@v4.2.0-2A9459E4
Detected basecaller_model: dna_r10.4.1_e8.2_5khz_400bps_sup@v4.2.0-2B27370F
Detected basecaller_model: dna_r10.4.1_e8.2_5khz_400bps_sup@v4.2.0-2CBFAB44
Detected basecaller_model: dna_r10.4.1_e8.2_5khz_400bps_sup@v4.2.0-2CF90782
Detected basecaller_model: dna_r10.4.1_e8.2_5khz_400bps_sup@v4.2.0-2DD011D1
Detected basecaller_model: dna_r10.4.1_e8.2_5khz_400bps_sup@v4.2.0-2FED2EDA
Detected basecaller_model: dna_r10.4.1_e8.2_5khz_400bps_sup@v4.2.0-304026A7
Detected basecaller_model: dna_r10.4.1_e8.2_5khz_400bps_sup@v4.2.0-31CB6B3E
Detected basecaller_model: dna_r10.4.1_e8.2_5khz_400bps_sup@v4.2.0-31E76A99
Detected basecaller_model: dna_r10.4.1_e8.2_5khz_400bps_sup@v4.2.0-32AF006F
Detected basecaller_model: dna_r10.4.1_e8.2_5khz_400bps_sup@v4.2.0-34060027
Detected basecaller_model: dna_r10.4.1_e8.2_5khz_400bps_sup@v4.2.0-34BD0F20
Detected basecaller_model: dna_r10.4.1_e8.2_5khz_400bps_sup@v4.2.0-350011DD
Detected basecaller_model: dna_r10.4.1_e8.2_5khz_400bps_sup@v4.2.0-361882AC
Detected basecaller_model: dna_r10.4.1_e8.2_5khz_400bps_sup@v4.2.0-361D2590
Detected basecaller_model: dna_r10.4.1_e8.2_5khz_400bps_sup@v4.2.0-36830487
Detected basecaller_model: dna_r10.4.1_e8.2_5khz_400bps_sup@v4.2.0-370737CC
Detected basecaller_model: dna_r10.4.1_e8.2_5khz_400bps_sup@v4.2.0-372369F7
Detected basecaller_model: dna_r10.4.1_e8.2_5khz_400bps_sup@v4.2.0-3866F59F
Detected basecaller_model: dna_r10.4.1_e8.2_5khz_400bps_sup@v4.2.0-3AB0BCE
Detected basecaller_model: dna_r10.4.1_e8.2_5khz_400bps_sup@v4.2.0-3AC8242E
Detected basecaller_model: dna_r10.4.1_e8.2_5khz_400bps_sup@v4.2.0-3C5AD182
Detected basecaller_model: dna_r10.4.1_e8.2_5khz_400bps_sup@v4.2.0-3CF1FB14
Detected basecaller_model: dna_r10.4.1_e8.2_5khz_400bps_sup@v4.2.0-3DE59B0D
Detected basecaller_model: dna_r10.4.1_e8.2_5khz_400bps_sup@v4.2.0-3E806201
Detected basecaller_model: dna_r10.4.1_e8.2_5khz_400bps_sup@v4.2.0-3E8B7F32
Detected basecaller_model: dna_r10.4.1_e8.2_5khz_400bps_sup@v4.2.0-3E98B38D
Detected basecaller_model: dna_r10.4.1_e8.2_5khz_400bps_sup@v4.2.0-3F181B47
Detected basecaller_model: dna_r10.4.1_e8.2_5khz_400bps_sup@v4.2.0-3FA6F9AE
Detected basecaller_model: dna_r10.4.1_e8.2_5khz_400bps_sup@v4.2.0-3FEEC44F
Detected basecaller_model: dna_r10.4.1_e8.2_5khz_400bps_sup@v4.2.0-4022B4AF
Detected basecaller_model: dna_r10.4.1_e8.2_5khz_400bps_sup@v4.2.0-42332D80
Detected basecaller_model: dna_r10.4.1_e8.2_5khz_400bps_sup@v4.2.0-429427DF
Detected basecaller_model: dna_r10.4.1_e8.2_5khz_400bps_sup@v4.2.0-42E0AB84
Detected basecaller_model: dna_r10.4.1_e8.2_5khz_400bps_sup@v4.2.0-4407C840
Detected basecaller_model: dna_r10.4.1_e8.2_5khz_400bps_sup@v4.2.0-44BC271F
Detected basecaller_model: dna_r10.4.1_e8.2_5khz_400bps_sup@v4.2.0-44DA4F9F
Detected basecaller_model: dna_r10.4.1_e8.2_5khz_400bps_sup@v4.2.0-44F5C266
Detected basecaller_model: dna_r10.4.1_e8.2_5khz_400bps_sup@v4.2.0-454E4757
Detected basecaller_model: dna_r10.4.1_e8.2_5khz_400bps_sup@v4.2.0-4569C49F
Detected basecaller_model: dna_r10.4.1_e8.2_5khz_400bps_sup@v4.2.0-46F9DBC1
Detected basecaller_model: dna_r10.4.1_e8.2_5khz_400bps_sup@v4.2.0-4874544C
Detected basecaller_model: dna_r10.4.1_e8.2_5khz_400bps_sup@v4.2.0-4A51586E
Detected basecaller_model: dna_r10.4.1_e8.2_5khz_400bps_sup@v4.2.0-4AC61680
Detected basecaller_model: dna_r10.4.1_e8.2_5khz_400bps_sup@v4.2.0-4C551A09
Detected basecaller_model: dna_r10.4.1_e8.2_5khz_400bps_sup@v4.2.0-4C8C12F1
Detected basecaller_model: dna_r10.4.1_e8.2_5khz_400bps_sup@v4.2.0-4CCDFD50
Detected basecaller_model: dna_r10.4.1_e8.2_5khz_400bps_sup@v4.2.0-4DD10682
Detected basecaller_model: dna_r10.4.1_e8.2_5khz_400bps_sup@v4.2.0-4E4BEA9E
Detected basecaller_model: dna_r10.4.1_e8.2_5khz_400bps_sup@v4.2.0-4F3DEB47
Detected basecaller_model: dna_r10.4.1_e8.2_5khz_400bps_sup@v4.2.0-50AAC6E7
Detected basecaller_model: dna_r10.4.1_e8.2_5khz_400bps_sup@v4.2.0-5115A0A9
Detected basecaller_model: dna_r10.4.1_e8.2_5khz_400bps_sup@v4.2.0-512B0C0F
Detected basecaller_model: dna_r10.4.1_e8.2_5khz_400bps_sup@v4.2.0-516ADC7E
Detected basecaller_model: dna_r10.4.1_e8.2_5khz_400bps_sup@v4.2.0-51823556
Detected basecaller_model: dna_r10.4.1_e8.2_5khz_400bps_sup@v4.2.0-5220ECF9
Detected basecaller_model: dna_r10.4.1_e8.2_5khz_400bps_sup@v4.2.0-5244064B
Detected basecaller_model: dna_r10.4.1_e8.2_5khz_400bps_sup@v4.2.0-530A5959
Detected basecaller_model: dna_r10.4.1_e8.2_5khz_400bps_sup@v4.2.0-5488FE2D
Detected basecaller_model: dna_r10.4.1_e8.2_5khz_400bps_sup@v4.2.0-5655F021
Detected basecaller_model: dna_r10.4.1_e8.2_5khz_400bps_sup@v4.2.0-56784563
Detected basecaller_model: dna_r10.4.1_e8.2_5khz_400bps_sup@v4.2.0-56A37152
Detected basecaller_model: dna_r10.4.1_e8.2_5khz_400bps_sup@v4.2.0-56EAE238
Detected basecaller_model: dna_r10.4.1_e8.2_5khz_400bps_sup@v4.2.0-57AA9EC
Detected basecaller_model: dna_r10.4.1_e8.2_5khz_400bps_sup@v4.2.0-57ECDD5F
Detected basecaller_model: dna_r10.4.1_e8.2_5khz_400bps_sup@v4.2.0-5C1F42B4
Detected basecaller_model: dna_r10.4.1_e8.2_5khz_400bps_sup@v4.2.0-5C208065
Detected basecaller_model: dna_r10.4.1_e8.2_5khz_400bps_sup@v4.2.0-5CF57153
Detected basecaller_model: dna_r10.4.1_e8.2_5khz_400bps_sup@v4.2.0-5D3C88AF
Detected basecaller_model: dna_r10.4.1_e8.2_5khz_400bps_sup@v4.2.0-5DA2DA64
Detected basecaller_model: dna_r10.4.1_e8.2_5khz_400bps_sup@v4.2.0-5E644304
Detected basecaller_model: dna_r10.4.1_e8.2_5khz_400bps_sup@v4.2.0-5E941310
Detected basecaller_model: dna_r10.4.1_e8.2_5khz_400bps_sup@v4.2.0-5F4E658E
Detected basecaller_model: dna_r10.4.1_e8.2_5khz_400bps_sup@v4.2.0-6025DD57
Detected basecaller_model: dna_r10.4.1_e8.2_5khz_400bps_sup@v4.2.0-60400326
Detected basecaller_model: dna_r10.4.1_e8.2_5khz_400bps_sup@v4.2.0-60D74565
Detected basecaller_model: dna_r10.4.1_e8.2_5khz_400bps_sup@v4.2.0-60DA8170
Detected basecaller_model: dna_r10.4.1_e8.2_5khz_400bps_sup@v4.2.0-61BF7FE9
Detected basecaller_model: dna_r10.4.1_e8.2_5khz_400bps_sup@v4.2.0-625CA63F
Detected basecaller_model: dna_r10.4.1_e8.2_5khz_400bps_sup@v4.2.0-62F72577
WARN: Overriding basecaller_model: dna_r10.4.1_e8.2_400bps_sup@v4.2.0
WARN: Overriding basecaller_model: dna_r10.4.1_e8.2_400bps_sup@v4.2.0
WARN: Overriding basecaller_model: dna_r10.4.1_e8.2_400bps_sup@v4.2.0
WARN: Overriding basecaller_model: dna_r10.4.1_e8.2_400bps_sup@v4.2.0
WARN: Overriding basecaller_model: dna_r10.4.1_e8.2_400bps_sup@v4.2.0
WARN: Overriding basecaller_model: dna_r10.4.1_e8.2_400bps_sup@v4.2.0
WARN: Overriding basecaller_model: dna_r10.4.1_e8.2_400bps_sup@v4.2.0
WARN: Overriding basecaller_model: dna_r10.4.1_e8.2_400bps_sup@v4.2.0
WARN: Overriding basecaller_model: dna_r10.4.1_e8.2_400bps_sup@v4.2.0
WARN: Overriding basecaller_model: dna_r10.4.1_e8.2_400bps_sup@v4.2.0
WARN: Overriding basecaller_model: dna_r10.4.1_e8.2_400bps_sup@v4.2.0
WARN: Overriding basecaller_model: dna_r10.4.1_e8.2_400bps_sup@v4.2.0
WARN: Overriding basecaller_model: dna_r10.4.1_e8.2_400bps_sup@v4.2.0
WARN: Overriding basecaller_model: dna_r10.4.1_e8.2_400bps_sup@v4.2.0
WARN: Overriding basecaller_model: dna_r10.4.1_e8.2_400bps_sup@v4.2.0
WARN: Overriding basecaller_model: dna_r10.4.1_e8.2_400bps_sup@v4.2.0
WARN: Overriding basecaller_model: dna_r10.4.1_e8.2_400bps_sup@v4.2.0
WARN: Overriding basecaller_model: dna_r10.4.1_e8.2_400bps_sup@v4.2.0
WARN: Overriding basecaller_model: dna_r10.4.1_e8.2_400bps_sup@v4.2.0
WARN: Overriding basecaller_model: dna_r10.4.1_e8.2_400bps_sup@v4.2.0
WARN: Overriding basecaller_model: dna_r10.4.1_e8.2_400bps_sup@v4.2.0
WARN: Overriding basecaller_model: dna_r10.4.1_e8.2_400bps_sup@v4.2.0
WARN: Overriding basecaller_model: dna_r10.4.1_e8.2_400bps_sup@v4.2.0
WARN: Overriding basecaller_model: dna_r10.4.1_e8.2_400bps_sup@v4.2.0
WARN: Overriding basecaller_model: dna_r10.4.1_e8.2_400bps_sup@v4.2.0
WARN: Overriding basecaller_model: dna_r10.4.1_e8.2_400bps_sup@v4.2.0
WARN: Overriding basecaller_model: dna_r10.4.1_e8.2_400bps_sup@v4.2.0
WARN: Overriding basecaller_model: dna_r10.4.1_e8.2_400bps_sup@v4.2.0
WARN: Overriding basecaller_model: dna_r10.4.1_e8.2_400bps_sup@v4.2.0
WARN: Overriding basecaller_model: dna_r10.4.1_e8.2_400bps_sup@v4.2.0
WARN: Overriding basecaller_model: dna_r10.4.1_e8.2_400bps_sup@v4.2.0
WARN: Overriding basecaller_model: dna_r10.4.1_e8.2_400bps_sup@v4.2.0
WARN: Overriding basecaller_model: dna_r10.4.1_e8.2_400bps_sup@v4.2.0
WARN: Overriding basecaller_model: dna_r10.4.1_e8.2_400bps_sup@v4.2.0
WARN: Overriding basecaller_model: dna_r10.4.1_e8.2_400bps_sup@v4.2.0
WARN: Overriding basecaller_model: dna_r10.4.1_e8.2_400bps_sup@v4.2.0
WARN: Overriding basecaller_model: dna_r10.4.1_e8.2_400bps_sup@v4.2.0
WARN: Overriding basecaller_model: dna_r10.4.1_e8.2_400bps_sup@v4.2.0
WARN: Overriding basecaller_model: dna_r10.4.1_e8.2_400bps_sup@v4.2.0
WARN: Overriding basecaller_model: dna_r10.4.1_e8.2_400bps_sup@v4.2.0
WARN: Overriding basecaller_model: dna_r10.4.1_e8.2_400bps_sup@v4.2.0
WARN: Overriding basecaller_model: dna_r10.4.1_e8.2_400bps_sup@v4.2.0
WARN: Overriding basecaller_model: dna_r10.4.1_e8.2_400bps_sup@v4.2.0
WARN: Overriding basecaller_model: dna_r10.4.1_e8.2_400bps_sup@v4.2.0
WARN: Overriding basecaller_model: dna_r10.4.1_e8.2_400bps_sup@v4.2.0
WARN: Overriding basecaller_model: dna_r10.4.1_e8.2_400bps_sup@v4.2.0
WARN: Overriding basecaller_model: dna_r10.4.1_e8.2_400bps_sup@v4.2.0
WARN: Overriding basecaller_model: dna_r10.4.1_e8.2_400bps_sup@v4.2.0
WARN: Overriding basecaller_model: dna_r10.4.1_e8.2_400bps_sup@v4.2.0
WARN: Overriding basecaller_model: dna_r10.4.1_e8.2_400bps_sup@v4.2.0
WARN: Overriding basecaller_model: dna_r10.4.1_e8.2_400bps_sup@v4.2.0
WARN: Overriding basecaller_model: dna_r10.4.1_e8.2_400bps_sup@v4.2.0
WARN: Overriding basecaller_model: dna_r10.4.1_e8.2_400bps_sup@v4.2.0
WARN: Overriding basecaller_model: dna_r10.4.1_e8.2_400bps_sup@v4.2.0
WARN: Overriding basecaller_model: dna_r10.4.1_e8.2_400bps_sup@v4.2.0
WARN: Overriding basecaller_model: dna_r10.4.1_e8.2_400bps_sup@v4.2.0
WARN: Overriding basecaller_model: dna_r10.4.1_e8.2_400bps_sup@v4.2.0
WARN: Overriding basecaller_model: dna_r10.4.1_e8.2_400bps_sup@v4.2.0
WARN: Overriding basecaller_model: dna_r10.4.1_e8.2_400bps_sup@v4.2.0
WARN: Overriding basecaller_model: dna_r10.4.1_e8.2_400bps_sup@v4.2.0
WARN: Overriding basecaller_model: dna_r10.4.1_e8.2_400bps_sup@v4.2.0
WARN: Overriding basecaller_model: dna_r10.4.1_e8.2_400bps_sup@v4.2.0
WARN: Overriding basecaller_model: dna_r10.4.1_e8.2_400bps_sup@v4.2.0
WARN: Overriding basecaller_model: dna_r10.4.1_e8.2_400bps_sup@v4.2.0
WARN: Overriding basecaller_model: dna_r10.4.1_e8.2_400bps_sup@v4.2.0
WARN: Overriding basecaller_model: dna_r10.4.1_e8.2_400bps_sup@v4.2.0
WARN: Overriding basecaller_model: dna_r10.4.1_e8.2_400bps_sup@v4.2.0
WARN: Overriding basecaller_model: dna_r10.4.1_e8.2_400bps_sup@v4.2.0
WARN: Overriding basecaller_model: dna_r10.4.1_e8.2_400bps_sup@v4.2.0
WARN: Overriding basecaller_model: dna_r10.4.1_e8.2_400bps_sup@v4.2.0
WARN: Overriding basecaller_model: dna_r10.4.1_e8.2_400bps_sup@v4.2.0
WARN: Overriding basecaller_model: dna_r10.4.1_e8.2_400bps_sup@v4.2.0
WARN: Overriding basecaller_model: dna_r10.4.1_e8.2_400bps_sup@v4.2.0
WARN: Overriding basecaller_model: dna_r10.4.1_e8.2_400bps_sup@v4.2.0
WARN: Overriding basecaller_model: dna_r10.4.1_e8.2_400bps_sup@v4.2.0
WARN: Overriding basecaller_model: dna_r10.4.1_e8.2_400bps_sup@v4.2.0
WARN: Overriding basecaller_model: dna_r10.4.1_e8.2_400bps_sup@v4.2.0
WARN: Overriding basecaller_model: dna_r10.4.1_e8.2_400bps_sup@v4.2.0
WARN: Overriding basecaller_model: dna_r10.4.1_e8.2_400bps_sup@v4.2.0
WARN: Overriding basecaller_model: dna_r10.4.1_e8.2_400bps_sup@v4.2.0
WARN: Overriding basecaller_model: dna_r10.4.1_e8.2_400bps_sup@v4.2.0
WARN: Overriding basecaller_model: dna_r10.4.1_e8.2_400bps_sup@v4.2.0
WARN: Overriding basecaller_model: dna_r10.4.1_e8.2_400bps_sup@v4.2.0
WARN: Overriding basecaller_model: dna_r10.4.1_e8.2_400bps_sup@v4.2.0
WARN: Overriding basecaller_model: dna_r10.4.1_e8.2_400bps_sup@v4.2.0
WARN: Overriding basecaller_model: dna_r10.4.1_e8.2_400bps_sup@v4.2.0
WARN: Overriding basecaller_model: dna_r10.4.1_e8.2_400bps_sup@v4.2.0
WARN: Overriding basecaller_model: dna_r10.4.1_e8.2_400bps_sup@v4.2.0
WARN: Overriding basecaller_model: dna_r10.4.1_e8.2_400bps_sup@v4.2.0
WARN: Overriding basecaller_model: dna_r10.4.1_e8.2_400bps_sup@v4.2.0
WARN: Overriding basecaller_model: dna_r10.4.1_e8.2_400bps_sup@v4.2.0
WARN: Overriding basecaller_model: dna_r10.4.1_e8.2_400bps_sup@v4.2.0
WARN: Overriding basecaller_model: dna_r10.4.1_e8.2_400bps_sup@v4.2.0
WARN: Overriding basecaller_model: dna_r10.4.1_e8.2_400bps_sup@v4.2.0
WARN: Overriding basecaller_model: dna_r10.4.1_e8.2_400bps_sup@v4.2.0
WARN: Overriding basecaller_model: dna_r10.4.1_e8.2_400bps_sup@v4.2.0
WARN: Overriding basecaller_model: dna_r10.4.1_e8.2_400bps_sup@v4.2.0
WARN: Overriding basecaller_model: dna_r10.4.1_e8.2_400bps_sup@v4.2.0
WARN: Overriding basecaller_model: dna_r10.4.1_e8.2_400bps_sup@v4.2.0
WARN: Overriding basecaller_model: dna_r10.4.1_e8.2_400bps_sup@v4.2.0
WARN: Overriding basecaller_model: dna_r10.4.1_e8.2_400bps_sup@v4.2.0
WARN: Overriding basecaller_model: dna_r10.4.1_e8.2_400bps_sup@v4.2.0
WARN: Overriding basecaller_model: dna_r10.4.1_e8.2_400bps_sup@v4.2.0
WARN: Overriding basecaller_model: dna_r10.4.1_e8.2_400bps_sup@v4.2.0
WARN: Overriding basecaller_model: dna_r10.4.1_e8.2_400bps_sup@v4.2.0
WARN: Overriding basecaller_model: dna_r10.4.1_e8.2_400bps_sup@v4.2.0
WARN: Overriding basecaller_model: dna_r10.4.1_e8.2_400bps_sup@v4.2.0
WARN: Overriding basecaller_model: dna_r10.4.1_e8.2_400bps_sup@v4.2.0
WARN: Overriding basecaller_model: dna_r10.4.1_e8.2_400bps_sup@v4.2.0
WARN: Overriding basecaller_model: dna_r10.4.1_e8.2_400bps_sup@v4.2.0
WARN: Overriding basecaller_model: dna_r10.4.1_e8.2_400bps_sup@v4.2.0
WARN: Overriding basecaller_model: dna_r10.4.1_e8.2_400bps_sup@v4.2.0
WARN: Overriding basecaller_model: dna_r10.4.1_e8.2_400bps_sup@v4.2.0
WARN: Overriding basecaller_model: dna_r10.4.1_e8.2_400bps_sup@v4.2.0

[ ... ]

WARN: Killing running tasks (1)
executor >  slurm (61)
[68/f0e205] cram_cache (1)                 | 1 of 1 ✔
[fb/0b8134] getVersions                    | 1 of 1 ✔
[6b/f9da35] getVersions_somvar             | 1 of 1 ✔
[e3/2830b8] getParams                      | 1 of 1 ✔
[96/0b9941] ing…normal:checkBamHeaders (1) | 1 of 1 ✔
[-        ] ingress_normal:sortBam         -
[-        ] ingress_normal:mergeBams       -
[-        ] ingress_normal:catSortBams     -
[-        ] ingress_normal:validateIndex   -
[b4/66c426] ing…_normal:samtools_index (1) | 1 of 1 ✔
[14/cb2bfa] ing…al:check_for_alignment (1) | 1 of 1 ✔
[5f/b28c5d] ing…mal:minimap2_alignment (1) | 1 of 1 ✔
[e6/e2ed89] ing…_tumor:checkBamHeaders (1) | 1 of 1 ✔
[-        ] ingress_tumor:sortBam          -
[-        ] ingress_tumor:mergeBams        -
[-        ] ingress_tumor:catSortBams      -
[-        ] ingress_tumor:validateIndex    -
[65/eba577] ing…s_tumor:samtools_index (1) | 1 of 1 ✔
[75/f2502c] ing…or:check_for_alignment (1) | 1 of 1 ✔
[95/d994e2] ing…mor:minimap2_alignment (1) | 1 of 1 ✔
[13/770ff3] getGenome (2)                  | 2 of 2 ✔
[2c/604090] getAllChromosomesBed (1)       | 1 of 1 ✔
[22/4edbde] alignment_stats:bamstats (2)   | 2 of 2 ✔
[55/d6f57b] alignment_stats:mosdepth (2)   | 2 of 2 ✔
[-        ] alignment_stats:makeQCreport   | 0 of 1
[fe/dfc28b] publish_qc (2)                 | 1 of 14
[-        ] snv:wf_build_regions           | 0 of 1
[13/a9e3d7] snv:getVersions                | 1 of 1 ✔
[62/6e26d6] snv:getParams                  | 1 of 1 ✔
[-        ] sv:severus                     | 0 of 1
[99/1b8785] sv:getVersions                 | 1 of 1 ✔
[77/b1e656] sv:getParams                   | 1 of 1 ✔
[14/5ecf9f] mod:validate_modbam (1)        | 1 of 2
[0b/0f3191] mod:sample_probs (1)           | 1 of 1
[de/c73b91] mod:modkit (23)                | 25 of 25
[ff/252ad6] mod:summary (1)                | 1 of 2
[c7/7ba38d] mod:getVersions                | 1 of 1 ✔
[5d/3f2293] mod:rVersions                  | 1 of 1 ✔
[da/bbc665] mod:getParams                  | 1 of 1 ✔
[05/8eb124] mod:publish_modbase (3)        | 3 of 3
[cc/9cbf9d] publish (2)                    | 2 of 2
Plus 43 more processes waiting for tasks…
Detected basecaller_model: dna_r10.4.1_e8.2_400bps_sup@v4.2.0
WARN: Overriding basecaller_model: dna_r10.4.1_e8.2_400bps_sup@v4.2.0
################################################################################
# INPUT DATA PROBLEM
Your input tumor and/or normal BAM files indicate two different basecall models
in the header.

wf-somatic-variation requires the basecall model used to call the tumor and 
normal samples to be identical.

## Next steps
You must re-run the workflow specifying the same basecaller model with the
--override_basecaller_cfg option.
################################################################################

Application activity log entry

No response

Were you able to successfully run the latest version of the workflow with the demo data?

yes

Other demo data information

Did not try, but successful runs with other data.
RenzoTale88 commented 3 months ago

Hi @sabackman thank you for opening a new issue with this and providing all the details. I can see you are running with a custom profile, would you mind confirming whether you are running in docker/singularity mode, and with which containers?

Andrea

sabackman commented 3 months ago

I'm running in singularity mode.

The container images were built using the following commands:

singularity build modkit.img docker://ontresearch/modkit:shae137452eb41f5f12b790774cafe15bc97f48d4d0
singularity build wf-common.img docker://ontresearch/wf-common:shad399cf22079b5b153920ac39ee40095a677933f1
singularity build wf-somatic-variation.img docker://ontresearch/wf-somatic-variation:sha18cc2ea1fae27fc772e7b728957996119c1ec81a
singularity build clairs-to.img docker://ontresearch/clairs-to:sha0bed006dfad374401e9410b44d67e6609f7936e9
singularity build wf-somatic-snp.img docker://ontresearch/wf-somatic-snp:sha86599ecdb46472d01e0b7c2704683f6aee8a82dd
singularity build dss.img docker://ontresearch/dss:sha28cc44cacad092f1b074aca3548ff8ba9b483fc0
singularity build wf-somatic-sv.img docker://ontresearch/wf-somatic-sv:shaf318a490c2625b8b8e4f2bc70b4dc546ce2d5ffc
singularity build snpeff.img docker://ontresearch/snpeff:sha313729d816872d70b410752001a0e2654eb88807
RenzoTale88 commented 3 months ago

I can see where the issue is, apologies for not spotting this immediately. The two models, while being the same, they have slightly different naming:

dna_r10.4.1_e8.2_5khz_400bps_sup@v4.2.0
dna_r10.4.1_e8.2_400bps_sup@v4.2.0

Which tricks the workflow into thinking something is off. We had a similar issue in the past in wf-human-variation.

I'll open a ticket to add this model to the list of acceptable models for the next release, apologies for the inconvenience.

RenzoTale88 commented 2 months ago

@sabackman we just released v1.4.0 that should address your issue. Now --override_basecaller_cfg will enforce the same basecaller model even if the two BAM have been basecalled with different models. To update simply run:

nextflow pull epi2me-labs/wf-somatic-variation -r v1.4.0

I'll close this, but please re-open if you still come across issues.

Andrea