PureCN.Rcheck/tests_i386/testthat.Rout
R version 3.5.3 (2019-03-11) -- "Great Truth"
Copyright (C) 2019 The R Foundation for Statistical Computing
Platform: i386-w64-mingw32/i386 (32-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
Loading required package: parallel
Attaching package: 'BiocGenerics'
The following objects are masked from 'package:parallel':
clusterApply, clusterApplyLB, clusterCall, clusterEvalQ,
clusterExport, clusterMap, parApply, parCapply, parLapply,
parLapplyLB, parRapply, parSapply, parSapplyLB
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, append,
as.data.frame, basename, cbind, colMeans, colSums, colnames,
dirname, do.call, duplicated, eval, evalq, get, grep, grepl,
intersect, is.unsorted, lapply, lengths, mapply, match, mget,
order, paste, pmax, pmax.int, pmin, pmin.int, rank, rbind,
rowMeans, rowSums, rownames, sapply, setdiff, sort, table, tapply,
union, unique, unsplit, which, which.max, which.min
Loading required package: GenomeInfoDb
Loading required package: S4Vectors
Loading required package: stats4
Attaching package: 'S4Vectors'
The following object is masked from 'package:base':
expand.grid
Loading required package: IRanges
Attaching package: 'IRanges'
The following object is masked from 'package:grDevices':
windows
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")'.
Loading required package: DelayedArray
Loading required package: matrixStats
Attaching package: 'matrixStats'
The following objects are masked from 'package:Biobase':
anyMissing, rowMedians
Loading required package: BiocParallel
Attaching package: 'DelayedArray'
The following objects are masked from 'package:matrixStats':
colMaxs, colMins, colRanges, rowMaxs, rowMins, rowRanges
The following objects are masked from 'package:base':
aperm, apply
Loading required package: Rsamtools
Loading required package: Biostrings
Loading required package: XVector
Attaching package: 'Biostrings'
The following object is masked from 'package:DelayedArray':
type
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 [2019-04-13 04:54:21] Attempted adding gene symbols to intervals. Heuristics have been used to pick symbols for overlapping genes.
WARN [2019-04-13 04:54:21] Chromosome naming style of txdb file (UCSC) was different from interval file (NCBI).
WARN [2019-04-13 04:54:25] Attempted adding gene symbols to intervals. Heuristics have been used to pick symbols for overlapping genes.
FATAL [2019-04-13 04:54:27] tumor.coverage.file and interval.file do not align.
FATAL [2019-04-13 04:54:27]
FATAL [2019-04-13 04:54:27] This is most likely a user error due to invalid input data or
FATAL [2019-04-13 04:54:27] parameters (PureCN 1.12.2).
FATAL [2019-04-13 04:54:27] Need either f or purity and ploidy.
FATAL [2019-04-13 04:54:27]
FATAL [2019-04-13 04:54:27] This is most likely a user error due to invalid input data or
FATAL [2019-04-13 04:54:27] parameters (PureCN 1.12.2).
FATAL [2019-04-13 04:54:27] f not in expected range.
FATAL [2019-04-13 04:54:27]
FATAL [2019-04-13 04:54:27] This is most likely a user error due to invalid input data or
FATAL [2019-04-13 04:54:27] parameters (PureCN 1.12.2).
FATAL [2019-04-13 04:54:27] coverage not in expected range (>=2)
FATAL [2019-04-13 04:54:27]
FATAL [2019-04-13 04:54:27] This is most likely a user error due to invalid input data or
FATAL [2019-04-13 04:54:27] parameters (PureCN 1.12.2).
FATAL [2019-04-13 04:54:27] purity not in expected range.
FATAL [2019-04-13 04:54:27]
FATAL [2019-04-13 04:54:27] This is most likely a user error due to invalid input data or
FATAL [2019-04-13 04:54:27] parameters (PureCN 1.12.2).
FATAL [2019-04-13 04:54:27] ploidy not in expected range.
FATAL [2019-04-13 04:54:27]
FATAL [2019-04-13 04:54:27] This is most likely a user error due to invalid input data or
FATAL [2019-04-13 04:54:27] parameters (PureCN 1.12.2).
FATAL [2019-04-13 04:54:27] cell.fraction not in expected range.
FATAL [2019-04-13 04:54:27]
FATAL [2019-04-13 04:54:27] This is most likely a user error due to invalid input data or
FATAL [2019-04-13 04:54:27] parameters (PureCN 1.12.2).
INFO [2019-04-13 04:54:30] ------------------------------------------------------------
INFO [2019-04-13 04:54:30] PureCN 1.12.2
INFO [2019-04-13 04:54:30] ------------------------------------------------------------
INFO [2019-04-13 04:54:30] Loading coverage files...
INFO [2019-04-13 04:54:30] Mean target coverages: 112X (tumor) 99X (normal).
WARN [2019-04-13 04:54:30] Allosome coverage missing, cannot determine sex.
WARN [2019-04-13 04:54:30] Allosome coverage missing, cannot determine sex.
INFO [2019-04-13 04:54:30] Removing 228 intervals with missing log.ratio.
INFO [2019-04-13 04:54:30] Removing 22 small (< 5bp) intervals.
INFO [2019-04-13 04:54:30] Removing 19 intervals with low total coverage in normal (< 150.00 reads).
WARN [2019-04-13 04:54:30] No normalDB provided. Provide one for better results.
INFO [2019-04-13 04:54:30] Removing 233 low coverage (< 15.0000X) intervals.
INFO [2019-04-13 04:54:30] Using 9547 intervals (9547 on-target, 0 off-target).
INFO [2019-04-13 04:54:30] No off-target intervals. If this is hybrid-capture data, consider adding them.
INFO [2019-04-13 04:54:30] No interval.file provided. Cannot check if data was GC-normalized. Was it?
INFO [2019-04-13 04:54:30] Loading VCF...
INFO [2019-04-13 04:54:30] Found 127 variants in VCF file.
INFO [2019-04-13 04:54:30] 127 (100.0%) variants annotated as likely germline (DB INFO flag).
INFO [2019-04-13 04:54:31] LIB-02240e4 is tumor in VCF file.
INFO [2019-04-13 04:54:31] No homozygous variants in VCF, provide unfiltered VCF.
INFO [2019-04-13 04:54:31] Removing 1 non heterozygous (in matched normal) germline SNPs.
INFO [2019-04-13 04:54:31] Removing 1 variants with AF < 0.030 or AF >= 1.000 or less than 4 supporting reads or depth < 15.
INFO [2019-04-13 04:54:31] Removing 0 low quality variants with BQ < 25.
INFO [2019-04-13 04:54:31] Total size of targeted genomic region: 1.63Mb (2.58Mb with 50bp padding).
INFO [2019-04-13 04:54:32] 1.0% of targets contain variants.
INFO [2019-04-13 04:54:32] Removing 2 variants outside intervals.
INFO [2019-04-13 04:54:32] Found SOMATIC annotation in VCF.
INFO [2019-04-13 04:54:32] Setting somatic prior probabilities for somatic variants to 0.999000 or to 0.000100 otherwise.
INFO [2019-04-13 04:54:32] Found SOMATIC annotation in VCF. Setting mapping bias to 0.973.
INFO [2019-04-13 04:54:32] Sample sex: ?
INFO [2019-04-13 04:54:32] Segmenting data...
INFO [2019-04-13 04:54:32] Loading pre-computed boundaries for DNAcopy...
INFO [2019-04-13 04:54:32] Setting undo.SD parameter to 1.000000.
INFO [2019-04-13 04:54:33] Setting prune.hclust.h parameter to 0.200000.
INFO [2019-04-13 04:54:34] Found 54 segments with median size of 24.88Mb.
INFO [2019-04-13 04:54:34] Using 123 variants.
INFO [2019-04-13 04:54:34] Mean standard deviation of log-ratios: 0.37
INFO [2019-04-13 04:54:34] 2D-grid search of purity and ploidy...
INFO [2019-04-13 04:54:36] Local optima: 0.63/1.9, 0.4/2.2, 0.5/2
INFO [2019-04-13 04:54:36] Testing local optimum 1/3 at purity 0.63 and total ploidy 1.90...
INFO [2019-04-13 04:54:37] Fitting variants for purity 0.65, tumor ploidy 1.74 and contamination 0.01.
INFO [2019-04-13 04:54:38] Optimized purity: 0.65
INFO [2019-04-13 04:54:38] Testing local optimum 2/3 at purity 0.40 and total ploidy 2.20...
INFO [2019-04-13 04:54:39] Fitting variants for purity 0.40, tumor ploidy 2.34 and contamination 0.01.
INFO [2019-04-13 04:54:40] Optimized purity: 0.40
INFO [2019-04-13 04:54:40] Testing local optimum 3/3 at purity 0.50 and total ploidy 2.00...
INFO [2019-04-13 04:54:43] Fitting variants for purity 0.65, tumor ploidy 1.74 and contamination 0.01.
INFO [2019-04-13 04:54:44] Optimized purity: 0.65
INFO [2019-04-13 04:54:44] Done.
INFO [2019-04-13 04:54:44] ------------------------------------------------------------
FATAL [2019-04-13 04:54:45] exclude not a GRanges object.
FATAL [2019-04-13 04:54:45]
FATAL [2019-04-13 04:54:45] This is most likely a user error due to invalid input data or
FATAL [2019-04-13 04:54:45] parameters (PureCN 1.12.2).
FATAL [2019-04-13 04:54:45] callable not a GRanges object.
FATAL [2019-04-13 04:54:45]
FATAL [2019-04-13 04:54:45] This is most likely a user error due to invalid input data or
FATAL [2019-04-13 04:54:45] parameters (PureCN 1.12.2).
WARN [2019-04-13 04:54:57] tumor.coverage.file and interval.file do not align.
FATAL [2019-04-13 04:54:58] No gc_bias column in interval.file.
FATAL [2019-04-13 04:54:58]
FATAL [2019-04-13 04:54:58] This is most likely a user error due to invalid input data or
FATAL [2019-04-13 04:54:58] parameters (PureCN 1.12.2).
FATAL [2019-04-13 04:54:59] Provided coverage is zero, most likely due to a corrupt BAM file.
FATAL [2019-04-13 04:54:59]
FATAL [2019-04-13 04:54:59] This is most likely a user error due to invalid input data or
FATAL [2019-04-13 04:54:59] parameters (PureCN 1.12.2).
INFO [2019-04-13 04:55:05] Reading C:\Users\biocbuild\bbs-3.8-bioc\tmpdir\RtmpK8Yhi6\file39fc49c83a95.rds...
INFO [2019-04-13 04:55:05] Reading C:\Users\biocbuild\bbs-3.8-bioc\tmpdir\RtmpK8Yhi6\file39fc49c83a95.rds...
INFO [2019-04-13 04:55:05] Reading C:\Users\biocbuild\bbs-3.8-bioc\tmpdir\RtmpK8Yhi6\file39fc49c83a95.rds...
INFO [2019-04-13 04:55:05] Reading C:\Users\biocbuild\bbs-3.8-bioc\tmpdir\RtmpK8Yhi6\file39fc49c83a95.rds...
INFO [2019-04-13 04:55:06] Reading C:\Users\biocbuild\bbs-3.8-bioc\tmpdir\RtmpK8Yhi6\file39fc49c83a95.rds...
INFO [2019-04-13 04:55:06] Reading C:\Users\biocbuild\bbs-3.8-bioc\tmpdir\RtmpK8Yhi6\file39fc49c83a95.rds...
FATAL [2019-04-13 04:55:06] Purity or Ploidy not numeric or in expected range.
FATAL [2019-04-13 04:55:06]
FATAL [2019-04-13 04:55:06] This is most likely a user error due to invalid input data or
FATAL [2019-04-13 04:55:06] parameters (PureCN 1.12.2).
INFO [2019-04-13 04:55:06] Reading C:\Users\biocbuild\bbs-3.8-bioc\tmpdir\RtmpK8Yhi6\file39fc49c83a95.rds...
INFO [2019-04-13 04:55:06] Reading C:\Users\biocbuild\bbs-3.8-bioc\tmpdir\RtmpK8Yhi6\file39fc49c83a95.rds...
FATAL [2019-04-13 04:55:06] 'Failed' column in
FATAL [2019-04-13 04:55:06] C:\Users\biocbuild\bbs-3.8-bioc\tmpdir\RtmpK8Yhi6\file39fc49c83a95.csv
FATAL [2019-04-13 04:55:06] not logical(1).
FATAL [2019-04-13 04:55:06]
FATAL [2019-04-13 04:55:06] This is most likely a user error due to invalid input data or
FATAL [2019-04-13 04:55:06] parameters (PureCN 1.12.2).
INFO [2019-04-13 04:55:07] 576 on-target bins with low coverage in all samples.
WARN [2019-04-13 04:55:07] You are likely not using the correct baits file!
WARN [2019-04-13 04:55:07] Allosome coverage missing, cannot determine sex.
WARN [2019-04-13 04:55:07] Allosome coverage missing, cannot determine sex.
INFO [2019-04-13 04:55:07] Processing on-target regions...
INFO [2019-04-13 04:55:08] Removing 930 intervals with low coverage in normalDB.
INFO [2019-04-13 04:55:08] Removing 1 intervals with zero coverage in more than 3% of normalDB.
INFO [2019-04-13 04:55:10] 576 on-target bins with low coverage in all samples.
WARN [2019-04-13 04:55:10] You are likely not using the correct baits file!
WARN [2019-04-13 04:55:10] Allosome coverage missing, cannot determine sex.
WARN [2019-04-13 04:55:10] Allosome coverage missing, cannot determine sex.
INFO [2019-04-13 04:55:10] Processing on-target regions...
INFO [2019-04-13 04:55:11] Removing 930 intervals with low coverage in normalDB.
INFO [2019-04-13 04:55:11] Removing 1 intervals with zero coverage in more than 3% of normalDB.
INFO [2019-04-13 04:55:13] 576 on-target bins with low coverage in all samples.
WARN [2019-04-13 04:55:13] You are likely not using the correct baits file!
WARN [2019-04-13 04:55:13] Allosome coverage missing, cannot determine sex.
WARN [2019-04-13 04:55:13] Allosome coverage missing, cannot determine sex.
INFO [2019-04-13 04:55:13] Processing on-target regions...
INFO [2019-04-13 04:55:13] Removing 930 intervals with low coverage in normalDB.
INFO [2019-04-13 04:55:13] Removing 1 intervals with zero coverage in more than 3% of normalDB.
INFO [2019-04-13 04:55:15] 576 on-target bins with low coverage in all samples.
WARN [2019-04-13 04:55:15] You are likely not using the correct baits file!
WARN [2019-04-13 04:55:15] Allosome coverage missing, cannot determine sex.
WARN [2019-04-13 04:55:15] Allosome coverage missing, cannot determine sex.
FATAL [2019-04-13 04:55:15] Length of normal.coverage.files and sex different
FATAL [2019-04-13 04:55:15]
FATAL [2019-04-13 04:55:15] This is most likely a user error due to invalid input data or
FATAL [2019-04-13 04:55:15] parameters (PureCN 1.12.2).
WARN [2019-04-13 04:55:19] Target intervals were not sorted.
INFO [2019-04-13 04:55:19] 560 on-target bins with low coverage in all samples.
WARN [2019-04-13 04:55:19] You are likely not using the correct baits file!
WARN [2019-04-13 04:55:19] Allosome coverage missing, cannot determine sex.
WARN [2019-04-13 04:55:19] Allosome coverage missing, cannot determine sex.
WARN [2019-04-13 04:55:19] Allosome coverage missing, cannot determine sex.
INFO [2019-04-13 04:55:19] Processing on-target regions...
INFO [2019-04-13 04:55:20] Removing 978 intervals with low coverage in normalDB.
INFO [2019-04-13 04:55:20] Removing 11 intervals with zero coverage in more than 3% of normalDB.
FATAL [2019-04-13 04:55:21] tumor.coverage.file and normalDB do not align.
FATAL [2019-04-13 04:55:21]
FATAL [2019-04-13 04:55:21] This is most likely a user error due to invalid input data or
FATAL [2019-04-13 04:55:21] parameters (PureCN 1.12.2).
INFO [2019-04-13 04:55:21] ------------------------------------------------------------
INFO [2019-04-13 04:55:21] PureCN 1.12.2
INFO [2019-04-13 04:55:21] ------------------------------------------------------------
INFO [2019-04-13 04:55:21] Loading coverage files...
INFO [2019-04-13 04:55:22] Mean target coverages: 112X (tumor) 99X (normal).
WARN [2019-04-13 04:55:22] Allosome coverage missing, cannot determine sex.
WARN [2019-04-13 04:55:22] Allosome coverage missing, cannot determine sex.
INFO [2019-04-13 04:55:22] Removing 228 intervals with missing log.ratio.
FATAL [2019-04-13 04:55:22] normalDB incompatible with this PureCN version. Please re-run
FATAL [2019-04-13 04:55:22] NormalDB.R.
FATAL [2019-04-13 04:55:22]
FATAL [2019-04-13 04:55:22] This is most likely a user error due to invalid input data or
FATAL [2019-04-13 04:55:22] parameters (PureCN 1.12.2).
FATAL [2019-04-13 04:55:22] normalDB incompatible with this PureCN version. Please re-run
FATAL [2019-04-13 04:55:22] NormalDB.R.
FATAL [2019-04-13 04:55:22]
FATAL [2019-04-13 04:55:22] This is most likely a user error due to invalid input data or
FATAL [2019-04-13 04:55:22] parameters (PureCN 1.12.2).
INFO [2019-04-13 04:55:23] Removing 22 non heterozygous (in matched normal) germline SNPs.
INFO [2019-04-13 04:55:23] Initial testing for significant sample cross-contamination: unlikely
INFO [2019-04-13 04:55:23] Removing 48 variants with AF < 0.030 or AF >= 1.000 or less than 4 supporting reads or depth < 15.
INFO [2019-04-13 04:55:23] Removing 16 low quality variants with BQ < 25.
INFO [2019-04-13 04:55:24] Removing 22 non heterozygous (in matched normal) germline SNPs.
INFO [2019-04-13 04:55:24] Initial testing for significant sample cross-contamination: unlikely
INFO [2019-04-13 04:55:24] Removing 48 variants with AF < 0.030 or AF >= 1.000 or less than 4 supporting reads or depth < 15.
INFO [2019-04-13 04:55:24] Removing 6 blacklisted variants.
INFO [2019-04-13 04:55:24] Removing 16 low quality variants with BQ < 25.
INFO [2019-04-13 04:55:25] Removing 22 non heterozygous (in matched normal) germline SNPs.
INFO [2019-04-13 04:55:25] Initial testing for significant sample cross-contamination: unlikely
INFO [2019-04-13 04:55:25] Removing 48 variants with AF < 0.030 or AF >= 1.000 or less than 4 supporting reads or depth < 15.
INFO [2019-04-13 04:55:25] Removing 16 low quality variants with BQ < 25.
WARN [2019-04-13 04:55:25] MuTect stats file lacks contig and position columns.
INFO [2019-04-13 04:55:25] Removing 22 non heterozygous (in matched normal) germline SNPs.
INFO [2019-04-13 04:55:26] Initial testing for significant sample cross-contamination: unlikely
INFO [2019-04-13 04:55:26] Removing 48 variants with AF < 0.030 or AF >= 1.000 or less than 4 supporting reads or depth < 15.
INFO [2019-04-13 04:55:26] Removing 16 low quality variants with BQ < 25.
WARN [2019-04-13 04:55:26] MuTect stats file and VCF file do not align perfectly. Will remove 2325 unmatched variants.
WARN [2019-04-13 04:55:26] MuTect stats file lacks failure_reasons column. Keeping all variants listed in stats file.
INFO [2019-04-13 04:55:26] Removing 0 non heterozygous (in matched normal) germline SNPs.
INFO [2019-04-13 04:55:27] Removing 0 variants with AF < 0.030 or AF >= 1.000 or less than 4 supporting reads or depth < 15.
INFO [2019-04-13 04:55:27] Removing 0 low quality variants with BQ < 25.
INFO [2019-04-13 04:55:27] Removing 22 non heterozygous (in matched normal) germline SNPs.
INFO [2019-04-13 04:55:27] Initial testing for significant sample cross-contamination: unlikely
INFO [2019-04-13 04:55:27] Removing 48 variants with AF < 0.030 or AF >= 1.000 or less than 4 supporting reads or depth < 15.
INFO [2019-04-13 04:55:28] Removing 22 non heterozygous (in matched normal) germline SNPs.
INFO [2019-04-13 04:55:28] Initial testing for significant sample cross-contamination: unlikely
INFO [2019-04-13 04:55:28] Removing 48 variants with AF < 0.030 or AF >= 1.000 or less than 4 supporting reads or depth < 15.
INFO [2019-04-13 04:55:28] Found 11 variants in VCF file.
WARN [2019-04-13 04:55:29] vcf.file has no DB info field for membership in germline databases. Found and used valid population allele frequency > 0.001000 instead.
INFO [2019-04-13 04:55:29] 1 (9.1%) variants annotated as likely germline (DB INFO flag).
INFO [2019-04-13 04:55:31] Found 11 variants in VCF file.
WARN [2019-04-13 04:55:32] vcf.file has no DB info field for membership in germline databases. Found and used valid population allele frequency > 0.001000 instead.
INFO [2019-04-13 04:55:32] 1 (9.1%) variants annotated as likely germline (DB INFO flag).
INFO [2019-04-13 04:55:33] Found 11 variants in VCF file.
WARN [2019-04-13 04:55:33] vcf.file has no DB or POP_AF info field for membership in germline databases. Guessing it based on available dbSNP ID.
INFO [2019-04-13 04:55:33] 1 (9.1%) variants annotated as likely germline (DB INFO flag).
WARN [2019-04-13 04:55:34] Allosome coverage missing, cannot determine sex.
WARN [2019-04-13 04:55:34] Allosome coverage missing, cannot determine sex.
INFO [2019-04-13 04:55:34] Mean coverages: chrX: 108.58, chrY: 108.58, chr1-22: 104.87.
INFO [2019-04-13 04:55:34] Mean coverages: chrX: 108.58, chrY: 2.17, chr1-22: 104.87.
INFO [2019-04-13 04:55:34] Mean coverages: chrX: 108.58, chrY: 5.17, chr1-22: 104.87.
INFO [2019-04-13 04:55:34] No homozygous variants in VCF, provide unfiltered VCF.
INFO [2019-04-13 04:55:35] No germline variants in VCF.
FATAL [2019-04-13 04:55:35] No solution with id hello
FATAL [2019-04-13 04:55:35]
FATAL [2019-04-13 04:55:35] This is most likely a user error due to invalid input data or
FATAL [2019-04-13 04:55:35] parameters (PureCN 1.12.2).
FATAL [2019-04-13 04:55:35] No solution with id 10
FATAL [2019-04-13 04:55:35]
FATAL [2019-04-13 04:55:35] This is most likely a user error due to invalid input data or
FATAL [2019-04-13 04:55:35] parameters (PureCN 1.12.2).
FATAL [2019-04-13 04:55:36] all.data and w have different lengths.
FATAL [2019-04-13 04:55:36]
FATAL [2019-04-13 04:55:36] This is most likely a user error due to invalid input data or
FATAL [2019-04-13 04:55:36] parameters (PureCN 1.12.2).
WARN [2019-04-13 04:55:36] No mappability scores provided.
WARN [2019-04-13 04:55:36] No reptiming scores provided.
INFO [2019-04-13 04:55:36] Calculating GC-content...
WARN [2019-04-13 04:55:37] No mappability scores provided.
WARN [2019-04-13 04:55:37] No reptiming scores provided.
INFO [2019-04-13 04:55:37] Calculating GC-content...
INFO [2019-04-13 04:55:37] Removing 1 targets overlapping with exclude.
WARN [2019-04-13 04:55:37] No mappability scores provided.
WARN [2019-04-13 04:55:37] No reptiming scores provided.
INFO [2019-04-13 04:55:37] Calculating GC-content...
FATAL [2019-04-13 04:55:37] off.target.padding must be negative.
FATAL [2019-04-13 04:55:37]
FATAL [2019-04-13 04:55:37] This is most likely a user error due to invalid input data or
FATAL [2019-04-13 04:55:37] parameters (PureCN 1.12.2).
FATAL [2019-04-13 04:55:37] Interval coordinates should start at 1, not at 0
FATAL [2019-04-13 04:55:37]
FATAL [2019-04-13 04:55:37] This is most likely a user error due to invalid input data or
FATAL [2019-04-13 04:55:37] parameters (PureCN 1.12.2).
FATAL [2019-04-13 04:55:37] No off-target regions after filtering for mappability and
FATAL [2019-04-13 04:55:37] off.target.padding
FATAL [2019-04-13 04:55:37]
FATAL [2019-04-13 04:55:37] This is most likely a user error due to invalid input data or
FATAL [2019-04-13 04:55:37] parameters (PureCN 1.12.2).
WARN [2019-04-13 04:55:37] No mappability scores provided.
INFO [2019-04-13 04:55:37] Calculating GC-content...
INFO [2019-04-13 04:55:38] Averaging reptiming into bins of size 200...
WARN [2019-04-13 04:55:38] No mappability scores provided.
INFO [2019-04-13 04:55:38] Calculating GC-content...
INFO [2019-04-13 04:55:38] Splitting 1 large targets to an average width of 200.
WARN [2019-04-13 04:55:38] No mappability scores provided.
WARN [2019-04-13 04:55:38] No reptiming scores provided.
INFO [2019-04-13 04:55:38] Calculating GC-content...
WARN [2019-04-13 04:55:38] No mappability scores provided.
WARN [2019-04-13 04:55:38] No reptiming scores provided.
INFO [2019-04-13 04:55:38] Calculating GC-content...
INFO [2019-04-13 04:55:38] Tiling off-target regions to an average width of 200000.
WARN [2019-04-13 04:55:38] No mappability scores provided.
WARN [2019-04-13 04:55:38] No reptiming scores provided.
INFO [2019-04-13 04:55:38] Calculating GC-content...
WARN [2019-04-13 04:55:38] Intervals contain off-target regions. Will not change intervals.
WARN [2019-04-13 04:55:38] No mappability scores provided.
WARN [2019-04-13 04:55:38] No reptiming scores provided.
INFO [2019-04-13 04:55:38] Calculating GC-content...
WARN [2019-04-13 04:55:39] No reptiming scores provided.
INFO [2019-04-13 04:55:39] Calculating GC-content...
WARN [2019-04-13 04:55:39] 1 intervals without mappability score (1 on-target).
INFO [2019-04-13 04:55:39] Removing 1 intervals with low mappability score (<0.50).
WARN [2019-04-13 04:55:39] No reptiming scores provided.
INFO [2019-04-13 04:55:39] Calculating GC-content...
INFO [2019-04-13 04:55:40] Tiling off-target regions to an average width of 200000.
WARN [2019-04-13 04:55:40] No reptiming scores provided.
INFO [2019-04-13 04:55:40] Calculating GC-content...
WARN [2019-04-13 04:55:40] No mappability scores provided.
WARN [2019-04-13 04:55:40] No reptiming scores provided.
INFO [2019-04-13 04:55:40] Calculating GC-content...
WARN [2019-04-13 04:55:40] Chromosome naming style of interval file (NCBI) was different from reference (UCSC).
WARN [2019-04-13 04:55:40] No mappability scores provided.
WARN [2019-04-13 04:55:40] No reptiming scores provided.
INFO [2019-04-13 04:55:40] Calculating GC-content...
FATAL [2019-04-13 04:55:40] Chromosome naming style of interval file unknown, should be UCSC.
FATAL [2019-04-13 04:55:40]
FATAL [2019-04-13 04:55:40] This is most likely a user error due to invalid input data or
FATAL [2019-04-13 04:55:40] parameters (PureCN 1.12.2).
WARN [2019-04-13 04:55:41] Chromosome naming style of interval file (NCBI) was different from reference (UCSC).
WARN [2019-04-13 04:55:41] Chromosome naming style of mappability file (NCBI) was different from reference (UCSC).
WARN [2019-04-13 04:55:41] No reptiming scores provided.
INFO [2019-04-13 04:55:41] Calculating GC-content...
INFO [2019-04-13 04:55:42] Loading coverage data...
INFO [2019-04-13 04:55:42] Mean target coverages: 71X (tumor) 99X (normal).
INFO [2019-04-13 04:55:42] Mean target coverages: 71X (tumor) 43X (normal).
INFO [2019-04-13 04:55:43] ------------------------------------------------------------
INFO [2019-04-13 04:55:43] PureCN 1.12.2
INFO [2019-04-13 04:55:43] ------------------------------------------------------------
INFO [2019-04-13 04:55:43] Using BiocParallel for parallel optimization.
INFO [2019-04-13 04:55:43] Loading coverage files...
INFO [2019-04-13 04:55:44] Mean target coverages: 112X (tumor) 99X (normal).
WARN [2019-04-13 04:55:44] Allosome coverage missing, cannot determine sex.
WARN [2019-04-13 04:55:44] Allosome coverage missing, cannot determine sex.
INFO [2019-04-13 04:55:44] Removing 228 intervals with missing log.ratio.
INFO [2019-04-13 04:55:44] Removing 22 small (< 5bp) intervals.
INFO [2019-04-13 04:55:44] Removing 19 intervals with low total coverage in normal (< 150.00 reads).
WARN [2019-04-13 04:55:44] No normalDB provided. Provide one for better results.
INFO [2019-04-13 04:55:44] Removing 233 low coverage (< 15.0000X) intervals.
INFO [2019-04-13 04:55:44] Using 9547 intervals (9547 on-target, 0 off-target).
INFO [2019-04-13 04:55:44] No off-target intervals. If this is hybrid-capture data, consider adding them.
INFO [2019-04-13 04:55:44] No interval.file provided. Cannot check if data was GC-normalized. Was it?
INFO [2019-04-13 04:55:44] Sample sex: ?
INFO [2019-04-13 04:55:44] Segmenting data...
INFO [2019-04-13 04:55:44] Interval weights found, will use weighted CBS.
INFO [2019-04-13 04:55:44] Loading pre-computed boundaries for DNAcopy...
INFO [2019-04-13 04:55:44] Setting undo.SD parameter to 1.000000.
Setting multi-figure configuration
INFO [2019-04-13 04:55:46] Found 54 segments with median size of 25.38Mb.
INFO [2019-04-13 04:55:46] Mean standard deviation of log-ratios: 0.37
INFO [2019-04-13 04:55:46] 2D-grid search of purity and ploidy...
INFO [2019-04-13 04:55:46] Local optima: 0.65/1.9, 0.33/1.9, 0.38/2.2, 0.27/2.1, 0.29/2, 0.78/2
INFO [2019-04-13 04:56:32] Done.
INFO [2019-04-13 04:56:32] ------------------------------------------------------------
INFO [2019-04-13 04:56:32] Reading C:\Users\biocbuild\bbs-3.8-bioc\tmpdir\RtmpK8Yhi6\file39fc6e0941ee.rds...
FATAL [2019-04-13 04:56:32] runAbsoluteCN was run without a VCF file.
FATAL [2019-04-13 04:56:32]
FATAL [2019-04-13 04:56:32] This is most likely a user error due to invalid input data or
FATAL [2019-04-13 04:56:32] parameters (PureCN 1.12.2).
FATAL [2019-04-13 04:56:32] runAbsoluteCN was run without a VCF file.
FATAL [2019-04-13 04:56:32]
FATAL [2019-04-13 04:56:32] This is most likely a user error due to invalid input data or
FATAL [2019-04-13 04:56:32] parameters (PureCN 1.12.2).
INFO [2019-04-13 04:56:32] ------------------------------------------------------------
INFO [2019-04-13 04:56:32] PureCN 1.12.2
INFO [2019-04-13 04:56:32] ------------------------------------------------------------
INFO [2019-04-13 04:56:32] Loading coverage files...
FATAL [2019-04-13 04:56:33] Need a normal coverage file if log.ratio and seg.file are not provided.
FATAL [2019-04-13 04:56:33]
FATAL [2019-04-13 04:56:33] This is most likely a user error due to invalid input data or
FATAL [2019-04-13 04:56:33] parameters (PureCN 1.12.2).
INFO [2019-04-13 04:56:33] ------------------------------------------------------------
INFO [2019-04-13 04:56:33] PureCN 1.12.2
INFO [2019-04-13 04:56:33] ------------------------------------------------------------
FATAL [2019-04-13 04:56:33] min.ploidy or max.ploidy not within expected range.
FATAL [2019-04-13 04:56:33]
FATAL [2019-04-13 04:56:33] This is most likely a user error due to invalid input data or
FATAL [2019-04-13 04:56:33] parameters (PureCN 1.12.2).
INFO [2019-04-13 04:56:33] ------------------------------------------------------------
INFO [2019-04-13 04:56:33] PureCN 1.12.2
INFO [2019-04-13 04:56:33] ------------------------------------------------------------
FATAL [2019-04-13 04:56:33] min.ploidy or max.ploidy not within expected range.
FATAL [2019-04-13 04:56:33]
FATAL [2019-04-13 04:56:33] This is most likely a user error due to invalid input data or
FATAL [2019-04-13 04:56:33] parameters (PureCN 1.12.2).
INFO [2019-04-13 04:56:33] ------------------------------------------------------------
INFO [2019-04-13 04:56:33] PureCN 1.12.2
INFO [2019-04-13 04:56:33] ------------------------------------------------------------
INFO [2019-04-13 04:56:33] ------------------------------------------------------------
INFO [2019-04-13 04:56:33] PureCN 1.12.2
INFO [2019-04-13 04:56:33] ------------------------------------------------------------
INFO [2019-04-13 04:56:33] ------------------------------------------------------------
INFO [2019-04-13 04:56:33] PureCN 1.12.2
INFO [2019-04-13 04:56:33] ------------------------------------------------------------
FATAL [2019-04-13 04:56:33] test.num.copy not within expected range.
FATAL [2019-04-13 04:56:33]
FATAL [2019-04-13 04:56:33] This is most likely a user error due to invalid input data or
FATAL [2019-04-13 04:56:33] parameters (PureCN 1.12.2).
INFO [2019-04-13 04:56:33] ------------------------------------------------------------
INFO [2019-04-13 04:56:33] PureCN 1.12.2
INFO [2019-04-13 04:56:33] ------------------------------------------------------------
WARN [2019-04-13 04:56:33] test.num.copy outside recommended range.
FATAL [2019-04-13 04:56:33] max.non.clonal not within expected range or format.
FATAL [2019-04-13 04:56:33]
FATAL [2019-04-13 04:56:33] This is most likely a user error due to invalid input data or
FATAL [2019-04-13 04:56:33] parameters (PureCN 1.12.2).
INFO [2019-04-13 04:56:33] ------------------------------------------------------------
INFO [2019-04-13 04:56:33] PureCN 1.12.2
INFO [2019-04-13 04:56:33] ------------------------------------------------------------
WARN [2019-04-13 04:56:33] test.num.copy outside recommended range.
FATAL [2019-04-13 04:56:33] max.non.clonal not within expected range or format.
FATAL [2019-04-13 04:56:33]
FATAL [2019-04-13 04:56:33] This is most likely a user error due to invalid input data or
FATAL [2019-04-13 04:56:33] parameters (PureCN 1.12.2).
INFO [2019-04-13 04:56:33] ------------------------------------------------------------
INFO [2019-04-13 04:56:33] PureCN 1.12.2
INFO [2019-04-13 04:56:33] ------------------------------------------------------------
FATAL [2019-04-13 04:56:33] test.purity not within expected range.
FATAL [2019-04-13 04:56:33]
FATAL [2019-04-13 04:56:33] This is most likely a user error due to invalid input data or
FATAL [2019-04-13 04:56:33] parameters (PureCN 1.12.2).
INFO [2019-04-13 04:56:33] ------------------------------------------------------------
INFO [2019-04-13 04:56:33] PureCN 1.12.2
INFO [2019-04-13 04:56:33] ------------------------------------------------------------
INFO [2019-04-13 04:56:33] Loading coverage files...
FATAL [2019-04-13 04:56:33] Tumor and normal are identical. This won't give any meaningful results
FATAL [2019-04-13 04:56:33] and I'm stopping here.
FATAL [2019-04-13 04:56:33]
FATAL [2019-04-13 04:56:33] This is most likely a user error due to invalid input data or
FATAL [2019-04-13 04:56:33] parameters (PureCN 1.12.2).
INFO [2019-04-13 04:56:33] ------------------------------------------------------------
INFO [2019-04-13 04:56:33] PureCN 1.12.2
INFO [2019-04-13 04:56:33] ------------------------------------------------------------
INFO [2019-04-13 04:56:33] Loading coverage files...
FATAL [2019-04-13 04:56:34] Length of log.ratio different from tumor coverage.
FATAL [2019-04-13 04:56:34]
FATAL [2019-04-13 04:56:34] This is most likely a user error due to invalid input data or
FATAL [2019-04-13 04:56:34] parameters (PureCN 1.12.2).
INFO [2019-04-13 04:56:34] ------------------------------------------------------------
INFO [2019-04-13 04:56:34] PureCN 1.12.2
INFO [2019-04-13 04:56:34] ------------------------------------------------------------
FATAL [2019-04-13 04:56:34] prior.purity must have the same length as test.purity.
FATAL [2019-04-13 04:56:34]
FATAL [2019-04-13 04:56:34] This is most likely a user error due to invalid input data or
FATAL [2019-04-13 04:56:34] parameters (PureCN 1.12.2).
INFO [2019-04-13 04:56:34] ------------------------------------------------------------
INFO [2019-04-13 04:56:34] PureCN 1.12.2
INFO [2019-04-13 04:56:34] ------------------------------------------------------------
FATAL [2019-04-13 04:56:34] min.gof not within expected range or format.
FATAL [2019-04-13 04:56:34]
FATAL [2019-04-13 04:56:34] This is most likely a user error due to invalid input data or
FATAL [2019-04-13 04:56:34] parameters (PureCN 1.12.2).
INFO [2019-04-13 04:56:34] ------------------------------------------------------------
INFO [2019-04-13 04:56:34] PureCN 1.12.2
INFO [2019-04-13 04:56:34] ------------------------------------------------------------
FATAL [2019-04-13 04:56:34] prior.purity not within expected range or format.
FATAL [2019-04-13 04:56:34]
FATAL [2019-04-13 04:56:34] This is most likely a user error due to invalid input data or
FATAL [2019-04-13 04:56:34] parameters (PureCN 1.12.2).
INFO [2019-04-13 04:56:34] ------------------------------------------------------------
INFO [2019-04-13 04:56:34] PureCN 1.12.2
INFO [2019-04-13 04:56:34] ------------------------------------------------------------
FATAL [2019-04-13 04:56:34] prior.purity must add to 1. Sum is 1.5
FATAL [2019-04-13 04:56:34]
FATAL [2019-04-13 04:56:34] This is most likely a user error due to invalid input data or
FATAL [2019-04-13 04:56:34] parameters (PureCN 1.12.2).
INFO [2019-04-13 04:56:34] ------------------------------------------------------------
INFO [2019-04-13 04:56:34] PureCN 1.12.2
INFO [2019-04-13 04:56:34] ------------------------------------------------------------
FATAL [2019-04-13 04:56:34] max.homozygous.loss not within expected range or format.
FATAL [2019-04-13 04:56:34]
FATAL [2019-04-13 04:56:34] This is most likely a user error due to invalid input data or
FATAL [2019-04-13 04:56:34] parameters (PureCN 1.12.2).
INFO [2019-04-13 04:56:34] ------------------------------------------------------------
INFO [2019-04-13 04:56:34] PureCN 1.12.2
INFO [2019-04-13 04:56:34] ------------------------------------------------------------
FATAL [2019-04-13 04:56:34] prior.K not within expected range or format.
FATAL [2019-04-13 04:56:34]
FATAL [2019-04-13 04:56:34] This is most likely a user error due to invalid input data or
FATAL [2019-04-13 04:56:34] parameters (PureCN 1.12.2).
INFO [2019-04-13 04:56:34] ------------------------------------------------------------
INFO [2019-04-13 04:56:34] PureCN 1.12.2
INFO [2019-04-13 04:56:34] ------------------------------------------------------------
FATAL [2019-04-13 04:56:34] prior.contamination not within expected range or format.
FATAL [2019-04-13 04:56:34]
FATAL [2019-04-13 04:56:34] This is most likely a user error due to invalid input data or
FATAL [2019-04-13 04:56:34] parameters (PureCN 1.12.2).
INFO [2019-04-13 04:56:34] ------------------------------------------------------------
INFO [2019-04-13 04:56:34] PureCN 1.12.2
INFO [2019-04-13 04:56:34] ------------------------------------------------------------
FATAL [2019-04-13 04:56:34] Iterations not in the expected range from 10 to 250.
FATAL [2019-04-13 04:56:34]
FATAL [2019-04-13 04:56:34] This is most likely a user error due to invalid input data or
FATAL [2019-04-13 04:56:34] parameters (PureCN 1.12.2).
INFO [2019-04-13 04:56:34] ------------------------------------------------------------
INFO [2019-04-13 04:56:34] PureCN 1.12.2
INFO [2019-04-13 04:56:34] ------------------------------------------------------------
FATAL [2019-04-13 04:56:34] Iterations not in the expected range from 10 to 250.
FATAL [2019-04-13 04:56:34]
FATAL [2019-04-13 04:56:34] This is most likely a user error due to invalid input data or
FATAL [2019-04-13 04:56:34] parameters (PureCN 1.12.2).
INFO [2019-04-13 04:56:34] ------------------------------------------------------------
INFO [2019-04-13 04:56:34] PureCN 1.12.2
INFO [2019-04-13 04:56:34] ------------------------------------------------------------
INFO [2019-04-13 04:56:34] Loading coverage files...
FATAL [2019-04-13 04:56:34] Missing tumor.coverage.file requires seg.file or log.ratio and
FATAL [2019-04-13 04:56:34] interval.file.
FATAL [2019-04-13 04:56:34]
FATAL [2019-04-13 04:56:34] This is most likely a user error due to invalid input data or
FATAL [2019-04-13 04:56:34] parameters (PureCN 1.12.2).
INFO [2019-04-13 04:56:34] ------------------------------------------------------------
INFO [2019-04-13 04:56:34] PureCN 1.12.2
INFO [2019-04-13 04:56:34] ------------------------------------------------------------
INFO [2019-04-13 04:56:34] ------------------------------------------------------------
INFO [2019-04-13 04:56:34] PureCN 1.12.2
INFO [2019-04-13 04:56:34] ------------------------------------------------------------
INFO [2019-04-13 04:56:34] Loading coverage files...
FATAL [2019-04-13 04:56:35] Interval files in normal and tumor different.
FATAL [2019-04-13 04:56:35]
FATAL [2019-04-13 04:56:35] This is most likely a user error due to invalid input data or
FATAL [2019-04-13 04:56:35] parameters (PureCN 1.12.2).
INFO [2019-04-13 04:56:35] ------------------------------------------------------------
INFO [2019-04-13 04:56:35] PureCN 1.12.2
INFO [2019-04-13 04:56:35] ------------------------------------------------------------
INFO [2019-04-13 04:56:35] Loading coverage files...
INFO [2019-04-13 04:56:35] Mean target coverages: 0X (tumor) 99X (normal).
WARN [2019-04-13 04:56:35] Large difference in coverage of tumor and normal.
WARN [2019-04-13 04:56:35] Allosome coverage missing, cannot determine sex.
WARN [2019-04-13 04:56:35] Allosome coverage missing, cannot determine sex.
INFO [2019-04-13 04:56:35] Removing 10049 intervals with missing log.ratio.
WARN [2019-04-13 04:56:35] No normalDB provided. Provide one for better results.
INFO [2019-04-13 04:56:35] Using 0 intervals (0 on-target, 0 off-target).
FATAL [2019-04-13 04:56:35] No intervals passing filters.
FATAL [2019-04-13 04:56:35]
FATAL [2019-04-13 04:56:35] This is most likely a user error due to invalid input data or
FATAL [2019-04-13 04:56:35] parameters (PureCN 1.12.2).
INFO [2019-04-13 04:56:35] ------------------------------------------------------------
INFO [2019-04-13 04:56:35] PureCN 1.12.2
INFO [2019-04-13 04:56:35] ------------------------------------------------------------
INFO [2019-04-13 04:56:35] Loading coverage files...
INFO [2019-04-13 04:56:36] Mean target coverages: 112X (tumor) 99X (normal).
WARN [2019-04-13 04:56:36] Allosome coverage missing, cannot determine sex.
WARN [2019-04-13 04:56:36] Allosome coverage missing, cannot determine sex.
INFO [2019-04-13 04:56:36] Removing 228 intervals with missing log.ratio.
INFO [2019-04-13 04:56:36] Removing 22 small (< 5bp) intervals.
INFO [2019-04-13 04:56:36] Removing 19 intervals with low total coverage in normal (< 150.00 reads).
WARN [2019-04-13 04:56:36] No normalDB provided. Provide one for better results.
INFO [2019-04-13 04:56:36] Removing 233 low coverage (< 15.0000X) intervals.
INFO [2019-04-13 04:56:36] Using 9547 intervals (9547 on-target, 0 off-target).
INFO [2019-04-13 04:56:36] No off-target intervals. If this is hybrid-capture data, consider adding them.
INFO [2019-04-13 04:56:36] No interval.file provided. Cannot check if data was GC-normalized. Was it?
INFO [2019-04-13 04:56:36] Loading VCF...
INFO [2019-04-13 04:56:36] Found 127 variants in VCF file.
INFO [2019-04-13 04:56:36] 127 (100.0%) variants annotated as likely germline (DB INFO flag).
INFO [2019-04-13 04:56:37] LIB-02240e4 is tumor in VCF file.
INFO [2019-04-13 04:56:37] No homozygous variants in VCF, provide unfiltered VCF.
INFO [2019-04-13 04:56:37] Removing 1 non heterozygous (in matched normal) germline SNPs.
INFO [2019-04-13 04:56:37] Removing 1 variants with AF < 0.030 or AF >= 1.000 or less than 4 supporting reads or depth < 15.
INFO [2019-04-13 04:56:37] Removing 0 low quality variants with BQ < 25.
INFO [2019-04-13 04:56:37] Total size of targeted genomic region: 1.63Mb (2.58Mb with 50bp padding).
INFO [2019-04-13 04:56:37] 1.0% of targets contain variants.
INFO [2019-04-13 04:56:37] Removing 2 variants outside intervals.
INFO [2019-04-13 04:56:37] Found SOMATIC annotation in VCF.
INFO [2019-04-13 04:56:37] Setting somatic prior probabilities for somatic variants to 0.999000 or to 0.000100 otherwise.
INFO [2019-04-13 04:56:37] Found SOMATIC annotation in VCF. Setting mapping bias to 0.973.
INFO [2019-04-13 04:56:37] Sample sex: ?
INFO [2019-04-13 04:56:37] Segmenting data...
INFO [2019-04-13 04:56:38] Loading pre-computed boundaries for DNAcopy...
INFO [2019-04-13 04:56:38] Setting undo.SD parameter to 1.000000.
INFO [2019-04-13 04:56:39] Setting prune.hclust.h parameter to 0.200000.
INFO [2019-04-13 04:56:39] Found 54 segments with median size of 24.88Mb.
INFO [2019-04-13 04:56:39] Using 123 variants.
INFO [2019-04-13 04:56:39] Mean standard deviation of log-ratios: 0.37
INFO [2019-04-13 04:56:39] 2D-grid search of purity and ploidy...
INFO [2019-04-13 04:56:40] Local optima: 0.62/1.9, 0.3/1.9, 0.32/2
INFO [2019-04-13 04:56:40] Testing local optimum 1/3 at purity 0.62 and total ploidy 1.90...
INFO [2019-04-13 04:56:42] Fitting variants for purity 0.65, tumor ploidy 1.74 and contamination 0.01.
INFO [2019-04-13 04:56:43] Optimized purity: 0.65
INFO [2019-04-13 04:56:43] Testing local optimum 2/3 at purity 0.30 and total ploidy 1.90...
INFO [2019-04-13 04:56:43] Recalibrating log-ratios...
INFO [2019-04-13 04:56:43] Testing local optimum 2/3 at purity 0.30 and total ploidy 1.90...
INFO [2019-04-13 04:56:43] Recalibrating log-ratios...
INFO [2019-04-13 04:56:43] Testing local optimum 2/3 at purity 0.30 and total ploidy 1.90...
INFO [2019-04-13 04:56:44] Recalibrating log-ratios...
INFO [2019-04-13 04:56:44] Testing local optimum 2/3 at purity 0.30 and total ploidy 1.90...
INFO [2019-04-13 04:56:44] Testing local optimum 3/3 at purity 0.32 and total ploidy 2.00...
INFO [2019-04-13 04:56:45] Recalibrating log-ratios...
INFO [2019-04-13 04:56:45] Testing local optimum 3/3 at purity 0.32 and total ploidy 2.00...
INFO [2019-04-13 04:56:45] Recalibrating log-ratios...
INFO [2019-04-13 04:56:45] Testing local optimum 3/3 at purity 0.32 and total ploidy 2.00...
INFO [2019-04-13 04:56:45] Recalibrating log-ratios...
INFO [2019-04-13 04:56:45] Testing local optimum 3/3 at purity 0.32 and total ploidy 2.00...
INFO [2019-04-13 04:56:46] Done.
INFO [2019-04-13 04:56:46] ------------------------------------------------------------
INFO [2019-04-13 04:56:46] ------------------------------------------------------------
INFO [2019-04-13 04:56:46] PureCN 1.12.2
INFO [2019-04-13 04:56:46] ------------------------------------------------------------
INFO [2019-04-13 04:56:46] Loading coverage files...
INFO [2019-04-13 04:56:47] Mean target coverages: 112X (tumor) 99X (normal).
WARN [2019-04-13 04:56:47] Allosome coverage missing, cannot determine sex.
WARN [2019-04-13 04:56:47] Allosome coverage missing, cannot determine sex.
INFO [2019-04-13 04:56:47] Removing 228 intervals with missing log.ratio.
INFO [2019-04-13 04:56:47] Removing 22 small (< 5bp) intervals.
INFO [2019-04-13 04:56:47] Removing 19 intervals with low total coverage in normal (< 150.00 reads).
WARN [2019-04-13 04:56:47] No normalDB provided. Provide one for better results.
INFO [2019-04-13 04:56:47] Removing 233 low coverage (< 15.0000X) intervals.
INFO [2019-04-13 04:56:47] Using 9547 intervals (9547 on-target, 0 off-target).
INFO [2019-04-13 04:56:47] No off-target intervals. If this is hybrid-capture data, consider adding them.
INFO [2019-04-13 04:56:47] No interval.file provided. Cannot check if data was GC-normalized. Was it?
INFO [2019-04-13 04:56:47] Loading VCF...
INFO [2019-04-13 04:56:47] Found 127 variants in VCF file.
WARN [2019-04-13 04:56:47] DP FORMAT field contains NAs. Removing 3 variants.
WARN [2019-04-13 04:56:47] DB INFO flag contains NAs
INFO [2019-04-13 04:56:47] 119 (96.0%) variants annotated as likely germline (DB INFO flag).
INFO [2019-04-13 04:56:47] LIB-02240e4 is tumor in VCF file.
INFO [2019-04-13 04:56:47] No homozygous variants in VCF, provide unfiltered VCF.
INFO [2019-04-13 04:56:48] Removing 1 variants with AF < 0.030 or AF >= 0.970 or less than 4 supporting reads or depth < 15.
INFO [2019-04-13 04:56:48] Removing 0 low quality variants with BQ < 25.
INFO [2019-04-13 04:56:48] Total size of targeted genomic region: 1.63Mb (2.58Mb with 50bp padding).
INFO [2019-04-13 04:56:48] 1.0% of targets contain variants.
INFO [2019-04-13 04:56:48] Setting somatic prior probabilities for dbSNP hits to 0.000500 or to 0.500000 otherwise.
INFO [2019-04-13 04:56:48] Sample sex: ?
INFO [2019-04-13 04:56:48] Segmenting data...
INFO [2019-04-13 04:56:48] Loading pre-computed boundaries for DNAcopy...
INFO [2019-04-13 04:56:48] Setting undo.SD parameter to 1.000000.
Setting multi-figure configuration
INFO [2019-04-13 04:56:50] Setting prune.hclust.h parameter to 0.200000.
INFO [2019-04-13 04:56:50] Found 54 segments with median size of 24.88Mb.
INFO [2019-04-13 04:56:50] Using 123 variants.
INFO [2019-04-13 04:56:50] Mean standard deviation of log-ratios: 0.37
INFO [2019-04-13 04:56:50] 2D-grid search of purity and ploidy...
INFO [2019-04-13 04:56:51] Local optima: 0.62/1.9, 0.3/1.9, 0.32/2
INFO [2019-04-13 04:56:51] Testing local optimum 1/3 at purity 0.62 and total ploidy 1.90...
INFO [2019-04-13 04:56:53] Fitting variants for purity 0.65, tumor ploidy 1.74 and contamination 0.01.
INFO [2019-04-13 04:56:54] Optimized purity: 0.65
INFO [2019-04-13 04:56:54] Testing local optimum 2/3 at purity 0.30 and total ploidy 1.90...
INFO [2019-04-13 04:56:54] Recalibrating log-ratios...
INFO [2019-04-13 04:56:54] Testing local optimum 2/3 at purity 0.30 and total ploidy 1.90...
INFO [2019-04-13 04:56:55] Recalibrating log-ratios...
INFO [2019-04-13 04:56:55] Testing local optimum 2/3 at purity 0.30 and total ploidy 1.90...
INFO [2019-04-13 04:56:55] Recalibrating log-ratios...
INFO [2019-04-13 04:56:55] Testing local optimum 2/3 at purity 0.30 and total ploidy 1.90...
INFO [2019-04-13 04:56:55] Testing local optimum 3/3 at purity 0.32 and total ploidy 2.00...
INFO [2019-04-13 04:56:56] Recalibrating log-ratios...
INFO [2019-04-13 04:56:56] Testing local optimum 3/3 at purity 0.32 and total ploidy 2.00...
INFO [2019-04-13 04:56:56] Recalibrating log-ratios...
INFO [2019-04-13 04:56:56] Testing local optimum 3/3 at purity 0.32 and total ploidy 2.00...
INFO [2019-04-13 04:56:57] Recalibrating log-ratios...
INFO [2019-04-13 04:56:57] Testing local optimum 3/3 at purity 0.32 and total ploidy 2.00...
INFO [2019-04-13 04:56:57] Done.
INFO [2019-04-13 04:56:57] ------------------------------------------------------------
INFO [2019-04-13 04:56:57] ------------------------------------------------------------
INFO [2019-04-13 04:56:57] PureCN 1.12.2
INFO [2019-04-13 04:56:57] ------------------------------------------------------------
INFO [2019-04-13 04:56:57] Loading coverage files...
INFO [2019-04-13 04:56:58] Mean target coverages: 112X (tumor) 99X (normal).
WARN [2019-04-13 04:56:58] Allosome coverage missing, cannot determine sex.
WARN [2019-04-13 04:56:58] Allosome coverage missing, cannot determine sex.
INFO [2019-04-13 04:56:58] No Gene column in interval.file. You won't get gene-level calls.
INFO [2019-04-13 04:56:58] Removing 228 intervals with missing log.ratio.
INFO [2019-04-13 04:56:58] Removing 15 low/high GC targets.
INFO [2019-04-13 04:56:58] Removing 21 small (< 5bp) intervals.
INFO [2019-04-13 04:56:58] Removing 19 intervals with low total coverage in normal (< 150.00 reads).
WARN [2019-04-13 04:56:58] No normalDB provided. Provide one for better results.
INFO [2019-04-13 04:56:58] Removing 232 low coverage (< 15.0000X) intervals.
INFO [2019-04-13 04:56:58] Using 9534 intervals (9534 on-target, 0 off-target).
INFO [2019-04-13 04:56:58] No off-target intervals. If this is hybrid-capture data, consider adding them.
INFO [2019-04-13 04:56:58] AT/GC dropout: 1.03 (tumor), 1.03 (normal).
INFO [2019-04-13 04:56:58] Loading VCF...
INFO [2019-04-13 04:56:59] Found 127 variants in VCF file.
INFO [2019-04-13 04:56:59] 127 (100.0%) variants annotated as likely germline (DB INFO flag).
INFO [2019-04-13 04:56:59] LIB-02240e4 is tumor in VCF file.
INFO [2019-04-13 04:56:59] No homozygous variants in VCF, provide unfiltered VCF.
INFO [2019-04-13 04:56:59] Removing 1 non heterozygous (in matched normal) germline SNPs.
INFO [2019-04-13 04:57:00] Removing 1 variants with AF < 0.030 or AF >= Inf or less than 4 supporting reads or depth < 15.
INFO [2019-04-13 04:57:00] Removing 0 low quality variants with BQ < 25.
INFO [2019-04-13 04:57:00] Total size of targeted genomic region: 1.63Mb (2.58Mb with 50bp padding).
INFO [2019-04-13 04:57:00] 1.0% of targets contain variants.
INFO [2019-04-13 04:57:00] Removing 2 variants outside intervals.
INFO [2019-04-13 04:57:00] Found SOMATIC annotation in VCF.
INFO [2019-04-13 04:57:00] Setting somatic prior probabilities for somatic variants to 0.999000 or to 0.000100 otherwise.
INFO [2019-04-13 04:57:00] Found SOMATIC annotation in VCF. Setting mapping bias to 0.973.
INFO [2019-04-13 04:57:00] Sample sex: ?
INFO [2019-04-13 04:57:00] Segmenting data...
INFO [2019-04-13 04:57:00] Loading pre-computed boundaries for DNAcopy...
INFO [2019-04-13 04:57:00] Setting undo.SD parameter to 1.000000.
Setting multi-figure configuration
INFO [2019-04-13 04:57:02] Setting prune.hclust.h parameter to 0.200000.
INFO [2019-04-13 04:57:02] Found 54 segments with median size of 24.88Mb.
INFO [2019-04-13 04:57:02] Using 123 variants.
INFO [2019-04-13 04:57:02] Mean standard deviation of log-ratios: 0.37
INFO [2019-04-13 04:57:02] 2D-grid search of purity and ploidy...
INFO [2019-04-13 04:57:02] Local optima: 0.65/1.9, 0.33/1.9, 0.38/2.2, 0.27/2.1, 0.29/2, 0.78/2
INFO [2019-04-13 04:57:02] Testing local optimum 1/6 at purity 0.65 and total ploidy 1.90...
INFO [2019-04-13 04:57:03] Recalibrating log-ratios...
INFO [2019-04-13 04:57:03] Testing local optimum 1/6 at purity 0.65 and total ploidy 1.90...
INFO [2019-04-13 04:57:03] Recalibrating log-ratios...
INFO [2019-04-13 04:57:03] Testing local optimum 1/6 at purity 0.65 and total ploidy 1.90...
INFO [2019-04-13 04:57:03] Recalibrating log-ratios...
INFO [2019-04-13 04:57:03] Testing local optimum 1/6 at purity 0.65 and total ploidy 1.90...
INFO [2019-04-13 04:57:04] Testing local optimum 2/6 at purity 0.33 and total ploidy 1.90...
INFO [2019-04-13 04:57:04] Recalibrating log-ratios...
INFO [2019-04-13 04:57:04] Testing local optimum 2/6 at purity 0.33 and total ploidy 1.90...
INFO [2019-04-13 04:57:05] Recalibrating log-ratios...
INFO [2019-04-13 04:57:05] Testing local optimum 2/6 at purity 0.33 and total ploidy 1.90...
INFO [2019-04-13 04:57:05] Recalibrating log-ratios...
INFO [2019-04-13 04:57:05] Testing local optimum 2/6 at purity 0.33 and total ploidy 1.90...
INFO [2019-04-13 04:57:05] Testing local optimum 3/6 at purity 0.38 and total ploidy 2.20...
INFO [2019-04-13 04:57:06] Fitting variants for purity 0.40, tumor ploidy 2.44 and contamination 0.01.
INFO [2019-04-13 04:57:08] Optimized purity: 0.40
INFO [2019-04-13 04:57:08] Testing local optimum 4/6 at purity 0.27 and total ploidy 2.10...
INFO [2019-04-13 04:57:08] Recalibrating log-ratios...
INFO [2019-04-13 04:57:08] Testing local optimum 4/6 at purity 0.27 and total ploidy 2.10...
INFO [2019-04-13 04:57:08] Recalibrating log-ratios...
INFO [2019-04-13 04:57:08] Testing local optimum 4/6 at purity 0.27 and total ploidy 2.10...
INFO [2019-04-13 04:57:09] Recalibrating log-ratios...
INFO [2019-04-13 04:57:09] Testing local optimum 4/6 at purity 0.27 and total ploidy 2.10...
INFO [2019-04-13 04:57:09] Testing local optimum 5/6 at purity 0.29 and total ploidy 2.00...
INFO [2019-04-13 04:57:09] Recalibrating log-ratios...
INFO [2019-04-13 04:57:09] Testing local optimum 5/6 at purity 0.29 and total ploidy 2.00...
INFO [2019-04-13 04:57:10] Recalibrating log-ratios...
INFO [2019-04-13 04:57:10] Testing local optimum 5/6 at purity 0.29 and total ploidy 2.00...
INFO [2019-04-13 04:57:10] Recalibrating log-ratios...
INFO [2019-04-13 04:57:10] Testing local optimum 5/6 at purity 0.29 and total ploidy 2.00...
INFO [2019-04-13 04:57:11] Testing local optimum 6/6 at purity 0.78 and total ploidy 2.00...
INFO [2019-04-13 04:57:11] Recalibrating log-ratios...
INFO [2019-04-13 04:57:11] Testing local optimum 6/6 at purity 0.78 and total ploidy 2.00...
INFO [2019-04-13 04:57:11] Recalibrating log-ratios...
INFO [2019-04-13 04:57:11] Testing local optimum 6/6 at purity 0.78 and total ploidy 2.00...
INFO [2019-04-13 04:57:12] Recalibrating log-ratios...
INFO [2019-04-13 04:57:12] Testing local optimum 6/6 at purity 0.78 and total ploidy 2.00...
INFO [2019-04-13 04:57:12] Done.
INFO [2019-04-13 04:57:12] ------------------------------------------------------------
FATAL [2019-04-13 04:57:12] This function requires gene-level calls. Please add a column 'Gene'
FATAL [2019-04-13 04:57:12] containing gene symbols to the interval.file.
FATAL [2019-04-13 04:57:12]
FATAL [2019-04-13 04:57:12] This is most likely a user error due to invalid input data or
FATAL [2019-04-13 04:57:12] parameters (PureCN 1.12.2).
INFO [2019-04-13 04:57:13] ------------------------------------------------------------
INFO [2019-04-13 04:57:13] PureCN 1.12.2
INFO [2019-04-13 04:57:13] ------------------------------------------------------------
INFO [2019-04-13 04:57:13] Loading coverage files...
INFO [2019-04-13 04:57:14] Mean target coverages: 112X (tumor) 99X (normal).
WARN [2019-04-13 04:57:14] Allosome coverage missing, cannot determine sex.
WARN [2019-04-13 04:57:14] Allosome coverage missing, cannot determine sex.
INFO [2019-04-13 04:57:14] Removing 228 intervals with missing log.ratio.
INFO [2019-04-13 04:57:14] Removing 22 small (< 5bp) intervals.
INFO [2019-04-13 04:57:14] Removing 19 intervals with low total coverage in normal (< 150.00 reads).
WARN [2019-04-13 04:57:14] No normalDB provided. Provide one for better results.
INFO [2019-04-13 04:57:14] Removing 233 low coverage (< 15.0000X) intervals.
INFO [2019-04-13 04:57:14] Using 9547 intervals (9547 on-target, 0 off-target).
INFO [2019-04-13 04:57:14] No off-target intervals. If this is hybrid-capture data, consider adding them.
INFO [2019-04-13 04:57:14] No interval.file provided. Cannot check if data was GC-normalized. Was it?
INFO [2019-04-13 04:57:14] Loading VCF...
INFO [2019-04-13 04:57:14] Found 127 variants in VCF file.
INFO [2019-04-13 04:57:14] 127 (100.0%) variants annotated as likely germline (DB INFO flag).
FATAL [2019-04-13 04:57:14] Different chromosome names in coverage and VCF.
FATAL [2019-04-13 04:57:14]
FATAL [2019-04-13 04:57:14] This is most likely a user error due to invalid input data or
FATAL [2019-04-13 04:57:14] parameters (PureCN 1.12.2).
INFO [2019-04-13 04:57:14] ------------------------------------------------------------
INFO [2019-04-13 04:57:14] PureCN 1.12.2
INFO [2019-04-13 04:57:14] ------------------------------------------------------------
INFO [2019-04-13 04:57:14] Loading coverage files...
INFO [2019-04-13 04:57:14] Mean target coverages: 112X (tumor) 99X (normal).
WARN [2019-04-13 04:57:14] Allosome coverage missing, cannot determine sex.
WARN [2019-04-13 04:57:14] Allosome coverage missing, cannot determine sex.
FATAL [2019-04-13 04:57:15] tumor.coverage.file and interval.file do not align.
FATAL [2019-04-13 04:57:15]
FATAL [2019-04-13 04:57:15] This is most likely a user error due to invalid input data or
FATAL [2019-04-13 04:57:15] parameters (PureCN 1.12.2).
INFO [2019-04-13 04:57:15] ------------------------------------------------------------
INFO [2019-04-13 04:57:15] PureCN 1.12.2
INFO [2019-04-13 04:57:15] ------------------------------------------------------------
INFO [2019-04-13 04:57:15] Loading coverage files...
INFO [2019-04-13 04:57:15] Mean target coverages: 112X (tumor) 99X (normal).
WARN [2019-04-13 04:57:15] Allosome coverage missing, cannot determine sex.
WARN [2019-04-13 04:57:15] Allosome coverage missing, cannot determine sex.
INFO [2019-04-13 04:57:16] Removing 228 intervals with missing log.ratio.
INFO [2019-04-13 04:57:16] Removing 15 low/high GC targets.
INFO [2019-04-13 04:57:16] Removing 21 small (< 5bp) intervals.
INFO [2019-04-13 04:57:16] Removing 19 intervals with low total coverage in normal (< 150.00 reads).
WARN [2019-04-13 04:57:16] No normalDB provided. Provide one for better results.
INFO [2019-04-13 04:57:16] Removing 232 low coverage (< 15.0000X) intervals.
INFO [2019-04-13 04:57:16] Using 9534 intervals (9534 on-target, 0 off-target).
INFO [2019-04-13 04:57:16] No off-target intervals. If this is hybrid-capture data, consider adding them.
INFO [2019-04-13 04:57:16] AT/GC dropout: 1.03 (tumor), 1.03 (normal).
INFO [2019-04-13 04:57:16] Loading VCF...
INFO [2019-04-13 04:57:16] Found 127 variants in VCF file.
INFO [2019-04-13 04:57:16] 127 (100.0%) variants annotated as likely germline (DB INFO flag).
INFO [2019-04-13 04:57:17] LIB-02240e4 is tumor in VCF file.
INFO [2019-04-13 04:57:17] No homozygous variants in VCF, provide unfiltered VCF.
INFO [2019-04-13 04:57:17] Removing 1 non heterozygous (in matched normal) germline SNPs.
INFO [2019-04-13 04:57:17] Removing 1 variants with AF < 0.030 or AF >= 1.000 or less than 4 supporting reads or depth < 15.
INFO [2019-04-13 04:57:17] Removing 0 low quality variants with BQ < 25.
INFO [2019-04-13 04:57:17] Total size of targeted genomic region: 1.63Mb (2.58Mb with 50bp padding).
INFO [2019-04-13 04:57:17] 1.0% of targets contain variants.
INFO [2019-04-13 04:57:17] Removing 2 variants outside intervals.
INFO [2019-04-13 04:57:17] Found SOMATIC annotation in VCF.
INFO [2019-04-13 04:57:17] Setting somatic prior probabilities for somatic variants to 0.999000 or to 0.000100 otherwise.
INFO [2019-04-13 04:57:18] Found SOMATIC annotation in VCF. Setting mapping bias to 0.973.
INFO [2019-04-13 04:57:18] Sample sex: ?
INFO [2019-04-13 04:57:18] Segmenting data...
INFO [2019-04-13 04:57:18] Loading pre-computed boundaries for DNAcopy...
INFO [2019-04-13 04:57:18] Setting undo.SD parameter to 1.000000.
INFO [2019-04-13 04:57:19] Setting prune.hclust.h parameter to 0.200000.
INFO [2019-04-13 04:57:19] Found 54 segments with median size of 24.88Mb.
INFO [2019-04-13 04:57:19] Using 123 variants.
INFO [2019-04-13 04:57:19] Mean standard deviation of log-ratios: 0.37
INFO [2019-04-13 04:57:19] 2D-grid search of purity and ploidy...
INFO [2019-04-13 04:57:20] Local optima: 0.63/1.9, 0.4/2.1
INFO [2019-04-13 04:57:20] Testing local optimum 1/2 at purity 0.63 and total ploidy 1.90...
INFO [2019-04-13 04:57:22] Fitting variants for purity 0.65, tumor ploidy 1.74 and contamination 0.01.
INFO [2019-04-13 04:57:23] Optimized purity: 0.65
INFO [2019-04-13 04:57:23] Testing local optimum 2/2 at purity 0.40 and total ploidy 2.10...
INFO [2019-04-13 04:57:23] Recalibrating log-ratios...
INFO [2019-04-13 04:57:23] Testing local optimum 2/2 at purity 0.40 and total ploidy 2.10...
INFO [2019-04-13 04:57:24] Recalibrating log-ratios...
INFO [2019-04-13 04:57:24] Testing local optimum 2/2 at purity 0.40 and total ploidy 2.10...
INFO [2019-04-13 04:57:24] Recalibrating log-ratios...
INFO [2019-04-13 04:57:24] Testing local optimum 2/2 at purity 0.40 and total ploidy 2.10...
INFO [2019-04-13 04:57:25] Done.
INFO [2019-04-13 04:57:25] ------------------------------------------------------------
INFO [2019-04-13 04:57:25] ------------------------------------------------------------
INFO [2019-04-13 04:57:25] PureCN 1.12.2
INFO [2019-04-13 04:57:25] ------------------------------------------------------------
INFO [2019-04-13 04:57:25] Loading coverage files...
INFO [2019-04-13 04:57:25] Mean target coverages: 112X (tumor) 99X (normal).
WARN [2019-04-13 04:57:25] Allosome coverage missing, cannot determine sex.
WARN [2019-04-13 04:57:25] Allosome coverage missing, cannot determine sex.
INFO [2019-04-13 04:57:26] Removing 228 intervals with missing log.ratio.
INFO [2019-04-13 04:57:26] Removing 15 low/high GC targets.
INFO [2019-04-13 04:57:26] Removing 21 small (< 5bp) intervals.
INFO [2019-04-13 04:57:26] Removing 19 intervals with low total coverage in normal (< 150.00 reads).
WARN [2019-04-13 04:57:26] No normalDB provided. Provide one for better results.
INFO [2019-04-13 04:57:26] Removing 232 low coverage (< 15.0000X) intervals.
INFO [2019-04-13 04:57:26] Using 9534 intervals (9534 on-target, 0 off-target).
INFO [2019-04-13 04:57:26] No off-target intervals. If this is hybrid-capture data, consider adding them.
INFO [2019-04-13 04:57:26] AT/GC dropout: 1.03 (tumor), 1.03 (normal).
INFO [2019-04-13 04:57:26] Loading VCF...
INFO [2019-04-13 04:57:26] Found 127 variants in VCF file.
INFO [2019-04-13 04:57:26] 127 (100.0%) variants annotated as likely germline (DB INFO flag).
INFO [2019-04-13 04:57:26] LIB-02240e4 is tumor in VCF file.
INFO [2019-04-13 04:57:26] No homozygous variants in VCF, provide unfiltered VCF.
INFO [2019-04-13 04:57:26] Removing 1 non heterozygous (in matched normal) germline SNPs.
INFO [2019-04-13 04:57:27] Removing 1 variants with AF < 0.030 or AF >= 1.000 or less than 4 supporting reads or depth < 15.
INFO [2019-04-13 04:57:27] Removing 0 low quality variants with BQ < 25.
INFO [2019-04-13 04:57:27] Total size of targeted genomic region: 1.63Mb (2.58Mb with 50bp padding).
INFO [2019-04-13 04:57:27] 1.0% of targets contain variants.
INFO [2019-04-13 04:57:27] Removing 2 variants outside intervals.
INFO [2019-04-13 04:57:27] Found SOMATIC annotation in VCF.
INFO [2019-04-13 04:57:27] Setting somatic prior probabilities for somatic variants to 0.999000 or to 0.000100 otherwise.
INFO [2019-04-13 04:57:27] Found SOMATIC annotation in VCF. Setting mapping bias to 0.973.
INFO [2019-04-13 04:57:27] Sample sex: ?
INFO [2019-04-13 04:57:27] Segmenting data...
INFO [2019-04-13 04:57:27] Loading pre-computed boundaries for DNAcopy...
INFO [2019-04-13 04:57:27] Setting undo.SD parameter to 1.000000.
INFO [2019-04-13 04:57:29] Setting prune.hclust.h parameter to 0.200000.
INFO [2019-04-13 04:57:29] Found 54 segments with median size of 24.88Mb.
INFO [2019-04-13 04:57:29] Using 123 variants.
INFO [2019-04-13 04:57:29] Mean standard deviation of log-ratios: 0.37
INFO [2019-04-13 04:57:29] 2D-grid search of purity and ploidy...
INFO [2019-04-13 04:57:30] Local optima: 0.63/1.9, 0.4/2.1
INFO [2019-04-13 04:57:30] Testing local optimum 1/2 at purity 0.63 and total ploidy 1.90...
INFO [2019-04-13 04:57:31] Fitting variants for purity 0.65, tumor ploidy 1.74 and contamination 0.01.
INFO [2019-04-13 04:57:32] Optimized purity: 0.65
INFO [2019-04-13 04:57:32] Testing local optimum 2/2 at purity 0.40 and total ploidy 2.10...
INFO [2019-04-13 04:57:33] Recalibrating log-ratios...
INFO [2019-04-13 04:57:33] Testing local optimum 2/2 at purity 0.40 and total ploidy 2.10...
INFO [2019-04-13 04:57:33] Recalibrating log-ratios...
INFO [2019-04-13 04:57:33] Testing local optimum 2/2 at purity 0.40 and total ploidy 2.10...
INFO [2019-04-13 04:57:33] Recalibrating log-ratios...
INFO [2019-04-13 04:57:33] Testing local optimum 2/2 at purity 0.40 and total ploidy 2.10...
INFO [2019-04-13 04:57:34] Done.
INFO [2019-04-13 04:57:34] ------------------------------------------------------------
FATAL [2019-04-13 04:57:34] chr1 not valid chromosome name(s). Valid names are:
FATAL [2019-04-13 04:57:34] 1,2,3,4,5,6,7,8,9,10,11,12,13,14,15,16,17,18,19,20,21,22
FATAL [2019-04-13 04:57:34]
FATAL [2019-04-13 04:57:34] This is most likely a user error due to invalid input data or
FATAL [2019-04-13 04:57:34] parameters (PureCN 1.12.2).
INFO [2019-04-13 04:57:34] ------------------------------------------------------------
INFO [2019-04-13 04:57:34] PureCN 1.12.2
INFO [2019-04-13 04:57:34] ------------------------------------------------------------
INFO [2019-04-13 04:57:34] Loading coverage files...
WARN [2019-04-13 04:57:34] Provided sampleid (Sample2) does not match Sample1 found in segmentation.
WARN [2019-04-13 04:57:34] Allosome coverage missing, cannot determine sex.
WARN [2019-04-13 04:57:34] Allosome coverage missing, cannot determine sex.
INFO [2019-04-13 04:57:34] Removing 10 intervals with missing log.ratio.
INFO [2019-04-13 04:57:34] Using 10039 intervals (10039 on-target, 0 off-target).
INFO [2019-04-13 04:57:34] No off-target intervals. If this is hybrid-capture data, consider adding them.
INFO [2019-04-13 04:57:34] Loading VCF...
INFO [2019-04-13 04:57:34] Found 127 variants in VCF file.
INFO [2019-04-13 04:57:35] 127 (100.0%) variants annotated as likely germline (DB INFO flag).
INFO [2019-04-13 04:57:35] LIB-02240e4 is tumor in VCF file.
INFO [2019-04-13 04:57:35] No homozygous variants in VCF, provide unfiltered VCF.
INFO [2019-04-13 04:57:35] Removing 1 non heterozygous (in matched normal) germline SNPs.
INFO [2019-04-13 04:57:35] Removing 1 variants with AF < 0.030 or AF >= 1.000 or less than 4 supporting reads or depth < 15.
INFO [2019-04-13 04:57:35] Removing 0 low quality variants with BQ < 25.
INFO [2019-04-13 04:57:35] Total size of targeted genomic region: 1.71Mb (2.70Mb with 50bp padding).
INFO [2019-04-13 04:57:36] 1.0% of targets contain variants.
INFO [2019-04-13 04:57:36] Removing 0 variants outside intervals.
INFO [2019-04-13 04:57:36] Found SOMATIC annotation in VCF.
INFO [2019-04-13 04:57:36] Setting somatic prior probabilities for somatic variants to 0.999000 or to 0.000100 otherwise.
INFO [2019-04-13 04:57:36] Found SOMATIC annotation in VCF. Setting mapping bias to 0.976.
INFO [2019-04-13 04:57:36] Sample sex: ?
INFO [2019-04-13 04:57:36] Segmenting data...
WARN [2019-04-13 04:57:36] Provided sampleid (Sample2) does not match Sample1 found in segmentation.
INFO [2019-04-13 04:57:36] Loading pre-computed boundaries for DNAcopy...
INFO [2019-04-13 04:57:36] Setting undo.SD parameter to 0.000000.
Setting multi-figure configuration
INFO [2019-04-13 04:57:36] Setting prune.hclust.h parameter to 0.200000.
INFO [2019-04-13 04:57:36] Found 54 segments with median size of 24.88Mb.
INFO [2019-04-13 04:57:36] Using 125 variants.
INFO [2019-04-13 04:57:36] Mean standard deviation of log-ratios: 0.40
INFO [2019-04-13 04:57:36] 2D-grid search of purity and ploidy...
INFO [2019-04-13 04:57:37] Local optima: 0.6/1.9, 0.4/2.1
INFO [2019-04-13 04:57:37] Testing local optimum 1/2 at purity 0.60 and total ploidy 1.90...
INFO [2019-04-13 04:57:38] Fitting variants for purity 0.65, tumor ploidy 1.74 and contamination 0.01.
INFO [2019-04-13 04:57:40] Optimized purity: 0.65
INFO [2019-04-13 04:57:40] Testing local optimum 2/2 at purity 0.40 and total ploidy 2.10...
INFO [2019-04-13 04:57:41] Recalibrating log-ratios...
INFO [2019-04-13 04:57:41] Testing local optimum 2/2 at purity 0.40 and total ploidy 2.10...
INFO [2019-04-13 04:57:41] Recalibrating log-ratios...
INFO [2019-04-13 04:57:41] Testing local optimum 2/2 at purity 0.40 and total ploidy 2.10...
INFO [2019-04-13 04:57:41] Recalibrating log-ratios...
INFO [2019-04-13 04:57:41] Testing local optimum 2/2 at purity 0.40 and total ploidy 2.10...
INFO [2019-04-13 04:57:42] Done.
INFO [2019-04-13 04:57:42] ------------------------------------------------------------
WARN [2019-04-13 04:57:42] Allosome coverage missing, cannot determine sex.
WARN [2019-04-13 04:57:42] Allosome coverage missing, cannot determine sex.
Setting multi-figure configuration
FATAL [2019-04-13 04:57:50] Segmentation file expected with colnames ID, chrom, loc.start, loc.end,
FATAL [2019-04-13 04:57:50] num.mark, seg.mean
FATAL [2019-04-13 04:57:50]
FATAL [2019-04-13 04:57:50] This is most likely a user error due to invalid input data or
FATAL [2019-04-13 04:57:50] parameters (PureCN 1.12.2).
FATAL [2019-04-13 04:57:51] seg.file contains multiple samples and sampleid missing.
FATAL [2019-04-13 04:57:51]
FATAL [2019-04-13 04:57:51] This is most likely a user error due to invalid input data or
FATAL [2019-04-13 04:57:51] parameters (PureCN 1.12.2).
FATAL [2019-04-13 04:57:51] seg.file contains multiple samples and sampleid does not match any.
FATAL [2019-04-13 04:57:51]
FATAL [2019-04-13 04:57:51] This is most likely a user error due to invalid input data or
FATAL [2019-04-13 04:57:51] parameters (PureCN 1.12.2).
WARN [2019-04-13 04:57:51] Allosome coverage missing, cannot determine sex.
WARN [2019-04-13 04:57:51] Allosome coverage missing, cannot determine sex.
Setting multi-figure configuration
WARN [2019-04-13 04:58:00] Provided sampleid (example_tumor.txt) does not match Sample1 found in segmentation.
WARN [2019-04-13 04:58:00] Allosome coverage missing, cannot determine sex.
WARN [2019-04-13 04:58:00] Allosome coverage missing, cannot determine sex.
WARN [2019-04-13 04:58:00] Provided sampleid (example_tumor.txt) does not match Sample1 found in segmentation.
WARN [2019-04-13 04:58:05] Allosome coverage missing, cannot determine sex.
WARN [2019-04-13 04:58:05] Allosome coverage missing, cannot determine sex.
WARN [2019-04-13 04:58:06] No normalDB provided. Provide one for better results.
WARN [2019-04-13 04:58:14] Allosome coverage missing, cannot determine sex.
WARN [2019-04-13 04:58:14] Allosome coverage missing, cannot determine sex.
WARN [2019-04-13 04:58:16] Provided sampleid (Sample.1) does not match Sample1 found in segmentation.
Setting multi-figure configuration
WARN [2019-04-13 04:58:26] Allosome coverage missing, cannot determine sex.
WARN [2019-04-13 04:58:26] Allosome coverage missing, cannot determine sex.
WARN [2019-04-13 04:58:26] No normalDB provided. Provide one for better results.
WARN [2019-04-13 04:58:26] Sampleid looks like a normal in VCF, not like a tumor.
WARN [2019-04-13 04:58:38] Allosome coverage missing, cannot determine sex.
WARN [2019-04-13 04:58:38] Allosome coverage missing, cannot determine sex.
FATAL [2019-04-13 04:58:38] normalDB not a valid normalDB object. Use createNormalDatabase to
FATAL [2019-04-13 04:58:38] create one.
FATAL [2019-04-13 04:58:38]
FATAL [2019-04-13 04:58:38] This is most likely a user error due to invalid input data or
FATAL [2019-04-13 04:58:38] parameters (PureCN 1.12.2).
WARN [2019-04-13 04:58:39] You are likely not using the correct baits file!
WARN [2019-04-13 04:58:39] Allosome coverage missing, cannot determine sex.
WARN [2019-04-13 04:58:39] Allosome coverage missing, cannot determine sex.
WARN [2019-04-13 04:58:42] Allosome coverage missing, cannot determine sex.
WARN [2019-04-13 04:58:42] Allosome coverage missing, cannot determine sex.
FATAL [2019-04-13 04:58:42] normalDB appears to be empty.
FATAL [2019-04-13 04:58:42]
FATAL [2019-04-13 04:58:42] This is most likely a user error due to invalid input data or
FATAL [2019-04-13 04:58:42] parameters (PureCN 1.12.2).
WARN [2019-04-13 04:58:42] Allosome coverage missing, cannot determine sex.
WARN [2019-04-13 04:58:42] Allosome coverage missing, cannot determine sex.
WARN [2019-04-13 04:58:43] Intervals in coverage and interval.file have conflicting on/off-target annotation.
FATAL [2019-04-13 04:59:15] min.normals must be >=2.
FATAL [2019-04-13 04:59:15]
FATAL [2019-04-13 04:59:15] This is most likely a user error due to invalid input data or
FATAL [2019-04-13 04:59:15] parameters (PureCN 1.12.2).
FATAL [2019-04-13 04:59:27] The normal.panel.vcf.file contains only a single sample.
FATAL [2019-04-13 04:59:27]
FATAL [2019-04-13 04:59:27] This is most likely a user error due to invalid input data or
FATAL [2019-04-13 04:59:27] parameters (PureCN 1.12.2).
== testthat results ===========================================================
OK: 305 SKIPPED: 0 FAILED: 0
>
> proc.time()
user system elapsed
272.79 3.23 322.96
|
PureCN.Rcheck/tests_x64/testthat.Rout
R version 3.5.3 (2019-03-11) -- "Great Truth"
Copyright (C) 2019 The R Foundation for Statistical Computing
Platform: x86_64-w64-mingw32/x64 (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
Loading required package: parallel
Attaching package: 'BiocGenerics'
The following objects are masked from 'package:parallel':
clusterApply, clusterApplyLB, clusterCall, clusterEvalQ,
clusterExport, clusterMap, parApply, parCapply, parLapply,
parLapplyLB, parRapply, parSapply, parSapplyLB
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, append,
as.data.frame, basename, cbind, colMeans, colSums, colnames,
dirname, do.call, duplicated, eval, evalq, get, grep, grepl,
intersect, is.unsorted, lapply, lengths, mapply, match, mget,
order, paste, pmax, pmax.int, pmin, pmin.int, rank, rbind,
rowMeans, rowSums, rownames, sapply, setdiff, sort, table, tapply,
union, unique, unsplit, which, which.max, which.min
Loading required package: GenomeInfoDb
Loading required package: S4Vectors
Loading required package: stats4
Attaching package: 'S4Vectors'
The following object is masked from 'package:base':
expand.grid
Loading required package: IRanges
Attaching package: 'IRanges'
The following object is masked from 'package:grDevices':
windows
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")'.
Loading required package: DelayedArray
Loading required package: matrixStats
Attaching package: 'matrixStats'
The following objects are masked from 'package:Biobase':
anyMissing, rowMedians
Loading required package: BiocParallel
Attaching package: 'DelayedArray'
The following objects are masked from 'package:matrixStats':
colMaxs, colMins, colRanges, rowMaxs, rowMins, rowRanges
The following objects are masked from 'package:base':
aperm, apply
Loading required package: Rsamtools
Loading required package: Biostrings
Loading required package: XVector
Attaching package: 'Biostrings'
The following object is masked from 'package:DelayedArray':
type
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 [2019-04-13 04:59:47] Attempted adding gene symbols to intervals. Heuristics have been used to pick symbols for overlapping genes.
WARN [2019-04-13 04:59:47] Chromosome naming style of txdb file (UCSC) was different from interval file (NCBI).
WARN [2019-04-13 04:59:50] Attempted adding gene symbols to intervals. Heuristics have been used to pick symbols for overlapping genes.
FATAL [2019-04-13 04:59:52] tumor.coverage.file and interval.file do not align.
FATAL [2019-04-13 04:59:52]
FATAL [2019-04-13 04:59:52] This is most likely a user error due to invalid input data or
FATAL [2019-04-13 04:59:52] parameters (PureCN 1.12.2).
FATAL [2019-04-13 04:59:52] Need either f or purity and ploidy.
FATAL [2019-04-13 04:59:52]
FATAL [2019-04-13 04:59:52] This is most likely a user error due to invalid input data or
FATAL [2019-04-13 04:59:52] parameters (PureCN 1.12.2).
FATAL [2019-04-13 04:59:52] f not in expected range.
FATAL [2019-04-13 04:59:52]
FATAL [2019-04-13 04:59:52] This is most likely a user error due to invalid input data or
FATAL [2019-04-13 04:59:52] parameters (PureCN 1.12.2).
FATAL [2019-04-13 04:59:53] coverage not in expected range (>=2)
FATAL [2019-04-13 04:59:53]
FATAL [2019-04-13 04:59:53] This is most likely a user error due to invalid input data or
FATAL [2019-04-13 04:59:53] parameters (PureCN 1.12.2).
FATAL [2019-04-13 04:59:53] purity not in expected range.
FATAL [2019-04-13 04:59:53]
FATAL [2019-04-13 04:59:53] This is most likely a user error due to invalid input data or
FATAL [2019-04-13 04:59:53] parameters (PureCN 1.12.2).
FATAL [2019-04-13 04:59:53] ploidy not in expected range.
FATAL [2019-04-13 04:59:53]
FATAL [2019-04-13 04:59:53] This is most likely a user error due to invalid input data or
FATAL [2019-04-13 04:59:53] parameters (PureCN 1.12.2).
FATAL [2019-04-13 04:59:53] cell.fraction not in expected range.
FATAL [2019-04-13 04:59:53]
FATAL [2019-04-13 04:59:53] This is most likely a user error due to invalid input data or
FATAL [2019-04-13 04:59:53] parameters (PureCN 1.12.2).
INFO [2019-04-13 04:59:56] ------------------------------------------------------------
INFO [2019-04-13 04:59:56] PureCN 1.12.2
INFO [2019-04-13 04:59:56] ------------------------------------------------------------
INFO [2019-04-13 04:59:56] Loading coverage files...
INFO [2019-04-13 04:59:56] Mean target coverages: 112X (tumor) 99X (normal).
WARN [2019-04-13 04:59:56] Allosome coverage missing, cannot determine sex.
WARN [2019-04-13 04:59:56] Allosome coverage missing, cannot determine sex.
INFO [2019-04-13 04:59:56] Removing 228 intervals with missing log.ratio.
INFO [2019-04-13 04:59:56] Removing 22 small (< 5bp) intervals.
INFO [2019-04-13 04:59:56] Removing 19 intervals with low total coverage in normal (< 150.00 reads).
WARN [2019-04-13 04:59:56] No normalDB provided. Provide one for better results.
INFO [2019-04-13 04:59:56] Removing 233 low coverage (< 15.0000X) intervals.
INFO [2019-04-13 04:59:56] Using 9547 intervals (9547 on-target, 0 off-target).
INFO [2019-04-13 04:59:56] No off-target intervals. If this is hybrid-capture data, consider adding them.
INFO [2019-04-13 04:59:56] No interval.file provided. Cannot check if data was GC-normalized. Was it?
INFO [2019-04-13 04:59:56] Loading VCF...
INFO [2019-04-13 04:59:56] Found 127 variants in VCF file.
INFO [2019-04-13 04:59:57] 127 (100.0%) variants annotated as likely germline (DB INFO flag).
INFO [2019-04-13 04:59:57] LIB-02240e4 is tumor in VCF file.
INFO [2019-04-13 04:59:57] No homozygous variants in VCF, provide unfiltered VCF.
INFO [2019-04-13 04:59:57] Removing 1 non heterozygous (in matched normal) germline SNPs.
INFO [2019-04-13 04:59:58] Removing 1 variants with AF < 0.030 or AF >= 1.000 or less than 4 supporting reads or depth < 15.
INFO [2019-04-13 04:59:58] Removing 0 low quality variants with BQ < 25.
INFO [2019-04-13 04:59:58] Total size of targeted genomic region: 1.63Mb (2.58Mb with 50bp padding).
INFO [2019-04-13 04:59:58] 1.0% of targets contain variants.
INFO [2019-04-13 04:59:58] Removing 2 variants outside intervals.
INFO [2019-04-13 04:59:58] Found SOMATIC annotation in VCF.
INFO [2019-04-13 04:59:58] Setting somatic prior probabilities for somatic variants to 0.999000 or to 0.000100 otherwise.
INFO [2019-04-13 04:59:58] Found SOMATIC annotation in VCF. Setting mapping bias to 0.973.
INFO [2019-04-13 04:59:58] Sample sex: ?
INFO [2019-04-13 04:59:58] Segmenting data...
INFO [2019-04-13 04:59:58] Loading pre-computed boundaries for DNAcopy...
INFO [2019-04-13 04:59:58] Setting undo.SD parameter to 1.000000.
INFO [2019-04-13 04:59:59] Setting prune.hclust.h parameter to 0.200000.
INFO [2019-04-13 05:00:00] Found 54 segments with median size of 24.88Mb.
INFO [2019-04-13 05:00:00] Using 123 variants.
INFO [2019-04-13 05:00:00] Mean standard deviation of log-ratios: 0.37
INFO [2019-04-13 05:00:00] 2D-grid search of purity and ploidy...
INFO [2019-04-13 05:00:02] Local optima: 0.63/1.9, 0.4/2.2, 0.5/2
INFO [2019-04-13 05:00:02] Testing local optimum 1/3 at purity 0.63 and total ploidy 1.90...
INFO [2019-04-13 05:00:03] Fitting variants for purity 0.65, tumor ploidy 1.74 and contamination 0.01.
INFO [2019-04-13 05:00:05] Optimized purity: 0.65
INFO [2019-04-13 05:00:05] Testing local optimum 2/3 at purity 0.40 and total ploidy 2.20...
INFO [2019-04-13 05:00:06] Fitting variants for purity 0.40, tumor ploidy 2.34 and contamination 0.01.
INFO [2019-04-13 05:00:07] Optimized purity: 0.40
INFO [2019-04-13 05:00:07] Testing local optimum 3/3 at purity 0.50 and total ploidy 2.00...
INFO [2019-04-13 05:00:08] Fitting variants for purity 0.65, tumor ploidy 1.74 and contamination 0.01.
INFO [2019-04-13 05:00:09] Optimized purity: 0.65
INFO [2019-04-13 05:00:09] Done.
INFO [2019-04-13 05:00:09] ------------------------------------------------------------
FATAL [2019-04-13 05:00:11] exclude not a GRanges object.
FATAL [2019-04-13 05:00:11]
FATAL [2019-04-13 05:00:11] This is most likely a user error due to invalid input data or
FATAL [2019-04-13 05:00:11] parameters (PureCN 1.12.2).
FATAL [2019-04-13 05:00:11] callable not a GRanges object.
FATAL [2019-04-13 05:00:11]
FATAL [2019-04-13 05:00:11] This is most likely a user error due to invalid input data or
FATAL [2019-04-13 05:00:11] parameters (PureCN 1.12.2).
WARN [2019-04-13 05:00:23] tumor.coverage.file and interval.file do not align.
FATAL [2019-04-13 05:00:25] No gc_bias column in interval.file.
FATAL [2019-04-13 05:00:25]
FATAL [2019-04-13 05:00:25] This is most likely a user error due to invalid input data or
FATAL [2019-04-13 05:00:25] parameters (PureCN 1.12.2).
FATAL [2019-04-13 05:00:25] Provided coverage is zero, most likely due to a corrupt BAM file.
FATAL [2019-04-13 05:00:25]
FATAL [2019-04-13 05:00:25] This is most likely a user error due to invalid input data or
FATAL [2019-04-13 05:00:25] parameters (PureCN 1.12.2).
INFO [2019-04-13 05:00:32] Reading C:\Users\biocbuild\bbs-3.8-bioc\tmpdir\RtmpSuRWdk\file2b603d941ada.rds...
INFO [2019-04-13 05:00:32] Reading C:\Users\biocbuild\bbs-3.8-bioc\tmpdir\RtmpSuRWdk\file2b603d941ada.rds...
INFO [2019-04-13 05:00:32] Reading C:\Users\biocbuild\bbs-3.8-bioc\tmpdir\RtmpSuRWdk\file2b603d941ada.rds...
INFO [2019-04-13 05:00:32] Reading C:\Users\biocbuild\bbs-3.8-bioc\tmpdir\RtmpSuRWdk\file2b603d941ada.rds...
INFO [2019-04-13 05:00:33] Reading C:\Users\biocbuild\bbs-3.8-bioc\tmpdir\RtmpSuRWdk\file2b603d941ada.rds...
INFO [2019-04-13 05:00:33] Reading C:\Users\biocbuild\bbs-3.8-bioc\tmpdir\RtmpSuRWdk\file2b603d941ada.rds...
FATAL [2019-04-13 05:00:33] Purity or Ploidy not numeric or in expected range.
FATAL [2019-04-13 05:00:33]
FATAL [2019-04-13 05:00:33] This is most likely a user error due to invalid input data or
FATAL [2019-04-13 05:00:33] parameters (PureCN 1.12.2).
INFO [2019-04-13 05:00:33] Reading C:\Users\biocbuild\bbs-3.8-bioc\tmpdir\RtmpSuRWdk\file2b603d941ada.rds...
INFO [2019-04-13 05:00:33] Reading C:\Users\biocbuild\bbs-3.8-bioc\tmpdir\RtmpSuRWdk\file2b603d941ada.rds...
FATAL [2019-04-13 05:00:33] 'Failed' column in
FATAL [2019-04-13 05:00:33] C:\Users\biocbuild\bbs-3.8-bioc\tmpdir\RtmpSuRWdk\file2b603d941ada.csv
FATAL [2019-04-13 05:00:33] not logical(1).
FATAL [2019-04-13 05:00:33]
FATAL [2019-04-13 05:00:33] This is most likely a user error due to invalid input data or
FATAL [2019-04-13 05:00:33] parameters (PureCN 1.12.2).
INFO [2019-04-13 05:00:34] 576 on-target bins with low coverage in all samples.
WARN [2019-04-13 05:00:34] You are likely not using the correct baits file!
WARN [2019-04-13 05:00:34] Allosome coverage missing, cannot determine sex.
WARN [2019-04-13 05:00:34] Allosome coverage missing, cannot determine sex.
INFO [2019-04-13 05:00:34] Processing on-target regions...
INFO [2019-04-13 05:00:35] Removing 930 intervals with low coverage in normalDB.
INFO [2019-04-13 05:00:35] Removing 1 intervals with zero coverage in more than 3% of normalDB.
INFO [2019-04-13 05:00:38] 576 on-target bins with low coverage in all samples.
WARN [2019-04-13 05:00:38] You are likely not using the correct baits file!
WARN [2019-04-13 05:00:38] Allosome coverage missing, cannot determine sex.
WARN [2019-04-13 05:00:38] Allosome coverage missing, cannot determine sex.
INFO [2019-04-13 05:00:38] Processing on-target regions...
INFO [2019-04-13 05:00:38] Removing 930 intervals with low coverage in normalDB.
INFO [2019-04-13 05:00:38] Removing 1 intervals with zero coverage in more than 3% of normalDB.
INFO [2019-04-13 05:00:40] 576 on-target bins with low coverage in all samples.
WARN [2019-04-13 05:00:40] You are likely not using the correct baits file!
WARN [2019-04-13 05:00:40] Allosome coverage missing, cannot determine sex.
WARN [2019-04-13 05:00:40] Allosome coverage missing, cannot determine sex.
INFO [2019-04-13 05:00:41] Processing on-target regions...
INFO [2019-04-13 05:00:41] Removing 930 intervals with low coverage in normalDB.
INFO [2019-04-13 05:00:41] Removing 1 intervals with zero coverage in more than 3% of normalDB.
INFO [2019-04-13 05:00:43] 576 on-target bins with low coverage in all samples.
WARN [2019-04-13 05:00:43] You are likely not using the correct baits file!
WARN [2019-04-13 05:00:43] Allosome coverage missing, cannot determine sex.
WARN [2019-04-13 05:00:43] Allosome coverage missing, cannot determine sex.
FATAL [2019-04-13 05:00:43] Length of normal.coverage.files and sex different
FATAL [2019-04-13 05:00:43]
FATAL [2019-04-13 05:00:43] This is most likely a user error due to invalid input data or
FATAL [2019-04-13 05:00:43] parameters (PureCN 1.12.2).
WARN [2019-04-13 05:00:47] Target intervals were not sorted.
INFO [2019-04-13 05:00:48] 560 on-target bins with low coverage in all samples.
WARN [2019-04-13 05:00:48] You are likely not using the correct baits file!
WARN [2019-04-13 05:00:48] Allosome coverage missing, cannot determine sex.
WARN [2019-04-13 05:00:48] Allosome coverage missing, cannot determine sex.
WARN [2019-04-13 05:00:48] Allosome coverage missing, cannot determine sex.
INFO [2019-04-13 05:00:48] Processing on-target regions...
INFO [2019-04-13 05:00:48] Removing 978 intervals with low coverage in normalDB.
INFO [2019-04-13 05:00:48] Removing 11 intervals with zero coverage in more than 3% of normalDB.
FATAL [2019-04-13 05:00:51] tumor.coverage.file and normalDB do not align.
FATAL [2019-04-13 05:00:51]
FATAL [2019-04-13 05:00:51] This is most likely a user error due to invalid input data or
FATAL [2019-04-13 05:00:51] parameters (PureCN 1.12.2).
INFO [2019-04-13 05:00:51] ------------------------------------------------------------
INFO [2019-04-13 05:00:51] PureCN 1.12.2
INFO [2019-04-13 05:00:51] ------------------------------------------------------------
INFO [2019-04-13 05:00:51] Loading coverage files...
INFO [2019-04-13 05:00:52] Mean target coverages: 112X (tumor) 99X (normal).
WARN [2019-04-13 05:00:52] Allosome coverage missing, cannot determine sex.
WARN [2019-04-13 05:00:52] Allosome coverage missing, cannot determine sex.
INFO [2019-04-13 05:00:52] Removing 228 intervals with missing log.ratio.
FATAL [2019-04-13 05:00:52] normalDB incompatible with this PureCN version. Please re-run
FATAL [2019-04-13 05:00:52] NormalDB.R.
FATAL [2019-04-13 05:00:52]
FATAL [2019-04-13 05:00:52] This is most likely a user error due to invalid input data or
FATAL [2019-04-13 05:00:52] parameters (PureCN 1.12.2).
FATAL [2019-04-13 05:00:52] normalDB incompatible with this PureCN version. Please re-run
FATAL [2019-04-13 05:00:52] NormalDB.R.
FATAL [2019-04-13 05:00:52]
FATAL [2019-04-13 05:00:52] This is most likely a user error due to invalid input data or
FATAL [2019-04-13 05:00:52] parameters (PureCN 1.12.2).
INFO [2019-04-13 05:00:53] Removing 22 non heterozygous (in matched normal) germline SNPs.
INFO [2019-04-13 05:00:54] Initial testing for significant sample cross-contamination: unlikely
INFO [2019-04-13 05:00:54] Removing 48 variants with AF < 0.030 or AF >= 1.000 or less than 4 supporting reads or depth < 15.
INFO [2019-04-13 05:00:54] Removing 16 low quality variants with BQ < 25.
INFO [2019-04-13 05:00:54] Removing 22 non heterozygous (in matched normal) germline SNPs.
INFO [2019-04-13 05:00:55] Initial testing for significant sample cross-contamination: unlikely
INFO [2019-04-13 05:00:56] Removing 48 variants with AF < 0.030 or AF >= 1.000 or less than 4 supporting reads or depth < 15.
INFO [2019-04-13 05:00:56] Removing 6 blacklisted variants.
INFO [2019-04-13 05:00:56] Removing 16 low quality variants with BQ < 25.
INFO [2019-04-13 05:00:57] Removing 22 non heterozygous (in matched normal) germline SNPs.
INFO [2019-04-13 05:00:58] Initial testing for significant sample cross-contamination: unlikely
INFO [2019-04-13 05:00:58] Removing 48 variants with AF < 0.030 or AF >= 1.000 or less than 4 supporting reads or depth < 15.
INFO [2019-04-13 05:00:58] Removing 16 low quality variants with BQ < 25.
WARN [2019-04-13 05:00:58] MuTect stats file lacks contig and position columns.
INFO [2019-04-13 05:00:59] Removing 22 non heterozygous (in matched normal) germline SNPs.
INFO [2019-04-13 05:01:00] Initial testing for significant sample cross-contamination: unlikely
INFO [2019-04-13 05:01:00] Removing 48 variants with AF < 0.030 or AF >= 1.000 or less than 4 supporting reads or depth < 15.
INFO [2019-04-13 05:01:00] Removing 16 low quality variants with BQ < 25.
WARN [2019-04-13 05:01:00] MuTect stats file and VCF file do not align perfectly. Will remove 2325 unmatched variants.
WARN [2019-04-13 05:01:00] MuTect stats file lacks failure_reasons column. Keeping all variants listed in stats file.
INFO [2019-04-13 05:01:00] Removing 0 non heterozygous (in matched normal) germline SNPs.
INFO [2019-04-13 05:01:01] Removing 0 variants with AF < 0.030 or AF >= 1.000 or less than 4 supporting reads or depth < 15.
INFO [2019-04-13 05:01:01] Removing 0 low quality variants with BQ < 25.
INFO [2019-04-13 05:01:02] Removing 22 non heterozygous (in matched normal) germline SNPs.
INFO [2019-04-13 05:01:02] Initial testing for significant sample cross-contamination: unlikely
INFO [2019-04-13 05:01:02] Removing 48 variants with AF < 0.030 or AF >= 1.000 or less than 4 supporting reads or depth < 15.
INFO [2019-04-13 05:01:03] Removing 22 non heterozygous (in matched normal) germline SNPs.
INFO [2019-04-13 05:01:03] Initial testing for significant sample cross-contamination: unlikely
INFO [2019-04-13 05:01:03] Removing 48 variants with AF < 0.030 or AF >= 1.000 or less than 4 supporting reads or depth < 15.
INFO [2019-04-13 05:01:04] Found 11 variants in VCF file.
WARN [2019-04-13 05:01:04] vcf.file has no DB info field for membership in germline databases. Found and used valid population allele frequency > 0.001000 instead.
INFO [2019-04-13 05:01:05] 1 (9.1%) variants annotated as likely germline (DB INFO flag).
INFO [2019-04-13 05:01:09] Found 11 variants in VCF file.
WARN [2019-04-13 05:01:10] vcf.file has no DB info field for membership in germline databases. Found and used valid population allele frequency > 0.001000 instead.
INFO [2019-04-13 05:01:10] 1 (9.1%) variants annotated as likely germline (DB INFO flag).
INFO [2019-04-13 05:01:11] Found 11 variants in VCF file.
WARN [2019-04-13 05:01:11] vcf.file has no DB or POP_AF info field for membership in germline databases. Guessing it based on available dbSNP ID.
INFO [2019-04-13 05:01:11] 1 (9.1%) variants annotated as likely germline (DB INFO flag).
WARN [2019-04-13 05:01:12] Allosome coverage missing, cannot determine sex.
WARN [2019-04-13 05:01:13] Allosome coverage missing, cannot determine sex.
INFO [2019-04-13 05:01:13] Mean coverages: chrX: 108.58, chrY: 108.58, chr1-22: 104.87.
INFO [2019-04-13 05:01:13] Mean coverages: chrX: 108.58, chrY: 2.17, chr1-22: 104.87.
INFO [2019-04-13 05:01:13] Mean coverages: chrX: 108.58, chrY: 5.17, chr1-22: 104.87.
INFO [2019-04-13 05:01:13] No homozygous variants in VCF, provide unfiltered VCF.
INFO [2019-04-13 05:01:13] No germline variants in VCF.
FATAL [2019-04-13 05:01:13] No solution with id hello
FATAL [2019-04-13 05:01:13]
FATAL [2019-04-13 05:01:13] This is most likely a user error due to invalid input data or
FATAL [2019-04-13 05:01:13] parameters (PureCN 1.12.2).
FATAL [2019-04-13 05:01:13] No solution with id 10
FATAL [2019-04-13 05:01:13]
FATAL [2019-04-13 05:01:13] This is most likely a user error due to invalid input data or
FATAL [2019-04-13 05:01:13] parameters (PureCN 1.12.2).
FATAL [2019-04-13 05:01:15] all.data and w have different lengths.
FATAL [2019-04-13 05:01:15]
FATAL [2019-04-13 05:01:15] This is most likely a user error due to invalid input data or
FATAL [2019-04-13 05:01:15] parameters (PureCN 1.12.2).
WARN [2019-04-13 05:01:16] No mappability scores provided.
WARN [2019-04-13 05:01:16] No reptiming scores provided.
INFO [2019-04-13 05:01:16] Calculating GC-content...
WARN [2019-04-13 05:01:16] No mappability scores provided.
WARN [2019-04-13 05:01:16] No reptiming scores provided.
INFO [2019-04-13 05:01:16] Calculating GC-content...
INFO [2019-04-13 05:01:17] Removing 1 targets overlapping with exclude.
WARN [2019-04-13 05:01:17] No mappability scores provided.
WARN [2019-04-13 05:01:17] No reptiming scores provided.
INFO [2019-04-13 05:01:17] Calculating GC-content...
FATAL [2019-04-13 05:01:17] off.target.padding must be negative.
FATAL [2019-04-13 05:01:17]
FATAL [2019-04-13 05:01:17] This is most likely a user error due to invalid input data or
FATAL [2019-04-13 05:01:17] parameters (PureCN 1.12.2).
FATAL [2019-04-13 05:01:17] Interval coordinates should start at 1, not at 0
FATAL [2019-04-13 05:01:17]
FATAL [2019-04-13 05:01:17] This is most likely a user error due to invalid input data or
FATAL [2019-04-13 05:01:17] parameters (PureCN 1.12.2).
FATAL [2019-04-13 05:01:18] No off-target regions after filtering for mappability and
FATAL [2019-04-13 05:01:18] off.target.padding
FATAL [2019-04-13 05:01:18]
FATAL [2019-04-13 05:01:18] This is most likely a user error due to invalid input data or
FATAL [2019-04-13 05:01:18] parameters (PureCN 1.12.2).
WARN [2019-04-13 05:01:18] No mappability scores provided.
INFO [2019-04-13 05:01:18] Calculating GC-content...
INFO [2019-04-13 05:01:18] Averaging reptiming into bins of size 200...
WARN [2019-04-13 05:01:18] No mappability scores provided.
INFO [2019-04-13 05:01:18] Calculating GC-content...
INFO [2019-04-13 05:01:19] Splitting 1 large targets to an average width of 200.
WARN [2019-04-13 05:01:19] No mappability scores provided.
WARN [2019-04-13 05:01:19] No reptiming scores provided.
INFO [2019-04-13 05:01:19] Calculating GC-content...
WARN [2019-04-13 05:01:19] No mappability scores provided.
WARN [2019-04-13 05:01:19] No reptiming scores provided.
INFO [2019-04-13 05:01:19] Calculating GC-content...
INFO [2019-04-13 05:01:19] Tiling off-target regions to an average width of 200000.
WARN [2019-04-13 05:01:19] No mappability scores provided.
WARN [2019-04-13 05:01:19] No reptiming scores provided.
INFO [2019-04-13 05:01:19] Calculating GC-content...
WARN [2019-04-13 05:01:19] Intervals contain off-target regions. Will not change intervals.
WARN [2019-04-13 05:01:19] No mappability scores provided.
WARN [2019-04-13 05:01:19] No reptiming scores provided.
INFO [2019-04-13 05:01:19] Calculating GC-content...
WARN [2019-04-13 05:01:20] No reptiming scores provided.
INFO [2019-04-13 05:01:20] Calculating GC-content...
WARN [2019-04-13 05:01:20] 1 intervals without mappability score (1 on-target).
INFO [2019-04-13 05:01:21] Removing 1 intervals with low mappability score (<0.50).
WARN [2019-04-13 05:01:21] No reptiming scores provided.
INFO [2019-04-13 05:01:21] Calculating GC-content...
INFO [2019-04-13 05:01:21] Tiling off-target regions to an average width of 200000.
WARN [2019-04-13 05:01:21] No reptiming scores provided.
INFO [2019-04-13 05:01:21] Calculating GC-content...
WARN [2019-04-13 05:01:22] No mappability scores provided.
WARN [2019-04-13 05:01:22] No reptiming scores provided.
INFO [2019-04-13 05:01:22] Calculating GC-content...
WARN [2019-04-13 05:01:22] Chromosome naming style of interval file (NCBI) was different from reference (UCSC).
WARN [2019-04-13 05:01:22] No mappability scores provided.
WARN [2019-04-13 05:01:22] No reptiming scores provided.
INFO [2019-04-13 05:01:22] Calculating GC-content...
FATAL [2019-04-13 05:01:22] Chromosome naming style of interval file unknown, should be UCSC.
FATAL [2019-04-13 05:01:22]
FATAL [2019-04-13 05:01:22] This is most likely a user error due to invalid input data or
FATAL [2019-04-13 05:01:22] parameters (PureCN 1.12.2).
WARN [2019-04-13 05:01:22] Chromosome naming style of interval file (NCBI) was different from reference (UCSC).
WARN [2019-04-13 05:01:24] Chromosome naming style of mappability file (NCBI) was different from reference (UCSC).
WARN [2019-04-13 05:01:24] No reptiming scores provided.
INFO [2019-04-13 05:01:24] Calculating GC-content...
INFO [2019-04-13 05:01:25] Loading coverage data...
INFO [2019-04-13 05:01:26] Mean target coverages: 71X (tumor) 99X (normal).
INFO [2019-04-13 05:01:26] Mean target coverages: 71X (tumor) 43X (normal).
INFO [2019-04-13 05:01:27] ------------------------------------------------------------
INFO [2019-04-13 05:01:27] PureCN 1.12.2
INFO [2019-04-13 05:01:27] ------------------------------------------------------------
INFO [2019-04-13 05:01:27] Using BiocParallel for parallel optimization.
INFO [2019-04-13 05:01:27] Loading coverage files...
INFO [2019-04-13 05:01:27] Mean target coverages: 112X (tumor) 99X (normal).
WARN [2019-04-13 05:01:27] Allosome coverage missing, cannot determine sex.
WARN [2019-04-13 05:01:27] Allosome coverage missing, cannot determine sex.
INFO [2019-04-13 05:01:27] Removing 228 intervals with missing log.ratio.
INFO [2019-04-13 05:01:27] Removing 22 small (< 5bp) intervals.
INFO [2019-04-13 05:01:27] Removing 19 intervals with low total coverage in normal (< 150.00 reads).
WARN [2019-04-13 05:01:27] No normalDB provided. Provide one for better results.
INFO [2019-04-13 05:01:27] Removing 233 low coverage (< 15.0000X) intervals.
INFO [2019-04-13 05:01:27] Using 9547 intervals (9547 on-target, 0 off-target).
INFO [2019-04-13 05:01:27] No off-target intervals. If this is hybrid-capture data, consider adding them.
INFO [2019-04-13 05:01:27] No interval.file provided. Cannot check if data was GC-normalized. Was it?
INFO [2019-04-13 05:01:27] Sample sex: ?
INFO [2019-04-13 05:01:27] Segmenting data...
INFO [2019-04-13 05:01:27] Interval weights found, will use weighted CBS.
INFO [2019-04-13 05:01:27] Loading pre-computed boundaries for DNAcopy...
INFO [2019-04-13 05:01:28] Setting undo.SD parameter to 1.000000.
Setting multi-figure configuration
INFO [2019-04-13 05:01:29] Found 54 segments with median size of 25.38Mb.
INFO [2019-04-13 05:01:29] Mean standard deviation of log-ratios: 0.37
INFO [2019-04-13 05:01:29] 2D-grid search of purity and ploidy...
INFO [2019-04-13 05:01:29] Local optima: 0.65/1.9, 0.33/1.9, 0.38/2.2, 0.27/2.1, 0.29/2, 0.78/2
INFO [2019-04-13 05:02:24] Done.
INFO [2019-04-13 05:02:24] ------------------------------------------------------------
INFO [2019-04-13 05:02:24] Reading C:\Users\biocbuild\bbs-3.8-bioc\tmpdir\RtmpSuRWdk\file2b6076101ccc.rds...
FATAL [2019-04-13 05:02:24] runAbsoluteCN was run without a VCF file.
FATAL [2019-04-13 05:02:24]
FATAL [2019-04-13 05:02:24] This is most likely a user error due to invalid input data or
FATAL [2019-04-13 05:02:24] parameters (PureCN 1.12.2).
FATAL [2019-04-13 05:02:24] runAbsoluteCN was run without a VCF file.
FATAL [2019-04-13 05:02:24]
FATAL [2019-04-13 05:02:24] This is most likely a user error due to invalid input data or
FATAL [2019-04-13 05:02:24] parameters (PureCN 1.12.2).
INFO [2019-04-13 05:02:24] ------------------------------------------------------------
INFO [2019-04-13 05:02:24] PureCN 1.12.2
INFO [2019-04-13 05:02:24] ------------------------------------------------------------
INFO [2019-04-13 05:02:24] Loading coverage files...
FATAL [2019-04-13 05:02:25] Need a normal coverage file if log.ratio and seg.file are not provided.
FATAL [2019-04-13 05:02:25]
FATAL [2019-04-13 05:02:25] This is most likely a user error due to invalid input data or
FATAL [2019-04-13 05:02:25] parameters (PureCN 1.12.2).
INFO [2019-04-13 05:02:25] ------------------------------------------------------------
INFO [2019-04-13 05:02:25] PureCN 1.12.2
INFO [2019-04-13 05:02:25] ------------------------------------------------------------
FATAL [2019-04-13 05:02:25] min.ploidy or max.ploidy not within expected range.
FATAL [2019-04-13 05:02:25]
FATAL [2019-04-13 05:02:25] This is most likely a user error due to invalid input data or
FATAL [2019-04-13 05:02:25] parameters (PureCN 1.12.2).
INFO [2019-04-13 05:02:25] ------------------------------------------------------------
INFO [2019-04-13 05:02:25] PureCN 1.12.2
INFO [2019-04-13 05:02:25] ------------------------------------------------------------
FATAL [2019-04-13 05:02:25] min.ploidy or max.ploidy not within expected range.
FATAL [2019-04-13 05:02:25]
FATAL [2019-04-13 05:02:25] This is most likely a user error due to invalid input data or
FATAL [2019-04-13 05:02:25] parameters (PureCN 1.12.2).
INFO [2019-04-13 05:02:25] ------------------------------------------------------------
INFO [2019-04-13 05:02:25] PureCN 1.12.2
INFO [2019-04-13 05:02:25] ------------------------------------------------------------
INFO [2019-04-13 05:02:25] ------------------------------------------------------------
INFO [2019-04-13 05:02:25] PureCN 1.12.2
INFO [2019-04-13 05:02:25] ------------------------------------------------------------
INFO [2019-04-13 05:02:25] ------------------------------------------------------------
INFO [2019-04-13 05:02:25] PureCN 1.12.2
INFO [2019-04-13 05:02:25] ------------------------------------------------------------
FATAL [2019-04-13 05:02:25] test.num.copy not within expected range.
FATAL [2019-04-13 05:02:25]
FATAL [2019-04-13 05:02:25] This is most likely a user error due to invalid input data or
FATAL [2019-04-13 05:02:25] parameters (PureCN 1.12.2).
INFO [2019-04-13 05:02:25] ------------------------------------------------------------
INFO [2019-04-13 05:02:25] PureCN 1.12.2
INFO [2019-04-13 05:02:25] ------------------------------------------------------------
WARN [2019-04-13 05:02:25] test.num.copy outside recommended range.
FATAL [2019-04-13 05:02:25] max.non.clonal not within expected range or format.
FATAL [2019-04-13 05:02:25]
FATAL [2019-04-13 05:02:25] This is most likely a user error due to invalid input data or
FATAL [2019-04-13 05:02:25] parameters (PureCN 1.12.2).
INFO [2019-04-13 05:02:25] ------------------------------------------------------------
INFO [2019-04-13 05:02:25] PureCN 1.12.2
INFO [2019-04-13 05:02:25] ------------------------------------------------------------
WARN [2019-04-13 05:02:25] test.num.copy outside recommended range.
FATAL [2019-04-13 05:02:25] max.non.clonal not within expected range or format.
FATAL [2019-04-13 05:02:25]
FATAL [2019-04-13 05:02:25] This is most likely a user error due to invalid input data or
FATAL [2019-04-13 05:02:25] parameters (PureCN 1.12.2).
INFO [2019-04-13 05:02:25] ------------------------------------------------------------
INFO [2019-04-13 05:02:25] PureCN 1.12.2
INFO [2019-04-13 05:02:25] ------------------------------------------------------------
FATAL [2019-04-13 05:02:25] test.purity not within expected range.
FATAL [2019-04-13 05:02:25]
FATAL [2019-04-13 05:02:25] This is most likely a user error due to invalid input data or
FATAL [2019-04-13 05:02:25] parameters (PureCN 1.12.2).
INFO [2019-04-13 05:02:25] ------------------------------------------------------------
INFO [2019-04-13 05:02:25] PureCN 1.12.2
INFO [2019-04-13 05:02:25] ------------------------------------------------------------
INFO [2019-04-13 05:02:25] Loading coverage files...
FATAL [2019-04-13 05:02:26] Tumor and normal are identical. This won't give any meaningful results
FATAL [2019-04-13 05:02:26] and I'm stopping here.
FATAL [2019-04-13 05:02:26]
FATAL [2019-04-13 05:02:26] This is most likely a user error due to invalid input data or
FATAL [2019-04-13 05:02:26] parameters (PureCN 1.12.2).
INFO [2019-04-13 05:02:26] ------------------------------------------------------------
INFO [2019-04-13 05:02:26] PureCN 1.12.2
INFO [2019-04-13 05:02:26] ------------------------------------------------------------
INFO [2019-04-13 05:02:26] Loading coverage files...
FATAL [2019-04-13 05:02:26] Length of log.ratio different from tumor coverage.
FATAL [2019-04-13 05:02:26]
FATAL [2019-04-13 05:02:26] This is most likely a user error due to invalid input data or
FATAL [2019-04-13 05:02:26] parameters (PureCN 1.12.2).
INFO [2019-04-13 05:02:26] ------------------------------------------------------------
INFO [2019-04-13 05:02:26] PureCN 1.12.2
INFO [2019-04-13 05:02:26] ------------------------------------------------------------
FATAL [2019-04-13 05:02:26] prior.purity must have the same length as test.purity.
FATAL [2019-04-13 05:02:26]
FATAL [2019-04-13 05:02:26] This is most likely a user error due to invalid input data or
FATAL [2019-04-13 05:02:26] parameters (PureCN 1.12.2).
INFO [2019-04-13 05:02:26] ------------------------------------------------------------
INFO [2019-04-13 05:02:26] PureCN 1.12.2
INFO [2019-04-13 05:02:26] ------------------------------------------------------------
FATAL [2019-04-13 05:02:26] min.gof not within expected range or format.
FATAL [2019-04-13 05:02:26]
FATAL [2019-04-13 05:02:26] This is most likely a user error due to invalid input data or
FATAL [2019-04-13 05:02:26] parameters (PureCN 1.12.2).
INFO [2019-04-13 05:02:26] ------------------------------------------------------------
INFO [2019-04-13 05:02:26] PureCN 1.12.2
INFO [2019-04-13 05:02:26] ------------------------------------------------------------
FATAL [2019-04-13 05:02:26] prior.purity not within expected range or format.
FATAL [2019-04-13 05:02:26]
FATAL [2019-04-13 05:02:26] This is most likely a user error due to invalid input data or
FATAL [2019-04-13 05:02:26] parameters (PureCN 1.12.2).
INFO [2019-04-13 05:02:26] ------------------------------------------------------------
INFO [2019-04-13 05:02:26] PureCN 1.12.2
INFO [2019-04-13 05:02:26] ------------------------------------------------------------
FATAL [2019-04-13 05:02:26] prior.purity must add to 1. Sum is 1.5
FATAL [2019-04-13 05:02:26]
FATAL [2019-04-13 05:02:26] This is most likely a user error due to invalid input data or
FATAL [2019-04-13 05:02:26] parameters (PureCN 1.12.2).
INFO [2019-04-13 05:02:26] ------------------------------------------------------------
INFO [2019-04-13 05:02:26] PureCN 1.12.2
INFO [2019-04-13 05:02:26] ------------------------------------------------------------
FATAL [2019-04-13 05:02:26] max.homozygous.loss not within expected range or format.
FATAL [2019-04-13 05:02:26]
FATAL [2019-04-13 05:02:26] This is most likely a user error due to invalid input data or
FATAL [2019-04-13 05:02:26] parameters (PureCN 1.12.2).
INFO [2019-04-13 05:02:26] ------------------------------------------------------------
INFO [2019-04-13 05:02:26] PureCN 1.12.2
INFO [2019-04-13 05:02:26] ------------------------------------------------------------
FATAL [2019-04-13 05:02:26] prior.K not within expected range or format.
FATAL [2019-04-13 05:02:26]
FATAL [2019-04-13 05:02:26] This is most likely a user error due to invalid input data or
FATAL [2019-04-13 05:02:26] parameters (PureCN 1.12.2).
INFO [2019-04-13 05:02:26] ------------------------------------------------------------
INFO [2019-04-13 05:02:26] PureCN 1.12.2
INFO [2019-04-13 05:02:26] ------------------------------------------------------------
FATAL [2019-04-13 05:02:26] prior.contamination not within expected range or format.
FATAL [2019-04-13 05:02:26]
FATAL [2019-04-13 05:02:26] This is most likely a user error due to invalid input data or
FATAL [2019-04-13 05:02:26] parameters (PureCN 1.12.2).
INFO [2019-04-13 05:02:26] ------------------------------------------------------------
INFO [2019-04-13 05:02:26] PureCN 1.12.2
INFO [2019-04-13 05:02:26] ------------------------------------------------------------
FATAL [2019-04-13 05:02:26] Iterations not in the expected range from 10 to 250.
FATAL [2019-04-13 05:02:26]
FATAL [2019-04-13 05:02:26] This is most likely a user error due to invalid input data or
FATAL [2019-04-13 05:02:26] parameters (PureCN 1.12.2).
INFO [2019-04-13 05:02:26] ------------------------------------------------------------
INFO [2019-04-13 05:02:26] PureCN 1.12.2
INFO [2019-04-13 05:02:26] ------------------------------------------------------------
FATAL [2019-04-13 05:02:26] Iterations not in the expected range from 10 to 250.
FATAL [2019-04-13 05:02:26]
FATAL [2019-04-13 05:02:26] This is most likely a user error due to invalid input data or
FATAL [2019-04-13 05:02:26] parameters (PureCN 1.12.2).
INFO [2019-04-13 05:02:26] ------------------------------------------------------------
INFO [2019-04-13 05:02:26] PureCN 1.12.2
INFO [2019-04-13 05:02:26] ------------------------------------------------------------
INFO [2019-04-13 05:02:26] Loading coverage files...
FATAL [2019-04-13 05:02:27] Missing tumor.coverage.file requires seg.file or log.ratio and
FATAL [2019-04-13 05:02:27] interval.file.
FATAL [2019-04-13 05:02:27]
FATAL [2019-04-13 05:02:27] This is most likely a user error due to invalid input data or
FATAL [2019-04-13 05:02:27] parameters (PureCN 1.12.2).
INFO [2019-04-13 05:02:27] ------------------------------------------------------------
INFO [2019-04-13 05:02:27] PureCN 1.12.2
INFO [2019-04-13 05:02:27] ------------------------------------------------------------
INFO [2019-04-13 05:02:27] ------------------------------------------------------------
INFO [2019-04-13 05:02:27] PureCN 1.12.2
INFO [2019-04-13 05:02:27] ------------------------------------------------------------
INFO [2019-04-13 05:02:27] Loading coverage files...
FATAL [2019-04-13 05:02:27] Interval files in normal and tumor different.
FATAL [2019-04-13 05:02:27]
FATAL [2019-04-13 05:02:27] This is most likely a user error due to invalid input data or
FATAL [2019-04-13 05:02:27] parameters (PureCN 1.12.2).
INFO [2019-04-13 05:02:28] ------------------------------------------------------------
INFO [2019-04-13 05:02:28] PureCN 1.12.2
INFO [2019-04-13 05:02:28] ------------------------------------------------------------
INFO [2019-04-13 05:02:28] Loading coverage files...
INFO [2019-04-13 05:02:28] Mean target coverages: 0X (tumor) 99X (normal).
WARN [2019-04-13 05:02:28] Large difference in coverage of tumor and normal.
WARN [2019-04-13 05:02:28] Allosome coverage missing, cannot determine sex.
WARN [2019-04-13 05:02:28] Allosome coverage missing, cannot determine sex.
INFO [2019-04-13 05:02:28] Removing 10049 intervals with missing log.ratio.
WARN [2019-04-13 05:02:28] No normalDB provided. Provide one for better results.
INFO [2019-04-13 05:02:28] Using 0 intervals (0 on-target, 0 off-target).
FATAL [2019-04-13 05:02:28] No intervals passing filters.
FATAL [2019-04-13 05:02:28]
FATAL [2019-04-13 05:02:28] This is most likely a user error due to invalid input data or
FATAL [2019-04-13 05:02:28] parameters (PureCN 1.12.2).
INFO [2019-04-13 05:02:28] ------------------------------------------------------------
INFO [2019-04-13 05:02:28] PureCN 1.12.2
INFO [2019-04-13 05:02:28] ------------------------------------------------------------
INFO [2019-04-13 05:02:28] Loading coverage files...
INFO [2019-04-13 05:02:28] Mean target coverages: 112X (tumor) 99X (normal).
WARN [2019-04-13 05:02:29] Allosome coverage missing, cannot determine sex.
WARN [2019-04-13 05:02:29] Allosome coverage missing, cannot determine sex.
INFO [2019-04-13 05:02:29] Removing 228 intervals with missing log.ratio.
INFO [2019-04-13 05:02:29] Removing 22 small (< 5bp) intervals.
INFO [2019-04-13 05:02:29] Removing 19 intervals with low total coverage in normal (< 150.00 reads).
WARN [2019-04-13 05:02:29] No normalDB provided. Provide one for better results.
INFO [2019-04-13 05:02:29] Removing 233 low coverage (< 15.0000X) intervals.
INFO [2019-04-13 05:02:29] Using 9547 intervals (9547 on-target, 0 off-target).
INFO [2019-04-13 05:02:29] No off-target intervals. If this is hybrid-capture data, consider adding them.
INFO [2019-04-13 05:02:29] No interval.file provided. Cannot check if data was GC-normalized. Was it?
INFO [2019-04-13 05:02:29] Loading VCF...
INFO [2019-04-13 05:02:29] Found 127 variants in VCF file.
INFO [2019-04-13 05:02:29] 127 (100.0%) variants annotated as likely germline (DB INFO flag).
INFO [2019-04-13 05:02:29] LIB-02240e4 is tumor in VCF file.
INFO [2019-04-13 05:02:29] No homozygous variants in VCF, provide unfiltered VCF.
INFO [2019-04-13 05:02:29] Removing 1 non heterozygous (in matched normal) germline SNPs.
INFO [2019-04-13 05:02:30] Removing 1 variants with AF < 0.030 or AF >= 1.000 or less than 4 supporting reads or depth < 15.
INFO [2019-04-13 05:02:30] Removing 0 low quality variants with BQ < 25.
INFO [2019-04-13 05:02:30] Total size of targeted genomic region: 1.63Mb (2.58Mb with 50bp padding).
INFO [2019-04-13 05:02:30] 1.0% of targets contain variants.
INFO [2019-04-13 05:02:30] Removing 2 variants outside intervals.
INFO [2019-04-13 05:02:30] Found SOMATIC annotation in VCF.
INFO [2019-04-13 05:02:30] Setting somatic prior probabilities for somatic variants to 0.999000 or to 0.000100 otherwise.
INFO [2019-04-13 05:02:30] Found SOMATIC annotation in VCF. Setting mapping bias to 0.973.
INFO [2019-04-13 05:02:30] Sample sex: ?
INFO [2019-04-13 05:02:30] Segmenting data...
INFO [2019-04-13 05:02:30] Loading pre-computed boundaries for DNAcopy...
INFO [2019-04-13 05:02:30] Setting undo.SD parameter to 1.000000.
INFO [2019-04-13 05:02:31] Setting prune.hclust.h parameter to 0.200000.
INFO [2019-04-13 05:02:32] Found 54 segments with median size of 24.88Mb.
INFO [2019-04-13 05:02:32] Using 123 variants.
INFO [2019-04-13 05:02:32] Mean standard deviation of log-ratios: 0.37
INFO [2019-04-13 05:02:32] 2D-grid search of purity and ploidy...
INFO [2019-04-13 05:02:33] Local optima: 0.62/1.9, 0.3/1.9, 0.32/2
INFO [2019-04-13 05:02:33] Testing local optimum 1/3 at purity 0.62 and total ploidy 1.90...
INFO [2019-04-13 05:02:34] Fitting variants for purity 0.65, tumor ploidy 1.74 and contamination 0.01.
INFO [2019-04-13 05:02:35] Optimized purity: 0.65
INFO [2019-04-13 05:02:35] Testing local optimum 2/3 at purity 0.30 and total ploidy 1.90...
INFO [2019-04-13 05:02:36] Recalibrating log-ratios...
INFO [2019-04-13 05:02:36] Testing local optimum 2/3 at purity 0.30 and total ploidy 1.90...
INFO [2019-04-13 05:02:36] Recalibrating log-ratios...
INFO [2019-04-13 05:02:36] Testing local optimum 2/3 at purity 0.30 and total ploidy 1.90...
INFO [2019-04-13 05:02:37] Recalibrating log-ratios...
INFO [2019-04-13 05:02:37] Testing local optimum 2/3 at purity 0.30 and total ploidy 1.90...
INFO [2019-04-13 05:02:38] Testing local optimum 3/3 at purity 0.32 and total ploidy 2.00...
INFO [2019-04-13 05:02:38] Recalibrating log-ratios...
INFO [2019-04-13 05:02:38] Testing local optimum 3/3 at purity 0.32 and total ploidy 2.00...
INFO [2019-04-13 05:02:38] Recalibrating log-ratios...
INFO [2019-04-13 05:02:38] Testing local optimum 3/3 at purity 0.32 and total ploidy 2.00...
INFO [2019-04-13 05:02:39] Recalibrating log-ratios...
INFO [2019-04-13 05:02:39] Testing local optimum 3/3 at purity 0.32 and total ploidy 2.00...
INFO [2019-04-13 05:02:39] Done.
INFO [2019-04-13 05:02:39] ------------------------------------------------------------
INFO [2019-04-13 05:02:40] ------------------------------------------------------------
INFO [2019-04-13 05:02:40] PureCN 1.12.2
INFO [2019-04-13 05:02:40] ------------------------------------------------------------
INFO [2019-04-13 05:02:40] Loading coverage files...
INFO [2019-04-13 05:02:40] Mean target coverages: 112X (tumor) 99X (normal).
WARN [2019-04-13 05:02:40] Allosome coverage missing, cannot determine sex.
WARN [2019-04-13 05:02:40] Allosome coverage missing, cannot determine sex.
INFO [2019-04-13 05:02:40] Removing 228 intervals with missing log.ratio.
INFO [2019-04-13 05:02:40] Removing 22 small (< 5bp) intervals.
INFO [2019-04-13 05:02:40] Removing 19 intervals with low total coverage in normal (< 150.00 reads).
WARN [2019-04-13 05:02:40] No normalDB provided. Provide one for better results.
INFO [2019-04-13 05:02:40] Removing 233 low coverage (< 15.0000X) intervals.
INFO [2019-04-13 05:02:40] Using 9547 intervals (9547 on-target, 0 off-target).
INFO [2019-04-13 05:02:40] No off-target intervals. If this is hybrid-capture data, consider adding them.
INFO [2019-04-13 05:02:40] No interval.file provided. Cannot check if data was GC-normalized. Was it?
INFO [2019-04-13 05:02:40] Loading VCF...
INFO [2019-04-13 05:02:40] Found 127 variants in VCF file.
WARN [2019-04-13 05:02:40] DP FORMAT field contains NAs. Removing 3 variants.
WARN [2019-04-13 05:02:40] DB INFO flag contains NAs
INFO [2019-04-13 05:02:40] 119 (96.0%) variants annotated as likely germline (DB INFO flag).
INFO [2019-04-13 05:02:40] LIB-02240e4 is tumor in VCF file.
INFO [2019-04-13 05:02:41] No homozygous variants in VCF, provide unfiltered VCF.
INFO [2019-04-13 05:02:41] Removing 1 variants with AF < 0.030 or AF >= 0.970 or less than 4 supporting reads or depth < 15.
INFO [2019-04-13 05:02:41] Removing 0 low quality variants with BQ < 25.
INFO [2019-04-13 05:02:41] Total size of targeted genomic region: 1.63Mb (2.58Mb with 50bp padding).
INFO [2019-04-13 05:02:41] 1.0% of targets contain variants.
INFO [2019-04-13 05:02:41] Setting somatic prior probabilities for dbSNP hits to 0.000500 or to 0.500000 otherwise.
INFO [2019-04-13 05:02:41] Sample sex: ?
INFO [2019-04-13 05:02:42] Segmenting data...
INFO [2019-04-13 05:02:42] Loading pre-computed boundaries for DNAcopy...
INFO [2019-04-13 05:02:42] Setting undo.SD parameter to 1.000000.
Setting multi-figure configuration
INFO [2019-04-13 05:02:43] Setting prune.hclust.h parameter to 0.200000.
INFO [2019-04-13 05:02:43] Found 54 segments with median size of 24.88Mb.
INFO [2019-04-13 05:02:43] Using 123 variants.
INFO [2019-04-13 05:02:43] Mean standard deviation of log-ratios: 0.37
INFO [2019-04-13 05:02:43] 2D-grid search of purity and ploidy...
INFO [2019-04-13 05:02:45] Local optima: 0.62/1.9, 0.3/1.9, 0.32/2
INFO [2019-04-13 05:02:45] Testing local optimum 1/3 at purity 0.62 and total ploidy 1.90...
INFO [2019-04-13 05:02:46] Fitting variants for purity 0.65, tumor ploidy 1.74 and contamination 0.01.
INFO [2019-04-13 05:02:47] Optimized purity: 0.65
INFO [2019-04-13 05:02:47] Testing local optimum 2/3 at purity 0.30 and total ploidy 1.90...
INFO [2019-04-13 05:02:47] Recalibrating log-ratios...
INFO [2019-04-13 05:02:47] Testing local optimum 2/3 at purity 0.30 and total ploidy 1.90...
INFO [2019-04-13 05:02:48] Recalibrating log-ratios...
INFO [2019-04-13 05:02:48] Testing local optimum 2/3 at purity 0.30 and total ploidy 1.90...
INFO [2019-04-13 05:02:48] Recalibrating log-ratios...
INFO [2019-04-13 05:02:48] Testing local optimum 2/3 at purity 0.30 and total ploidy 1.90...
INFO [2019-04-13 05:02:49] Testing local optimum 3/3 at purity 0.32 and total ploidy 2.00...
INFO [2019-04-13 05:02:49] Recalibrating log-ratios...
INFO [2019-04-13 05:02:49] Testing local optimum 3/3 at purity 0.32 and total ploidy 2.00...
INFO [2019-04-13 05:02:49] Recalibrating log-ratios...
INFO [2019-04-13 05:02:49] Testing local optimum 3/3 at purity 0.32 and total ploidy 2.00...
INFO [2019-04-13 05:02:50] Recalibrating log-ratios...
INFO [2019-04-13 05:02:50] Testing local optimum 3/3 at purity 0.32 and total ploidy 2.00...
INFO [2019-04-13 05:02:50] Done.
INFO [2019-04-13 05:02:50] ------------------------------------------------------------
INFO [2019-04-13 05:02:50] ------------------------------------------------------------
INFO [2019-04-13 05:02:50] PureCN 1.12.2
INFO [2019-04-13 05:02:50] ------------------------------------------------------------
INFO [2019-04-13 05:02:50] Loading coverage files...
INFO [2019-04-13 05:02:51] Mean target coverages: 112X (tumor) 99X (normal).
WARN [2019-04-13 05:02:51] Allosome coverage missing, cannot determine sex.
WARN [2019-04-13 05:02:51] Allosome coverage missing, cannot determine sex.
INFO [2019-04-13 05:02:51] No Gene column in interval.file. You won't get gene-level calls.
INFO [2019-04-13 05:02:51] Removing 228 intervals with missing log.ratio.
INFO [2019-04-13 05:02:51] Removing 15 low/high GC targets.
INFO [2019-04-13 05:02:51] Removing 21 small (< 5bp) intervals.
INFO [2019-04-13 05:02:51] Removing 19 intervals with low total coverage in normal (< 150.00 reads).
WARN [2019-04-13 05:02:51] No normalDB provided. Provide one for better results.
INFO [2019-04-13 05:02:51] Removing 232 low coverage (< 15.0000X) intervals.
INFO [2019-04-13 05:02:51] Using 9534 intervals (9534 on-target, 0 off-target).
INFO [2019-04-13 05:02:51] No off-target intervals. If this is hybrid-capture data, consider adding them.
INFO [2019-04-13 05:02:51] AT/GC dropout: 1.03 (tumor), 1.03 (normal).
INFO [2019-04-13 05:02:51] Loading VCF...
INFO [2019-04-13 05:02:51] Found 127 variants in VCF file.
INFO [2019-04-13 05:02:52] 127 (100.0%) variants annotated as likely germline (DB INFO flag).
INFO [2019-04-13 05:02:52] LIB-02240e4 is tumor in VCF file.
INFO [2019-04-13 05:02:52] No homozygous variants in VCF, provide unfiltered VCF.
INFO [2019-04-13 05:02:52] Removing 1 non heterozygous (in matched normal) germline SNPs.
INFO [2019-04-13 05:02:53] Removing 1 variants with AF < 0.030 or AF >= Inf or less than 4 supporting reads or depth < 15.
INFO [2019-04-13 05:02:53] Removing 0 low quality variants with BQ < 25.
INFO [2019-04-13 05:02:53] Total size of targeted genomic region: 1.63Mb (2.58Mb with 50bp padding).
INFO [2019-04-13 05:02:53] 1.0% of targets contain variants.
INFO [2019-04-13 05:02:53] Removing 2 variants outside intervals.
INFO [2019-04-13 05:02:53] Found SOMATIC annotation in VCF.
INFO [2019-04-13 05:02:53] Setting somatic prior probabilities for somatic variants to 0.999000 or to 0.000100 otherwise.
INFO [2019-04-13 05:02:53] Found SOMATIC annotation in VCF. Setting mapping bias to 0.973.
INFO [2019-04-13 05:02:53] Sample sex: ?
INFO [2019-04-13 05:02:53] Segmenting data...
INFO [2019-04-13 05:02:53] Loading pre-computed boundaries for DNAcopy...
INFO [2019-04-13 05:02:53] Setting undo.SD parameter to 1.000000.
Setting multi-figure configuration
INFO [2019-04-13 05:02:54] Setting prune.hclust.h parameter to 0.200000.
INFO [2019-04-13 05:02:55] Found 54 segments with median size of 24.88Mb.
INFO [2019-04-13 05:02:55] Using 123 variants.
INFO [2019-04-13 05:02:55] Mean standard deviation of log-ratios: 0.37
INFO [2019-04-13 05:02:55] 2D-grid search of purity and ploidy...
INFO [2019-04-13 05:02:55] Local optima: 0.65/1.9, 0.33/1.9, 0.38/2.2, 0.27/2.1, 0.29/2, 0.78/2
INFO [2019-04-13 05:02:55] Testing local optimum 1/6 at purity 0.65 and total ploidy 1.90...
INFO [2019-04-13 05:02:55] Recalibrating log-ratios...
INFO [2019-04-13 05:02:55] Testing local optimum 1/6 at purity 0.65 and total ploidy 1.90...
INFO [2019-04-13 05:02:56] Recalibrating log-ratios...
INFO [2019-04-13 05:02:56] Testing local optimum 1/6 at purity 0.65 and total ploidy 1.90...
INFO [2019-04-13 05:02:56] Recalibrating log-ratios...
INFO [2019-04-13 05:02:56] Testing local optimum 1/6 at purity 0.65 and total ploidy 1.90...
INFO [2019-04-13 05:02:57] Testing local optimum 2/6 at purity 0.33 and total ploidy 1.90...
INFO [2019-04-13 05:02:57] Recalibrating log-ratios...
INFO [2019-04-13 05:02:57] Testing local optimum 2/6 at purity 0.33 and total ploidy 1.90...
INFO [2019-04-13 05:02:57] Recalibrating log-ratios...
INFO [2019-04-13 05:02:57] Testing local optimum 2/6 at purity 0.33 and total ploidy 1.90...
INFO [2019-04-13 05:02:58] Recalibrating log-ratios...
INFO [2019-04-13 05:02:58] Testing local optimum 2/6 at purity 0.33 and total ploidy 1.90...
INFO [2019-04-13 05:02:58] Testing local optimum 3/6 at purity 0.38 and total ploidy 2.20...
INFO [2019-04-13 05:02:59] Fitting variants for purity 0.40, tumor ploidy 2.44 and contamination 0.01.
INFO [2019-04-13 05:03:00] Optimized purity: 0.40
INFO [2019-04-13 05:03:00] Testing local optimum 4/6 at purity 0.27 and total ploidy 2.10...
INFO [2019-04-13 05:03:01] Recalibrating log-ratios...
INFO [2019-04-13 05:03:01] Testing local optimum 4/6 at purity 0.27 and total ploidy 2.10...
INFO [2019-04-13 05:03:01] Recalibrating log-ratios...
INFO [2019-04-13 05:03:01] Testing local optimum 4/6 at purity 0.27 and total ploidy 2.10...
INFO [2019-04-13 05:03:02] Recalibrating log-ratios...
INFO [2019-04-13 05:03:02] Testing local optimum 4/6 at purity 0.27 and total ploidy 2.10...
INFO [2019-04-13 05:03:02] Testing local optimum 5/6 at purity 0.29 and total ploidy 2.00...
INFO [2019-04-13 05:03:02] Recalibrating log-ratios...
INFO [2019-04-13 05:03:02] Testing local optimum 5/6 at purity 0.29 and total ploidy 2.00...
INFO [2019-04-13 05:03:03] Recalibrating log-ratios...
INFO [2019-04-13 05:03:03] Testing local optimum 5/6 at purity 0.29 and total ploidy 2.00...
INFO [2019-04-13 05:03:04] Recalibrating log-ratios...
INFO [2019-04-13 05:03:04] Testing local optimum 5/6 at purity 0.29 and total ploidy 2.00...
INFO [2019-04-13 05:03:04] Testing local optimum 6/6 at purity 0.78 and total ploidy 2.00...
INFO [2019-04-13 05:03:05] Recalibrating log-ratios...
INFO [2019-04-13 05:03:05] Testing local optimum 6/6 at purity 0.78 and total ploidy 2.00...
INFO [2019-04-13 05:03:05] Recalibrating log-ratios...
INFO [2019-04-13 05:03:05] Testing local optimum 6/6 at purity 0.78 and total ploidy 2.00...
INFO [2019-04-13 05:03:06] Recalibrating log-ratios...
INFO [2019-04-13 05:03:06] Testing local optimum 6/6 at purity 0.78 and total ploidy 2.00...
INFO [2019-04-13 05:03:06] Done.
INFO [2019-04-13 05:03:06] ------------------------------------------------------------
FATAL [2019-04-13 05:03:06] This function requires gene-level calls. Please add a column 'Gene'
FATAL [2019-04-13 05:03:06] containing gene symbols to the interval.file.
FATAL [2019-04-13 05:03:06]
FATAL [2019-04-13 05:03:06] This is most likely a user error due to invalid input data or
FATAL [2019-04-13 05:03:06] parameters (PureCN 1.12.2).
INFO [2019-04-13 05:03:08] ------------------------------------------------------------
INFO [2019-04-13 05:03:08] PureCN 1.12.2
INFO [2019-04-13 05:03:08] ------------------------------------------------------------
INFO [2019-04-13 05:03:08] Loading coverage files...
INFO [2019-04-13 05:03:08] Mean target coverages: 112X (tumor) 99X (normal).
WARN [2019-04-13 05:03:08] Allosome coverage missing, cannot determine sex.
WARN [2019-04-13 05:03:08] Allosome coverage missing, cannot determine sex.
INFO [2019-04-13 05:03:08] Removing 228 intervals with missing log.ratio.
INFO [2019-04-13 05:03:08] Removing 22 small (< 5bp) intervals.
INFO [2019-04-13 05:03:08] Removing 19 intervals with low total coverage in normal (< 150.00 reads).
WARN [2019-04-13 05:03:08] No normalDB provided. Provide one for better results.
INFO [2019-04-13 05:03:08] Removing 233 low coverage (< 15.0000X) intervals.
INFO [2019-04-13 05:03:08] Using 9547 intervals (9547 on-target, 0 off-target).
INFO [2019-04-13 05:03:08] No off-target intervals. If this is hybrid-capture data, consider adding them.
INFO [2019-04-13 05:03:08] No interval.file provided. Cannot check if data was GC-normalized. Was it?
INFO [2019-04-13 05:03:08] Loading VCF...
INFO [2019-04-13 05:03:08] Found 127 variants in VCF file.
INFO [2019-04-13 05:03:09] 127 (100.0%) variants annotated as likely germline (DB INFO flag).
FATAL [2019-04-13 05:03:09] Different chromosome names in coverage and VCF.
FATAL [2019-04-13 05:03:09]
FATAL [2019-04-13 05:03:09] This is most likely a user error due to invalid input data or
FATAL [2019-04-13 05:03:09] parameters (PureCN 1.12.2).
INFO [2019-04-13 05:03:09] ------------------------------------------------------------
INFO [2019-04-13 05:03:09] PureCN 1.12.2
INFO [2019-04-13 05:03:09] ------------------------------------------------------------
INFO [2019-04-13 05:03:09] Loading coverage files...
INFO [2019-04-13 05:03:09] Mean target coverages: 112X (tumor) 99X (normal).
WARN [2019-04-13 05:03:09] Allosome coverage missing, cannot determine sex.
WARN [2019-04-13 05:03:09] Allosome coverage missing, cannot determine sex.
FATAL [2019-04-13 05:03:10] tumor.coverage.file and interval.file do not align.
FATAL [2019-04-13 05:03:10]
FATAL [2019-04-13 05:03:10] This is most likely a user error due to invalid input data or
FATAL [2019-04-13 05:03:10] parameters (PureCN 1.12.2).
INFO [2019-04-13 05:03:10] ------------------------------------------------------------
INFO [2019-04-13 05:03:10] PureCN 1.12.2
INFO [2019-04-13 05:03:10] ------------------------------------------------------------
INFO [2019-04-13 05:03:10] Loading coverage files...
INFO [2019-04-13 05:03:10] Mean target coverages: 112X (tumor) 99X (normal).
WARN [2019-04-13 05:03:10] Allosome coverage missing, cannot determine sex.
WARN [2019-04-13 05:03:10] Allosome coverage missing, cannot determine sex.
INFO [2019-04-13 05:03:12] Removing 228 intervals with missing log.ratio.
INFO [2019-04-13 05:03:12] Removing 15 low/high GC targets.
INFO [2019-04-13 05:03:12] Removing 21 small (< 5bp) intervals.
INFO [2019-04-13 05:03:12] Removing 19 intervals with low total coverage in normal (< 150.00 reads).
WARN [2019-04-13 05:03:12] No normalDB provided. Provide one for better results.
INFO [2019-04-13 05:03:12] Removing 232 low coverage (< 15.0000X) intervals.
INFO [2019-04-13 05:03:12] Using 9534 intervals (9534 on-target, 0 off-target).
INFO [2019-04-13 05:03:12] No off-target intervals. If this is hybrid-capture data, consider adding them.
INFO [2019-04-13 05:03:12] AT/GC dropout: 1.03 (tumor), 1.03 (normal).
INFO [2019-04-13 05:03:12] Loading VCF...
INFO [2019-04-13 05:03:12] Found 127 variants in VCF file.
INFO [2019-04-13 05:03:12] 127 (100.0%) variants annotated as likely germline (DB INFO flag).
INFO [2019-04-13 05:03:12] LIB-02240e4 is tumor in VCF file.
INFO [2019-04-13 05:03:13] No homozygous variants in VCF, provide unfiltered VCF.
INFO [2019-04-13 05:03:13] Removing 1 non heterozygous (in matched normal) germline SNPs.
INFO [2019-04-13 05:03:13] Removing 1 variants with AF < 0.030 or AF >= 1.000 or less than 4 supporting reads or depth < 15.
INFO [2019-04-13 05:03:13] Removing 0 low quality variants with BQ < 25.
INFO [2019-04-13 05:03:13] Total size of targeted genomic region: 1.63Mb (2.58Mb with 50bp padding).
INFO [2019-04-13 05:03:13] 1.0% of targets contain variants.
INFO [2019-04-13 05:03:14] Removing 2 variants outside intervals.
INFO [2019-04-13 05:03:14] Found SOMATIC annotation in VCF.
INFO [2019-04-13 05:03:14] Setting somatic prior probabilities for somatic variants to 0.999000 or to 0.000100 otherwise.
INFO [2019-04-13 05:03:14] Found SOMATIC annotation in VCF. Setting mapping bias to 0.973.
INFO [2019-04-13 05:03:14] Sample sex: ?
INFO [2019-04-13 05:03:14] Segmenting data...
INFO [2019-04-13 05:03:14] Loading pre-computed boundaries for DNAcopy...
INFO [2019-04-13 05:03:14] Setting undo.SD parameter to 1.000000.
INFO [2019-04-13 05:03:15] Setting prune.hclust.h parameter to 0.200000.
INFO [2019-04-13 05:03:15] Found 54 segments with median size of 24.88Mb.
INFO [2019-04-13 05:03:15] Using 123 variants.
INFO [2019-04-13 05:03:15] Mean standard deviation of log-ratios: 0.37
INFO [2019-04-13 05:03:15] 2D-grid search of purity and ploidy...
INFO [2019-04-13 05:03:16] Local optima: 0.63/1.9, 0.4/2.1
INFO [2019-04-13 05:03:16] Testing local optimum 1/2 at purity 0.63 and total ploidy 1.90...
INFO [2019-04-13 05:03:18] Fitting variants for purity 0.65, tumor ploidy 1.74 and contamination 0.01.
INFO [2019-04-13 05:03:19] Optimized purity: 0.65
INFO [2019-04-13 05:03:19] Testing local optimum 2/2 at purity 0.40 and total ploidy 2.10...
INFO [2019-04-13 05:03:19] Recalibrating log-ratios...
INFO [2019-04-13 05:03:19] Testing local optimum 2/2 at purity 0.40 and total ploidy 2.10...
INFO [2019-04-13 05:03:19] Recalibrating log-ratios...
INFO [2019-04-13 05:03:19] Testing local optimum 2/2 at purity 0.40 and total ploidy 2.10...
INFO [2019-04-13 05:03:20] Recalibrating log-ratios...
INFO [2019-04-13 05:03:20] Testing local optimum 2/2 at purity 0.40 and total ploidy 2.10...
INFO [2019-04-13 05:03:21] Done.
INFO [2019-04-13 05:03:21] ------------------------------------------------------------
INFO [2019-04-13 05:03:21] ------------------------------------------------------------
INFO [2019-04-13 05:03:21] PureCN 1.12.2
INFO [2019-04-13 05:03:21] ------------------------------------------------------------
INFO [2019-04-13 05:03:21] Loading coverage files...
INFO [2019-04-13 05:03:21] Mean target coverages: 112X (tumor) 99X (normal).
WARN [2019-04-13 05:03:21] Allosome coverage missing, cannot determine sex.
WARN [2019-04-13 05:03:21] Allosome coverage missing, cannot determine sex.
INFO [2019-04-13 05:03:23] Removing 228 intervals with missing log.ratio.
INFO [2019-04-13 05:03:23] Removing 15 low/high GC targets.
INFO [2019-04-13 05:03:23] Removing 21 small (< 5bp) intervals.
INFO [2019-04-13 05:03:23] Removing 19 intervals with low total coverage in normal (< 150.00 reads).
WARN [2019-04-13 05:03:23] No normalDB provided. Provide one for better results.
INFO [2019-04-13 05:03:23] Removing 232 low coverage (< 15.0000X) intervals.
INFO [2019-04-13 05:03:23] Using 9534 intervals (9534 on-target, 0 off-target).
INFO [2019-04-13 05:03:23] No off-target intervals. If this is hybrid-capture data, consider adding them.
INFO [2019-04-13 05:03:23] AT/GC dropout: 1.03 (tumor), 1.03 (normal).
INFO [2019-04-13 05:03:23] Loading VCF...
INFO [2019-04-13 05:03:23] Found 127 variants in VCF file.
INFO [2019-04-13 05:03:23] 127 (100.0%) variants annotated as likely germline (DB INFO flag).
INFO [2019-04-13 05:03:23] LIB-02240e4 is tumor in VCF file.
INFO [2019-04-13 05:03:23] No homozygous variants in VCF, provide unfiltered VCF.
INFO [2019-04-13 05:03:24] Removing 1 non heterozygous (in matched normal) germline SNPs.
INFO [2019-04-13 05:03:24] Removing 1 variants with AF < 0.030 or AF >= 1.000 or less than 4 supporting reads or depth < 15.
INFO [2019-04-13 05:03:24] Removing 0 low quality variants with BQ < 25.
INFO [2019-04-13 05:03:25] Total size of targeted genomic region: 1.63Mb (2.58Mb with 50bp padding).
INFO [2019-04-13 05:03:25] 1.0% of targets contain variants.
INFO [2019-04-13 05:03:25] Removing 2 variants outside intervals.
INFO [2019-04-13 05:03:25] Found SOMATIC annotation in VCF.
INFO [2019-04-13 05:03:25] Setting somatic prior probabilities for somatic variants to 0.999000 or to 0.000100 otherwise.
INFO [2019-04-13 05:03:25] Found SOMATIC annotation in VCF. Setting mapping bias to 0.973.
INFO [2019-04-13 05:03:25] Sample sex: ?
INFO [2019-04-13 05:03:25] Segmenting data...
INFO [2019-04-13 05:03:25] Loading pre-computed boundaries for DNAcopy...
INFO [2019-04-13 05:03:25] Setting undo.SD parameter to 1.000000.
INFO [2019-04-13 05:03:26] Setting prune.hclust.h parameter to 0.200000.
INFO [2019-04-13 05:03:26] Found 54 segments with median size of 24.88Mb.
INFO [2019-04-13 05:03:27] Using 123 variants.
INFO [2019-04-13 05:03:27] Mean standard deviation of log-ratios: 0.37
INFO [2019-04-13 05:03:27] 2D-grid search of purity and ploidy...
INFO [2019-04-13 05:03:28] Local optima: 0.63/1.9, 0.4/2.1
INFO [2019-04-13 05:03:28] Testing local optimum 1/2 at purity 0.63 and total ploidy 1.90...
INFO [2019-04-13 05:03:30] Fitting variants for purity 0.65, tumor ploidy 1.74 and contamination 0.01.
INFO [2019-04-13 05:03:32] Optimized purity: 0.65
INFO [2019-04-13 05:03:32] Testing local optimum 2/2 at purity 0.40 and total ploidy 2.10...
INFO [2019-04-13 05:03:32] Recalibrating log-ratios...
INFO [2019-04-13 05:03:32] Testing local optimum 2/2 at purity 0.40 and total ploidy 2.10...
INFO [2019-04-13 05:03:33] Recalibrating log-ratios...
INFO [2019-04-13 05:03:33] Testing local optimum 2/2 at purity 0.40 and total ploidy 2.10...
INFO [2019-04-13 05:03:33] Recalibrating log-ratios...
INFO [2019-04-13 05:03:33] Testing local optimum 2/2 at purity 0.40 and total ploidy 2.10...
INFO [2019-04-13 05:03:34] Done.
INFO [2019-04-13 05:03:34] ------------------------------------------------------------
FATAL [2019-04-13 05:03:34] chr1 not valid chromosome name(s). Valid names are:
FATAL [2019-04-13 05:03:34] 1,2,3,4,5,6,7,8,9,10,11,12,13,14,15,16,17,18,19,20,21,22
FATAL [2019-04-13 05:03:34]
FATAL [2019-04-13 05:03:34] This is most likely a user error due to invalid input data or
FATAL [2019-04-13 05:03:34] parameters (PureCN 1.12.2).
INFO [2019-04-13 05:03:34] ------------------------------------------------------------
INFO [2019-04-13 05:03:34] PureCN 1.12.2
INFO [2019-04-13 05:03:34] ------------------------------------------------------------
INFO [2019-04-13 05:03:34] Loading coverage files...
WARN [2019-04-13 05:03:35] Provided sampleid (Sample2) does not match Sample1 found in segmentation.
WARN [2019-04-13 05:03:35] Allosome coverage missing, cannot determine sex.
WARN [2019-04-13 05:03:35] Allosome coverage missing, cannot determine sex.
INFO [2019-04-13 05:03:35] Removing 10 intervals with missing log.ratio.
INFO [2019-04-13 05:03:35] Using 10039 intervals (10039 on-target, 0 off-target).
INFO [2019-04-13 05:03:35] No off-target intervals. If this is hybrid-capture data, consider adding them.
INFO [2019-04-13 05:03:35] Loading VCF...
INFO [2019-04-13 05:03:35] Found 127 variants in VCF file.
INFO [2019-04-13 05:03:35] 127 (100.0%) variants annotated as likely germline (DB INFO flag).
INFO [2019-04-13 05:03:36] LIB-02240e4 is tumor in VCF file.
INFO [2019-04-13 05:03:36] No homozygous variants in VCF, provide unfiltered VCF.
INFO [2019-04-13 05:03:36] Removing 1 non heterozygous (in matched normal) germline SNPs.
INFO [2019-04-13 05:03:37] Removing 1 variants with AF < 0.030 or AF >= 1.000 or less than 4 supporting reads or depth < 15.
INFO [2019-04-13 05:03:37] Removing 0 low quality variants with BQ < 25.
INFO [2019-04-13 05:03:37] Total size of targeted genomic region: 1.71Mb (2.70Mb with 50bp padding).
INFO [2019-04-13 05:03:37] 1.0% of targets contain variants.
INFO [2019-04-13 05:03:37] Removing 0 variants outside intervals.
INFO [2019-04-13 05:03:37] Found SOMATIC annotation in VCF.
INFO [2019-04-13 05:03:37] Setting somatic prior probabilities for somatic variants to 0.999000 or to 0.000100 otherwise.
INFO [2019-04-13 05:03:38] Found SOMATIC annotation in VCF. Setting mapping bias to 0.976.
INFO [2019-04-13 05:03:38] Sample sex: ?
INFO [2019-04-13 05:03:38] Segmenting data...
WARN [2019-04-13 05:03:38] Provided sampleid (Sample2) does not match Sample1 found in segmentation.
INFO [2019-04-13 05:03:38] Loading pre-computed boundaries for DNAcopy...
INFO [2019-04-13 05:03:38] Setting undo.SD parameter to 0.000000.
Setting multi-figure configuration
INFO [2019-04-13 05:03:38] Setting prune.hclust.h parameter to 0.200000.
INFO [2019-04-13 05:03:38] Found 54 segments with median size of 24.88Mb.
INFO [2019-04-13 05:03:38] Using 125 variants.
INFO [2019-04-13 05:03:38] Mean standard deviation of log-ratios: 0.40
INFO [2019-04-13 05:03:38] 2D-grid search of purity and ploidy...
INFO [2019-04-13 05:03:39] Local optima: 0.6/1.9, 0.4/2.1
INFO [2019-04-13 05:03:39] Testing local optimum 1/2 at purity 0.60 and total ploidy 1.90...
INFO [2019-04-13 05:03:41] Fitting variants for purity 0.65, tumor ploidy 1.74 and contamination 0.01.
INFO [2019-04-13 05:03:42] Optimized purity: 0.65
INFO [2019-04-13 05:03:42] Testing local optimum 2/2 at purity 0.40 and total ploidy 2.10...
INFO [2019-04-13 05:03:43] Recalibrating log-ratios...
INFO [2019-04-13 05:03:43] Testing local optimum 2/2 at purity 0.40 and total ploidy 2.10...
INFO [2019-04-13 05:03:43] Recalibrating log-ratios...
INFO [2019-04-13 05:03:43] Testing local optimum 2/2 at purity 0.40 and total ploidy 2.10...
INFO [2019-04-13 05:03:43] Recalibrating log-ratios...
INFO [2019-04-13 05:03:43] Testing local optimum 2/2 at purity 0.40 and total ploidy 2.10...
INFO [2019-04-13 05:03:44] Done.
INFO [2019-04-13 05:03:44] ------------------------------------------------------------
WARN [2019-04-13 05:03:44] Allosome coverage missing, cannot determine sex.
WARN [2019-04-13 05:03:44] Allosome coverage missing, cannot determine sex.
Setting multi-figure configuration
FATAL [2019-04-13 05:03:55] Segmentation file expected with colnames ID, chrom, loc.start, loc.end,
FATAL [2019-04-13 05:03:55] num.mark, seg.mean
FATAL [2019-04-13 05:03:55]
FATAL [2019-04-13 05:03:55] This is most likely a user error due to invalid input data or
FATAL [2019-04-13 05:03:55] parameters (PureCN 1.12.2).
FATAL [2019-04-13 05:03:55] seg.file contains multiple samples and sampleid missing.
FATAL [2019-04-13 05:03:55]
FATAL [2019-04-13 05:03:55] This is most likely a user error due to invalid input data or
FATAL [2019-04-13 05:03:55] parameters (PureCN 1.12.2).
FATAL [2019-04-13 05:03:55] seg.file contains multiple samples and sampleid does not match any.
FATAL [2019-04-13 05:03:55]
FATAL [2019-04-13 05:03:55] This is most likely a user error due to invalid input data or
FATAL [2019-04-13 05:03:55] parameters (PureCN 1.12.2).
WARN [2019-04-13 05:03:56] Allosome coverage missing, cannot determine sex.
WARN [2019-04-13 05:03:56] Allosome coverage missing, cannot determine sex.
Setting multi-figure configuration
WARN [2019-04-13 05:04:07] Provided sampleid (example_tumor.txt) does not match Sample1 found in segmentation.
WARN [2019-04-13 05:04:07] Allosome coverage missing, cannot determine sex.
WARN [2019-04-13 05:04:07] Allosome coverage missing, cannot determine sex.
WARN [2019-04-13 05:04:07] Provided sampleid (example_tumor.txt) does not match Sample1 found in segmentation.
WARN [2019-04-13 05:04:13] Allosome coverage missing, cannot determine sex.
WARN [2019-04-13 05:04:13] Allosome coverage missing, cannot determine sex.
WARN [2019-04-13 05:04:14] No normalDB provided. Provide one for better results.
WARN [2019-04-13 05:04:24] Allosome coverage missing, cannot determine sex.
WARN [2019-04-13 05:04:24] Allosome coverage missing, cannot determine sex.
WARN [2019-04-13 05:04:27] Provided sampleid (Sample.1) does not match Sample1 found in segmentation.
Setting multi-figure configuration
WARN [2019-04-13 05:04:36] Allosome coverage missing, cannot determine sex.
WARN [2019-04-13 05:04:36] Allosome coverage missing, cannot determine sex.
WARN [2019-04-13 05:04:36] No normalDB provided. Provide one for better results.
WARN [2019-04-13 05:04:36] Sampleid looks like a normal in VCF, not like a tumor.
WARN [2019-04-13 05:04:46] Allosome coverage missing, cannot determine sex.
WARN [2019-04-13 05:04:46] Allosome coverage missing, cannot determine sex.
FATAL [2019-04-13 05:04:46] normalDB not a valid normalDB object. Use createNormalDatabase to
FATAL [2019-04-13 05:04:46] create one.
FATAL [2019-04-13 05:04:46]
FATAL [2019-04-13 05:04:46] This is most likely a user error due to invalid input data or
FATAL [2019-04-13 05:04:46] parameters (PureCN 1.12.2).
WARN [2019-04-13 05:04:46] You are likely not using the correct baits file!
WARN [2019-04-13 05:04:46] Allosome coverage missing, cannot determine sex.
WARN [2019-04-13 05:04:46] Allosome coverage missing, cannot determine sex.
WARN [2019-04-13 05:04:49] Allosome coverage missing, cannot determine sex.
WARN [2019-04-13 05:04:49] Allosome coverage missing, cannot determine sex.
FATAL [2019-04-13 05:04:49] normalDB appears to be empty.
FATAL [2019-04-13 05:04:49]
FATAL [2019-04-13 05:04:49] This is most likely a user error due to invalid input data or
FATAL [2019-04-13 05:04:49] parameters (PureCN 1.12.2).
WARN [2019-04-13 05:04:49] Allosome coverage missing, cannot determine sex.
WARN [2019-04-13 05:04:49] Allosome coverage missing, cannot determine sex.
WARN [2019-04-13 05:04:50] Intervals in coverage and interval.file have conflicting on/off-target annotation.
FATAL [2019-04-13 05:05:17] min.normals must be >=2.
FATAL [2019-04-13 05:05:17]
FATAL [2019-04-13 05:05:17] This is most likely a user error due to invalid input data or
FATAL [2019-04-13 05:05:17] parameters (PureCN 1.12.2).
FATAL [2019-04-13 05:05:31] The normal.panel.vcf.file contains only a single sample.
FATAL [2019-04-13 05:05:31]
FATAL [2019-04-13 05:05:31] This is most likely a user error due to invalid input data or
FATAL [2019-04-13 05:05:31] parameters (PureCN 1.12.2).
== testthat results ===========================================================
OK: 307 SKIPPED: 0 FAILED: 0
>
> proc.time()
user system elapsed
305.95 2.15 363.01
|