| 123456789101112131415161718192021222324252627282930313233343536373839404142434445464748495051525354555657585960616263646566676869707172737475767778798081828384858687888990919293949596979899100101102103104105106107108109110111112113114115116117118119120121122123124125126127128129130131132133134135136137138 |
- WARNING: While bind mounting '/mnt/beegfs02/scratch/t_steimle/data/wgs/CHAHA/diag/ClairS/part4:/mnt/beegfs02/scratch/t_steimle/data/wgs/CHAHA/diag/ClairS/part4': destination is already in the mount point list
- [COMMAND] /opt/bin/run_clairs --tumor_bam_fn /mnt/beegfs02/scratch/t_steimle/data/wgs/CHAHA/diag/CHAHA_diag_hs1.bam --normal_bam_fn /mnt/beegfs02/scratch/t_steimle/data/wgs/CHAHA/norm/CHAHA_norm_hs1.bam --ref_fn /home/t_steimle/ref/hs1/chm13v2.0.fa --threads 40 --platform ont_r10_dorado_sup_5khz_ssrs --output_dir /mnt/beegfs02/scratch/t_steimle/data/wgs/CHAHA/diag/ClairS/part4 --region chr10:63341881-134758134 --sample_name CHAHA_diag --include_all_ctgs --print_germline_calls --enable_indel_calling --enable_clair3_germline_output --use_longphase_for_intermediate_haplotagging True
- [INFO] Create folder /mnt/beegfs02/scratch/t_steimle/data/wgs/CHAHA/diag/ClairS/part4/logs
- [INFO] Create folder /mnt/beegfs02/scratch/t_steimle/data/wgs/CHAHA/diag/ClairS/part4/tmp
- [INFO] Create folder /mnt/beegfs02/scratch/t_steimle/data/wgs/CHAHA/diag/ClairS/part4/tmp/split_beds
- [INFO] Create folder /mnt/beegfs02/scratch/t_steimle/data/wgs/CHAHA/diag/ClairS/part4/tmp/candidates
- [INFO] Create folder /mnt/beegfs02/scratch/t_steimle/data/wgs/CHAHA/diag/ClairS/part4/tmp/pileup_tensor_can
- [INFO] Create folder /mnt/beegfs02/scratch/t_steimle/data/wgs/CHAHA/diag/ClairS/part4/tmp/fa_tensor_can
- [INFO] Create folder /mnt/beegfs02/scratch/t_steimle/data/wgs/CHAHA/diag/ClairS/part4/tmp/vcf_output
- [INFO] Create folder /mnt/beegfs02/scratch/t_steimle/data/wgs/CHAHA/diag/ClairS/part4/tmp/tmp_vcf_output
- [INFO] Create folder /mnt/beegfs02/scratch/t_steimle/data/wgs/CHAHA/diag/ClairS/part4/logs/clair3_log
- [INFO] Create folder /mnt/beegfs02/scratch/t_steimle/data/wgs/CHAHA/diag/ClairS/part4/tmp/clair3_output/phased_output
- [INFO] Create folder /mnt/beegfs02/scratch/t_steimle/data/wgs/CHAHA/diag/ClairS/part4/tmp/clair3_output/vcf
- [INFO] --include_all_ctgs enabled
- [INFO] Call variants in contigs: chr10
- [INFO] Number of chunks for each contig: 27
- [INFO] CALLER VERSION: 0.4.4
- [INFO] NORMAL BAM FILE PATH: /mnt/beegfs02/scratch/t_steimle/data/wgs/CHAHA/norm/CHAHA_norm_hs1.bam
- [INFO] TUMOR BAM FILE PATH: /mnt/beegfs02/scratch/t_steimle/data/wgs/CHAHA/diag/CHAHA_diag_hs1.bam
- [INFO] REFERENCE FILE PATH: /home/t_steimle/ref/hs1/chm13v2.0.fa
- [INFO] PLATFORM: ont_r10_dorado_sup_5khz_ssrs
- [INFO] THREADS: 40
- [INFO] OUTPUT FOLDER: /mnt/beegfs02/scratch/t_steimle/data/wgs/CHAHA/diag/ClairS/part4
- [INFO] OUTPUT VCF PATH: /mnt/beegfs02/scratch/t_steimle/data/wgs/CHAHA/diag/ClairS/part4/output.vcf.gz
- [INFO] PILEUP MODEL PATH: /opt/conda/envs/clairs/bin/clairs_models/ont_r10_dorado_sup_5khz_ssrs/pileup.pkl
- [INFO] FULL-ALIGNMENT MODEL PATH: /opt/conda/envs/clairs/bin/clairs_models/ont_r10_dorado_sup_5khz_ssrs/full_alignment.pkl
- [INFO] BED FILE PATH: /mnt/beegfs02/scratch/t_steimle/data/wgs/CHAHA/diag/ClairS/part4/tmp/region.bed
- [INFO] GENOTYPING MODE VCF FILE PATH: None
- [INFO] HYBRID MODE VCF FILE PATH: None
- [INFO] REGION FOR CALLING: chr10:63341881-134758134
- [INFO] CONTIGS FOR CALLING: None
- [INFO] CONDA BINARY PREFIX: /opt/conda/envs/clairs
- [INFO] SAMTOOLS BINARY PATH: samtools
- [INFO] PYTHON BINARY PATH: python3
- [INFO] PYPY BINARY PATH: pypy3
- [INFO] PARALLEL BINARY PATH: parallel
- [INFO] LONGPHASE BINARY PATH: /opt/conda/envs/clairs/bin/longphase
- [INFO] CHUNK SIZE: 5000000
- [INFO] SNV MINIMUM AF: 0.05
- [INFO] SNV MINIMUM QUAL: 8
- [INFO] INDEL MINIMUM AF: 0.1
- [INFO] INDEL PILEUP MODEL PATH: /opt/conda/envs/clairs/bin/clairs_models/ont_r10_dorado_sup_5khz_ssrs/indel/pileup.pkl
- [INFO] INDEL FULL-ALIGNMENT MODEL PATH: /opt/conda/envs/clairs/bin/clairs_models/ont_r10_dorado_sup_5khz_ssrs/indel/full_alignment.pkl
- [INFO] INDEL MINIMUM QUAL: 8
- [INFO] NORMAL VCF FILE PATH: None
- [INFO] DISABLE PHASING: False
- [INFO] ENABLE DRY RUN: False
- [INFO] ENABLE INDEL CALLING: True
- [INFO] ENABLE PRINTING REFERENCE CALLS: False
- [INFO] ENABLE PRINTING GERMLINE CALLS: True
- [INFO] ENABLE INCLUDING ALL CTGS FOR CALLING: True
- [INFO] ENABLE REMOVING INTERMEDIATE FILES: False
- [INFO] Call Germline Variants in Normal BAM using Clair3
- [INFO] RUN THE FOLLOWING COMMAND:
- ( /opt/conda/envs/clairs/bin/run_clair3.sh --bam_fn /mnt/beegfs02/scratch/t_steimle/data/wgs/CHAHA/norm/CHAHA_norm_hs1.bam --ref_fn /home/t_steimle/ref/hs1/chm13v2.0.fa --model_path /opt/conda/envs/clairs/bin/clairs_models/clair3_models/r1041_e82_400bps_sup_v420 --platform ont --threads 40 --output /mnt/beegfs02/scratch/t_steimle/data/wgs/CHAHA/diag/ClairS/part4/tmp/clair3_output/clair3_normal_output --ctg_name=chr10 --samtools=samtools --pypy=pypy3 --python=python3 --min_coverage=4 --snp_min_af=0.08 --indel_min_af=0.15 --longphase_for_phasing ) 2>&1 | tee /mnt/beegfs02/scratch/t_steimle/data/wgs/CHAHA/diag/ClairS/part4/logs/clair3_log/1_CLAIR3_NORMAL.log
- [INFO] CLAIR3 VERSION: v1.0.8
- [INFO] BAM FILE PATH: /mnt/beegfs02/scratch/t_steimle/data/wgs/CHAHA/norm/CHAHA_norm_hs1.bam
- [INFO] REFERENCE FILE PATH: /home/t_steimle/ref/hs1/chm13v2.0.fa
- [INFO] MODEL PATH: /opt/conda/envs/clairs/bin/clairs_models/clair3_models/r1041_e82_400bps_sup_v420
- [INFO] OUTPUT FOLDER: /mnt/beegfs02/scratch/t_steimle/data/wgs/CHAHA/diag/ClairS/part4/tmp/clair3_output/clair3_normal_output
- [INFO] PLATFORM: ont
- [INFO] THREADS: 40
- [INFO] BED FILE PATH: EMPTY
- [INFO] VCF FILE PATH: EMPTY
- [INFO] CONTIGS: chr10
- [INFO] CONDA PREFIX:
- [INFO] SAMTOOLS PATH: samtools
- [INFO] PYTHON PATH: python3
- [INFO] PYPY PATH: pypy3
- [INFO] PARALLEL PATH: parallel
- [INFO] WHATSHAP PATH: whatshap
- [INFO] LONGPHASE PATH: /opt/conda/envs/clairs/bin/longphase
- [INFO] CHUNK SIZE: 5000000
- [INFO] FULL ALIGN PROPORTION: 0.7
- [INFO] FULL ALIGN REFERENCE PROPORTION: 0.1
- [INFO] PHASING PROPORTION: 0.7
- [INFO] MINIMUM MQ: 5
- [INFO] MINIMUM COVERAGE: 4
- [INFO] SNP AF THRESHOLD: 0.08
- [INFO] INDEL AF THRESHOLD: 0.15
- [INFO] BASE ERROR IN GVCF: 0.001
- [INFO] GQ BIN SIZE IN GVCF: 5
- [INFO] ENABLE FILEUP ONLY CALLING: False
- [INFO] ENABLE FAST MODE CALLING: False
- [INFO] ENABLE CALLING SNP CANDIDATES ONLY: False
- [INFO] ENABLE PRINTING REFERENCE CALLS: False
- [INFO] ENABLE OUTPUT GVCF: False
- [INFO] ENABLE HAPLOID PRECISE MODE: False
- [INFO] ENABLE HAPLOID SENSITIVE MODE: False
- [INFO] ENABLE INCLUDE ALL CTGS CALLING: False
- [INFO] ENABLE NO PHASING FOR FULL ALIGNMENT: False
- [INFO] ENABLE REMOVING INTERMEDIATE FILES: False
- [INFO] ENABLE LONGPHASE FOR INTERMEDIATE VCF PHASING: True
- [INFO] ENABLE PHASING FINAL VCF OUTPUT USING WHATSHAP: False
- [INFO] ENABLE PHASING FINAL VCF OUTPUT USING LONGPHASE: False
- [INFO] ENABLE HAPLOTAGGING FINAL BAM: False
- [INFO] ENABLE LONG INDEL CALLING: False
- [INFO] ENABLE C_IMPLEMENT: True
- + /opt/conda/envs/clairs/bin/scripts/clair3_c_impl.sh --bam_fn /mnt/beegfs02/scratch/t_steimle/data/wgs/CHAHA/norm/CHAHA_norm_hs1.bam --ref_fn /home/t_steimle/ref/hs1/chm13v2.0.fa --threads 40 --model_path /opt/conda/envs/clairs/bin/clairs_models/clair3_models/r1041_e82_400bps_sup_v420 --platform ont --output /mnt/beegfs02/scratch/t_steimle/data/wgs/CHAHA/diag/ClairS/part4/tmp/clair3_output/clair3_normal_output --bed_fn=EMPTY --vcf_fn=EMPTY --ctg_name=chr10 --sample_name=SAMPLE --chunk_num=0 --chunk_size=5000000 --samtools=samtools --python=python3 --pypy=pypy3 --parallel=parallel --whatshap=whatshap --qual=2 --var_pct_full=0.7 --ref_pct_full=0.1 --var_pct_phasing=0.7 --snp_min_af=0.08 --indel_min_af=0.15 --min_mq=5 --min_coverage=4 --min_contig_size=0 --pileup_only=False --gvcf=False --base_err=0.001 --gq_bin_size=5 --fast_mode=False --call_snp_only=False --print_ref_calls=False --haploid_precise=False --haploid_sensitive=False --include_all_ctgs=False --no_phasing_for_fa=False --pileup_model_prefix=pileup --fa_model_prefix=full_alignment --remove_intermediate_dir=False --enable_phasing=False --enable_long_indel=False --keep_iupac_bases=False --use_gpu=False --longphase_for_phasing=True --longphase=/opt/conda/envs/clairs/bin/longphase --use_whatshap_for_intermediate_phasing=True --use_longphase_for_intermediate_phasing=True --use_whatshap_for_final_output_phasing=False --use_longphase_for_final_output_phasing=False --use_whatshap_for_final_output_haplotagging=False
- [INFO] Check environment variables
- [INFO] Create folder /mnt/beegfs02/scratch/t_steimle/data/wgs/CHAHA/diag/ClairS/part4/tmp/clair3_output/clair3_normal_output/log
- [INFO] Create folder /mnt/beegfs02/scratch/t_steimle/data/wgs/CHAHA/diag/ClairS/part4/tmp/clair3_output/clair3_normal_output/tmp/pileup_output
- [INFO] Create folder /mnt/beegfs02/scratch/t_steimle/data/wgs/CHAHA/diag/ClairS/part4/tmp/clair3_output/clair3_normal_output/tmp/merge_output
- [INFO] Create folder /mnt/beegfs02/scratch/t_steimle/data/wgs/CHAHA/diag/ClairS/part4/tmp/clair3_output/clair3_normal_output/tmp/phase_output
- [INFO] Create folder /mnt/beegfs02/scratch/t_steimle/data/wgs/CHAHA/diag/ClairS/part4/tmp/clair3_output/clair3_normal_output/tmp/gvcf_tmp_output
- [INFO] Create folder /mnt/beegfs02/scratch/t_steimle/data/wgs/CHAHA/diag/ClairS/part4/tmp/clair3_output/clair3_normal_output/tmp/full_alignment_output
- [INFO] Create folder /mnt/beegfs02/scratch/t_steimle/data/wgs/CHAHA/diag/ClairS/part4/tmp/clair3_output/clair3_normal_output/tmp/phase_output/phase_vcf
- [INFO] Create folder /mnt/beegfs02/scratch/t_steimle/data/wgs/CHAHA/diag/ClairS/part4/tmp/clair3_output/clair3_normal_output/tmp/phase_output/phase_bam
- [INFO] Create folder /mnt/beegfs02/scratch/t_steimle/data/wgs/CHAHA/diag/ClairS/part4/tmp/clair3_output/clair3_normal_output/tmp/full_alignment_output/candidate_bed
- [INFO] Call variant in contigs: chr10
- [INFO] Chunk number for each contig: 27
- [INFO] 1/7 Call variants using pileup model
- slurmstepd: error: *** JOB 2556296 ON n05 CANCELLED AT 2025-12-08T16:49:54 ***
- slurmstepd: error: Detected 8 oom_kill events in StepId=2556296.batch. Some of the step tasks have been OOM Killed.
- ----------------------------------------------
- | SLURM EPILOG |
- ----------------------------------------------
- Job ID: 2556296
- Cluster: flamingo
- User/Group: t_steimle/gs_hpc_biopath
- Nodelist: n05
- Nodes: 1
- Cores per node: 40
- Job started at: 2025-12-08 16:47:24
- Job ended at: 2025-12-08 16:49:54
- Job Wall-clock time: 00:02:30
- CPU Utilized: 00:59:54
- CPU Efficiency: 59.90% of 01:40:00 core-walltime
- Memory Utilized: 39.58 GB
- Memory Efficiency: 98.96% of 40.00 GB
|