File Info

Filename
.command.log
Full Path
s3://natera-rnd-pltf-dev-nextflow-scratch-01/work/db/4e59882341fa7865f432f440dc1e3a/.command.log
Size
3.6 KB
Attempt
  Downloading: s3://natera-platform-sandbox/pipeline-resources/ngi-igenomes/igenomes/Homo_sapiens/GATK/GRCh38/Sequence/WholeGenomeFasta/Homo_sapiens_assembly38.fasta
  Downloading: s3://natera-platform-sandbox/pipeline-inputs/test_sarek/end_to_end_regression/01_preprocess_regression_fastqs/regression_preprocess__60492e9--20260310-034416/alignments/recalibrated/HCC1395_tumor/HCC1395_tumor.recalibrated.bam.bai
  Downloading: s3://natera-rnd-pltf-dev-nextflow-scratch-01/work/db/4e59882341fa7865f432f440dc1e3a/.command.sh
  Downloading: s3://natera-platform-sandbox/pipeline-inputs/test_sarek/end_to_end_regression/01_preprocess_regression_fastqs/regression_preprocess__60492e9--20260310-034416/alignments/recalibrated/HCC1395_tumor/HCC1395_tumor.recalibrated.bam
  Downloading: s3://natera-rnd-pltf-dev-nextflow-scratch-01/work/e0/fb262875a3360e87179755a222def4/HCC1395_tumor_vs_HCC1395_BL.tnseq.filtered.vcf.gz
  Downloading: s3://natera-platform-sandbox/pipeline-resources/ngi-igenomes/igenomes/Homo_sapiens/GATK/GRCh38/Sequence/WholeGenomeFasta/Homo_sapiens_assembly38.fasta.fai
  Downloading: s3://natera-rnd-pltf-dev-nextflow-scratch-01/work/e0/fb262875a3360e87179755a222def4/HCC1395_tumor_vs_HCC1395_BL.tnseq.filtered.vcf.gz.tbi
  Downloading: s3://natera-rnd-pltf-dev-nextflow-scratch-01/work/db/4e59882341fa7865f432f440dc1e3a/.command.run
==> STAGING COMPLETE (8 inputs)

This is WhatsHap 2.8 running under Python 3.9.23
[W::hts_idx_load3] The index file is older than the data file: HCC1395_tumor.recalibrated.bam.bai
[W::hts_idx_load3] The index file is older than the data file: HCC1395_tumor.recalibrated.bam.bai
[E::idx_find_and_load] Could not retrieve index file for 'input_filtered.vcf.gz'
[E::idx_find_and_load] Could not retrieve index file for 'input_filtered.vcf.gz'
[E::idx_find_and_load] Could not retrieve index file for 'input_filtered.vcf.gz'
Working on 2 samples from 2 families

# Working on contig chr21 in individual HCC1395_BL
Found 0 usable heterozygous variants (0 skipped due to missing genotypes)
WARNING: Sample 'HCC1395_BL' not found in any BAM/CRAM file.
Found 0 reads covering 0 variants
Kept 0 reads that cover at least two variants each
Selected 0 most phase-informative reads covering 0 variants
Best-case phasing would result in 0 non-singleton phased blocks (0 singletons). 
Phasing 1 sample by solving the MEC problem ...
No. of phased blocks: 0

# Working on contig chr21 in individual HCC1395_tumor
Found 204 usable heterozygous variants (0 skipped due to missing genotypes)
[W::hts_idx_load3] The index file is older than the data file: HCC1395_tumor.recalibrated.bam.bai
Number of supplementary alignments: 0
Number of non-singleton groups: 12264
Skipped 0 groups
Found 45727 reads covering 204 variants
Kept 12991 reads that cover at least two variants each
Selected 619 most phase-informative reads covering 121 variants
Best-case phasing would result in 40 non-singleton phased blocks (0 singletons). 
Phasing 1 sample by solving the MEC problem ...
Largest block contains 8 variants (6.6% of accessible variants) between position 44573789 and 44592359
WARNING: Ignoring existing phasing information found in input VCF (PS tag exists).

# Resource usage
Maximum memory usage: 0.161 GB
Time spent reading BAM/CRAM:                   40.2 s
Time spent parsing VCF:                         0.0 s
Time spent selecting reads:                     1.0 s
Time spent phasing:                             0.6 s
Time spent writing VCF:                         0.0 s
Time spent finding components:                  0.0 s
Time spent on rest:                             0.2 s
Total elapsed time:                            42.0 s