Back to Multiple platform build/check report for BioC 3.18:   simplified   long
ABCDEFGHIJKLMNO[P]QRSTUVWXYZ

This page was generated on 2024-03-27 11:37:58 -0400 (Wed, 27 Mar 2024).

HostnameOSArch (*)R versionInstalled pkgs
nebbiolo2Linux (Ubuntu 22.04.3 LTS)x86_644.3.3 (2024-02-29) -- "Angel Food Cake" 4667
palomino4Windows Server 2022 Datacenterx644.3.3 (2024-02-29 ucrt) -- "Angel Food Cake" 4403
merida1macOS 12.7.1 Montereyx86_644.3.3 (2024-02-29) -- "Angel Food Cake" 4426
Click on any hostname to see more info about the system (e.g. compilers)      (*) as reported by 'uname -p', except on Windows and Mac OS X

Package 1613/2266HostnameOS / ArchINSTALLBUILDCHECKBUILD BIN
PureCN 2.8.1  (landing page)
Markus Riester
Snapshot Date: 2024-03-25 14:05:07 -0400 (Mon, 25 Mar 2024)
git_url: https://git.bioconductor.org/packages/PureCN
git_branch: RELEASE_3_18
git_last_commit: 43f6a49
git_last_commit_date: 2023-11-20 16:23:17 -0400 (Mon, 20 Nov 2023)
nebbiolo2Linux (Ubuntu 22.04.3 LTS) / x86_64  OK    OK    OK  UNNEEDED, same version is already published
palomino4Windows Server 2022 Datacenter / x64  OK    OK    OK    OK  UNNEEDED, same version is already published
merida1macOS 12.7.1 Monterey / x86_64  OK    OK    OK    OK  UNNEEDED, same version is already published
kjohnson1macOS 13.6.1 Ventura / arm64see weekly results here

CHECK results for PureCN on merida1


To the developers/maintainers of the PureCN package:
- Allow up to 24 hours (and sometimes 48 hours) for your latest push to git@git.bioconductor.org:packages/PureCN.git to reflect on this report. See Troubleshooting Build Report for more information.
- Use the following Renviron settings to reproduce errors and warnings.
- If 'R CMD check' started to fail recently on the Linux builder(s) over a missing dependency, add the missing dependency to 'Suggests:' in your DESCRIPTION file. See Renviron.bioc for more information.

raw results


Summary

Package: PureCN
Version: 2.8.1
Command: /Library/Frameworks/R.framework/Resources/bin/R CMD check --install=check:PureCN.install-out.txt --library=/Library/Frameworks/R.framework/Resources/library --no-vignettes --timings PureCN_2.8.1.tar.gz
StartedAt: 2024-03-26 07:01:50 -0400 (Tue, 26 Mar 2024)
EndedAt: 2024-03-26 07:35:54 -0400 (Tue, 26 Mar 2024)
EllapsedTime: 2043.2 seconds
RetCode: 0
Status:   OK  
CheckDir: PureCN.Rcheck
Warnings: 0

Command output

##############################################################################
##############################################################################
###
### Running command:
###
###   /Library/Frameworks/R.framework/Resources/bin/R CMD check --install=check:PureCN.install-out.txt --library=/Library/Frameworks/R.framework/Resources/library --no-vignettes --timings PureCN_2.8.1.tar.gz
###
##############################################################################
##############################################################################


* using log directory ‘/Users/biocbuild/bbs-3.18-bioc/meat/PureCN.Rcheck’
* using R version 4.3.3 (2024-02-29)
* using platform: x86_64-apple-darwin20 (64-bit)
* R was compiled by
    Apple clang version 14.0.0 (clang-1400.0.29.202)
    GNU Fortran (GCC) 12.2.0
* running under: macOS Monterey 12.7.1
* using session charset: UTF-8
* using option ‘--no-vignettes’
* checking for file ‘PureCN/DESCRIPTION’ ... OK
* checking extension type ... Package
* this is package ‘PureCN’ version ‘2.8.1’
* checking package namespace information ... OK
* checking package dependencies ... NOTE
Package which this enhances but not available for checking: ‘genomicsdb’
* checking if this is a source package ... OK
* checking if there is a namespace ... OK
* checking for hidden files and directories ... OK
* checking for portable file names ... OK
* checking for sufficient/correct file permissions ... OK
* checking whether package ‘PureCN’ can be installed ... OK
* checking installed package size ... OK
* checking package directory ... OK
* checking ‘build’ directory ... OK
* checking DESCRIPTION meta-information ... OK
* checking top-level files ... OK
* checking for left-over files ... OK
* checking index information ... OK
* checking package subdirectories ... OK
* checking R files for non-ASCII characters ... OK
* checking R files for syntax errors ... OK
* checking whether the package can be loaded ... OK
* checking whether the package can be loaded with stated dependencies ... OK
* checking whether the package can be unloaded cleanly ... OK
* checking whether the namespace can be loaded with stated dependencies ... OK
* checking whether the namespace can be unloaded cleanly ... OK
* checking dependencies in R code ... OK
* checking S3 generic/method consistency ... OK
* checking replacement functions ... OK
* checking foreign function calls ... OK
* checking R code for possible problems ... OK
* checking Rd files ... OK
* checking Rd metadata ... OK
* checking Rd cross-references ... OK
* checking for missing documentation entries ... OK
* checking for code/documentation mismatches ... OK
* checking Rd \usage sections ... OK
* checking Rd contents ... OK
* checking for unstated dependencies in examples ... OK
* checking contents of ‘data’ directory ... OK
* checking data for non-ASCII characters ... OK
* checking data for ASCII and uncompressed saves ... OK
* checking sizes of PDF files under ‘inst/doc’ ... OK
* checking files in ‘vignettes’ ... OK
* checking examples ... OK
Examples with CPU (user + system) or elapsed time > 5s
                                 user system elapsed
callAmplificationsInLowPurity 100.793  0.390 101.712
segmentationPSCBS              65.927  0.947  74.123
filterIntervals                52.531  0.606  55.007
runAbsoluteCN                  48.072  0.776  52.927
segmentationHclust             45.153  0.664  50.356
findFocal                      30.141  0.180  31.394
segmentationCBS                23.282  0.185  25.931
annotateTargets                14.673  0.635  15.362
calculateTangentNormal          7.829  0.084   8.016
processMultipleSamples          6.962  0.077   7.415
createNormalDatabase            6.830  0.068   7.264
* checking for unstated dependencies in ‘tests’ ... OK
* checking tests ...
  Running ‘testthat.R’
 OK
* checking for unstated dependencies in vignettes ... OK
* checking package vignettes in ‘inst/doc’ ... OK
* checking running R code from vignettes ... SKIPPED
* checking re-building of vignette outputs ... SKIPPED
* checking PDF version of manual ... OK
* DONE

Status: 1 NOTE
See
  ‘/Users/biocbuild/bbs-3.18-bioc/meat/PureCN.Rcheck/00check.log’
for details.



Installation output

PureCN.Rcheck/00install.out

##############################################################################
##############################################################################
###
### Running command:
###
###   /Library/Frameworks/R.framework/Resources/bin/R CMD INSTALL PureCN
###
##############################################################################
##############################################################################


* installing to library ‘/Library/Frameworks/R.framework/Versions/4.3-x86_64/Resources/library’
* installing *source* package ‘PureCN’ ...
** using staged installation
** R
** data
** inst
** byte-compile and prepare package for lazy loading
** help
*** installing help indices
** building package indices
** installing vignettes
** testing if installed package can be loaded from temporary location
** testing if installed package can be loaded from final location
** testing if installed package keeps a record of temporary installation path
* DONE (PureCN)

Tests output

PureCN.Rcheck/tests/testthat.Rout


R version 4.3.3 (2024-02-29) -- "Angel Food Cake"
Copyright (C) 2024 The R Foundation for Statistical Computing
Platform: x86_64-apple-darwin20 (64-bit)

R is free software and comes with ABSOLUTELY NO WARRANTY.
You are welcome to redistribute it under certain conditions.
Type 'license()' or 'licence()' for distribution details.

R is a collaborative project with many contributors.
Type 'contributors()' for more information and
'citation()' on how to cite R or R packages in publications.

Type 'demo()' for some demos, 'help()' for on-line help, or
'help.start()' for an HTML browser interface to help.
Type 'q()' to quit R.

> library(testthat)
> library(PureCN)
Loading required package: DNAcopy
Loading required package: VariantAnnotation
Loading required package: BiocGenerics

Attaching package: 'BiocGenerics'

The following objects are masked from 'package:stats':

    IQR, mad, sd, var, xtabs

The following objects are masked from 'package:base':

    Filter, Find, Map, Position, Reduce, anyDuplicated, aperm, append,
    as.data.frame, basename, cbind, colnames, dirname, do.call,
    duplicated, eval, evalq, get, grep, grepl, intersect, is.unsorted,
    lapply, mapply, match, mget, order, paste, pmax, pmax.int, pmin,
    pmin.int, rank, rbind, rownames, sapply, setdiff, sort, table,
    tapply, union, unique, unsplit, which.max, which.min

Loading required package: MatrixGenerics
Loading required package: matrixStats

Attaching package: 'MatrixGenerics'

The following objects are masked from 'package:matrixStats':

    colAlls, colAnyNAs, colAnys, colAvgsPerRowSet, colCollapse,
    colCounts, colCummaxs, colCummins, colCumprods, colCumsums,
    colDiffs, colIQRDiffs, colIQRs, colLogSumExps, colMadDiffs,
    colMads, colMaxs, colMeans2, colMedians, colMins, colOrderStats,
    colProds, colQuantiles, colRanges, colRanks, colSdDiffs, colSds,
    colSums2, colTabulates, colVarDiffs, colVars, colWeightedMads,
    colWeightedMeans, colWeightedMedians, colWeightedSds,
    colWeightedVars, rowAlls, rowAnyNAs, rowAnys, rowAvgsPerColSet,
    rowCollapse, rowCounts, rowCummaxs, rowCummins, rowCumprods,
    rowCumsums, rowDiffs, rowIQRDiffs, rowIQRs, rowLogSumExps,
    rowMadDiffs, rowMads, rowMaxs, rowMeans2, rowMedians, rowMins,
    rowOrderStats, rowProds, rowQuantiles, rowRanges, rowRanks,
    rowSdDiffs, rowSds, rowSums2, rowTabulates, rowVarDiffs, rowVars,
    rowWeightedMads, rowWeightedMeans, rowWeightedMedians,
    rowWeightedSds, rowWeightedVars

Loading required package: GenomeInfoDb
Loading required package: S4Vectors
Loading required package: stats4

Attaching package: 'S4Vectors'

The following object is masked from 'package:utils':

    findMatches

The following objects are masked from 'package:base':

    I, expand.grid, unname

Loading required package: IRanges
Loading required package: GenomicRanges
Loading required package: SummarizedExperiment
Loading required package: Biobase
Welcome to Bioconductor

    Vignettes contain introductory material; view with
    'browseVignettes()'. To cite Bioconductor, see
    'citation("Biobase")', and for packages 'citation("pkgname")'.


Attaching package: 'Biobase'

The following object is masked from 'package:MatrixGenerics':

    rowMedians

The following objects are masked from 'package:matrixStats':

    anyMissing, rowMedians

Loading required package: Rsamtools
Loading required package: Biostrings
Loading required package: XVector

Attaching package: 'Biostrings'

The following object is masked from 'package:base':

    strsplit


Attaching package: 'VariantAnnotation'

The following object is masked from 'package:base':

    tabulate

> 
> test_check("PureCN")
WARN [2024-03-26 07:20:27] Attempted adding gene symbols to intervals. Heuristics have been used to pick symbols for overlapping genes.
WARN [2024-03-26 07:20:27] Chromosome naming style of txdb file (UCSC) was different from interval file (Ensembl).
WARN [2024-03-26 07:20:37] Attempted adding gene symbols to intervals. Heuristics have been used to pick symbols for overlapping genes.
FATAL [2024-03-26 07:20:43] tumor.coverage.file and interval.file do not align. 

FATAL [2024-03-26 07:20:43]  

FATAL [2024-03-26 07:20:43] This is most likely a user error due to invalid input data or 

FATAL [2024-03-26 07:20:43] parameters (PureCN 2.8.1). 

WARN [2024-03-26 07:20:44] Cannot find all contig lengths while exporting interval file.
INFO [2024-03-26 07:20:47] Processing seq1:1-21 (1/3)...
INFO [2024-03-26 07:20:48] Processing seq1:1227-1247 (2/3)...
INFO [2024-03-26 07:20:49] Processing seq2:594-614 (3/3)...
WARN [2024-03-26 07:20:51] Large potential mis-calibration of on- and off-target log2 ratios: 0.26
FATAL [2024-03-26 07:20:51] Need either f or purity and ploidy. 

FATAL [2024-03-26 07:20:51]  

FATAL [2024-03-26 07:20:51] This is most likely a user error due to invalid input data or 

FATAL [2024-03-26 07:20:51] parameters (PureCN 2.8.1). 

FATAL [2024-03-26 07:20:51] f not in expected range. 

FATAL [2024-03-26 07:20:51]  

FATAL [2024-03-26 07:20:51] This is most likely a user error due to invalid input data or 

FATAL [2024-03-26 07:20:51] parameters (PureCN 2.8.1). 

FATAL [2024-03-26 07:20:51] coverage not in expected range (>=2) 

FATAL [2024-03-26 07:20:51]  

FATAL [2024-03-26 07:20:51] This is most likely a user error due to invalid input data or 

FATAL [2024-03-26 07:20:51] parameters (PureCN 2.8.1). 

FATAL [2024-03-26 07:20:51] purity not in expected range. 

FATAL [2024-03-26 07:20:51]  

FATAL [2024-03-26 07:20:51] This is most likely a user error due to invalid input data or 

FATAL [2024-03-26 07:20:51] parameters (PureCN 2.8.1). 

FATAL [2024-03-26 07:20:51] ploidy not in expected range. 

FATAL [2024-03-26 07:20:51]  

FATAL [2024-03-26 07:20:51] This is most likely a user error due to invalid input data or 

FATAL [2024-03-26 07:20:51] parameters (PureCN 2.8.1). 

FATAL [2024-03-26 07:20:51] cell.fraction not in expected range. 

FATAL [2024-03-26 07:20:51]  

FATAL [2024-03-26 07:20:51] This is most likely a user error due to invalid input data or 

FATAL [2024-03-26 07:20:51] parameters (PureCN 2.8.1). 

INFO [2024-03-26 07:20:57] 576 on-target bins with low coverage in all samples.
WARN [2024-03-26 07:20:57] You are likely not using the correct baits file!
WARN [2024-03-26 07:20:57] Allosome coverage missing, cannot determine sex.
WARN [2024-03-26 07:20:57] Allosome coverage missing, cannot determine sex.
INFO [2024-03-26 07:20:57] Processing on-target regions...
INFO [2024-03-26 07:20:59] Removing 930 intervals with low coverage in normalDB.
INFO [2024-03-26 07:20:59] Removing 1 intervals with zero coverage in more than 3% of normalDB.
INFO [2024-03-26 07:21:03] Tumor/normal noise ratio: 19.041
WARN [2024-03-26 07:21:03] Extensive noise in tumor compared to normals.
INFO [2024-03-26 07:23:06] Tumor/normal noise ratio: 19.041
WARN [2024-03-26 07:23:06] Extensive noise in tumor compared to normals.
INFO [2024-03-26 07:23:11] Using BiocParallel for parallel optimization.
FATAL [2024-03-26 07:23:47] pvalue.cutoff not within expected range or format. 

FATAL [2024-03-26 07:23:47]  

FATAL [2024-03-26 07:23:47] This is most likely a user error due to invalid input data or 

FATAL [2024-03-26 07:23:47] parameters (PureCN 2.8.1). 

FATAL [2024-03-26 07:23:47] pvalue.cutoff not within expected range or format. 

FATAL [2024-03-26 07:23:47]  

FATAL [2024-03-26 07:23:47] This is most likely a user error due to invalid input data or 

FATAL [2024-03-26 07:23:47] parameters (PureCN 2.8.1). 

FATAL [2024-03-26 07:23:47] percentile.cutoff not in expected range (0 to 100). 

FATAL [2024-03-26 07:23:47]  

FATAL [2024-03-26 07:23:47] This is most likely a user error due to invalid input data or 

FATAL [2024-03-26 07:23:47] parameters (PureCN 2.8.1). 

FATAL [2024-03-26 07:23:47] percentile.cutoff not in expected range (0 to 100). 

FATAL [2024-03-26 07:23:47]  

FATAL [2024-03-26 07:23:47] This is most likely a user error due to invalid input data or 

FATAL [2024-03-26 07:23:47] parameters (PureCN 2.8.1). 

FATAL [2024-03-26 07:23:47] purity not within expected range or format. 

FATAL [2024-03-26 07:23:47]  

FATAL [2024-03-26 07:23:47] This is most likely a user error due to invalid input data or 

FATAL [2024-03-26 07:23:47] parameters (PureCN 2.8.1). 

FATAL [2024-03-26 07:23:47] purity not within expected range or format. 

FATAL [2024-03-26 07:23:47]  

FATAL [2024-03-26 07:23:47] This is most likely a user error due to invalid input data or 

FATAL [2024-03-26 07:23:47] parameters (PureCN 2.8.1). 

INFO [2024-03-26 07:23:52] ------------------------------------------------------------
INFO [2024-03-26 07:23:52] PureCN 2.8.1
INFO [2024-03-26 07:23:52] ------------------------------------------------------------
INFO [2024-03-26 07:23:52] Loading coverage files...
INFO [2024-03-26 07:23:52] Mean target coverages: 112X (tumor) 99X (normal).
WARN [2024-03-26 07:23:53] Allosome coverage missing, cannot determine sex.
WARN [2024-03-26 07:23:53] Allosome coverage missing, cannot determine sex.
INFO [2024-03-26 07:23:53] Removing 228 intervals with missing log.ratio.
INFO [2024-03-26 07:23:53] Removing 22 small (< 5bp) intervals.
INFO [2024-03-26 07:23:53] Removing 19 intervals with low total coverage in normal (< 150.00 reads).
WARN [2024-03-26 07:23:53] No normalDB provided. Provide one for better results.
INFO [2024-03-26 07:23:53] Removing 1724 low count (< 100 total reads) intervals.
INFO [2024-03-26 07:23:53] Removing 6 low coverage (< 15.0000X) intervals.
INFO [2024-03-26 07:23:53] Using 8050 intervals (8050 on-target, 0 off-target).
INFO [2024-03-26 07:23:53] No off-target intervals. If this is hybrid-capture data, consider adding them.
INFO [2024-03-26 07:23:53] No interval.file provided. Cannot check for any GC-biases.
INFO [2024-03-26 07:23:53] Loading VCF...
INFO [2024-03-26 07:23:53] Found 127 variants in VCF file.
INFO [2024-03-26 07:23:53] 127 (100.0%) variants annotated as likely germline (DB INFO flag).
INFO [2024-03-26 07:23:54] LIB-02240e4 is tumor in VCF file.
INFO [2024-03-26 07:23:54] No homozygous variants in VCF, provide unfiltered VCF.
INFO [2024-03-26 07:23:54] Removing 1 non heterozygous (in matched normal) germline SNPs.
INFO [2024-03-26 07:23:54] Removing 0 low quality variants with non-offset BQ < 25.
INFO [2024-03-26 07:23:54] Base quality scores range from 29 to 35 (offset by 1)
INFO [2024-03-26 07:23:54] Minimum number of supporting reads ranges from 2 to 5, depending on coverage and BQS.
INFO [2024-03-26 07:23:55] Removing 1 variants with AF < 0.030 or AF >= 1.000 or insufficient supporting reads or depth < 15.
INFO [2024-03-26 07:23:55] Total size of targeted genomic region: 1.53Mb (2.33Mb with 50bp padding).
INFO [2024-03-26 07:23:56] 1.2% of targets contain variants.
INFO [2024-03-26 07:23:56] Removing 4 variants outside intervals.
INFO [2024-03-26 07:23:56] Found SOMATIC annotation in VCF.
INFO [2024-03-26 07:23:56] Setting somatic prior probabilities for somatic variants to 0.999000 or to 0.000100 otherwise.
INFO [2024-03-26 07:23:57] Found SOMATIC annotation in VCF. Setting mapping bias to 0.972.
INFO [2024-03-26 07:23:57] Excluding 0 novel or poor quality variants from segmentation.
INFO [2024-03-26 07:23:57] Sample sex: ?
INFO [2024-03-26 07:23:57] Segmenting data...
INFO [2024-03-26 07:23:57] Loading pre-computed boundaries for DNAcopy...
INFO [2024-03-26 07:23:57] Setting undo.SD parameter to 1.000000.
INFO [2024-03-26 07:23:59] Setting prune.hclust.h parameter to 0.200000.
INFO [2024-03-26 07:23:59] Found 59 segments with median size of 17.67Mb.
INFO [2024-03-26 07:23:59] Using 121 variants.
INFO [2024-03-26 07:23:59] Mean standard deviation of log-ratios: 0.34 (MAPD: 0.30)
INFO [2024-03-26 07:23:59] 2D-grid search of purity and ploidy...
INFO [2024-03-26 07:24:03] Local optima: 0.63/1.9, 0.5/2
INFO [2024-03-26 07:24:03] Testing local optimum 1/2 at purity 0.63 and total ploidy 1.90...
INFO [2024-03-26 07:24:05] Testing local optimum 2/2 at purity 0.50 and total ploidy 2.00...
INFO [2024-03-26 07:24:08] Skipping 1 solutions that converged to the same optima.
INFO [2024-03-26 07:24:08] Fitting variants with beta model for local optimum 1/2...
INFO [2024-03-26 07:24:08] Fitting variants for purity 0.65, tumor ploidy 1.75 and contamination 0.01.
INFO [2024-03-26 07:24:17] Optimized purity: 0.65
INFO [2024-03-26 07:24:17] Done.
INFO [2024-03-26 07:24:17] ------------------------------------------------------------
INFO [2024-03-26 07:24:18] Estimating callable regions.
FATAL [2024-03-26 07:24:22] exclude not a GRanges object. 

FATAL [2024-03-26 07:24:22]  

FATAL [2024-03-26 07:24:22] This is most likely a user error due to invalid input data or 

FATAL [2024-03-26 07:24:22] parameters (PureCN 2.8.1). 

FATAL [2024-03-26 07:24:22] callable not a GRanges object. 

FATAL [2024-03-26 07:24:22]  

FATAL [2024-03-26 07:24:22] This is most likely a user error due to invalid input data or 

FATAL [2024-03-26 07:24:22] parameters (PureCN 2.8.1). 

WARN [2024-03-26 07:24:42] tumor.coverage.file and interval.file do not align.
INFO [2024-03-26 07:24:44] No Gene column in interval.file. You won't get gene-level calls.
FATAL [2024-03-26 07:24:47] No gc_bias column in interval.file. 

FATAL [2024-03-26 07:24:47]  

FATAL [2024-03-26 07:24:47] This is most likely a user error due to invalid input data or 

FATAL [2024-03-26 07:24:47] parameters (PureCN 2.8.1). 

FATAL [2024-03-26 07:24:48] Provided coverage is zero, most likely due to a corrupt BAM file. 

FATAL [2024-03-26 07:24:48]  

FATAL [2024-03-26 07:24:48] This is most likely a user error due to invalid input data or 

FATAL [2024-03-26 07:24:48] parameters (PureCN 2.8.1). 

INFO [2024-03-26 07:24:54] No reptiming column in interval.file.
INFO [2024-03-26 07:25:02] Reading /tmp/RtmpP77MMO/file91854563f224.rds...
INFO [2024-03-26 07:25:02] Reading /tmp/RtmpP77MMO/file91854563f224.rds...
INFO [2024-03-26 07:25:02] Reading /tmp/RtmpP77MMO/file91854563f224.rds...
INFO [2024-03-26 07:25:02] Reading /tmp/RtmpP77MMO/file91854563f224.rds...
INFO [2024-03-26 07:25:03] Reading /tmp/RtmpP77MMO/file91854563f224.rds...
INFO [2024-03-26 07:25:03] Reading /tmp/RtmpP77MMO/file91854563f224.rds...
FATAL [2024-03-26 07:25:03] Purity or Ploidy not numeric or in expected range. 

FATAL [2024-03-26 07:25:03]  

FATAL [2024-03-26 07:25:03] This is most likely a user error due to invalid input data or 

FATAL [2024-03-26 07:25:03] parameters (PureCN 2.8.1). 

INFO [2024-03-26 07:25:03] Reading /tmp/RtmpP77MMO/file91854563f224.rds...
INFO [2024-03-26 07:25:03] Reading /tmp/RtmpP77MMO/file91854563f224.rds...
FATAL [2024-03-26 07:25:03] 'Failed' column in /tmp/RtmpP77MMO/file91854563f224.csv not logical(1). 

FATAL [2024-03-26 07:25:03]  

FATAL [2024-03-26 07:25:03] This is most likely a user error due to invalid input data or 

FATAL [2024-03-26 07:25:03] parameters (PureCN 2.8.1). 

INFO [2024-03-26 07:25:05] 576 on-target bins with low coverage in all samples.
WARN [2024-03-26 07:25:05] You are likely not using the correct baits file!
WARN [2024-03-26 07:25:05] Allosome coverage missing, cannot determine sex.
WARN [2024-03-26 07:25:05] Allosome coverage missing, cannot determine sex.
INFO [2024-03-26 07:25:05] Processing on-target regions...
INFO [2024-03-26 07:25:07] Removing 930 intervals with low coverage in normalDB.
INFO [2024-03-26 07:25:07] Removing 1 intervals with zero coverage in more than 3% of normalDB.
INFO [2024-03-26 07:25:14] 576 on-target bins with low coverage in all samples.
WARN [2024-03-26 07:25:14] You are likely not using the correct baits file!
WARN [2024-03-26 07:25:14] Allosome coverage missing, cannot determine sex.
WARN [2024-03-26 07:25:14] Allosome coverage missing, cannot determine sex.
INFO [2024-03-26 07:25:14] Processing on-target regions...
INFO [2024-03-26 07:25:16] Removing 930 intervals with low coverage in normalDB.
INFO [2024-03-26 07:25:16] Removing 1 intervals with zero coverage in more than 3% of normalDB.
INFO [2024-03-26 07:25:21] 576 on-target bins with low coverage in all samples.
WARN [2024-03-26 07:25:21] You are likely not using the correct baits file!
WARN [2024-03-26 07:25:21] Allosome coverage missing, cannot determine sex.
WARN [2024-03-26 07:25:21] Allosome coverage missing, cannot determine sex.
INFO [2024-03-26 07:25:21] Processing on-target regions...
INFO [2024-03-26 07:25:23] Removing 930 intervals with low coverage in normalDB.
INFO [2024-03-26 07:25:23] Removing 1 intervals with zero coverage in more than 3% of normalDB.
WARN [2024-03-26 07:25:25] Allosome coverage missing, cannot determine sex.
INFO [2024-03-26 07:25:25] Sample sex: NA
WARN [2024-03-26 07:25:26] Allosome coverage missing, cannot determine sex.
INFO [2024-03-26 07:25:26] Sample sex: NA
INFO [2024-03-26 07:25:28] 576 on-target bins with low coverage in all samples.
WARN [2024-03-26 07:25:28] You are likely not using the correct baits file!
WARN [2024-03-26 07:25:28] Allosome coverage missing, cannot determine sex.
WARN [2024-03-26 07:25:28] Allosome coverage missing, cannot determine sex.
FATAL [2024-03-26 07:25:28] Length of normal.coverage.files and sex different 

FATAL [2024-03-26 07:25:28]  

FATAL [2024-03-26 07:25:28] This is most likely a user error due to invalid input data or 

FATAL [2024-03-26 07:25:28] parameters (PureCN 2.8.1). 

WARN [2024-03-26 07:25:38] Target intervals were not sorted.
INFO [2024-03-26 07:25:38] 560 on-target bins with low coverage in all samples.
WARN [2024-03-26 07:25:38] You are likely not using the correct baits file!
WARN [2024-03-26 07:25:38] Allosome coverage missing, cannot determine sex.
WARN [2024-03-26 07:25:38] Allosome coverage missing, cannot determine sex.
WARN [2024-03-26 07:25:38] Allosome coverage missing, cannot determine sex.
INFO [2024-03-26 07:25:38] Processing on-target regions...
INFO [2024-03-26 07:25:40] Removing 978 intervals with low coverage in normalDB.
INFO [2024-03-26 07:25:40] Removing 11 intervals with zero coverage in more than 3% of normalDB.
FATAL [2024-03-26 07:25:45] tumor.coverage.file and normalDB do not align. 

FATAL [2024-03-26 07:25:45]  

FATAL [2024-03-26 07:25:45] This is most likely a user error due to invalid input data or 

FATAL [2024-03-26 07:25:45] parameters (PureCN 2.8.1). 

FATAL [2024-03-26 07:25:45] At least 2 normal.coverage.files required. 

FATAL [2024-03-26 07:25:45]  

FATAL [2024-03-26 07:25:45] This is most likely a user error due to invalid input data or 

FATAL [2024-03-26 07:25:45] parameters (PureCN 2.8.1). 

INFO [2024-03-26 07:25:53] ------------------------------------------------------------
INFO [2024-03-26 07:25:53] PureCN 2.8.1
INFO [2024-03-26 07:25:53] ------------------------------------------------------------
INFO [2024-03-26 07:25:53] Loading coverage files...
INFO [2024-03-26 07:25:54] Mean target coverages: 112X (tumor) 99X (normal).
WARN [2024-03-26 07:25:54] Allosome coverage missing, cannot determine sex.
WARN [2024-03-26 07:25:54] Allosome coverage missing, cannot determine sex.
INFO [2024-03-26 07:25:54] Removing 228 intervals with missing log.ratio.
FATAL [2024-03-26 07:25:54] normalDB incompatible with this PureCN version. Please re-run 

FATAL [2024-03-26 07:25:54] NormalDB.R. 

FATAL [2024-03-26 07:25:54]  

FATAL [2024-03-26 07:25:54] This is most likely a user error due to invalid input data or 

FATAL [2024-03-26 07:25:54] parameters (PureCN 2.8.1). 

FATAL [2024-03-26 07:25:55] normalDB incompatible with this PureCN version. Please re-run 

FATAL [2024-03-26 07:25:55] NormalDB.R. 

FATAL [2024-03-26 07:25:55]  

FATAL [2024-03-26 07:25:55] This is most likely a user error due to invalid input data or 

FATAL [2024-03-26 07:25:55] parameters (PureCN 2.8.1). 

FATAL [2024-03-26 07:25:55] normal.coverage.files with _coverage.txt and _loess.txt suffix 

FATAL [2024-03-26 07:25:55] provided. Provide either only GC-normalized or raw coverage files! 

FATAL [2024-03-26 07:25:55]  

FATAL [2024-03-26 07:25:55] This is most likely a user error due to invalid input data or 

FATAL [2024-03-26 07:25:55] parameters (PureCN 2.8.1). 

INFO [2024-03-26 07:25:56] Removing 22 non heterozygous (in matched normal) germline SNPs.
INFO [2024-03-26 07:25:56] Removing 19 low quality variants with non-offset BQ < 25.
INFO [2024-03-26 07:25:56] Base quality scores range from 24 to 35 (offset by 1)
INFO [2024-03-26 07:25:56] Minimum number of supporting reads ranges from 2 to 9, depending on coverage and BQS.
INFO [2024-03-26 07:25:57] Initial testing for significant sample cross-contamination: unlikely
INFO [2024-03-26 07:25:57] Removing 45 variants with AF < 0.030 or AF >= 1.000 or insufficient supporting reads or depth < 15.
INFO [2024-03-26 07:25:58] Removing 22 non heterozygous (in matched normal) germline SNPs.
INFO [2024-03-26 07:25:58] Removing 19 low quality variants with non-offset BQ < 25.
INFO [2024-03-26 07:25:58] Base quality scores range from 24 to 35 (offset by 1)
INFO [2024-03-26 07:25:58] Minimum number of supporting reads ranges from 2 to 9, depending on coverage and BQS.
INFO [2024-03-26 07:25:59] Initial testing for significant sample cross-contamination: unlikely
INFO [2024-03-26 07:25:59] Removing 45 variants with AF < 0.030 or AF >= 1.000 or insufficient supporting reads or depth < 15.
INFO [2024-03-26 07:26:00] Removing 6 blacklisted variants.
INFO [2024-03-26 07:26:00] Removing 22 non heterozygous (in matched normal) germline SNPs.
INFO [2024-03-26 07:26:00] Removing 19 low quality variants with non-offset BQ < 25.
INFO [2024-03-26 07:26:00] Base quality scores range from 24 to 35 (offset by 1)
INFO [2024-03-26 07:26:00] Minimum number of supporting reads ranges from 2 to 9, depending on coverage and BQS.
INFO [2024-03-26 07:26:02] Initial testing for significant sample cross-contamination: unlikely
INFO [2024-03-26 07:26:02] Removing 45 variants with AF < 0.030 or AF >= 1.000 or insufficient supporting reads or depth < 15.
WARN [2024-03-26 07:26:02] MuTect stats file lacks contig and position columns.
INFO [2024-03-26 07:26:02] Removing 22 non heterozygous (in matched normal) germline SNPs.
INFO [2024-03-26 07:26:02] Removing 19 low quality variants with non-offset BQ < 25.
INFO [2024-03-26 07:26:02] Base quality scores range from 24 to 35 (offset by 1)
INFO [2024-03-26 07:26:02] Minimum number of supporting reads ranges from 2 to 9, depending on coverage and BQS.
INFO [2024-03-26 07:26:03] Initial testing for significant sample cross-contamination: unlikely
INFO [2024-03-26 07:26:04] Removing 45 variants with AF < 0.030 or AF >= 1.000 or insufficient supporting reads or depth < 15.
WARN [2024-03-26 07:26:04] MuTect stats file and VCF file do not align perfectly. Will remove 2325 unmatched variants.
WARN [2024-03-26 07:26:04] MuTect stats file lacks failure_reasons column. Keeping all variants listed in stats file.
INFO [2024-03-26 07:26:04] Removing 0 non heterozygous (in matched normal) germline SNPs.
INFO [2024-03-26 07:26:04] Removing 0 low quality variants with non-offset BQ < 25.
INFO [2024-03-26 07:26:04] Base quality scores range from 31 to 33 (offset by 1)
INFO [2024-03-26 07:26:05] Minimum number of supporting reads ranges from 3 to 4, depending on coverage and BQS.
INFO [2024-03-26 07:26:05] Removing 0 variants with AF < 0.030 or AF >= 1.000 or insufficient supporting reads or depth < 15.
INFO [2024-03-26 07:26:06] Removing 22 non heterozygous (in matched normal) germline SNPs.
INFO [2024-03-26 07:26:06] Base quality scores range from 7 to 35 (offset by 1)
INFO [2024-03-26 07:26:06] Minimum number of supporting reads ranges from 1 to 41, depending on coverage and BQS.
INFO [2024-03-26 07:26:07] Initial testing for significant sample cross-contamination: unlikely
INFO [2024-03-26 07:26:07] Removing 49 variants with AF < 0.030 or AF >= 1.000 or insufficient supporting reads or depth < 15.
INFO [2024-03-26 07:26:07] Removing 22 non heterozygous (in matched normal) germline SNPs.
INFO [2024-03-26 07:26:07] Base quality scores range from 7 to 35 (offset by 1)
INFO [2024-03-26 07:26:08] Minimum number of supporting reads ranges from 1 to 41, depending on coverage and BQS.
INFO [2024-03-26 07:26:09] Initial testing for significant sample cross-contamination: unlikely
INFO [2024-03-26 07:26:09] Removing 49 variants with AF < 0.030 or AF >= 1.000 or insufficient supporting reads or depth < 15.
INFO [2024-03-26 07:26:09] Removing 22 non heterozygous (in matched normal) germline SNPs.
WARN [2024-03-26 07:26:09] Many variants removed by min.base.quality (50) filter. You might want to lower the cutoff.
FATAL [2024-03-26 07:26:09] No variants passed filter BQ. 

FATAL [2024-03-26 07:26:09]  

FATAL [2024-03-26 07:26:09] This is most likely a user error due to invalid input data or 

FATAL [2024-03-26 07:26:09] parameters (PureCN 2.8.1). 

INFO [2024-03-26 07:26:10] Found 11 variants in VCF file.
WARN [2024-03-26 07:26:11] vcf.file has no DB info field for membership in germline databases. Found and used valid population allele frequency > 0.001000 instead.
INFO [2024-03-26 07:26:11] 1 (9.1%) variants annotated as likely germline (DB INFO flag).
INFO [2024-03-26 07:26:16] Found 11 variants in VCF file.
WARN [2024-03-26 07:26:16] vcf.file has no DB info field for membership in germline databases. Found and used valid population allele frequency > 0.001000 instead.
INFO [2024-03-26 07:26:17] 1 (9.1%) variants annotated as likely germline (DB INFO flag).
INFO [2024-03-26 07:26:19] Found 11 variants in VCF file.
WARN [2024-03-26 07:26:19] vcf.file has no DB or POP_AF info field for membership in germline databases. Guessing it based on available dbSNP ID.
INFO [2024-03-26 07:26:19] 1 (9.1%) variants annotated as likely germline (DB INFO flag).
INFO [2024-03-26 07:26:20] Found 1000 variants in VCF file.
INFO [2024-03-26 07:26:20] Removing 2 triallelic sites.
WARN [2024-03-26 07:26:21] Having trouble guessing SOMATIC status...
WARN [2024-03-26 07:26:21] DP FORMAT field contains NAs. Removing 44 variants.
INFO [2024-03-26 07:26:21] 954 (100.0%) variants annotated as likely germline (DB INFO flag).
INFO [2024-03-26 07:26:26] Found 12 variants in VCF file.
INFO [2024-03-26 07:26:26] Removing 1 triallelic sites.
WARN [2024-03-26 07:26:26] DP FORMAT field contains NAs. Removing 1 variants.
WARN [2024-03-26 07:26:26] vcf.file has no DB info field for membership in germline databases. Found and used somatic status instead.
INFO [2024-03-26 07:26:26] 8 (80.0%) variants annotated as likely germline (DB INFO flag).
WARN [2024-03-26 07:26:27] BQ FORMAT field contains NAs.
INFO [2024-03-26 07:26:28] Found 2331 variants in VCF file.
INFO [2024-03-26 07:26:28] 2282 (97.9%) variants annotated as likely germline (DB INFO flag).
INFO [2024-03-26 07:26:29] Found 2331 variants in VCF file.
INFO [2024-03-26 07:26:30] 2282 (97.9%) variants annotated as likely germline (DB INFO flag).
WARN [2024-03-26 07:26:30] Did not find base quality scores, will use global error rate of 0.0010 instead.
INFO [2024-03-26 07:26:31] Removing 21 low quality variants with non-offset BQ < 25.
INFO [2024-03-26 07:26:31] Base quality scores range from 24 to 35 (offset by 1)
INFO [2024-03-26 07:26:31] Initial testing for significant sample cross-contamination: unlikely
INFO [2024-03-26 07:26:31] Removing 89 variants with AF < 0.030 or AF >= 0.970 or insufficient supporting reads or depth < 15.
INFO [2024-03-26 07:26:32] Found 2331 variants in VCF file.
INFO [2024-03-26 07:26:32] 2282 (97.9%) variants annotated as likely germline (DB INFO flag).
WARN [2024-03-26 07:26:32] BQ FORMAT field contains NAs.
WARN [2024-03-26 07:26:35] Allosome coverage missing, cannot determine sex.
WARN [2024-03-26 07:26:35] Allosome coverage missing, cannot determine sex.
INFO [2024-03-26 07:26:35] Mean coverages: chrX: 108.58, chrY: 108.58, chr1-22: 104.87.
INFO [2024-03-26 07:26:36] Mean coverages: chrX: 108.58, chrY: 2.17, chr1-22: 104.87.
INFO [2024-03-26 07:26:36] Mean coverages: chrX: 108.58, chrY: 5.17, chr1-22: 104.87.
INFO [2024-03-26 07:26:37] No homozygous variants in VCF, provide unfiltered VCF.
INFO [2024-03-26 07:26:37] No germline variants in VCF.
FATAL [2024-03-26 07:26:37] No solution with id hello 

FATAL [2024-03-26 07:26:37]  

FATAL [2024-03-26 07:26:37] This is most likely a user error due to invalid input data or 

FATAL [2024-03-26 07:26:37] parameters (PureCN 2.8.1). 

FATAL [2024-03-26 07:26:37] No solution with id 100 

FATAL [2024-03-26 07:26:37]  

FATAL [2024-03-26 07:26:37] This is most likely a user error due to invalid input data or 

FATAL [2024-03-26 07:26:37] parameters (PureCN 2.8.1). 

FATAL [2024-03-26 07:26:40] all.data and w have different lengths. 

FATAL [2024-03-26 07:26:40]  

FATAL [2024-03-26 07:26:40] This is most likely a user error due to invalid input data or 

FATAL [2024-03-26 07:26:40] parameters (PureCN 2.8.1). 

INFO [2024-03-26 07:26:42] Splitting 5 large targets to an average width of 400.
WARN [2024-03-26 07:26:42] No mappability scores provided.
WARN [2024-03-26 07:26:42] No reptiming scores provided.
INFO [2024-03-26 07:26:42] Calculating GC-content...
INFO [2024-03-26 07:26:43] Splitting 5 large targets to an average width of 400.
WARN [2024-03-26 07:26:43] No mappability scores provided.
WARN [2024-03-26 07:26:43] No reptiming scores provided.
INFO [2024-03-26 07:26:43] Calculating GC-content...
INFO [2024-03-26 07:26:44] Splitting 5 large targets to an average width of 400.
INFO [2024-03-26 07:26:44] Removing 1 targets overlapping with exclude.
WARN [2024-03-26 07:26:44] No mappability scores provided.
WARN [2024-03-26 07:26:44] No reptiming scores provided.
INFO [2024-03-26 07:26:44] Calculating GC-content...
WARN [2024-03-26 07:26:44] Found small target regions (< 100bp). Will resize them.
INFO [2024-03-26 07:26:45] Splitting 5 large targets to an average width of 400.
FATAL [2024-03-26 07:26:45] off.target.padding must be negative. 

FATAL [2024-03-26 07:26:45]  

FATAL [2024-03-26 07:26:45] This is most likely a user error due to invalid input data or 

FATAL [2024-03-26 07:26:45] parameters (PureCN 2.8.1). 

FATAL [2024-03-26 07:26:45] Interval coordinates should start at 1, not at 0 

FATAL [2024-03-26 07:26:45]  

FATAL [2024-03-26 07:26:45] This is most likely a user error due to invalid input data or 

FATAL [2024-03-26 07:26:45] parameters (PureCN 2.8.1). 

WARN [2024-03-26 07:26:45] Found small target regions (< 100bp). Will resize them.
INFO [2024-03-26 07:26:46] Splitting 5 large targets to an average width of 400.
FATAL [2024-03-26 07:26:46] No off-target regions after filtering for mappability and 

FATAL [2024-03-26 07:26:46] off.target.padding 

FATAL [2024-03-26 07:26:46]  

FATAL [2024-03-26 07:26:46] This is most likely a user error due to invalid input data or 

FATAL [2024-03-26 07:26:46] parameters (PureCN 2.8.1). 

WARN [2024-03-26 07:26:46] Found small target regions (< 100bp). Will resize them.
INFO [2024-03-26 07:26:47] Splitting 5 large targets to an average width of 400.
WARN [2024-03-26 07:26:47] No mappability scores provided.
INFO [2024-03-26 07:26:47] Calculating GC-content...
WARN [2024-03-26 07:26:47] Found small target regions (< 100bp). Will resize them.
INFO [2024-03-26 07:26:47] Averaging reptiming into bins of size 200...
INFO [2024-03-26 07:26:48] Splitting 5 large targets to an average width of 400.
WARN [2024-03-26 07:26:48] No mappability scores provided.
INFO [2024-03-26 07:26:48] Calculating GC-content...
WARN [2024-03-26 07:26:48] Found small target regions (< 100bp). Will resize them.
INFO [2024-03-26 07:26:49] Splitting 6 large targets to an average width of 200.
WARN [2024-03-26 07:26:49] No mappability scores provided.
WARN [2024-03-26 07:26:49] No reptiming scores provided.
INFO [2024-03-26 07:26:49] Calculating GC-content...
WARN [2024-03-26 07:26:49] Found small target regions (< 100bp). Will resize them.
INFO [2024-03-26 07:26:49] Splitting 5 large targets to an average width of 400.
WARN [2024-03-26 07:26:49] No mappability scores provided.
WARN [2024-03-26 07:26:49] No reptiming scores provided.
INFO [2024-03-26 07:26:50] Calculating GC-content...
INFO [2024-03-26 07:26:50] Splitting 5 large targets to an average width of 400.
INFO [2024-03-26 07:26:50] Tiling off-target regions to an average width of 200000.
WARN [2024-03-26 07:26:51] No mappability scores provided.
WARN [2024-03-26 07:26:51] No reptiming scores provided.
INFO [2024-03-26 07:26:51] Calculating GC-content...
WARN [2024-03-26 07:26:51] Intervals contain off-target regions. Will not change intervals.
WARN [2024-03-26 07:26:51] No mappability scores provided.
WARN [2024-03-26 07:26:51] No reptiming scores provided.
INFO [2024-03-26 07:26:51] Calculating GC-content...
INFO [2024-03-26 07:26:52] Splitting 5 large targets to an average width of 400.
WARN [2024-03-26 07:26:52] No reptiming scores provided.
INFO [2024-03-26 07:26:52] Calculating GC-content...
WARN [2024-03-26 07:26:52] Found small target regions (< 100bp). Will resize them.
INFO [2024-03-26 07:26:53] Splitting 5 large targets to an average width of 400.
WARN [2024-03-26 07:26:53] No reptiming scores provided.
INFO [2024-03-26 07:26:53] Calculating GC-content...
WARN [2024-03-26 07:26:55] Found small target regions (< 100bp). Will resize them.
INFO [2024-03-26 07:26:55] Splitting 5 large targets to an average width of 400.
WARN [2024-03-26 07:26:55] 1 intervals without mappability score (1 on-target).
INFO [2024-03-26 07:26:55] Removing 1 intervals with low mappability score (<0.60).
WARN [2024-03-26 07:26:55] No reptiming scores provided.
INFO [2024-03-26 07:26:55] Calculating GC-content...
WARN [2024-03-26 07:26:56] Found small target regions (< 100bp). Will resize them.
INFO [2024-03-26 07:26:56] Splitting 5 large targets to an average width of 400.
INFO [2024-03-26 07:26:57] Tiling off-target regions to an average width of 200000.
WARN [2024-03-26 07:26:57] No reptiming scores provided.
INFO [2024-03-26 07:26:57] Calculating GC-content...
INFO [2024-03-26 07:26:59] Splitting 5 large targets to an average width of 400.
WARN [2024-03-26 07:26:59] No mappability scores provided.
WARN [2024-03-26 07:26:59] No reptiming scores provided.
INFO [2024-03-26 07:26:59] Calculating GC-content...
WARN [2024-03-26 07:26:59] Chromosome naming style of interval file (Ensembl) was different from reference (UCSC).
INFO [2024-03-26 07:27:00] Splitting 5 large targets to an average width of 400.
WARN [2024-03-26 07:27:00] No mappability scores provided.
WARN [2024-03-26 07:27:00] No reptiming scores provided.
INFO [2024-03-26 07:27:00] Calculating GC-content...
FATAL [2024-03-26 07:27:00] Chromosome naming style of interval file unknown, should be UCSC. 

FATAL [2024-03-26 07:27:00]  

FATAL [2024-03-26 07:27:00] This is most likely a user error due to invalid input data or 

FATAL [2024-03-26 07:27:00] parameters (PureCN 2.8.1). 

WARN [2024-03-26 07:27:00] Chromosome naming style of interval file (Ensembl) was different from reference (UCSC).
WARN [2024-03-26 07:27:01] Chromosome naming style of mappability file (Ensembl) was different from reference (UCSC).
INFO [2024-03-26 07:27:01] Splitting 5 large targets to an average width of 400.
WARN [2024-03-26 07:27:01] No reptiming scores provided.
INFO [2024-03-26 07:27:01] Calculating GC-content...
WARN [2024-03-26 07:27:02] Found small target regions (< 60bp). Will resize them.
INFO [2024-03-26 07:27:02] Splitting 5 large targets to an average width of 400.
WARN [2024-03-26 07:27:02] No mappability scores provided.
WARN [2024-03-26 07:27:02] No reptiming scores provided.
INFO [2024-03-26 07:27:02] Calculating GC-content...
WARN [2024-03-26 07:27:02] Found small target regions (< 60bp). Will drop them.
INFO [2024-03-26 07:27:02] Splitting 4 large targets to an average width of 400.
WARN [2024-03-26 07:27:02] No mappability scores provided.
WARN [2024-03-26 07:27:02] No reptiming scores provided.
INFO [2024-03-26 07:27:02] Calculating GC-content...
WARN [2024-03-26 07:27:03] Found small target regions (< 200bp). Will resize them.
INFO [2024-03-26 07:27:03] Splitting 5 large targets to an average width of 400.
WARN [2024-03-26 07:27:03] No mappability scores provided.
WARN [2024-03-26 07:27:03] No reptiming scores provided.
INFO [2024-03-26 07:27:03] Calculating GC-content...
INFO [2024-03-26 07:27:05] Found 20 variants in VCF file.
INFO [2024-03-26 07:27:06] 20 (100.0%) variants annotated as likely germline (DB INFO flag).
WARN [2024-03-26 07:27:06] Did not find base quality scores, will use global error rate of 0.0010 instead.
INFO [2024-03-26 07:27:09] Found 127 variants in VCF file.
INFO [2024-03-26 07:27:09] 127 (100.0%) variants annotated as likely germline (DB INFO flag).
WARN [2024-03-26 07:27:09] Did not find base quality scores, will use global error rate of 0.0010 instead.
INFO [2024-03-26 07:27:10] ------------------------------------------------------------
INFO [2024-03-26 07:27:10] PureCN 2.8.1
INFO [2024-03-26 07:27:10] ------------------------------------------------------------
INFO [2024-03-26 07:27:10] Loading coverage files...
INFO [2024-03-26 07:27:11] Mean target coverages: 112X (tumor) 99X (normal).
WARN [2024-03-26 07:27:11] Allosome coverage missing, cannot determine sex.
WARN [2024-03-26 07:27:11] Allosome coverage missing, cannot determine sex.
INFO [2024-03-26 07:27:11] Removing 228 intervals with missing log.ratio.
INFO [2024-03-26 07:27:11] Removing 22 small (< 5bp) intervals.
INFO [2024-03-26 07:27:11] Removing 19 intervals with low total coverage in normal (< 150.00 reads).
WARN [2024-03-26 07:27:11] No normalDB provided. Provide one for better results.
INFO [2024-03-26 07:27:11] Removing 1724 low count (< 100 total reads) intervals.
INFO [2024-03-26 07:27:11] Removing 6 low coverage (< 15.0000X) intervals.
INFO [2024-03-26 07:27:11] Using 8050 intervals (8050 on-target, 0 off-target).
INFO [2024-03-26 07:27:11] No off-target intervals. If this is hybrid-capture data, consider adding them.
INFO [2024-03-26 07:27:11] No interval.file provided. Cannot check for any GC-biases.
INFO [2024-03-26 07:27:11] Loading VCF...
INFO [2024-03-26 07:27:12] Found 127 variants in VCF file.
INFO [2024-03-26 07:27:12] 127 (100.0%) variants annotated as likely germline (DB INFO flag).
INFO [2024-03-26 07:27:12] LIB-02240e4 is tumor in VCF file.
INFO [2024-03-26 07:27:12] No homozygous variants in VCF, provide unfiltered VCF.
INFO [2024-03-26 07:27:13] Removing 1 non heterozygous (in matched normal) germline SNPs.
INFO [2024-03-26 07:27:13] Removing 0 low quality variants with non-offset BQ < 25.
INFO [2024-03-26 07:27:13] Base quality scores range from 29 to 35 (offset by 1)
INFO [2024-03-26 07:27:13] Minimum number of supporting reads ranges from 2 to 5, depending on coverage and BQS.
INFO [2024-03-26 07:27:14] Removing 1 variants with AF < 0.030 or AF >= 1.000 or insufficient supporting reads or depth < 15.
INFO [2024-03-26 07:27:14] Total size of targeted genomic region: 1.53Mb (2.33Mb with 50bp padding).
INFO [2024-03-26 07:27:14] 1.2% of targets contain variants.
INFO [2024-03-26 07:27:15] Removing 4 variants outside intervals.
INFO [2024-03-26 07:27:15] Found SOMATIC annotation in VCF.
INFO [2024-03-26 07:27:15] Setting somatic prior probabilities for somatic variants to 0.999000 or to 0.000100 otherwise.
INFO [2024-03-26 07:27:15] Found SOMATIC annotation in VCF. Setting mapping bias to 0.972.
INFO [2024-03-26 07:27:15] Excluding 0 novel or poor quality variants from segmentation.
INFO [2024-03-26 07:27:15] Sample sex: ?
INFO [2024-03-26 07:27:15] Segmenting data...
INFO [2024-03-26 07:27:15] Loading pre-computed boundaries for DNAcopy...
INFO [2024-03-26 07:27:15] Setting undo.SD parameter to 1.000000.
INFO [2024-03-26 07:27:17] Setting prune.hclust.h parameter to 0.200000.
INFO [2024-03-26 07:27:17] Found 59 segments with median size of 17.67Mb.
INFO [2024-03-26 07:27:18] Using 121 variants.
INFO [2024-03-26 07:27:18] Mean standard deviation of log-ratios: 0.34 (MAPD: 0.30)
INFO [2024-03-26 07:27:18] 2D-grid search of purity and ploidy...
INFO [2024-03-26 07:27:18] Local optima: 0.65/1.8, 0.52/2
INFO [2024-03-26 07:27:18] Testing local optimum 1/2 at purity 0.65 and total ploidy 1.80...
INFO [2024-03-26 07:27:20] Testing local optimum 2/2 at purity 0.52 and total ploidy 2.00...
INFO [2024-03-26 07:27:22] Skipping 1 solutions that converged to the same optima.
INFO [2024-03-26 07:27:22] Fitting variants with beta model for local optimum 1/2...
INFO [2024-03-26 07:27:22] Fitting variants for purity 0.65, tumor ploidy 1.75 and contamination 0.01.
INFO [2024-03-26 07:27:30] Optimized purity: 0.65
INFO [2024-03-26 07:27:30] Done.
INFO [2024-03-26 07:27:30] ------------------------------------------------------------
FATAL [2024-03-26 07:27:34] log.ratio NULL in .writeLogRatioFileGATK4 

FATAL [2024-03-26 07:27:34]  

FATAL [2024-03-26 07:27:34] This runtime error might be caused by invalid input data or parameters. 

FATAL [2024-03-26 07:27:34] Please report bug (PureCN 2.8.1). 

INFO [2024-03-26 07:27:34] Loaded provided segmentation file example_seg.txt (format DNAcopy).
INFO [2024-03-26 07:27:34] Re-centering provided segment means (offset -0.0033).
INFO [2024-03-26 07:27:34] Loaded provided segmentation file example_gatk4_modelfinal.seg.gz (format GATK4).
WARN [2024-03-26 07:27:34] Expecting numeric chromosome names in seg.file, assuming file is properly sorted.
INFO [2024-03-26 07:27:34] Re-centering provided segment means (offset -0.0037).
INFO [2024-03-26 07:27:36] 576 on-target bins with low coverage in all samples.
WARN [2024-03-26 07:27:36] You are likely not using the correct baits file!
WARN [2024-03-26 07:27:36] Allosome coverage missing, cannot determine sex.
WARN [2024-03-26 07:27:36] Allosome coverage missing, cannot determine sex.
INFO [2024-03-26 07:27:36] Processing on-target regions...
INFO [2024-03-26 07:27:37] Removing 930 intervals with low coverage in normalDB.
INFO [2024-03-26 07:27:37] Removing 1 intervals with zero coverage in more than 3% of normalDB.
INFO [2024-03-26 07:27:42] ------------------------------------------------------------
INFO [2024-03-26 07:27:42] PureCN 2.8.1
INFO [2024-03-26 07:27:42] ------------------------------------------------------------
INFO [2024-03-26 07:27:42] Using BiocParallel for parallel optimization.
INFO [2024-03-26 07:27:42] Loading coverage files...
INFO [2024-03-26 07:27:43] Mean target coverages: 112X (tumor) 99X (normal).
WARN [2024-03-26 07:27:43] Allosome coverage missing, cannot determine sex.
WARN [2024-03-26 07:27:43] Allosome coverage missing, cannot determine sex.
INFO [2024-03-26 07:27:43] Removing 228 intervals with missing log.ratio.
INFO [2024-03-26 07:27:43] Removing 705 intervals excluded in normalDB.
INFO [2024-03-26 07:27:43] normalDB provided. Setting minimum coverage for segmentation to 0.0015X.
INFO [2024-03-26 07:27:44] Removing 1066 low count (< 100 total reads) intervals.
INFO [2024-03-26 07:27:44] Using 8050 intervals (8050 on-target, 0 off-target).
INFO [2024-03-26 07:27:44] No off-target intervals. If this is hybrid-capture data, consider adding them.
INFO [2024-03-26 07:27:44] No interval.file provided. Cannot check for any GC-biases.
INFO [2024-03-26 07:27:44] Sample sex: ?
INFO [2024-03-26 07:27:44] Segmenting data...
INFO [2024-03-26 07:27:44] Interval weights found, will use weighted CBS.
INFO [2024-03-26 07:27:44] Loading pre-computed boundaries for DNAcopy...
INFO [2024-03-26 07:27:44] Setting undo.SD parameter to 1.000000.
Setting multi-figure configuration
INFO [2024-03-26 07:27:47] Found 52 segments with median size of 29.35Mb.
INFO [2024-03-26 07:27:47] Mean standard deviation of log-ratios: 0.34 (MAPD: 0.30)
INFO [2024-03-26 07:27:47] 2D-grid search of purity and ploidy...
INFO [2024-03-26 07:27:47] Local optima: 0.65/1.8, 0.52/2
INFO [2024-03-26 07:27:52] Skipping 1 solutions that converged to the same optima.
INFO [2024-03-26 07:27:52] Done.
INFO [2024-03-26 07:27:52] ------------------------------------------------------------
INFO [2024-03-26 07:27:52] Reading /tmp/RtmpP77MMO/file91855fbb47fd.rds...
FATAL [2024-03-26 07:27:52] runAbsoluteCN was run without a VCF file. 

FATAL [2024-03-26 07:27:52]  

FATAL [2024-03-26 07:27:52] This is most likely a user error due to invalid input data or 

FATAL [2024-03-26 07:27:52] parameters (PureCN 2.8.1). 

FATAL [2024-03-26 07:27:52] runAbsoluteCN was run without a VCF file. 

FATAL [2024-03-26 07:27:52]  

FATAL [2024-03-26 07:27:52] This is most likely a user error due to invalid input data or 

FATAL [2024-03-26 07:27:52] parameters (PureCN 2.8.1). 

INFO [2024-03-26 07:27:52] ------------------------------------------------------------
INFO [2024-03-26 07:27:52] PureCN 2.8.1
INFO [2024-03-26 07:27:52] ------------------------------------------------------------
INFO [2024-03-26 07:27:52] Loading coverage files...
INFO [2024-03-26 07:27:54] Mean target coverages: 112X (tumor) 99X (normal).
WARN [2024-03-26 07:27:54] Allosome coverage missing, cannot determine sex.
WARN [2024-03-26 07:27:54] Allosome coverage missing, cannot determine sex.
INFO [2024-03-26 07:27:54] Removing 228 intervals with missing log.ratio.
INFO [2024-03-26 07:27:54] Removing 705 intervals excluded in normalDB.
INFO [2024-03-26 07:27:54] normalDB provided. Setting minimum coverage for segmentation to 0.0015X.
INFO [2024-03-26 07:27:54] Removing 1066 low count (< 100 total reads) intervals.
INFO [2024-03-26 07:27:55] Using 8050 intervals (8050 on-target, 0 off-target).
INFO [2024-03-26 07:27:55] No off-target intervals. If this is hybrid-capture data, consider adding them.
INFO [2024-03-26 07:27:55] No interval.file provided. Cannot check for any GC-biases.
INFO [2024-03-26 07:27:55] Sample sex: ?
INFO [2024-03-26 07:27:55] Segmenting data...
INFO [2024-03-26 07:27:55] Interval weights found, will use weighted PSCBS.
FATAL [2024-03-26 07:27:55] segmentationPSCBS requires VCF file. 

FATAL [2024-03-26 07:27:55]  

FATAL [2024-03-26 07:27:55] This is most likely a user error due to invalid input data or 

FATAL [2024-03-26 07:27:55] parameters (PureCN 2.8.1). 

INFO [2024-03-26 07:28:09] ------------------------------------------------------------
INFO [2024-03-26 07:28:09] PureCN 2.8.1
INFO [2024-03-26 07:28:09] ------------------------------------------------------------
INFO [2024-03-26 07:28:09] Loading coverage files...
FATAL [2024-03-26 07:28:10] Need a normal coverage file if log.ratio and seg.file are not provided. 

FATAL [2024-03-26 07:28:10]  

FATAL [2024-03-26 07:28:10] This is most likely a user error due to invalid input data or 

FATAL [2024-03-26 07:28:10] parameters (PureCN 2.8.1). 

INFO [2024-03-26 07:28:10] ------------------------------------------------------------
INFO [2024-03-26 07:28:10] PureCN 2.8.1
INFO [2024-03-26 07:28:10] ------------------------------------------------------------
FATAL [2024-03-26 07:28:10] min.ploidy or max.ploidy not within expected range. 

FATAL [2024-03-26 07:28:10]  

FATAL [2024-03-26 07:28:10] This is most likely a user error due to invalid input data or 

FATAL [2024-03-26 07:28:10] parameters (PureCN 2.8.1). 

INFO [2024-03-26 07:28:10] ------------------------------------------------------------
INFO [2024-03-26 07:28:10] PureCN 2.8.1
INFO [2024-03-26 07:28:10] ------------------------------------------------------------
FATAL [2024-03-26 07:28:10] min.ploidy or max.ploidy not within expected range. 

FATAL [2024-03-26 07:28:10]  

FATAL [2024-03-26 07:28:10] This is most likely a user error due to invalid input data or 

FATAL [2024-03-26 07:28:10] parameters (PureCN 2.8.1). 

INFO [2024-03-26 07:28:10] ------------------------------------------------------------
INFO [2024-03-26 07:28:10] PureCN 2.8.1
INFO [2024-03-26 07:28:10] ------------------------------------------------------------
INFO [2024-03-26 07:28:10] ------------------------------------------------------------
INFO [2024-03-26 07:28:10] PureCN 2.8.1
INFO [2024-03-26 07:28:10] ------------------------------------------------------------
INFO [2024-03-26 07:28:10] ------------------------------------------------------------
INFO [2024-03-26 07:28:10] PureCN 2.8.1
INFO [2024-03-26 07:28:10] ------------------------------------------------------------
FATAL [2024-03-26 07:28:10] test.num.copy not within expected range. 

FATAL [2024-03-26 07:28:10]  

FATAL [2024-03-26 07:28:10] This is most likely a user error due to invalid input data or 

FATAL [2024-03-26 07:28:10] parameters (PureCN 2.8.1). 

INFO [2024-03-26 07:28:10] ------------------------------------------------------------
INFO [2024-03-26 07:28:10] PureCN 2.8.1
INFO [2024-03-26 07:28:10] ------------------------------------------------------------
WARN [2024-03-26 07:28:10] test.num.copy outside recommended range.
FATAL [2024-03-26 07:28:10] max.non.clonal not within expected range or format. 

FATAL [2024-03-26 07:28:10]  

FATAL [2024-03-26 07:28:10] This is most likely a user error due to invalid input data or 

FATAL [2024-03-26 07:28:10] parameters (PureCN 2.8.1). 

INFO [2024-03-26 07:28:10] ------------------------------------------------------------
INFO [2024-03-26 07:28:10] PureCN 2.8.1
INFO [2024-03-26 07:28:10] ------------------------------------------------------------
WARN [2024-03-26 07:28:10] test.num.copy outside recommended range.
FATAL [2024-03-26 07:28:10] max.non.clonal not within expected range or format. 

FATAL [2024-03-26 07:28:10]  

FATAL [2024-03-26 07:28:10] This is most likely a user error due to invalid input data or 

FATAL [2024-03-26 07:28:10] parameters (PureCN 2.8.1). 

INFO [2024-03-26 07:28:10] ------------------------------------------------------------
INFO [2024-03-26 07:28:10] PureCN 2.8.1
INFO [2024-03-26 07:28:10] ------------------------------------------------------------
FATAL [2024-03-26 07:28:10] test.purity not within expected range. 

FATAL [2024-03-26 07:28:10]  

FATAL [2024-03-26 07:28:10] This is most likely a user error due to invalid input data or 

FATAL [2024-03-26 07:28:10] parameters (PureCN 2.8.1). 

INFO [2024-03-26 07:28:10] ------------------------------------------------------------
INFO [2024-03-26 07:28:10] PureCN 2.8.1
INFO [2024-03-26 07:28:10] ------------------------------------------------------------
INFO [2024-03-26 07:28:10] Loading coverage files...
FATAL [2024-03-26 07:28:12] Tumor and normal are identical. This won't give any meaningful results 

FATAL [2024-03-26 07:28:12] and I'm stopping here. 

FATAL [2024-03-26 07:28:12]  

FATAL [2024-03-26 07:28:12] This is most likely a user error due to invalid input data or 

FATAL [2024-03-26 07:28:12] parameters (PureCN 2.8.1). 

INFO [2024-03-26 07:28:12] ------------------------------------------------------------
INFO [2024-03-26 07:28:12] PureCN 2.8.1
INFO [2024-03-26 07:28:12] ------------------------------------------------------------
INFO [2024-03-26 07:28:12] Loading coverage files...
FATAL [2024-03-26 07:28:12] Length of log.ratio different from tumor coverage. 

FATAL [2024-03-26 07:28:12]  

FATAL [2024-03-26 07:28:12] This is most likely a user error due to invalid input data or 

FATAL [2024-03-26 07:28:12] parameters (PureCN 2.8.1). 

INFO [2024-03-26 07:28:12] ------------------------------------------------------------
INFO [2024-03-26 07:28:12] PureCN 2.8.1
INFO [2024-03-26 07:28:12] ------------------------------------------------------------
FATAL [2024-03-26 07:28:12] prior.purity must have the same length as test.purity. 

FATAL [2024-03-26 07:28:12]  

FATAL [2024-03-26 07:28:12] This is most likely a user error due to invalid input data or 

FATAL [2024-03-26 07:28:12] parameters (PureCN 2.8.1). 

INFO [2024-03-26 07:28:12] ------------------------------------------------------------
INFO [2024-03-26 07:28:12] PureCN 2.8.1
INFO [2024-03-26 07:28:13] ------------------------------------------------------------
FATAL [2024-03-26 07:28:13] min.gof not within expected range or format. 

FATAL [2024-03-26 07:28:13]  

FATAL [2024-03-26 07:28:13] This is most likely a user error due to invalid input data or 

FATAL [2024-03-26 07:28:13] parameters (PureCN 2.8.1). 

INFO [2024-03-26 07:28:13] ------------------------------------------------------------
INFO [2024-03-26 07:28:13] PureCN 2.8.1
INFO [2024-03-26 07:28:13] ------------------------------------------------------------
FATAL [2024-03-26 07:28:13] prior.purity not within expected range or format. 

FATAL [2024-03-26 07:28:13]  

FATAL [2024-03-26 07:28:13] This is most likely a user error due to invalid input data or 

FATAL [2024-03-26 07:28:13] parameters (PureCN 2.8.1). 

INFO [2024-03-26 07:28:13] ------------------------------------------------------------
INFO [2024-03-26 07:28:13] PureCN 2.8.1
INFO [2024-03-26 07:28:13] ------------------------------------------------------------
FATAL [2024-03-26 07:28:13] prior.purity must add to 1. Sum is 1.5 

FATAL [2024-03-26 07:28:13]  

FATAL [2024-03-26 07:28:13] This is most likely a user error due to invalid input data or 

FATAL [2024-03-26 07:28:13] parameters (PureCN 2.8.1). 

INFO [2024-03-26 07:28:13] ------------------------------------------------------------
INFO [2024-03-26 07:28:13] PureCN 2.8.1
INFO [2024-03-26 07:28:13] ------------------------------------------------------------
FATAL [2024-03-26 07:28:13] max.homozygous.loss not within expected range or format. 

FATAL [2024-03-26 07:28:13]  

FATAL [2024-03-26 07:28:13] This is most likely a user error due to invalid input data or 

FATAL [2024-03-26 07:28:13] parameters (PureCN 2.8.1). 

INFO [2024-03-26 07:28:13] ------------------------------------------------------------
INFO [2024-03-26 07:28:13] PureCN 2.8.1
INFO [2024-03-26 07:28:13] ------------------------------------------------------------
FATAL [2024-03-26 07:28:13] prior.K not within expected range or format. 

FATAL [2024-03-26 07:28:13]  

FATAL [2024-03-26 07:28:13] This is most likely a user error due to invalid input data or 

FATAL [2024-03-26 07:28:13] parameters (PureCN 2.8.1). 

INFO [2024-03-26 07:28:13] ------------------------------------------------------------
INFO [2024-03-26 07:28:13] PureCN 2.8.1
INFO [2024-03-26 07:28:13] ------------------------------------------------------------
FATAL [2024-03-26 07:28:13] prior.contamination not within expected range or format. 

FATAL [2024-03-26 07:28:13]  

FATAL [2024-03-26 07:28:13] This is most likely a user error due to invalid input data or 

FATAL [2024-03-26 07:28:13] parameters (PureCN 2.8.1). 

INFO [2024-03-26 07:28:13] ------------------------------------------------------------
INFO [2024-03-26 07:28:13] PureCN 2.8.1
INFO [2024-03-26 07:28:13] ------------------------------------------------------------
FATAL [2024-03-26 07:28:13] Iterations not in the expected range from 10 to 250. 

FATAL [2024-03-26 07:28:13]  

FATAL [2024-03-26 07:28:13] This is most likely a user error due to invalid input data or 

FATAL [2024-03-26 07:28:13] parameters (PureCN 2.8.1). 

INFO [2024-03-26 07:28:13] ------------------------------------------------------------
INFO [2024-03-26 07:28:13] PureCN 2.8.1
INFO [2024-03-26 07:28:13] ------------------------------------------------------------
FATAL [2024-03-26 07:28:13] Iterations not in the expected range from 10 to 250. 

FATAL [2024-03-26 07:28:13]  

FATAL [2024-03-26 07:28:13] This is most likely a user error due to invalid input data or 

FATAL [2024-03-26 07:28:13] parameters (PureCN 2.8.1). 

INFO [2024-03-26 07:28:13] ------------------------------------------------------------
INFO [2024-03-26 07:28:13] PureCN 2.8.1
INFO [2024-03-26 07:28:13] ------------------------------------------------------------
INFO [2024-03-26 07:28:13] Loading coverage files...
FATAL [2024-03-26 07:28:14] Missing tumor.coverage.file requires seg.file or log.ratio and 

FATAL [2024-03-26 07:28:14] interval.file. 

FATAL [2024-03-26 07:28:14]  

FATAL [2024-03-26 07:28:14] This is most likely a user error due to invalid input data or 

FATAL [2024-03-26 07:28:14] parameters (PureCN 2.8.1). 

INFO [2024-03-26 07:28:14] ------------------------------------------------------------
INFO [2024-03-26 07:28:14] PureCN 2.8.1
INFO [2024-03-26 07:28:14] ------------------------------------------------------------
INFO [2024-03-26 07:28:14] ------------------------------------------------------------
INFO [2024-03-26 07:28:14] PureCN 2.8.1
INFO [2024-03-26 07:28:14] ------------------------------------------------------------
INFO [2024-03-26 07:28:14] Loading coverage files...
FATAL [2024-03-26 07:28:15] Interval files in normal and tumor different. 

FATAL [2024-03-26 07:28:15]  

FATAL [2024-03-26 07:28:15] This is most likely a user error due to invalid input data or 

FATAL [2024-03-26 07:28:15] parameters (PureCN 2.8.1). 

INFO [2024-03-26 07:28:16] ------------------------------------------------------------
INFO [2024-03-26 07:28:16] PureCN 2.8.1
INFO [2024-03-26 07:28:16] ------------------------------------------------------------
INFO [2024-03-26 07:28:16] Loading coverage files...
INFO [2024-03-26 07:28:16] Mean target coverages: 0X (tumor) 99X (normal).
WARN [2024-03-26 07:28:16] Large difference in coverage of tumor and normal.
FATAL [2024-03-26 07:28:16] No finite intervals. 

FATAL [2024-03-26 07:28:16]  

FATAL [2024-03-26 07:28:16] This is most likely a user error due to invalid input data or 

FATAL [2024-03-26 07:28:16] parameters (PureCN 2.8.1). 

INFO [2024-03-26 07:28:17] ------------------------------------------------------------
INFO [2024-03-26 07:28:17] PureCN 2.8.1
INFO [2024-03-26 07:28:17] ------------------------------------------------------------
INFO [2024-03-26 07:28:17] Loading coverage files...
INFO [2024-03-26 07:28:18] Mean target coverages: 112X (tumor) 99X (normal).
WARN [2024-03-26 07:28:18] Allosome coverage missing, cannot determine sex.
WARN [2024-03-26 07:28:18] Allosome coverage missing, cannot determine sex.
INFO [2024-03-26 07:28:18] Removing 228 intervals with missing log.ratio.
INFO [2024-03-26 07:28:18] Removing 22 small (< 5bp) intervals.
INFO [2024-03-26 07:28:18] Removing 19 intervals with low total coverage in normal (< 150.00 reads).
WARN [2024-03-26 07:28:18] No normalDB provided. Provide one for better results.
INFO [2024-03-26 07:28:18] Removing 1724 low count (< 100 total reads) intervals.
INFO [2024-03-26 07:28:18] Removing 6 low coverage (< 15.0000X) intervals.
INFO [2024-03-26 07:28:18] Using 8050 intervals (8050 on-target, 0 off-target).
INFO [2024-03-26 07:28:18] No off-target intervals. If this is hybrid-capture data, consider adding them.
INFO [2024-03-26 07:28:18] No interval.file provided. Cannot check for any GC-biases.
INFO [2024-03-26 07:28:18] Loading VCF...
INFO [2024-03-26 07:28:19] Found 127 variants in VCF file.
INFO [2024-03-26 07:28:19] 127 (100.0%) variants annotated as likely germline (DB INFO flag).
INFO [2024-03-26 07:28:20] LIB-02240e4 is tumor in VCF file.
INFO [2024-03-26 07:28:20] No homozygous variants in VCF, provide unfiltered VCF.
INFO [2024-03-26 07:28:20] Removing 1 non heterozygous (in matched normal) germline SNPs.
INFO [2024-03-26 07:28:20] Removing 0 low quality variants with non-offset BQ < 25.
INFO [2024-03-26 07:28:20] Base quality scores range from 29 to 35 (offset by 1)
INFO [2024-03-26 07:28:20] Minimum number of supporting reads ranges from 2 to 5, depending on coverage and BQS.
INFO [2024-03-26 07:28:21] Removing 1 variants with AF < 0.030 or AF >= 1.000 or insufficient supporting reads or depth < 15.
Error in FileForFormat(resource(con), fileFormat(con)) : 
  Format 'txt' unsupported
In addition: Warning messages:
1: In for (getMethodName in getMethodNames) { :
  closing unused connection 5 (/tmp/RtmpP77MMO/file91857989a08d.tsv)
2: In for (getMethodName in getMethodNames) { :
  closing unused connection 4 (/Library/Frameworks/R.framework/Versions/4.3-x86_64/Resources/library/PureCN/extdata/example_gatk4_denoised_cr.tsv.gz)
FATAL [2024-03-26 07:28:21] Could not import snp.blacklist 

FATAL [2024-03-26 07:28:21] /Library/Frameworks/R.framework/Versions/4.3-x86_64/Resources/library/PureCN/extdata/example_normal.txt.gz:Error 

FATAL [2024-03-26 07:28:21] in FileForFormat(resource(con), fileFormat(con)) : Format 'txt' 

FATAL [2024-03-26 07:28:21] unsupported 

FATAL [2024-03-26 07:28:21]  

FATAL [2024-03-26 07:28:21] This is most likely a user error due to invalid input data or 

FATAL [2024-03-26 07:28:21] parameters (PureCN 2.8.1). 

INFO [2024-03-26 07:28:21] ------------------------------------------------------------
INFO [2024-03-26 07:28:21] PureCN 2.8.1
INFO [2024-03-26 07:28:21] ------------------------------------------------------------
INFO [2024-03-26 07:28:21] Loading coverage files...
INFO [2024-03-26 07:28:23] Mean target coverages: 112X (tumor) 99X (normal).
WARN [2024-03-26 07:28:23] Allosome coverage missing, cannot determine sex.
WARN [2024-03-26 07:28:23] Allosome coverage missing, cannot determine sex.
INFO [2024-03-26 07:28:23] Removing 228 intervals with missing log.ratio.
INFO [2024-03-26 07:28:23] Removing 22 small (< 5bp) intervals.
INFO [2024-03-26 07:28:23] Removing 19 intervals with low total coverage in normal (< 150.00 reads).
WARN [2024-03-26 07:28:23] No normalDB provided. Provide one for better results.
INFO [2024-03-26 07:28:23] Removing 233 low coverage (< 15.0000X) intervals.
INFO [2024-03-26 07:28:23] Using 9547 intervals (9547 on-target, 0 off-target).
INFO [2024-03-26 07:28:23] No off-target intervals. If this is hybrid-capture data, consider adding them.
INFO [2024-03-26 07:28:23] No interval.file provided. Cannot check for any GC-biases.
INFO [2024-03-26 07:28:23] Loading VCF...
INFO [2024-03-26 07:28:24] Found 127 variants in VCF file.
INFO [2024-03-26 07:28:24] 127 (100.0%) variants annotated as likely germline (DB INFO flag).
INFO [2024-03-26 07:28:24] LIB-02240e4 is tumor in VCF file.
INFO [2024-03-26 07:28:25] No homozygous variants in VCF, provide unfiltered VCF.
INFO [2024-03-26 07:28:25] Removing 1 non heterozygous (in matched normal) germline SNPs.
INFO [2024-03-26 07:28:25] Removing 0 low quality variants with non-offset BQ < 25.
INFO [2024-03-26 07:28:25] Base quality scores range from 29 to 35 (offset by 1)
INFO [2024-03-26 07:28:25] Minimum number of supporting reads ranges from 2 to 5, depending on coverage and BQS.
INFO [2024-03-26 07:28:26] Removing 1 variants with AF < 0.030 or AF >= 1.000 or insufficient supporting reads or depth < 15.
INFO [2024-03-26 07:28:26] Total size of targeted genomic region: 1.63Mb (2.58Mb with 50bp padding).
INFO [2024-03-26 07:28:27] 1.0% of targets contain variants.
INFO [2024-03-26 07:28:27] Removing 2 variants outside intervals.
INFO [2024-03-26 07:28:27] Found SOMATIC annotation in VCF.
INFO [2024-03-26 07:28:27] Setting somatic prior probabilities for somatic variants to 0.999000 or to 0.000100 otherwise.
INFO [2024-03-26 07:28:27] Found SOMATIC annotation in VCF. Setting mapping bias to 0.973.
INFO [2024-03-26 07:28:27] Excluding 0 novel or poor quality variants from segmentation.
INFO [2024-03-26 07:28:28] Sample sex: ?
INFO [2024-03-26 07:28:28] Segmenting data...
INFO [2024-03-26 07:28:28] Loading pre-computed boundaries for DNAcopy...
INFO [2024-03-26 07:28:28] Setting undo.SD parameter to 1.000000.
INFO [2024-03-26 07:28:29] Setting prune.hclust.h parameter to 0.200000.
INFO [2024-03-26 07:28:29] Found 54 segments with median size of 24.88Mb.
INFO [2024-03-26 07:28:30] Using 123 variants.
INFO [2024-03-26 07:28:30] Mean standard deviation of log-ratios: 0.37 (MAPD: 0.32)
INFO [2024-03-26 07:28:30] 2D-grid search of purity and ploidy...
INFO [2024-03-26 07:28:32] Local optima: 0.62/1.9, 0.3/1.9, 0.32/2
INFO [2024-03-26 07:28:32] Testing local optimum 1/3 at purity 0.62 and total ploidy 1.90...
INFO [2024-03-26 07:28:35] Testing local optimum 2/3 at purity 0.30 and total ploidy 1.90...
INFO [2024-03-26 07:28:36] Recalibrating log-ratios...
INFO [2024-03-26 07:28:36] Testing local optimum 2/3 at purity 0.30 and total ploidy 1.90...
INFO [2024-03-26 07:28:37] Recalibrating log-ratios...
INFO [2024-03-26 07:28:37] Testing local optimum 2/3 at purity 0.30 and total ploidy 1.90...
INFO [2024-03-26 07:28:37] Recalibrating log-ratios...
INFO [2024-03-26 07:28:37] Testing local optimum 2/3 at purity 0.30 and total ploidy 1.90...
INFO [2024-03-26 07:28:38] Testing local optimum 3/3 at purity 0.32 and total ploidy 2.00...
INFO [2024-03-26 07:28:39] Recalibrating log-ratios...
INFO [2024-03-26 07:28:39] Testing local optimum 3/3 at purity 0.32 and total ploidy 2.00...
INFO [2024-03-26 07:28:40] Recalibrating log-ratios...
INFO [2024-03-26 07:28:40] Testing local optimum 3/3 at purity 0.32 and total ploidy 2.00...
INFO [2024-03-26 07:28:41] Recalibrating log-ratios...
INFO [2024-03-26 07:28:41] Testing local optimum 3/3 at purity 0.32 and total ploidy 2.00...
INFO [2024-03-26 07:28:42] Skipping 1 solutions that converged to the same optima.
INFO [2024-03-26 07:28:42] Skipping 1 solutions exceeding max.non.clonal (0.20): 0.3/1.9 (purity/tumor ploidy)
INFO [2024-03-26 07:28:42] Fitting variants with beta model for local optimum 1/3...
INFO [2024-03-26 07:28:42] Fitting variants for purity 0.65, tumor ploidy 1.74 and contamination 0.01.
INFO [2024-03-26 07:28:49] Optimized purity: 0.65
INFO [2024-03-26 07:28:49] Done.
INFO [2024-03-26 07:28:49] ------------------------------------------------------------
INFO [2024-03-26 07:28:50] ------------------------------------------------------------
INFO [2024-03-26 07:28:50] PureCN 2.8.1
INFO [2024-03-26 07:28:50] ------------------------------------------------------------
INFO [2024-03-26 07:28:50] Loading coverage files...
INFO [2024-03-26 07:28:51] Mean target coverages: 112X (tumor) 99X (normal).
WARN [2024-03-26 07:28:51] Allosome coverage missing, cannot determine sex.
WARN [2024-03-26 07:28:52] Allosome coverage missing, cannot determine sex.
INFO [2024-03-26 07:28:52] Removing 228 intervals with missing log.ratio.
INFO [2024-03-26 07:28:52] Removing 22 small (< 5bp) intervals.
INFO [2024-03-26 07:28:52] Removing 19 intervals with low total coverage in normal (< 150.00 reads).
WARN [2024-03-26 07:28:52] No normalDB provided. Provide one for better results.
INFO [2024-03-26 07:28:52] Removing 1724 low count (< 100 total reads) intervals.
INFO [2024-03-26 07:28:52] Removing 6 low coverage (< 15.0000X) intervals.
INFO [2024-03-26 07:28:52] Using 8050 intervals (8050 on-target, 0 off-target).
INFO [2024-03-26 07:28:52] No off-target intervals. If this is hybrid-capture data, consider adding them.
INFO [2024-03-26 07:28:52] No interval.file provided. Cannot check for any GC-biases.
INFO [2024-03-26 07:28:52] Loading VCF...
INFO [2024-03-26 07:28:52] Found 127 variants in VCF file.
WARN [2024-03-26 07:28:52] DP FORMAT field contains NAs. Removing 3 variants.
WARN [2024-03-26 07:28:52] DB INFO flag contains NAs
INFO [2024-03-26 07:28:52] 119 (96.0%) variants annotated as likely germline (DB INFO flag).
INFO [2024-03-26 07:28:53] LIB-02240e4 is tumor in VCF file.
INFO [2024-03-26 07:28:53] No homozygous variants in VCF, provide unfiltered VCF.
INFO [2024-03-26 07:28:53] Removing 0 low quality variants with non-offset BQ < 25.
INFO [2024-03-26 07:28:53] Base quality scores range from 29 to 35 (offset by 1)
INFO [2024-03-26 07:28:53] Minimum number of supporting reads ranges from 2 to 5, depending on coverage and BQS.
INFO [2024-03-26 07:28:54] Removing 1 variants with AF < 0.030 or AF >= 0.970 or insufficient supporting reads or depth < 15.
INFO [2024-03-26 07:28:54] Total size of targeted genomic region: 1.53Mb (2.33Mb with 50bp padding).
INFO [2024-03-26 07:28:55] 1.2% of targets contain variants.
INFO [2024-03-26 07:28:55] Setting somatic prior probabilities for dbSNP hits to 0.000500 or to 0.500000 otherwise.
INFO [2024-03-26 07:28:55] Excluding 5 novel or poor quality variants from segmentation.
INFO [2024-03-26 07:28:55] Sample sex: ?
INFO [2024-03-26 07:28:55] Segmenting data...
INFO [2024-03-26 07:28:55] Loading pre-computed boundaries for DNAcopy...
INFO [2024-03-26 07:28:55] Setting undo.SD parameter to 1.000000.
Setting multi-figure configuration
INFO [2024-03-26 07:28:57] Setting prune.hclust.h parameter to 0.200000.
INFO [2024-03-26 07:28:58] Found 59 segments with median size of 17.67Mb.
INFO [2024-03-26 07:28:58] Using 123 variants.
INFO [2024-03-26 07:28:58] Mean standard deviation of log-ratios: 0.34 (MAPD: 0.30)
INFO [2024-03-26 07:28:58] 2D-grid search of purity and ploidy...
INFO [2024-03-26 07:29:00] Local optima: 0.62/1.9, 0.3/1.9, 0.32/2
INFO [2024-03-26 07:29:00] Testing local optimum 1/3 at purity 0.62 and total ploidy 1.90...
INFO [2024-03-26 07:29:02] Testing local optimum 2/3 at purity 0.30 and total ploidy 1.90...
INFO [2024-03-26 07:29:03] Recalibrating log-ratios...
INFO [2024-03-26 07:29:03] Testing local optimum 2/3 at purity 0.30 and total ploidy 1.90...
INFO [2024-03-26 07:29:04] Recalibrating log-ratios...
INFO [2024-03-26 07:29:04] Testing local optimum 2/3 at purity 0.30 and total ploidy 1.90...
INFO [2024-03-26 07:29:05] Recalibrating log-ratios...
INFO [2024-03-26 07:29:05] Testing local optimum 2/3 at purity 0.30 and total ploidy 1.90...
INFO [2024-03-26 07:29:05] Testing local optimum 3/3 at purity 0.32 and total ploidy 2.00...
INFO [2024-03-26 07:29:06] Recalibrating log-ratios...
INFO [2024-03-26 07:29:06] Testing local optimum 3/3 at purity 0.32 and total ploidy 2.00...
INFO [2024-03-26 07:29:07] Recalibrating log-ratios...
INFO [2024-03-26 07:29:07] Testing local optimum 3/3 at purity 0.32 and total ploidy 2.00...
INFO [2024-03-26 07:29:08] Recalibrating log-ratios...
INFO [2024-03-26 07:29:08] Testing local optimum 3/3 at purity 0.32 and total ploidy 2.00...
INFO [2024-03-26 07:29:09] Skipping 2 solutions exceeding max.non.clonal (0.20): 0.3/1.9, 0.35/1.61 (purity/tumor ploidy)
INFO [2024-03-26 07:29:09] Fitting variants with beta model for local optimum 1/3...
INFO [2024-03-26 07:29:09] Fitting variants for purity 0.65, tumor ploidy 1.75 and contamination 0.01.
INFO [2024-03-26 07:29:16] Optimized purity: 0.65
INFO [2024-03-26 07:29:16] Done.
INFO [2024-03-26 07:29:16] ------------------------------------------------------------
INFO [2024-03-26 07:29:17] ------------------------------------------------------------
INFO [2024-03-26 07:29:17] PureCN 2.8.1
INFO [2024-03-26 07:29:17] ------------------------------------------------------------
INFO [2024-03-26 07:29:17] Loading coverage files...
INFO [2024-03-26 07:29:18] Mean target coverages: 112X (tumor) 99X (normal).
WARN [2024-03-26 07:29:18] Allosome coverage missing, cannot determine sex.
WARN [2024-03-26 07:29:18] Allosome coverage missing, cannot determine sex.
INFO [2024-03-26 07:29:18] No Gene column in interval.file. You won't get gene-level calls.
INFO [2024-03-26 07:29:20] Removing 228 intervals with missing log.ratio.
INFO [2024-03-26 07:29:20] Removing 15 low/high GC targets.
INFO [2024-03-26 07:29:20] Removing 21 small (< 5bp) intervals.
INFO [2024-03-26 07:29:20] Removing 19 intervals with low total coverage in normal (< 150.00 reads).
WARN [2024-03-26 07:29:20] No normalDB provided. Provide one for better results.
INFO [2024-03-26 07:29:20] Removing 1718 low count (< 100 total reads) intervals.
INFO [2024-03-26 07:29:20] Removing 6 low coverage (< 15.0000X) intervals.
INFO [2024-03-26 07:29:20] Removing 36 low mappability intervals.
INFO [2024-03-26 07:29:20] Using 8006 intervals (8006 on-target, 0 off-target).
INFO [2024-03-26 07:29:20] No off-target intervals. If this is hybrid-capture data, consider adding them.
INFO [2024-03-26 07:29:20] AT/GC dropout: 1.01 (tumor), 1.03 (normal), 0.98 (coverage log-ratio). 
INFO [2024-03-26 07:29:20] Loading VCF...
INFO [2024-03-26 07:29:20] Found 127 variants in VCF file.
INFO [2024-03-26 07:29:21] 127 (100.0%) variants annotated as likely germline (DB INFO flag).
INFO [2024-03-26 07:29:21] LIB-02240e4 is tumor in VCF file.
INFO [2024-03-26 07:29:21] No homozygous variants in VCF, provide unfiltered VCF.
INFO [2024-03-26 07:29:21] Removing 1 non heterozygous (in matched normal) germline SNPs.
INFO [2024-03-26 07:29:22] Removing 0 low quality variants with non-offset BQ < 25.
INFO [2024-03-26 07:29:22] Base quality scores range from 29 to 35 (offset by 1)
INFO [2024-03-26 07:29:22] Minimum number of supporting reads ranges from 2 to 5, depending on coverage and BQS.
INFO [2024-03-26 07:29:23] Removing 1 variants with AF < 0.030 or AF >= Inf or insufficient supporting reads or depth < 15.
INFO [2024-03-26 07:29:23] Total size of targeted genomic region: 1.51Mb (2.31Mb with 50bp padding).
INFO [2024-03-26 07:29:23] 1.2% of targets contain variants.
INFO [2024-03-26 07:29:23] Removing 11 variants outside intervals.
INFO [2024-03-26 07:29:23] Found SOMATIC annotation in VCF.
INFO [2024-03-26 07:29:23] Setting somatic prior probabilities for somatic variants to 0.999000 or to 0.000100 otherwise.
INFO [2024-03-26 07:29:24] Found SOMATIC annotation in VCF. Setting mapping bias to 0.970.
INFO [2024-03-26 07:29:24] Excluding 0 novel or poor quality variants from segmentation.
INFO [2024-03-26 07:29:24] Sample sex: ?
INFO [2024-03-26 07:29:24] Segmenting data...
INFO [2024-03-26 07:29:24] MAPD of 114 allelic fractions: 0.14 (0.15 adjusted).
INFO [2024-03-26 07:29:25] Using unweighted PSCBS.
INFO [2024-03-26 07:29:25] Setting undo.SD parameter to 1.000000.
INFO [2024-03-26 07:30:01] Setting prune.hclust.h parameter to 0.200000.
INFO [2024-03-26 07:30:01] Found 71 segments with median size of 27.72Mb.
INFO [2024-03-26 07:30:01] Using 114 variants.
INFO [2024-03-26 07:30:01] Mean standard deviation of log-ratios: 0.34 (MAPD: 0.30)
INFO [2024-03-26 07:30:01] 2D-grid search of purity and ploidy...
INFO [2024-03-26 07:30:01] Local optima: 0.65/1.8, 0.38/2.2
INFO [2024-03-26 07:30:01] Testing local optimum 1/2 at purity 0.65 and total ploidy 1.80...
INFO [2024-03-26 07:30:02] Recalibrating log-ratios...
INFO [2024-03-26 07:30:02] Testing local optimum 1/2 at purity 0.65 and total ploidy 1.80...
INFO [2024-03-26 07:30:03] Recalibrating log-ratios...
INFO [2024-03-26 07:30:03] Testing local optimum 1/2 at purity 0.65 and total ploidy 1.80...
INFO [2024-03-26 07:30:04] Recalibrating log-ratios...
INFO [2024-03-26 07:30:04] Testing local optimum 1/2 at purity 0.65 and total ploidy 1.80...
INFO [2024-03-26 07:30:05] Testing local optimum 2/2 at purity 0.38 and total ploidy 2.20...
INFO [2024-03-26 07:30:08] Skipping 1 solutions exceeding max.non.clonal (0.20): 0.65/2.04 (purity/tumor ploidy)
INFO [2024-03-26 07:30:08] Fitting variants with beta model for local optimum 2/2...
INFO [2024-03-26 07:30:08] Fitting variants for purity 0.40, tumor ploidy 2.53 and contamination 0.01.
INFO [2024-03-26 07:30:19] Optimized purity: 0.40
INFO [2024-03-26 07:30:19] Done.
INFO [2024-03-26 07:30:19] ------------------------------------------------------------
FATAL [2024-03-26 07:30:19] This function requires gene-level calls. Please add a column 'Gene' 

FATAL [2024-03-26 07:30:19] containing gene symbols to the interval.file. 

FATAL [2024-03-26 07:30:19]  

FATAL [2024-03-26 07:30:19] This is most likely a user error due to invalid input data or 

FATAL [2024-03-26 07:30:19] parameters (PureCN 2.8.1). 

INFO [2024-03-26 07:30:23] ------------------------------------------------------------
INFO [2024-03-26 07:30:23] PureCN 2.8.1
INFO [2024-03-26 07:30:23] ------------------------------------------------------------
INFO [2024-03-26 07:30:23] Loading coverage files...
INFO [2024-03-26 07:30:24] Mean target coverages: 112X (tumor) 99X (normal).
WARN [2024-03-26 07:30:24] Allosome coverage missing, cannot determine sex.
WARN [2024-03-26 07:30:24] Allosome coverage missing, cannot determine sex.
INFO [2024-03-26 07:30:24] Removing 228 intervals with missing log.ratio.
INFO [2024-03-26 07:30:24] Removing 22 small (< 5bp) intervals.
INFO [2024-03-26 07:30:24] Removing 19 intervals with low total coverage in normal (< 150.00 reads).
WARN [2024-03-26 07:30:24] No normalDB provided. Provide one for better results.
INFO [2024-03-26 07:30:24] Removing 1724 low count (< 100 total reads) intervals.
INFO [2024-03-26 07:30:24] Removing 6 low coverage (< 15.0000X) intervals.
INFO [2024-03-26 07:30:24] Using 8050 intervals (8050 on-target, 0 off-target).
INFO [2024-03-26 07:30:24] No off-target intervals. If this is hybrid-capture data, consider adding them.
INFO [2024-03-26 07:30:24] No interval.file provided. Cannot check for any GC-biases.
INFO [2024-03-26 07:30:24] Loading VCF...
INFO [2024-03-26 07:30:25] Found 127 variants in VCF file.
INFO [2024-03-26 07:30:25] 127 (100.0%) variants annotated as likely germline (DB INFO flag).
FATAL [2024-03-26 07:30:25] Different chromosome names in coverage and VCF. 

FATAL [2024-03-26 07:30:25]  

FATAL [2024-03-26 07:30:25] This is most likely a user error due to invalid input data or 

FATAL [2024-03-26 07:30:25] parameters (PureCN 2.8.1). 

INFO [2024-03-26 07:30:25] ------------------------------------------------------------
INFO [2024-03-26 07:30:25] PureCN 2.8.1
INFO [2024-03-26 07:30:25] ------------------------------------------------------------
INFO [2024-03-26 07:30:25] Loading coverage files...
INFO [2024-03-26 07:30:25] Mean target coverages: 112X (tumor) 99X (normal).
WARN [2024-03-26 07:30:25] Allosome coverage missing, cannot determine sex.
WARN [2024-03-26 07:30:25] Allosome coverage missing, cannot determine sex.
FATAL [2024-03-26 07:30:26] tumor.coverage.file and interval.file do not align. 

FATAL [2024-03-26 07:30:26]  

FATAL [2024-03-26 07:30:26] This is most likely a user error due to invalid input data or 

FATAL [2024-03-26 07:30:26] parameters (PureCN 2.8.1). 

INFO [2024-03-26 07:30:27] ------------------------------------------------------------
INFO [2024-03-26 07:30:27] PureCN 2.8.1
INFO [2024-03-26 07:30:27] ------------------------------------------------------------
INFO [2024-03-26 07:30:27] Loading coverage files...
INFO [2024-03-26 07:30:27] Mean target coverages: 112X (tumor) 99X (normal).
WARN [2024-03-26 07:30:27] Allosome coverage missing, cannot determine sex.
WARN [2024-03-26 07:30:27] Allosome coverage missing, cannot determine sex.
INFO [2024-03-26 07:30:30] Removing 228 intervals with missing log.ratio.
INFO [2024-03-26 07:30:30] Removing 15 low/high GC targets.
INFO [2024-03-26 07:30:30] Removing 21 small (< 5bp) intervals.
INFO [2024-03-26 07:30:30] Removing 19 intervals with low total coverage in normal (< 150.00 reads).
WARN [2024-03-26 07:30:30] No normalDB provided. Provide one for better results.
INFO [2024-03-26 07:30:30] Removing 1718 low count (< 100 total reads) intervals.
INFO [2024-03-26 07:30:30] Removing 6 low coverage (< 15.0000X) intervals.
INFO [2024-03-26 07:30:30] Removing 36 low mappability intervals.
INFO [2024-03-26 07:30:30] Using 8006 intervals (8006 on-target, 0 off-target).
INFO [2024-03-26 07:30:30] No off-target intervals. If this is hybrid-capture data, consider adding them.
INFO [2024-03-26 07:30:30] AT/GC dropout: 1.01 (tumor), 1.03 (normal), 0.98 (coverage log-ratio). 
INFO [2024-03-26 07:30:30] Loading VCF...
INFO [2024-03-26 07:30:30] Found 127 variants in VCF file.
INFO [2024-03-26 07:30:31] 127 (100.0%) variants annotated as likely germline (DB INFO flag).
INFO [2024-03-26 07:30:31] LIB-02240e4 is tumor in VCF file.
INFO [2024-03-26 07:30:31] No homozygous variants in VCF, provide unfiltered VCF.
INFO [2024-03-26 07:30:31] Removing 1 non heterozygous (in matched normal) germline SNPs.
INFO [2024-03-26 07:30:31] Removing 0 low quality variants with non-offset BQ < 25.
INFO [2024-03-26 07:30:31] Base quality scores range from 29 to 35 (offset by 1)
INFO [2024-03-26 07:30:31] Minimum number of supporting reads ranges from 2 to 5, depending on coverage and BQS.
INFO [2024-03-26 07:30:32] Removing 1 variants with AF < 0.030 or AF >= 1.000 or insufficient supporting reads or depth < 15.
INFO [2024-03-26 07:30:32] Total size of targeted genomic region: 1.51Mb (2.31Mb with 50bp padding).
INFO [2024-03-26 07:30:33] 1.2% of targets contain variants.
INFO [2024-03-26 07:30:33] Removing 11 variants outside intervals.
INFO [2024-03-26 07:30:33] Found SOMATIC annotation in VCF.
INFO [2024-03-26 07:30:33] Setting somatic prior probabilities for somatic variants to 0.999000 or to 0.000100 otherwise.
INFO [2024-03-26 07:30:33] Found SOMATIC annotation in VCF. Setting mapping bias to 0.970.
INFO [2024-03-26 07:30:34] Excluding 0 novel or poor quality variants from segmentation.
INFO [2024-03-26 07:30:34] Sample sex: ?
INFO [2024-03-26 07:30:34] Segmenting data...
INFO [2024-03-26 07:30:34] Loading pre-computed boundaries for DNAcopy...
INFO [2024-03-26 07:30:34] Setting undo.SD parameter to 1.000000.
INFO [2024-03-26 07:30:36] Setting prune.hclust.h parameter to 0.200000.
INFO [2024-03-26 07:30:36] Found 59 segments with median size of 17.67Mb.
INFO [2024-03-26 07:30:36] Using 114 variants.
INFO [2024-03-26 07:30:36] Mean standard deviation of log-ratios: 0.34 (MAPD: 0.30)
INFO [2024-03-26 07:30:36] 2D-grid search of purity and ploidy...
INFO [2024-03-26 07:30:38] Local optima: 0.63/1.9
INFO [2024-03-26 07:30:38] Testing local optimum 1/1 at purity 0.63 and total ploidy 1.90...
INFO [2024-03-26 07:30:40] Fitting variants with beta model for local optimum 1/1...
INFO [2024-03-26 07:30:40] Fitting variants for purity 0.65, tumor ploidy 1.75 and contamination 0.01.
INFO [2024-03-26 07:30:48] Optimized purity: 0.65
INFO [2024-03-26 07:30:48] Done.
INFO [2024-03-26 07:30:48] ------------------------------------------------------------
INFO [2024-03-26 07:30:49] ------------------------------------------------------------
INFO [2024-03-26 07:30:49] PureCN 2.8.1
INFO [2024-03-26 07:30:49] ------------------------------------------------------------
INFO [2024-03-26 07:30:49] Loading coverage files...
INFO [2024-03-26 07:30:49] Mean target coverages: 112X (tumor) 99X (normal).
WARN [2024-03-26 07:30:49] Allosome coverage missing, cannot determine sex.
WARN [2024-03-26 07:30:49] Allosome coverage missing, cannot determine sex.
INFO [2024-03-26 07:30:52] Removing 228 intervals with missing log.ratio.
INFO [2024-03-26 07:30:52] Removing 15 low/high GC targets.
INFO [2024-03-26 07:30:52] Removing 21 small (< 5bp) intervals.
INFO [2024-03-26 07:30:52] Removing 19 intervals with low total coverage in normal (< 150.00 reads).
WARN [2024-03-26 07:30:52] No normalDB provided. Provide one for better results.
INFO [2024-03-26 07:30:52] Removing 1718 low count (< 100 total reads) intervals.
INFO [2024-03-26 07:30:52] Removing 6 low coverage (< 15.0000X) intervals.
INFO [2024-03-26 07:30:52] Removing 36 low mappability intervals.
INFO [2024-03-26 07:30:52] Using 8006 intervals (8006 on-target, 0 off-target).
INFO [2024-03-26 07:30:52] No off-target intervals. If this is hybrid-capture data, consider adding them.
INFO [2024-03-26 07:30:52] AT/GC dropout: 1.01 (tumor), 1.03 (normal), 0.98 (coverage log-ratio). 
INFO [2024-03-26 07:30:52] Loading VCF...
INFO [2024-03-26 07:30:52] Found 127 variants in VCF file.
INFO [2024-03-26 07:30:53] 127 (100.0%) variants annotated as likely germline (DB INFO flag).
INFO [2024-03-26 07:30:53] LIB-02240e4 is tumor in VCF file.
INFO [2024-03-26 07:30:53] No homozygous variants in VCF, provide unfiltered VCF.
INFO [2024-03-26 07:30:53] Removing 1 non heterozygous (in matched normal) germline SNPs.
INFO [2024-03-26 07:30:53] Removing 0 low quality variants with non-offset BQ < 25.
INFO [2024-03-26 07:30:53] Base quality scores range from 29 to 35 (offset by 1)
INFO [2024-03-26 07:30:53] Minimum number of supporting reads ranges from 2 to 5, depending on coverage and BQS.
INFO [2024-03-26 07:30:55] Removing 1 variants with AF < 0.030 or AF >= 1.000 or insufficient supporting reads or depth < 15.
INFO [2024-03-26 07:30:55] Total size of targeted genomic region: 1.51Mb (2.31Mb with 50bp padding).
INFO [2024-03-26 07:30:56] 1.2% of targets contain variants.
INFO [2024-03-26 07:30:56] Removing 11 variants outside intervals.
INFO [2024-03-26 07:30:56] Found SOMATIC annotation in VCF.
INFO [2024-03-26 07:30:56] Setting somatic prior probabilities for somatic variants to 0.999000 or to 0.000100 otherwise.
INFO [2024-03-26 07:30:56] Found SOMATIC annotation in VCF. Setting mapping bias to 0.970.
INFO [2024-03-26 07:30:56] Excluding 0 novel or poor quality variants from segmentation.
INFO [2024-03-26 07:30:57] Sample sex: ?
INFO [2024-03-26 07:30:57] Segmenting data...
INFO [2024-03-26 07:30:57] Loading pre-computed boundaries for DNAcopy...
INFO [2024-03-26 07:30:57] Setting undo.SD parameter to 1.000000.
INFO [2024-03-26 07:30:58] Setting prune.hclust.h parameter to 0.200000.
INFO [2024-03-26 07:30:59] Found 59 segments with median size of 17.67Mb.
INFO [2024-03-26 07:30:59] Using 114 variants.
INFO [2024-03-26 07:30:59] Mean standard deviation of log-ratios: 0.34 (MAPD: 0.30)
INFO [2024-03-26 07:30:59] 2D-grid search of purity and ploidy...
INFO [2024-03-26 07:31:00] Local optima: 0.63/1.9
INFO [2024-03-26 07:31:00] Testing local optimum 1/1 at purity 0.63 and total ploidy 1.90...
INFO [2024-03-26 07:31:02] Fitting variants with beta model for local optimum 1/1...
WARN [2024-03-26 07:31:03] Some gene symbols found on multiple chromosomes. Use of approved symbols suggested.
INFO [2024-03-26 07:31:03] Fitting variants for purity 0.65, tumor ploidy 1.75 and contamination 0.01.
INFO [2024-03-26 07:31:11] Optimized purity: 0.65
INFO [2024-03-26 07:31:11] Done.
INFO [2024-03-26 07:31:11] ------------------------------------------------------------
FATAL [2024-03-26 07:31:12] chr1 not valid chromosome name(s). Valid names are: 

FATAL [2024-03-26 07:31:12] 1,2,3,4,5,6,7,8,9,10,11,12,13,14,15,16,17,18,19,20,21,22 

FATAL [2024-03-26 07:31:12]  

FATAL [2024-03-26 07:31:12] This is most likely a user error due to invalid input data or 

FATAL [2024-03-26 07:31:12] parameters (PureCN 2.8.1). 

INFO [2024-03-26 07:31:12] ------------------------------------------------------------
INFO [2024-03-26 07:31:12] PureCN 2.8.1
INFO [2024-03-26 07:31:12] ------------------------------------------------------------
INFO [2024-03-26 07:31:12] Loading coverage files...
WARN [2024-03-26 07:31:13] Provided sampleid (Sample2) does not match Sample1 found in segmentation.
WARN [2024-03-26 07:31:13] Allosome coverage missing, cannot determine sex.
WARN [2024-03-26 07:31:13] Allosome coverage missing, cannot determine sex.
INFO [2024-03-26 07:31:13] Removing 10 intervals with missing log.ratio.
INFO [2024-03-26 07:31:13] Using 10039 intervals (10039 on-target, 0 off-target).
INFO [2024-03-26 07:31:13] No off-target intervals. If this is hybrid-capture data, consider adding them.
INFO [2024-03-26 07:31:13] Loading VCF...
INFO [2024-03-26 07:31:14] Found 127 variants in VCF file.
INFO [2024-03-26 07:31:14] 127 (100.0%) variants annotated as likely germline (DB INFO flag).
INFO [2024-03-26 07:31:14] LIB-02240e4 is tumor in VCF file.
INFO [2024-03-26 07:31:14] No homozygous variants in VCF, provide unfiltered VCF.
INFO [2024-03-26 07:31:15] Removing 1 non heterozygous (in matched normal) germline SNPs.
INFO [2024-03-26 07:31:15] Removing 0 low quality variants with non-offset BQ < 25.
INFO [2024-03-26 07:31:15] Base quality scores range from 29 to 35 (offset by 1)
INFO [2024-03-26 07:31:15] Minimum number of supporting reads ranges from 2 to 5, depending on coverage and BQS.
INFO [2024-03-26 07:31:16] Removing 1 variants with AF < 0.030 or AF >= 1.000 or insufficient supporting reads or depth < 15.
INFO [2024-03-26 07:31:16] Total size of targeted genomic region: 1.71Mb (2.70Mb with 50bp padding).
INFO [2024-03-26 07:31:16] 1.0% of targets contain variants.
INFO [2024-03-26 07:31:16] Removing 0 variants outside intervals.
INFO [2024-03-26 07:31:17] Found SOMATIC annotation in VCF.
INFO [2024-03-26 07:31:17] Setting somatic prior probabilities for somatic variants to 0.999000 or to 0.000100 otherwise.
INFO [2024-03-26 07:31:17] Found SOMATIC annotation in VCF. Setting mapping bias to 0.976.
INFO [2024-03-26 07:31:17] Excluding 0 novel or poor quality variants from segmentation.
INFO [2024-03-26 07:31:17] Sample sex: ?
INFO [2024-03-26 07:31:17] Segmenting data...
INFO [2024-03-26 07:31:17] Loaded provided segmentation file example_seg.txt (format DNAcopy).
WARN [2024-03-26 07:31:17] Provided sampleid (Sample2) does not match Sample1 found in segmentation.
INFO [2024-03-26 07:31:17] Re-centering provided segment means (offset -0.0033).
INFO [2024-03-26 07:31:17] Loading pre-computed boundaries for DNAcopy...
INFO [2024-03-26 07:31:17] Setting undo.SD parameter to 0.000000.
Setting multi-figure configuration
INFO [2024-03-26 07:31:18] Setting prune.hclust.h parameter to 0.200000.
INFO [2024-03-26 07:31:18] Found 54 segments with median size of 24.88Mb.
INFO [2024-03-26 07:31:18] Using 125 variants.
INFO [2024-03-26 07:31:19] Mean standard deviation of log-ratios: 0.40 (MAPD: 0.00)
INFO [2024-03-26 07:31:19] 2D-grid search of purity and ploidy...
INFO [2024-03-26 07:31:21] Local optima: 0.6/1.9
INFO [2024-03-26 07:31:21] Testing local optimum 1/1 at purity 0.60 and total ploidy 1.90...
INFO [2024-03-26 07:31:23] Fitting variants with beta model for local optimum 1/1...
INFO [2024-03-26 07:31:23] Fitting variants for purity 0.65, tumor ploidy 1.74 and contamination 0.01.
INFO [2024-03-26 07:31:32] Optimized purity: 0.65
INFO [2024-03-26 07:31:32] Done.
INFO [2024-03-26 07:31:32] ------------------------------------------------------------
WARN [2024-03-26 07:31:33] Allosome coverage missing, cannot determine sex.
WARN [2024-03-26 07:31:33] Allosome coverage missing, cannot determine sex.
Setting multi-figure configuration
FATAL [2024-03-26 07:31:53] Segmentation file expected with colnames ID, chrom, loc.start, loc.end, 

FATAL [2024-03-26 07:31:53] num.mark, seg.mean 

FATAL [2024-03-26 07:31:53]  

FATAL [2024-03-26 07:31:53] This is most likely a user error due to invalid input data or 

FATAL [2024-03-26 07:31:53] parameters (PureCN 2.8.1). 

FATAL [2024-03-26 07:31:54] seg.file contains multiple samples and sampleid missing. 

FATAL [2024-03-26 07:31:54]  

FATAL [2024-03-26 07:31:54] This is most likely a user error due to invalid input data or 

FATAL [2024-03-26 07:31:54] parameters (PureCN 2.8.1). 

FATAL [2024-03-26 07:31:55] seg.file contains multiple samples and sampleid does not match any. 

FATAL [2024-03-26 07:31:55]  

FATAL [2024-03-26 07:31:55] This is most likely a user error due to invalid input data or 

FATAL [2024-03-26 07:31:55] parameters (PureCN 2.8.1). 

WARN [2024-03-26 07:31:55] Allosome coverage missing, cannot determine sex.
WARN [2024-03-26 07:31:55] Allosome coverage missing, cannot determine sex.
WARN [2024-03-26 07:32:17] Allosome coverage missing, cannot determine sex.
WARN [2024-03-26 07:32:17] Allosome coverage missing, cannot determine sex.
WARN [2024-03-26 07:32:17] Provided sampleid (example_tumor.txt.gz) does not match Sample1 found in segmentation.
WARN [2024-03-26 07:32:24] Allosome coverage missing, cannot determine sex.
WARN [2024-03-26 07:32:24] Allosome coverage missing, cannot determine sex.
WARN [2024-03-26 07:32:25] No normalDB provided. Provide one for better results.
WARN [2024-03-26 07:32:43] Allosome coverage missing, cannot determine sex.
WARN [2024-03-26 07:32:43] Allosome coverage missing, cannot determine sex.
WARN [2024-03-26 07:32:50] Provided sampleid (Sample.1) does not match Sample1 found in segmentation.
Setting multi-figure configuration
WARN [2024-03-26 07:33:09] Allosome coverage missing, cannot determine sex.
WARN [2024-03-26 07:33:09] Allosome coverage missing, cannot determine sex.
WARN [2024-03-26 07:33:09] No normalDB provided. Provide one for better results.
WARN [2024-03-26 07:33:10] Sampleid looks like a normal in VCF, not like a tumor.
WARN [2024-03-26 07:33:26] Allosome coverage missing, cannot determine sex.
WARN [2024-03-26 07:33:26] Allosome coverage missing, cannot determine sex.
FATAL [2024-03-26 07:33:26] normalDB not a valid normalDB object. Use createNormalDatabase to 

FATAL [2024-03-26 07:33:26] create one. 

FATAL [2024-03-26 07:33:26]  

FATAL [2024-03-26 07:33:26] This is most likely a user error due to invalid input data or 

FATAL [2024-03-26 07:33:26] parameters (PureCN 2.8.1). 

WARN [2024-03-26 07:33:27] You are likely not using the correct baits file!
WARN [2024-03-26 07:33:27] Allosome coverage missing, cannot determine sex.
WARN [2024-03-26 07:33:27] Allosome coverage missing, cannot determine sex.
WARN [2024-03-26 07:33:35] Allosome coverage missing, cannot determine sex.
WARN [2024-03-26 07:33:35] Allosome coverage missing, cannot determine sex.
FATAL [2024-03-26 07:33:35] normalDB appears to be empty. 

FATAL [2024-03-26 07:33:35]  

FATAL [2024-03-26 07:33:35] This is most likely a user error due to invalid input data or 

FATAL [2024-03-26 07:33:35] parameters (PureCN 2.8.1). 

WARN [2024-03-26 07:33:37] Allosome coverage missing, cannot determine sex.
WARN [2024-03-26 07:33:37] Allosome coverage missing, cannot determine sex.
WARN [2024-03-26 07:33:38] Intervals in coverage and interval.file have conflicting on/off-target annotation.
WARN [2024-03-26 07:34:07] Allosome coverage missing, cannot determine sex.
WARN [2024-03-26 07:34:07] Allosome coverage missing, cannot determine sex.
FATAL [2024-03-26 07:34:07] Seqlevels missing in provided segmentation: 6 

FATAL [2024-03-26 07:34:07]  

FATAL [2024-03-26 07:34:07] This is most likely a user error due to invalid input data or 

FATAL [2024-03-26 07:34:07] parameters (PureCN 2.8.1). 

sh: gatk: command not found
WARN [2024-03-26 07:34:11] Cannot find gatk binary in path.
WARN [2024-03-26 07:34:12] Allosome coverage missing, cannot determine sex.
WARN [2024-03-26 07:34:12] Allosome coverage missing, cannot determine sex.
WARN [2024-03-26 07:34:12] No normalDB provided. Provide one for better results.
FATAL [2024-03-26 07:34:12] segmentationHclust requires an input segmentation. 

FATAL [2024-03-26 07:34:12]  

FATAL [2024-03-26 07:34:12] This is most likely a user error due to invalid input data or 

FATAL [2024-03-26 07:34:12] parameters (PureCN 2.8.1). 

WARN [2024-03-26 07:34:13] Allosome coverage missing, cannot determine sex.
WARN [2024-03-26 07:34:13] Allosome coverage missing, cannot determine sex.
WARN [2024-03-26 07:34:13] No normalDB provided. Provide one for better results.
FATAL [2024-03-26 07:35:28] The normal.panel.vcf.file contains only a single sample. 

FATAL [2024-03-26 07:35:28]  

FATAL [2024-03-26 07:35:28] This is most likely a user error due to invalid input data or 

FATAL [2024-03-26 07:35:28] parameters (PureCN 2.8.1). 

FATAL [2024-03-26 07:35:28] mapping.bias.file must be a file with *.rds suffix. 

FATAL [2024-03-26 07:35:28]  

FATAL [2024-03-26 07:35:28] This is most likely a user error due to invalid input data or 

FATAL [2024-03-26 07:35:28] parameters (PureCN 2.8.1). 

Failed with error:  'there is no package called 'genomicsdb''
FATAL [2024-03-26 07:35:31] min.normals (0) must be >= 1. 

FATAL [2024-03-26 07:35:31]  

FATAL [2024-03-26 07:35:31] This is most likely a user error due to invalid input data or 

FATAL [2024-03-26 07:35:31] parameters (PureCN 2.8.1). 

FATAL [2024-03-26 07:35:32] min.normals (10) cannot be larger than min.normals.assign.betafit (3). 

FATAL [2024-03-26 07:35:32]  

FATAL [2024-03-26 07:35:32] This is most likely a user error due to invalid input data or 

FATAL [2024-03-26 07:35:32] parameters (PureCN 2.8.1). 

FATAL [2024-03-26 07:35:32] min.normals.assign.betafit (10) cannot be larger than 

FATAL [2024-03-26 07:35:32] min.normals.betafit (7). 

FATAL [2024-03-26 07:35:32]  

FATAL [2024-03-26 07:35:32] This is most likely a user error due to invalid input data or 

FATAL [2024-03-26 07:35:32] parameters (PureCN 2.8.1). 

FATAL [2024-03-26 07:35:33] min.normals.betafit (20) cannot be larger than 

FATAL [2024-03-26 07:35:33] min.normals.position.specific.fit (10). 

FATAL [2024-03-26 07:35:33]  

FATAL [2024-03-26 07:35:33] This is most likely a user error due to invalid input data or 

FATAL [2024-03-26 07:35:33] parameters (PureCN 2.8.1). 

FATAL [2024-03-26 07:35:33] min.betafit.rho not within expected range or format. 

FATAL [2024-03-26 07:35:33]  

FATAL [2024-03-26 07:35:33] This is most likely a user error due to invalid input data or 

FATAL [2024-03-26 07:35:33] parameters (PureCN 2.8.1). 

FATAL [2024-03-26 07:35:34] max.betafit.rho not within expected range or format. 

FATAL [2024-03-26 07:35:34]  

FATAL [2024-03-26 07:35:34] This is most likely a user error due to invalid input data or 

FATAL [2024-03-26 07:35:34] parameters (PureCN 2.8.1). 

FATAL [2024-03-26 07:35:35] min.betafit.rho (0.9) cannot be larger than max.betafit.rho (0.2). 

FATAL [2024-03-26 07:35:35]  

FATAL [2024-03-26 07:35:35] This is most likely a user error due to invalid input data or 

FATAL [2024-03-26 07:35:35] parameters (PureCN 2.8.1). 

[ FAIL 0 | WARN 13 | SKIP 2 | PASS 402 ]

══ Skipped tests (2) ═══════════════════════════════════════════════════════════
• gatk binary > 4.1.7.0 required (1): 'test_segmentation.R:25:4'
• genomicsdb required (1): 'test_setMappingBiasVcf.R:45:5'

[ FAIL 0 | WARN 13 | SKIP 2 | PASS 402 ]
> 
> proc.time()
   user  system elapsed 
990.319  23.035 955.458 

Example timings

PureCN.Rcheck/PureCN-Ex.timings

nameusersystemelapsed
annotateTargets14.673 0.63515.362
bootstrapResults1.3430.1361.484
calculateBamCoverageByInterval0.5750.0310.609
calculateLogRatio1.4070.1091.524
calculateMappingBiasGatk40.0000.0010.000
calculateMappingBiasVcf4.6360.1674.840
calculatePowerDetectSomatic3.4970.0143.518
calculateTangentNormal7.8290.0848.016
callAlterations0.2060.0140.221
callAlterationsFromSegmentation3.7220.0383.781
callAmplificationsInLowPurity100.793 0.390101.712
callCIN0.3990.0100.410
callLOH0.3390.0110.351
callMutationBurden3.5910.0533.651
centromeres0.0020.0050.008
correctCoverageBias3.9640.0464.073
createCurationFile0.5110.0290.544
createNormalDatabase6.8300.0687.264
filterIntervals52.531 0.60655.007
filterVcfBasic1.9030.0222.183
filterVcfMuTect1.8600.0212.029
filterVcfMuTect21.9300.0252.121
findFocal30.141 0.18031.394
getSexFromCoverage0.5930.0190.681
getSexFromVcf0.7610.0100.817
plotAbs0.4320.0180.576
poolCoverage1.2800.0401.465
predictSomatic1.1510.0181.178
preprocessIntervals1.7670.0261.939
processMultipleSamples6.9620.0777.415
readAllelicCountsFile1.8540.0081.943
readCoverageFile0.5760.0190.623
readCurationFile0.3460.0110.375
readIntervalFile0.3590.0030.379
readLogRatioFile0.0710.0020.079
readSegmentationFile0.0150.0010.017
runAbsoluteCN48.072 0.77652.927
segmentationCBS23.282 0.18525.931
segmentationGATK40.0050.0020.008
segmentationHclust45.153 0.66450.356
segmentationPSCBS65.927 0.94774.123
setMappingBiasVcf0.8510.0120.960
setPriorVcf0.8100.0100.895