PureCN.Rcheck/tests_i386/testthat.Rout
R version 3.5.1 Patched (2018-07-24 r75005) -- "Feather Spray"
Copyright (C) 2018 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 [2018-10-17 04:20:44] Attempted adding gene symbols to intervals. Heuristics have been used to pick symbols for overlapping genes.
WARN [2018-10-17 04:20:44] Chromosome naming style of txdb file (UCSC) was different from interval file (NCBI).
WARN [2018-10-17 04:20:51] Attempted adding gene symbols to intervals. Heuristics have been used to pick symbols for overlapping genes.
FATAL [2018-10-17 04:20:53] tumor.coverage.file and interval.file do not align.
FATAL [2018-10-17 04:20:53]
FATAL [2018-10-17 04:20:53] This is most likely a user error due to invalid input data or
FATAL [2018-10-17 04:20:53] parameters (PureCN 1.10.0).
FATAL [2018-10-17 04:20:53] Need either f or purity and ploidy.
FATAL [2018-10-17 04:20:53]
FATAL [2018-10-17 04:20:53] This is most likely a user error due to invalid input data or
FATAL [2018-10-17 04:20:53] parameters (PureCN 1.10.0).
FATAL [2018-10-17 04:20:53] f not in expected range.
FATAL [2018-10-17 04:20:53]
FATAL [2018-10-17 04:20:53] This is most likely a user error due to invalid input data or
FATAL [2018-10-17 04:20:53] parameters (PureCN 1.10.0).
FATAL [2018-10-17 04:20:53] coverage not in expected range (>=2)
FATAL [2018-10-17 04:20:53]
FATAL [2018-10-17 04:20:53] This is most likely a user error due to invalid input data or
FATAL [2018-10-17 04:20:53] parameters (PureCN 1.10.0).
FATAL [2018-10-17 04:20:53] purity not in expected range.
FATAL [2018-10-17 04:20:53]
FATAL [2018-10-17 04:20:53] This is most likely a user error due to invalid input data or
FATAL [2018-10-17 04:20:53] parameters (PureCN 1.10.0).
FATAL [2018-10-17 04:20:53] ploidy not in expected range.
FATAL [2018-10-17 04:20:53]
FATAL [2018-10-17 04:20:53] This is most likely a user error due to invalid input data or
FATAL [2018-10-17 04:20:53] parameters (PureCN 1.10.0).
FATAL [2018-10-17 04:20:53] cell.fraction not in expected range.
FATAL [2018-10-17 04:20:53]
FATAL [2018-10-17 04:20:53] This is most likely a user error due to invalid input data or
FATAL [2018-10-17 04:20:53] parameters (PureCN 1.10.0).
INFO [2018-10-17 04:20:59] ------------------------------------------------------------
INFO [2018-10-17 04:20:59] PureCN 1.10.0
INFO [2018-10-17 04:20:59] ------------------------------------------------------------
INFO [2018-10-17 04:20:59] Loading coverage files...
INFO [2018-10-17 04:20:59] Mean target coverages: 112X (tumor) 99X (normal).
WARN [2018-10-17 04:20:59] Allosome coverage missing, cannot determine sex.
WARN [2018-10-17 04:20:59] Allosome coverage missing, cannot determine sex.
INFO [2018-10-17 04:20:59] Removing 228 intervals with missing log.ratio.
INFO [2018-10-17 04:20:59] Removing 22 small (< 5bp) targets.
INFO [2018-10-17 04:20:59] Removing 19 targets with low total coverage in normal (< 150.00 reads).
WARN [2018-10-17 04:20:59] No normalDB provided. Provide one for better results.
INFO [2018-10-17 04:20:59] Removing 233 low coverage (< 15.0000X) targets.
INFO [2018-10-17 04:20:59] Using 9547 intervals (9547 on-target, 0 off-target).
INFO [2018-10-17 04:20:59] No off-target intervals. If this is hybrid-capture data, consider adding them.
INFO [2018-10-17 04:20:59] No interval.file provided. Cannot check if data was GC-normalized. Was it?
INFO [2018-10-17 04:20:59] Loading VCF...
INFO [2018-10-17 04:20:59] Found 2331 variants in VCF file.
INFO [2018-10-17 04:20:59] 2282 (97.9%) variants annotated as likely germline (DB INFO flag).
INFO [2018-10-17 04:20:59] LIB-02240e4 is tumor in VCF file.
INFO [2018-10-17 04:20:59] No homozygous variants in VCF, provide unfiltered VCF.
INFO [2018-10-17 04:21:00] Removing 22 non heterozygous (in matched normal) germline SNPs.
INFO [2018-10-17 04:21:00] Initial testing for significant sample cross-contamination: unlikely
INFO [2018-10-17 04:21:01] Removing 48 variants with AF < 0.030 or AF >= 1.000 or less than 4 supporting reads or depth < 15.
INFO [2018-10-17 04:21:01] Removing 16 low quality variants with BQ < 25.
INFO [2018-10-17 04:21:01] Total size of targeted genomic region: 1.63Mb (2.58Mb with 50bp padding).
INFO [2018-10-17 04:21:01] 1.0% of targets contain variants.
INFO [2018-10-17 04:21:01] Removing 2122 variants outside intervals.
INFO [2018-10-17 04:21:01] Found SOMATIC annotation in VCF.
INFO [2018-10-17 04:21:01] Setting somatic prior probabilities for somatic variants to 0.999000 or to 0.000100 otherwise.
INFO [2018-10-17 04:21:01] Found SOMATIC annotation in VCF. Setting mapping bias to 0.973.
INFO [2018-10-17 04:21:01] Sample sex: ?
INFO [2018-10-17 04:21:01] Segmenting data...
INFO [2018-10-17 04:21:01] Loading pre-computed boundaries for DNAcopy...
INFO [2018-10-17 04:21:01] Setting undo.SD parameter to 1.000000.
INFO [2018-10-17 04:21:03] Setting prune.hclust.h parameter to 0.200000.
INFO [2018-10-17 04:21:03] Found 54 segments with median size of 24.88Mb.
INFO [2018-10-17 04:21:03] Using 123 variants.
INFO [2018-10-17 04:21:03] Mean standard deviation of log-ratios: 0.37
INFO [2018-10-17 04:21:03] 2D-grid search of purity and ploidy...
INFO [2018-10-17 04:21:06] Local optima: 0.63/1.9, 0.4/2.2, 0.5/2
INFO [2018-10-17 04:21:06] Testing local optimum 1/3 at purity 0.63 and total ploidy 1.90...
INFO [2018-10-17 04:21:08] Fitting variants for purity 0.65, tumor ploidy 1.74 and contamination 0.01.
INFO [2018-10-17 04:21:09] Optimized purity: 0.65
INFO [2018-10-17 04:21:09] Testing local optimum 2/3 at purity 0.40 and total ploidy 2.20...
INFO [2018-10-17 04:21:10] Fitting variants for purity 0.40, tumor ploidy 2.34 and contamination 0.01.
INFO [2018-10-17 04:21:12] Optimized purity: 0.40
INFO [2018-10-17 04:21:12] Testing local optimum 3/3 at purity 0.50 and total ploidy 2.00...
INFO [2018-10-17 04:21:13] Fitting variants for purity 0.65, tumor ploidy 1.74 and contamination 0.01.
INFO [2018-10-17 04:21:14] Optimized purity: 0.65
INFO [2018-10-17 04:21:14] Done.
INFO [2018-10-17 04:21:14] ------------------------------------------------------------
FATAL [2018-10-17 04:21:16] exclude not a GRanges object.
FATAL [2018-10-17 04:21:16]
FATAL [2018-10-17 04:21:16] This is most likely a user error due to invalid input data or
FATAL [2018-10-17 04:21:16] parameters (PureCN 1.10.0).
FATAL [2018-10-17 04:21:16] callable not a GRanges object.
FATAL [2018-10-17 04:21:16]
FATAL [2018-10-17 04:21:16] This is most likely a user error due to invalid input data or
FATAL [2018-10-17 04:21:16] parameters (PureCN 1.10.0).
WARN [2018-10-17 04:21:29] tumor.coverage.file and interval.file do not align.
FATAL [2018-10-17 04:21:33] No gc_bias column in interval.file.
FATAL [2018-10-17 04:21:33]
FATAL [2018-10-17 04:21:33] This is most likely a user error due to invalid input data or
FATAL [2018-10-17 04:21:33] parameters (PureCN 1.10.0).
FATAL [2018-10-17 04:21:34] Provided coverage is zero, most likely due to a corrupt BAM file.
FATAL [2018-10-17 04:21:34]
FATAL [2018-10-17 04:21:34] This is most likely a user error due to invalid input data or
FATAL [2018-10-17 04:21:34] parameters (PureCN 1.10.0).
INFO [2018-10-17 04:21:41] Reading C:\Users\biocbuild\bbs-3.7-bioc\tmpdir\RtmpCStehs\fileca016315d45.rds...
INFO [2018-10-17 04:21:42] Reading C:\Users\biocbuild\bbs-3.7-bioc\tmpdir\RtmpCStehs\fileca016315d45.rds...
INFO [2018-10-17 04:21:42] Reading C:\Users\biocbuild\bbs-3.7-bioc\tmpdir\RtmpCStehs\fileca016315d45.rds...
INFO [2018-10-17 04:21:42] Reading C:\Users\biocbuild\bbs-3.7-bioc\tmpdir\RtmpCStehs\fileca016315d45.rds...
INFO [2018-10-17 04:21:43] Reading C:\Users\biocbuild\bbs-3.7-bioc\tmpdir\RtmpCStehs\fileca016315d45.rds...
INFO [2018-10-17 04:21:43] Reading C:\Users\biocbuild\bbs-3.7-bioc\tmpdir\RtmpCStehs\fileca016315d45.rds...
FATAL [2018-10-17 04:21:43] Purity or Ploidy not numeric or in expected range.
FATAL [2018-10-17 04:21:43]
FATAL [2018-10-17 04:21:43] This is most likely a user error due to invalid input data or
FATAL [2018-10-17 04:21:43] parameters (PureCN 1.10.0).
INFO [2018-10-17 04:21:43] Reading C:\Users\biocbuild\bbs-3.7-bioc\tmpdir\RtmpCStehs\fileca016315d45.rds...
INFO [2018-10-17 04:21:43] Reading C:\Users\biocbuild\bbs-3.7-bioc\tmpdir\RtmpCStehs\fileca016315d45.rds...
FATAL [2018-10-17 04:21:43] 'Failed' column in
FATAL [2018-10-17 04:21:43] C:\Users\biocbuild\bbs-3.7-bioc\tmpdir\RtmpCStehs\fileca016315d45.csv
FATAL [2018-10-17 04:21:43] not logical(1).
FATAL [2018-10-17 04:21:43]
FATAL [2018-10-17 04:21:43] This is most likely a user error due to invalid input data or
FATAL [2018-10-17 04:21:43] parameters (PureCN 1.10.0).
INFO [2018-10-17 04:21:44] 576 on-target bins with low coverage in all samples.
WARN [2018-10-17 04:21:44] You are likely not using the correct baits file!
WARN [2018-10-17 04:21:44] Allosome coverage missing, cannot determine sex.
WARN [2018-10-17 04:21:44] Allosome coverage missing, cannot determine sex.
INFO [2018-10-17 04:21:44] Processing on-target regions...
INFO [2018-10-17 04:21:44] Removing 930 targets with low coverage in normalDB.
INFO [2018-10-17 04:21:44] Removing 1 targets with zero coverage in more than 3% of normalDB.
INFO [2018-10-17 04:21:48] 576 on-target bins with low coverage in all samples.
WARN [2018-10-17 04:21:48] You are likely not using the correct baits file!
WARN [2018-10-17 04:21:48] Allosome coverage missing, cannot determine sex.
WARN [2018-10-17 04:21:48] Allosome coverage missing, cannot determine sex.
INFO [2018-10-17 04:21:48] Processing on-target regions...
INFO [2018-10-17 04:21:49] Removing 930 targets with low coverage in normalDB.
INFO [2018-10-17 04:21:49] Removing 1 targets with zero coverage in more than 3% of normalDB.
INFO [2018-10-17 04:21:50] 576 on-target bins with low coverage in all samples.
WARN [2018-10-17 04:21:50] You are likely not using the correct baits file!
WARN [2018-10-17 04:21:50] Allosome coverage missing, cannot determine sex.
WARN [2018-10-17 04:21:50] Allosome coverage missing, cannot determine sex.
INFO [2018-10-17 04:21:50] Processing on-target regions...
INFO [2018-10-17 04:21:51] Removing 930 targets with low coverage in normalDB.
INFO [2018-10-17 04:21:51] Removing 1 targets with zero coverage in more than 3% of normalDB.
INFO [2018-10-17 04:21:53] 576 on-target bins with low coverage in all samples.
WARN [2018-10-17 04:21:53] You are likely not using the correct baits file!
WARN [2018-10-17 04:21:53] Allosome coverage missing, cannot determine sex.
WARN [2018-10-17 04:21:53] Allosome coverage missing, cannot determine sex.
FATAL [2018-10-17 04:21:53] Length of normal.coverage.files and sex different
FATAL [2018-10-17 04:21:53]
FATAL [2018-10-17 04:21:53] This is most likely a user error due to invalid input data or
FATAL [2018-10-17 04:21:53] parameters (PureCN 1.10.0).
WARN [2018-10-17 04:21:59] Target intervals were not sorted.
INFO [2018-10-17 04:21:59] 560 on-target bins with low coverage in all samples.
WARN [2018-10-17 04:21:59] You are likely not using the correct baits file!
WARN [2018-10-17 04:21:59] Allosome coverage missing, cannot determine sex.
WARN [2018-10-17 04:22:00] Allosome coverage missing, cannot determine sex.
WARN [2018-10-17 04:22:00] Allosome coverage missing, cannot determine sex.
INFO [2018-10-17 04:22:00] Processing on-target regions...
INFO [2018-10-17 04:22:00] Removing 978 targets with low coverage in normalDB.
INFO [2018-10-17 04:22:00] Removing 11 targets with zero coverage in more than 3% of normalDB.
FATAL [2018-10-17 04:22:03] tumor.coverage.file and normalDB do not align.
FATAL [2018-10-17 04:22:03]
FATAL [2018-10-17 04:22:03] This is most likely a user error due to invalid input data or
FATAL [2018-10-17 04:22:03] parameters (PureCN 1.10.0).
INFO [2018-10-17 04:22:03] ------------------------------------------------------------
INFO [2018-10-17 04:22:03] PureCN 1.10.0
INFO [2018-10-17 04:22:03] ------------------------------------------------------------
INFO [2018-10-17 04:22:03] Loading coverage files...
INFO [2018-10-17 04:22:03] Mean target coverages: 112X (tumor) 99X (normal).
WARN [2018-10-17 04:22:03] Allosome coverage missing, cannot determine sex.
WARN [2018-10-17 04:22:03] Allosome coverage missing, cannot determine sex.
INFO [2018-10-17 04:22:04] Removing 228 intervals with missing log.ratio.
FATAL [2018-10-17 04:22:04] normalDB incompatible with this PureCN version. Please re-run
FATAL [2018-10-17 04:22:04] NormalDB.R.
FATAL [2018-10-17 04:22:04]
FATAL [2018-10-17 04:22:04] This is most likely a user error due to invalid input data or
FATAL [2018-10-17 04:22:04] parameters (PureCN 1.10.0).
FATAL [2018-10-17 04:22:04] normalDB incompatible with this PureCN version. Please re-run
FATAL [2018-10-17 04:22:04] NormalDB.R.
FATAL [2018-10-17 04:22:04]
FATAL [2018-10-17 04:22:04] This is most likely a user error due to invalid input data or
FATAL [2018-10-17 04:22:04] parameters (PureCN 1.10.0).
INFO [2018-10-17 04:22:04] Removing 22 non heterozygous (in matched normal) germline SNPs.
INFO [2018-10-17 04:22:05] Initial testing for significant sample cross-contamination: unlikely
INFO [2018-10-17 04:22:05] Removing 48 variants with AF < 0.030 or AF >= 1.000 or less than 4 supporting reads or depth < 15.
INFO [2018-10-17 04:22:05] Removing 16 low quality variants with BQ < 25.
INFO [2018-10-17 04:22:06] Removing 22 non heterozygous (in matched normal) germline SNPs.
INFO [2018-10-17 04:22:07] Initial testing for significant sample cross-contamination: unlikely
INFO [2018-10-17 04:22:07] Removing 48 variants with AF < 0.030 or AF >= 1.000 or less than 4 supporting reads or depth < 15.
INFO [2018-10-17 04:22:07] Removing 6 blacklisted variants.
INFO [2018-10-17 04:22:07] Removing 16 low quality variants with BQ < 25.
INFO [2018-10-17 04:22:08] Removing 22 non heterozygous (in matched normal) germline SNPs.
INFO [2018-10-17 04:22:08] Initial testing for significant sample cross-contamination: unlikely
INFO [2018-10-17 04:22:09] Removing 48 variants with AF < 0.030 or AF >= 1.000 or less than 4 supporting reads or depth < 15.
INFO [2018-10-17 04:22:09] Removing 16 low quality variants with BQ < 25.
WARN [2018-10-17 04:22:09] MuTect stats file lacks contig and position columns.
INFO [2018-10-17 04:22:09] Removing 22 non heterozygous (in matched normal) germline SNPs.
INFO [2018-10-17 04:22:10] Initial testing for significant sample cross-contamination: unlikely
INFO [2018-10-17 04:22:10] Removing 48 variants with AF < 0.030 or AF >= 1.000 or less than 4 supporting reads or depth < 15.
INFO [2018-10-17 04:22:10] Removing 16 low quality variants with BQ < 25.
WARN [2018-10-17 04:22:10] MuTect stats file and VCF file do not align perfectly. Will remove 2325 unmatched variants.
WARN [2018-10-17 04:22:10] MuTect stats file lacks failure_reasons column. Keeping all variants listed in stats file.
INFO [2018-10-17 04:22:11] Removing 0 non heterozygous (in matched normal) germline SNPs.
INFO [2018-10-17 04:22:11] Removing 0 variants with AF < 0.030 or AF >= 1.000 or less than 4 supporting reads or depth < 15.
INFO [2018-10-17 04:22:12] Removing 0 low quality variants with BQ < 25.
INFO [2018-10-17 04:22:12] Found 11 variants in VCF file.
WARN [2018-10-17 04:22:12] vcf.file has no DB info field for membership in germline databases. Guessing it based on available dbSNP ID.
INFO [2018-10-17 04:22:12] 1 (9.1%) variants annotated as likely germline (DB INFO flag).
INFO [2018-10-17 04:22:16] Found 11 variants in VCF file.
WARN [2018-10-17 04:22:16] vcf.file has no DB info field for membership in germline databases. Guessing it based on available dbSNP ID.
INFO [2018-10-17 04:22:17] 1 (9.1%) variants annotated as likely germline (DB INFO flag).
WARN [2018-10-17 04:22:18] Allosome coverage missing, cannot determine sex.
WARN [2018-10-17 04:22:18] Allosome coverage missing, cannot determine sex.
INFO [2018-10-17 04:22:18] Mean coverages: chrX: 108.58, chrY: 108.58, chr1-22: 104.87.
INFO [2018-10-17 04:22:18] Mean coverages: chrX: 108.58, chrY: 2.17, chr1-22: 104.87.
INFO [2018-10-17 04:22:18] Mean coverages: chrX: 108.58, chrY: 5.17, chr1-22: 104.87.
INFO [2018-10-17 04:22:18] No homozygous variants in VCF, provide unfiltered VCF.
INFO [2018-10-17 04:22:19] No germline variants in VCF.
FATAL [2018-10-17 04:22:19] No solution with id hello
FATAL [2018-10-17 04:22:19]
FATAL [2018-10-17 04:22:19] This is most likely a user error due to invalid input data or
FATAL [2018-10-17 04:22:19] parameters (PureCN 1.10.0).
FATAL [2018-10-17 04:22:19] No solution with id 10
FATAL [2018-10-17 04:22:19]
FATAL [2018-10-17 04:22:19] This is most likely a user error due to invalid input data or
FATAL [2018-10-17 04:22:19] parameters (PureCN 1.10.0).
FATAL [2018-10-17 04:22:20] all.data and w have different lengths.
FATAL [2018-10-17 04:22:20]
FATAL [2018-10-17 04:22:20] This is most likely a user error due to invalid input data or
FATAL [2018-10-17 04:22:20] parameters (PureCN 1.10.0).
WARN [2018-10-17 04:22:21] No mappability scores provided.
WARN [2018-10-17 04:22:21] No reptiming scores provided.
INFO [2018-10-17 04:22:21] Calculating GC-content...
WARN [2018-10-17 04:22:21] No mappability scores provided.
WARN [2018-10-17 04:22:21] No reptiming scores provided.
INFO [2018-10-17 04:22:21] Calculating GC-content...
INFO [2018-10-17 04:22:21] Removing 1 targets overlapping with exclude.
WARN [2018-10-17 04:22:21] No mappability scores provided.
WARN [2018-10-17 04:22:21] No reptiming scores provided.
INFO [2018-10-17 04:22:21] Calculating GC-content...
FATAL [2018-10-17 04:22:21] off.target.padding must be negative.
FATAL [2018-10-17 04:22:21]
FATAL [2018-10-17 04:22:21] This is most likely a user error due to invalid input data or
FATAL [2018-10-17 04:22:21] parameters (PureCN 1.10.0).
FATAL [2018-10-17 04:22:21] Interval coordinates should start at 1, not at 0
FATAL [2018-10-17 04:22:21]
FATAL [2018-10-17 04:22:21] This is most likely a user error due to invalid input data or
FATAL [2018-10-17 04:22:21] parameters (PureCN 1.10.0).
FATAL [2018-10-17 04:22:22] No off-target regions after filtering for mappability and
FATAL [2018-10-17 04:22:22] off.target.padding
FATAL [2018-10-17 04:22:22]
FATAL [2018-10-17 04:22:22] This is most likely a user error due to invalid input data or
FATAL [2018-10-17 04:22:22] parameters (PureCN 1.10.0).
WARN [2018-10-17 04:22:22] No mappability scores provided.
INFO [2018-10-17 04:22:22] Calculating GC-content...
INFO [2018-10-17 04:22:22] Splitting 1 large targets to an average width of 200.
WARN [2018-10-17 04:22:22] No mappability scores provided.
WARN [2018-10-17 04:22:22] No reptiming scores provided.
INFO [2018-10-17 04:22:22] Calculating GC-content...
WARN [2018-10-17 04:22:22] No mappability scores provided.
WARN [2018-10-17 04:22:22] No reptiming scores provided.
INFO [2018-10-17 04:22:22] Calculating GC-content...
INFO [2018-10-17 04:22:23] Tiling off-target regions to an average width of 200000.
WARN [2018-10-17 04:22:23] No mappability scores provided.
WARN [2018-10-17 04:22:23] No reptiming scores provided.
INFO [2018-10-17 04:22:23] Calculating GC-content...
WARN [2018-10-17 04:22:23] Intervals contain off-target regions. Will not change intervals.
WARN [2018-10-17 04:22:23] No mappability scores provided.
WARN [2018-10-17 04:22:23] No reptiming scores provided.
INFO [2018-10-17 04:22:23] Calculating GC-content...
WARN [2018-10-17 04:22:23] No reptiming scores provided.
INFO [2018-10-17 04:22:23] Calculating GC-content...
WARN [2018-10-17 04:22:23] 1 intervals without mapping score.
INFO [2018-10-17 04:22:24] Removing 1 targets with low mappability score (<0.50).
WARN [2018-10-17 04:22:24] No reptiming scores provided.
INFO [2018-10-17 04:22:24] Calculating GC-content...
INFO [2018-10-17 04:22:24] Tiling off-target regions to an average width of 200000.
WARN [2018-10-17 04:22:24] No reptiming scores provided.
INFO [2018-10-17 04:22:24] Calculating GC-content...
WARN [2018-10-17 04:22:25] No mappability scores provided.
WARN [2018-10-17 04:22:25] No reptiming scores provided.
INFO [2018-10-17 04:22:25] Calculating GC-content...
WARN [2018-10-17 04:22:25] Chromosome naming style of interval file (NCBI) was different from reference (UCSC).
WARN [2018-10-17 04:22:25] No mappability scores provided.
WARN [2018-10-17 04:22:25] No reptiming scores provided.
INFO [2018-10-17 04:22:25] Calculating GC-content...
FATAL [2018-10-17 04:22:25] Chromosome naming style of interval file unknown, should be UCSC.
FATAL [2018-10-17 04:22:25]
FATAL [2018-10-17 04:22:25] This is most likely a user error due to invalid input data or
FATAL [2018-10-17 04:22:25] parameters (PureCN 1.10.0).
WARN [2018-10-17 04:22:25] Chromosome naming style of interval file (NCBI) was different from reference (UCSC).
WARN [2018-10-17 04:22:25] Chromosome naming style of mappability file (NCBI) was different from reference (UCSC).
WARN [2018-10-17 04:22:25] No reptiming scores provided.
INFO [2018-10-17 04:22:25] Calculating GC-content...
INFO [2018-10-17 04:22:27] Loading coverage data...
INFO [2018-10-17 04:22:28] Mean target coverages: 71X (tumor) 99X (normal).
INFO [2018-10-17 04:22:28] Mean target coverages: 71X (tumor) 43X (normal).
INFO [2018-10-17 04:22:29] ------------------------------------------------------------
INFO [2018-10-17 04:22:29] PureCN 1.10.0
INFO [2018-10-17 04:22:29] ------------------------------------------------------------
INFO [2018-10-17 04:22:29] Using BiocParallel for parallel optimization.
INFO [2018-10-17 04:22:29] Loading coverage files...
INFO [2018-10-17 04:22:30] Mean target coverages: 112X (tumor) 99X (normal).
WARN [2018-10-17 04:22:30] Allosome coverage missing, cannot determine sex.
WARN [2018-10-17 04:22:30] Allosome coverage missing, cannot determine sex.
INFO [2018-10-17 04:22:30] Removing 228 intervals with missing log.ratio.
INFO [2018-10-17 04:22:30] Removing 22 small (< 5bp) targets.
INFO [2018-10-17 04:22:30] Removing 19 targets with low total coverage in normal (< 150.00 reads).
WARN [2018-10-17 04:22:30] No normalDB provided. Provide one for better results.
INFO [2018-10-17 04:22:30] Removing 233 low coverage (< 15.0000X) targets.
INFO [2018-10-17 04:22:30] Using 9547 intervals (9547 on-target, 0 off-target).
INFO [2018-10-17 04:22:30] No off-target intervals. If this is hybrid-capture data, consider adding them.
INFO [2018-10-17 04:22:30] No interval.file provided. Cannot check if data was GC-normalized. Was it?
INFO [2018-10-17 04:22:30] Sample sex: ?
INFO [2018-10-17 04:22:30] Segmenting data...
INFO [2018-10-17 04:22:30] Target weights found, will use weighted CBS.
INFO [2018-10-17 04:22:30] Loading pre-computed boundaries for DNAcopy...
INFO [2018-10-17 04:22:30] Setting undo.SD parameter to 1.000000.
Setting multi-figure configuration
INFO [2018-10-17 04:22:32] Found 54 segments with median size of 25.38Mb.
INFO [2018-10-17 04:22:32] Mean standard deviation of log-ratios: 0.37
INFO [2018-10-17 04:22:32] 2D-grid search of purity and ploidy...
INFO [2018-10-17 04:22:32] 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 [2018-10-17 04:23:26] Done.
INFO [2018-10-17 04:23:26] ------------------------------------------------------------
INFO [2018-10-17 04:23:26] Reading C:\Users\biocbuild\bbs-3.7-bioc\tmpdir\RtmpCStehs\fileca06b324f0.rds...
FATAL [2018-10-17 04:23:26] runAbsoluteCN was run without a VCF file.
FATAL [2018-10-17 04:23:26]
FATAL [2018-10-17 04:23:26] This is most likely a user error due to invalid input data or
FATAL [2018-10-17 04:23:26] parameters (PureCN 1.10.0).
FATAL [2018-10-17 04:23:26] runAbsoluteCN was run without a VCF file.
FATAL [2018-10-17 04:23:26]
FATAL [2018-10-17 04:23:26] This is most likely a user error due to invalid input data or
FATAL [2018-10-17 04:23:26] parameters (PureCN 1.10.0).
INFO [2018-10-17 04:23:26] ------------------------------------------------------------
INFO [2018-10-17 04:23:26] PureCN 1.10.0
INFO [2018-10-17 04:23:26] ------------------------------------------------------------
INFO [2018-10-17 04:23:26] Loading coverage files...
FATAL [2018-10-17 04:23:26] Need a normal coverage file if log.ratio and seg.file are not provided.
FATAL [2018-10-17 04:23:26]
FATAL [2018-10-17 04:23:26] This is most likely a user error due to invalid input data or
FATAL [2018-10-17 04:23:26] parameters (PureCN 1.10.0).
INFO [2018-10-17 04:23:26] ------------------------------------------------------------
INFO [2018-10-17 04:23:26] PureCN 1.10.0
INFO [2018-10-17 04:23:26] ------------------------------------------------------------
FATAL [2018-10-17 04:23:26] min.ploidy or max.ploidy not within expected range.
FATAL [2018-10-17 04:23:26]
FATAL [2018-10-17 04:23:26] This is most likely a user error due to invalid input data or
FATAL [2018-10-17 04:23:26] parameters (PureCN 1.10.0).
INFO [2018-10-17 04:23:26] ------------------------------------------------------------
INFO [2018-10-17 04:23:26] PureCN 1.10.0
INFO [2018-10-17 04:23:26] ------------------------------------------------------------
FATAL [2018-10-17 04:23:26] min.ploidy or max.ploidy not within expected range.
FATAL [2018-10-17 04:23:26]
FATAL [2018-10-17 04:23:26] This is most likely a user error due to invalid input data or
FATAL [2018-10-17 04:23:26] parameters (PureCN 1.10.0).
INFO [2018-10-17 04:23:26] ------------------------------------------------------------
INFO [2018-10-17 04:23:26] PureCN 1.10.0
INFO [2018-10-17 04:23:26] ------------------------------------------------------------
INFO [2018-10-17 04:23:26] ------------------------------------------------------------
INFO [2018-10-17 04:23:26] PureCN 1.10.0
INFO [2018-10-17 04:23:26] ------------------------------------------------------------
INFO [2018-10-17 04:23:26] ------------------------------------------------------------
INFO [2018-10-17 04:23:26] PureCN 1.10.0
INFO [2018-10-17 04:23:26] ------------------------------------------------------------
FATAL [2018-10-17 04:23:26] test.num.copy not within expected range.
FATAL [2018-10-17 04:23:26]
FATAL [2018-10-17 04:23:26] This is most likely a user error due to invalid input data or
FATAL [2018-10-17 04:23:26] parameters (PureCN 1.10.0).
INFO [2018-10-17 04:23:26] ------------------------------------------------------------
INFO [2018-10-17 04:23:26] PureCN 1.10.0
INFO [2018-10-17 04:23:26] ------------------------------------------------------------
WARN [2018-10-17 04:23:27] test.num.copy outside recommended range.
FATAL [2018-10-17 04:23:27] max.non.clonal not within expected range or format.
FATAL [2018-10-17 04:23:27]
FATAL [2018-10-17 04:23:27] This is most likely a user error due to invalid input data or
FATAL [2018-10-17 04:23:27] parameters (PureCN 1.10.0).
INFO [2018-10-17 04:23:27] ------------------------------------------------------------
INFO [2018-10-17 04:23:27] PureCN 1.10.0
INFO [2018-10-17 04:23:27] ------------------------------------------------------------
WARN [2018-10-17 04:23:27] test.num.copy outside recommended range.
FATAL [2018-10-17 04:23:27] max.non.clonal not within expected range or format.
FATAL [2018-10-17 04:23:27]
FATAL [2018-10-17 04:23:27] This is most likely a user error due to invalid input data or
FATAL [2018-10-17 04:23:27] parameters (PureCN 1.10.0).
INFO [2018-10-17 04:23:27] ------------------------------------------------------------
INFO [2018-10-17 04:23:27] PureCN 1.10.0
INFO [2018-10-17 04:23:27] ------------------------------------------------------------
FATAL [2018-10-17 04:23:27] test.purity not within expected range.
FATAL [2018-10-17 04:23:27]
FATAL [2018-10-17 04:23:27] This is most likely a user error due to invalid input data or
FATAL [2018-10-17 04:23:27] parameters (PureCN 1.10.0).
INFO [2018-10-17 04:23:27] ------------------------------------------------------------
INFO [2018-10-17 04:23:27] PureCN 1.10.0
INFO [2018-10-17 04:23:27] ------------------------------------------------------------
INFO [2018-10-17 04:23:27] Loading coverage files...
FATAL [2018-10-17 04:23:27] Tumor and normal are identical. This won't give any meaningful results
FATAL [2018-10-17 04:23:27] and I'm stopping here.
FATAL [2018-10-17 04:23:27]
FATAL [2018-10-17 04:23:27] This is most likely a user error due to invalid input data or
FATAL [2018-10-17 04:23:27] parameters (PureCN 1.10.0).
INFO [2018-10-17 04:23:27] ------------------------------------------------------------
INFO [2018-10-17 04:23:27] PureCN 1.10.0
INFO [2018-10-17 04:23:27] ------------------------------------------------------------
INFO [2018-10-17 04:23:27] Loading coverage files...
FATAL [2018-10-17 04:23:27] Length of log.ratio different from tumor coverage.
FATAL [2018-10-17 04:23:27]
FATAL [2018-10-17 04:23:27] This is most likely a user error due to invalid input data or
FATAL [2018-10-17 04:23:27] parameters (PureCN 1.10.0).
INFO [2018-10-17 04:23:27] ------------------------------------------------------------
INFO [2018-10-17 04:23:27] PureCN 1.10.0
INFO [2018-10-17 04:23:27] ------------------------------------------------------------
FATAL [2018-10-17 04:23:27] prior.purity must have the same length as test.purity.
FATAL [2018-10-17 04:23:27]
FATAL [2018-10-17 04:23:27] This is most likely a user error due to invalid input data or
FATAL [2018-10-17 04:23:27] parameters (PureCN 1.10.0).
INFO [2018-10-17 04:23:27] ------------------------------------------------------------
INFO [2018-10-17 04:23:27] PureCN 1.10.0
INFO [2018-10-17 04:23:27] ------------------------------------------------------------
FATAL [2018-10-17 04:23:27] min.gof not within expected range or format.
FATAL [2018-10-17 04:23:27]
FATAL [2018-10-17 04:23:27] This is most likely a user error due to invalid input data or
FATAL [2018-10-17 04:23:27] parameters (PureCN 1.10.0).
INFO [2018-10-17 04:23:27] ------------------------------------------------------------
INFO [2018-10-17 04:23:27] PureCN 1.10.0
INFO [2018-10-17 04:23:27] ------------------------------------------------------------
FATAL [2018-10-17 04:23:27] prior.purity not within expected range or format.
FATAL [2018-10-17 04:23:27]
FATAL [2018-10-17 04:23:27] This is most likely a user error due to invalid input data or
FATAL [2018-10-17 04:23:27] parameters (PureCN 1.10.0).
INFO [2018-10-17 04:23:27] ------------------------------------------------------------
INFO [2018-10-17 04:23:27] PureCN 1.10.0
INFO [2018-10-17 04:23:27] ------------------------------------------------------------
FATAL [2018-10-17 04:23:27] prior.purity must add to 1. Sum is 1.5
FATAL [2018-10-17 04:23:27]
FATAL [2018-10-17 04:23:27] This is most likely a user error due to invalid input data or
FATAL [2018-10-17 04:23:27] parameters (PureCN 1.10.0).
INFO [2018-10-17 04:23:27] ------------------------------------------------------------
INFO [2018-10-17 04:23:27] PureCN 1.10.0
INFO [2018-10-17 04:23:27] ------------------------------------------------------------
FATAL [2018-10-17 04:23:27] max.homozygous.loss not within expected range or format.
FATAL [2018-10-17 04:23:27]
FATAL [2018-10-17 04:23:27] This is most likely a user error due to invalid input data or
FATAL [2018-10-17 04:23:27] parameters (PureCN 1.10.0).
INFO [2018-10-17 04:23:27] ------------------------------------------------------------
INFO [2018-10-17 04:23:27] PureCN 1.10.0
INFO [2018-10-17 04:23:27] ------------------------------------------------------------
FATAL [2018-10-17 04:23:28] prior.K not within expected range or format.
FATAL [2018-10-17 04:23:28]
FATAL [2018-10-17 04:23:28] This is most likely a user error due to invalid input data or
FATAL [2018-10-17 04:23:28] parameters (PureCN 1.10.0).
INFO [2018-10-17 04:23:28] ------------------------------------------------------------
INFO [2018-10-17 04:23:28] PureCN 1.10.0
INFO [2018-10-17 04:23:28] ------------------------------------------------------------
FATAL [2018-10-17 04:23:28] prior.contamination not within expected range or format.
FATAL [2018-10-17 04:23:28]
FATAL [2018-10-17 04:23:28] This is most likely a user error due to invalid input data or
FATAL [2018-10-17 04:23:28] parameters (PureCN 1.10.0).
INFO [2018-10-17 04:23:28] ------------------------------------------------------------
INFO [2018-10-17 04:23:28] PureCN 1.10.0
INFO [2018-10-17 04:23:28] ------------------------------------------------------------
FATAL [2018-10-17 04:23:28] Iterations not in the expected range from 10 to 250.
FATAL [2018-10-17 04:23:28]
FATAL [2018-10-17 04:23:28] This is most likely a user error due to invalid input data or
FATAL [2018-10-17 04:23:28] parameters (PureCN 1.10.0).
INFO [2018-10-17 04:23:28] ------------------------------------------------------------
INFO [2018-10-17 04:23:28] PureCN 1.10.0
INFO [2018-10-17 04:23:28] ------------------------------------------------------------
FATAL [2018-10-17 04:23:28] Iterations not in the expected range from 10 to 250.
FATAL [2018-10-17 04:23:28]
FATAL [2018-10-17 04:23:28] This is most likely a user error due to invalid input data or
FATAL [2018-10-17 04:23:28] parameters (PureCN 1.10.0).
INFO [2018-10-17 04:23:28] ------------------------------------------------------------
INFO [2018-10-17 04:23:28] PureCN 1.10.0
INFO [2018-10-17 04:23:28] ------------------------------------------------------------
INFO [2018-10-17 04:23:28] ------------------------------------------------------------
INFO [2018-10-17 04:23:28] PureCN 1.10.0
INFO [2018-10-17 04:23:28] ------------------------------------------------------------
INFO [2018-10-17 04:23:28] Loading coverage files...
FATAL [2018-10-17 04:23:28] Interval files in normal and tumor different.
FATAL [2018-10-17 04:23:28]
FATAL [2018-10-17 04:23:28] This is most likely a user error due to invalid input data or
FATAL [2018-10-17 04:23:28] parameters (PureCN 1.10.0).
INFO [2018-10-17 04:23:28] ------------------------------------------------------------
INFO [2018-10-17 04:23:28] PureCN 1.10.0
INFO [2018-10-17 04:23:28] ------------------------------------------------------------
INFO [2018-10-17 04:23:28] Loading coverage files...
INFO [2018-10-17 04:23:29] Mean target coverages: 112X (tumor) 99X (normal).
WARN [2018-10-17 04:23:29] Allosome coverage missing, cannot determine sex.
WARN [2018-10-17 04:23:29] Allosome coverage missing, cannot determine sex.
INFO [2018-10-17 04:23:29] Removing 228 intervals with missing log.ratio.
INFO [2018-10-17 04:23:29] Removing 22 small (< 5bp) targets.
INFO [2018-10-17 04:23:29] Removing 19 targets with low total coverage in normal (< 150.00 reads).
WARN [2018-10-17 04:23:29] No normalDB provided. Provide one for better results.
INFO [2018-10-17 04:23:29] Removing 233 low coverage (< 15.0000X) targets.
INFO [2018-10-17 04:23:29] Using 9547 intervals (9547 on-target, 0 off-target).
INFO [2018-10-17 04:23:29] No off-target intervals. If this is hybrid-capture data, consider adding them.
INFO [2018-10-17 04:23:29] No interval.file provided. Cannot check if data was GC-normalized. Was it?
INFO [2018-10-17 04:23:29] Loading VCF...
INFO [2018-10-17 04:23:29] Found 2331 variants in VCF file.
INFO [2018-10-17 04:23:29] 2282 (97.9%) variants annotated as likely germline (DB INFO flag).
INFO [2018-10-17 04:23:30] LIB-02240e4 is tumor in VCF file.
INFO [2018-10-17 04:23:30] No homozygous variants in VCF, provide unfiltered VCF.
INFO [2018-10-17 04:23:30] Removing 22 non heterozygous (in matched normal) germline SNPs.
INFO [2018-10-17 04:23:31] Initial testing for significant sample cross-contamination: unlikely
INFO [2018-10-17 04:23:31] Removing 48 variants with AF < 0.030 or AF >= 1.000 or less than 4 supporting reads or depth < 15.
INFO [2018-10-17 04:23:31] Removing 16 low quality variants with BQ < 25.
INFO [2018-10-17 04:23:31] Total size of targeted genomic region: 1.63Mb (2.58Mb with 50bp padding).
INFO [2018-10-17 04:23:31] 1.0% of targets contain variants.
INFO [2018-10-17 04:23:31] Removing 2122 variants outside intervals.
INFO [2018-10-17 04:23:31] Found SOMATIC annotation in VCF.
INFO [2018-10-17 04:23:31] Setting somatic prior probabilities for somatic variants to 0.999000 or to 0.000100 otherwise.
INFO [2018-10-17 04:23:31] Found SOMATIC annotation in VCF. Setting mapping bias to 0.973.
INFO [2018-10-17 04:23:31] Sample sex: ?
INFO [2018-10-17 04:23:31] Segmenting data...
INFO [2018-10-17 04:23:32] Loading pre-computed boundaries for DNAcopy...
INFO [2018-10-17 04:23:32] Setting undo.SD parameter to 1.000000.
INFO [2018-10-17 04:23:33] Setting prune.hclust.h parameter to 0.200000.
INFO [2018-10-17 04:23:33] Found 54 segments with median size of 24.88Mb.
INFO [2018-10-17 04:23:33] Using 123 variants.
INFO [2018-10-17 04:23:33] Mean standard deviation of log-ratios: 0.37
INFO [2018-10-17 04:23:33] 2D-grid search of purity and ploidy...
INFO [2018-10-17 04:23:35] Local optima: 0.62/1.9, 0.3/1.9, 0.32/2
INFO [2018-10-17 04:23:35] Testing local optimum 1/3 at purity 0.62 and total ploidy 1.90...
INFO [2018-10-17 04:23:36] Fitting variants for purity 0.65, tumor ploidy 1.74 and contamination 0.01.
INFO [2018-10-17 04:23:38] Optimized purity: 0.65
INFO [2018-10-17 04:23:38] Testing local optimum 2/3 at purity 0.30 and total ploidy 1.90...
INFO [2018-10-17 04:23:38] Recalibrating log-ratios...
INFO [2018-10-17 04:23:38] Testing local optimum 2/3 at purity 0.30 and total ploidy 1.90...
INFO [2018-10-17 04:23:39] Recalibrating log-ratios...
INFO [2018-10-17 04:23:39] Testing local optimum 2/3 at purity 0.30 and total ploidy 1.90...
INFO [2018-10-17 04:23:39] Recalibrating log-ratios...
INFO [2018-10-17 04:23:39] Testing local optimum 2/3 at purity 0.30 and total ploidy 1.90...
INFO [2018-10-17 04:23:39] Testing local optimum 3/3 at purity 0.32 and total ploidy 2.00...
INFO [2018-10-17 04:23:40] Recalibrating log-ratios...
INFO [2018-10-17 04:23:40] Testing local optimum 3/3 at purity 0.32 and total ploidy 2.00...
INFO [2018-10-17 04:23:40] Recalibrating log-ratios...
INFO [2018-10-17 04:23:40] Testing local optimum 3/3 at purity 0.32 and total ploidy 2.00...
INFO [2018-10-17 04:23:41] Recalibrating log-ratios...
INFO [2018-10-17 04:23:41] Testing local optimum 3/3 at purity 0.32 and total ploidy 2.00...
INFO [2018-10-17 04:23:41] Done.
INFO [2018-10-17 04:23:41] ------------------------------------------------------------
INFO [2018-10-17 04:23:41] ------------------------------------------------------------
INFO [2018-10-17 04:23:41] PureCN 1.10.0
INFO [2018-10-17 04:23:41] ------------------------------------------------------------
INFO [2018-10-17 04:23:41] Loading coverage files...
INFO [2018-10-17 04:23:42] Mean target coverages: 112X (tumor) 99X (normal).
WARN [2018-10-17 04:23:42] Allosome coverage missing, cannot determine sex.
WARN [2018-10-17 04:23:42] Allosome coverage missing, cannot determine sex.
INFO [2018-10-17 04:23:42] Removing 228 intervals with missing log.ratio.
INFO [2018-10-17 04:23:42] Removing 22 small (< 5bp) targets.
INFO [2018-10-17 04:23:42] Removing 19 targets with low total coverage in normal (< 150.00 reads).
WARN [2018-10-17 04:23:42] No normalDB provided. Provide one for better results.
INFO [2018-10-17 04:23:42] Removing 233 low coverage (< 15.0000X) targets.
INFO [2018-10-17 04:23:42] Using 9547 intervals (9547 on-target, 0 off-target).
INFO [2018-10-17 04:23:42] No off-target intervals. If this is hybrid-capture data, consider adding them.
INFO [2018-10-17 04:23:42] No interval.file provided. Cannot check if data was GC-normalized. Was it?
INFO [2018-10-17 04:23:42] Loading VCF...
INFO [2018-10-17 04:23:42] Found 2331 variants in VCF file.
WARN [2018-10-17 04:23:42] DB INFO flag contains NAs
INFO [2018-10-17 04:23:42] 2277 (97.7%) variants annotated as likely germline (DB INFO flag).
WARN [2018-10-17 04:23:42] DP FORMAT field contains NAs. Removing 3 variants.
INFO [2018-10-17 04:23:42] LIB-02240e4 is tumor in VCF file.
INFO [2018-10-17 04:23:43] No homozygous variants in VCF, provide unfiltered VCF.
INFO [2018-10-17 04:23:44] Initial testing for significant sample cross-contamination: unlikely
INFO [2018-10-17 04:23:44] Removing 62 variants with AF < 0.030 or AF >= 0.970 or less than 4 supporting reads or depth < 15.
INFO [2018-10-17 04:23:44] Removing 16 low quality variants with BQ < 25.
INFO [2018-10-17 04:23:44] Total size of targeted genomic region: 1.63Mb (2.58Mb with 50bp padding).
INFO [2018-10-17 04:23:44] 1.0% of targets contain variants.
INFO [2018-10-17 04:23:44] Setting somatic prior probabilities for dbSNP hits to 0.000500 or to 0.500000 otherwise.
INFO [2018-10-17 04:23:44] Sample sex: ?
INFO [2018-10-17 04:23:44] Segmenting data...
INFO [2018-10-17 04:23:45] Loading pre-computed boundaries for DNAcopy...
INFO [2018-10-17 04:23:45] Setting undo.SD parameter to 1.000000.
Setting multi-figure configuration
INFO [2018-10-17 04:23:47] Setting prune.hclust.h parameter to 0.200000.
INFO [2018-10-17 04:23:47] Found 59 segments with median size of 17.67Mb.
INFO [2018-10-17 04:23:47] Removing 280 variants outside segments.
INFO [2018-10-17 04:23:48] Using 1970 variants.
INFO [2018-10-17 04:23:48] Mean standard deviation of log-ratios: 0.37
INFO [2018-10-17 04:23:48] 2D-grid search of purity and ploidy...
INFO [2018-10-17 04:23:49] Local optima: 0.62/1.9, 0.3/1.9, 0.32/2
INFO [2018-10-17 04:23:49] Testing local optimum 1/3 at purity 0.62 and total ploidy 1.90...
INFO [2018-10-17 04:23:51] Fitting variants for purity 0.65, tumor ploidy 1.74 and contamination 0.01.
INFO [2018-10-17 04:24:04] Optimized purity: 0.65
INFO [2018-10-17 04:24:04] Testing local optimum 2/3 at purity 0.30 and total ploidy 1.90...
INFO [2018-10-17 04:24:05] Recalibrating log-ratios...
INFO [2018-10-17 04:24:05] Testing local optimum 2/3 at purity 0.30 and total ploidy 1.90...
INFO [2018-10-17 04:24:05] Recalibrating log-ratios...
INFO [2018-10-17 04:24:05] Testing local optimum 2/3 at purity 0.30 and total ploidy 1.90...
INFO [2018-10-17 04:24:06] Recalibrating log-ratios...
INFO [2018-10-17 04:24:06] Testing local optimum 2/3 at purity 0.30 and total ploidy 1.90...
INFO [2018-10-17 04:24:06] Testing local optimum 3/3 at purity 0.32 and total ploidy 2.00...
INFO [2018-10-17 04:24:06] Recalibrating log-ratios...
INFO [2018-10-17 04:24:06] Testing local optimum 3/3 at purity 0.32 and total ploidy 2.00...
INFO [2018-10-17 04:24:07] Recalibrating log-ratios...
INFO [2018-10-17 04:24:07] Testing local optimum 3/3 at purity 0.32 and total ploidy 2.00...
INFO [2018-10-17 04:24:07] Recalibrating log-ratios...
INFO [2018-10-17 04:24:07] Testing local optimum 3/3 at purity 0.32 and total ploidy 2.00...
INFO [2018-10-17 04:24:07] Done.
INFO [2018-10-17 04:24:07] ------------------------------------------------------------
INFO [2018-10-17 04:24:08] ------------------------------------------------------------
INFO [2018-10-17 04:24:08] PureCN 1.10.0
INFO [2018-10-17 04:24:08] ------------------------------------------------------------
INFO [2018-10-17 04:24:08] Loading coverage files...
INFO [2018-10-17 04:24:08] Mean target coverages: 112X (tumor) 99X (normal).
WARN [2018-10-17 04:24:08] Allosome coverage missing, cannot determine sex.
WARN [2018-10-17 04:24:08] Allosome coverage missing, cannot determine sex.
INFO [2018-10-17 04:24:08] No Gene column in interval.file. You won't get gene-level calls.
INFO [2018-10-17 04:24:09] Removing 228 intervals with missing log.ratio.
INFO [2018-10-17 04:24:09] Removing 15 low/high GC targets.
INFO [2018-10-17 04:24:09] Removing 21 small (< 5bp) targets.
INFO [2018-10-17 04:24:09] Removing 19 targets with low total coverage in normal (< 150.00 reads).
WARN [2018-10-17 04:24:09] No normalDB provided. Provide one for better results.
INFO [2018-10-17 04:24:09] Removing 232 low coverage (< 15.0000X) targets.
INFO [2018-10-17 04:24:09] Using 9534 intervals (9534 on-target, 0 off-target).
INFO [2018-10-17 04:24:09] No off-target intervals. If this is hybrid-capture data, consider adding them.
INFO [2018-10-17 04:24:09] AT/GC dropout: 1.03 (tumor), 1.03 (normal).
INFO [2018-10-17 04:24:09] Loading VCF...
INFO [2018-10-17 04:24:09] Found 2331 variants in VCF file.
INFO [2018-10-17 04:24:09] 2282 (97.9%) variants annotated as likely germline (DB INFO flag).
INFO [2018-10-17 04:24:09] LIB-02240e4 is tumor in VCF file.
INFO [2018-10-17 04:24:09] No homozygous variants in VCF, provide unfiltered VCF.
INFO [2018-10-17 04:24:10] Removing 22 non heterozygous (in matched normal) germline SNPs.
INFO [2018-10-17 04:24:10] Initial testing for significant sample cross-contamination: unlikely
INFO [2018-10-17 04:24:11] Removing 48 variants with AF < 0.030 or AF >= Inf or less than 4 supporting reads or depth < 15.
INFO [2018-10-17 04:24:11] Removing 16 low quality variants with BQ < 25.
INFO [2018-10-17 04:24:11] Total size of targeted genomic region: 1.63Mb (2.58Mb with 50bp padding).
INFO [2018-10-17 04:24:11] 1.0% of targets contain variants.
INFO [2018-10-17 04:24:11] Removing 2122 variants outside intervals.
INFO [2018-10-17 04:24:11] Found SOMATIC annotation in VCF.
INFO [2018-10-17 04:24:11] Setting somatic prior probabilities for somatic variants to 0.999000 or to 0.000100 otherwise.
INFO [2018-10-17 04:24:11] Found SOMATIC annotation in VCF. Setting mapping bias to 0.973.
INFO [2018-10-17 04:24:11] Sample sex: ?
INFO [2018-10-17 04:24:11] Segmenting data...
INFO [2018-10-17 04:24:11] Loading pre-computed boundaries for DNAcopy...
INFO [2018-10-17 04:24:11] Setting undo.SD parameter to 1.000000.
Setting multi-figure configuration
INFO [2018-10-17 04:24:13] Setting prune.hclust.h parameter to 0.200000.
INFO [2018-10-17 04:24:14] Found 54 segments with median size of 24.88Mb.
INFO [2018-10-17 04:24:14] Using 123 variants.
INFO [2018-10-17 04:24:14] Mean standard deviation of log-ratios: 0.37
INFO [2018-10-17 04:24:14] 2D-grid search of purity and ploidy...
INFO [2018-10-17 04:24:14] 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 [2018-10-17 04:24:14] Testing local optimum 1/6 at purity 0.65 and total ploidy 1.90...
INFO [2018-10-17 04:24:15] Fitting variants for purity 0.65, tumor ploidy 1.74 and contamination 0.01.
INFO [2018-10-17 04:24:16] Optimized purity: 0.65
INFO [2018-10-17 04:24:16] Testing local optimum 2/6 at purity 0.33 and total ploidy 1.90...
INFO [2018-10-17 04:24:16] Recalibrating log-ratios...
INFO [2018-10-17 04:24:16] Testing local optimum 2/6 at purity 0.33 and total ploidy 1.90...
INFO [2018-10-17 04:24:17] Recalibrating log-ratios...
INFO [2018-10-17 04:24:17] Testing local optimum 2/6 at purity 0.33 and total ploidy 1.90...
INFO [2018-10-17 04:24:17] Recalibrating log-ratios...
INFO [2018-10-17 04:24:17] Testing local optimum 2/6 at purity 0.33 and total ploidy 1.90...
INFO [2018-10-17 04:24:17] Testing local optimum 3/6 at purity 0.38 and total ploidy 2.20...
INFO [2018-10-17 04:24:19] Fitting variants for purity 0.40, tumor ploidy 2.34 and contamination 0.01.
INFO [2018-10-17 04:24:20] Optimized purity: 0.40
INFO [2018-10-17 04:24:20] Testing local optimum 4/6 at purity 0.27 and total ploidy 2.10...
INFO [2018-10-17 04:24:20] Recalibrating log-ratios...
INFO [2018-10-17 04:24:20] Testing local optimum 4/6 at purity 0.27 and total ploidy 2.10...
INFO [2018-10-17 04:24:21] Recalibrating log-ratios...
INFO [2018-10-17 04:24:21] Testing local optimum 4/6 at purity 0.27 and total ploidy 2.10...
INFO [2018-10-17 04:24:21] Recalibrating log-ratios...
INFO [2018-10-17 04:24:21] Testing local optimum 4/6 at purity 0.27 and total ploidy 2.10...
INFO [2018-10-17 04:24:21] Testing local optimum 5/6 at purity 0.29 and total ploidy 2.00...
INFO [2018-10-17 04:24:22] Recalibrating log-ratios...
INFO [2018-10-17 04:24:22] Testing local optimum 5/6 at purity 0.29 and total ploidy 2.00...
INFO [2018-10-17 04:24:22] Recalibrating log-ratios...
INFO [2018-10-17 04:24:22] Testing local optimum 5/6 at purity 0.29 and total ploidy 2.00...
INFO [2018-10-17 04:24:22] Recalibrating log-ratios...
INFO [2018-10-17 04:24:22] Testing local optimum 5/6 at purity 0.29 and total ploidy 2.00...
INFO [2018-10-17 04:24:23] Testing local optimum 6/6 at purity 0.78 and total ploidy 2.00...
INFO [2018-10-17 04:24:24] Fitting variants for purity 0.65, tumor ploidy 1.74 and contamination 0.01.
INFO [2018-10-17 04:24:25] Optimized purity: 0.65
INFO [2018-10-17 04:24:25] Done.
INFO [2018-10-17 04:24:25] ------------------------------------------------------------
FATAL [2018-10-17 04:24:25] This function requires gene-level calls. Please add a column 'Gene'
FATAL [2018-10-17 04:24:25] containing gene symbols to the interval.file.
FATAL [2018-10-17 04:24:25]
FATAL [2018-10-17 04:24:25] This is most likely a user error due to invalid input data or
FATAL [2018-10-17 04:24:25] parameters (PureCN 1.10.0).
INFO [2018-10-17 04:24:26] ------------------------------------------------------------
INFO [2018-10-17 04:24:26] PureCN 1.10.0
INFO [2018-10-17 04:24:26] ------------------------------------------------------------
INFO [2018-10-17 04:24:26] Loading coverage files...
INFO [2018-10-17 04:24:27] Mean target coverages: 112X (tumor) 99X (normal).
WARN [2018-10-17 04:24:27] Allosome coverage missing, cannot determine sex.
WARN [2018-10-17 04:24:27] Allosome coverage missing, cannot determine sex.
INFO [2018-10-17 04:24:27] Removing 228 intervals with missing log.ratio.
INFO [2018-10-17 04:24:27] Removing 22 small (< 5bp) targets.
INFO [2018-10-17 04:24:27] Removing 19 targets with low total coverage in normal (< 150.00 reads).
WARN [2018-10-17 04:24:27] No normalDB provided. Provide one for better results.
INFO [2018-10-17 04:24:27] Removing 233 low coverage (< 15.0000X) targets.
INFO [2018-10-17 04:24:27] Using 9547 intervals (9547 on-target, 0 off-target).
INFO [2018-10-17 04:24:27] No off-target intervals. If this is hybrid-capture data, consider adding them.
INFO [2018-10-17 04:24:27] No interval.file provided. Cannot check if data was GC-normalized. Was it?
INFO [2018-10-17 04:24:27] Loading VCF...
INFO [2018-10-17 04:24:27] Found 2331 variants in VCF file.
INFO [2018-10-17 04:24:27] 2282 (97.9%) variants annotated as likely germline (DB INFO flag).
FATAL [2018-10-17 04:24:27] Different chromosome names in coverage and VCF.
FATAL [2018-10-17 04:24:27]
FATAL [2018-10-17 04:24:27] This is most likely a user error due to invalid input data or
FATAL [2018-10-17 04:24:27] parameters (PureCN 1.10.0).
INFO [2018-10-17 04:24:27] ------------------------------------------------------------
INFO [2018-10-17 04:24:27] PureCN 1.10.0
INFO [2018-10-17 04:24:27] ------------------------------------------------------------
INFO [2018-10-17 04:24:27] Loading coverage files...
INFO [2018-10-17 04:24:27] Mean target coverages: 112X (tumor) 99X (normal).
WARN [2018-10-17 04:24:27] Allosome coverage missing, cannot determine sex.
WARN [2018-10-17 04:24:27] Allosome coverage missing, cannot determine sex.
FATAL [2018-10-17 04:24:27] tumor.coverage.file and interval.file do not align.
FATAL [2018-10-17 04:24:27]
FATAL [2018-10-17 04:24:27] This is most likely a user error due to invalid input data or
FATAL [2018-10-17 04:24:27] parameters (PureCN 1.10.0).
INFO [2018-10-17 04:24:28] ------------------------------------------------------------
INFO [2018-10-17 04:24:28] PureCN 1.10.0
INFO [2018-10-17 04:24:28] ------------------------------------------------------------
INFO [2018-10-17 04:24:28] Loading coverage files...
INFO [2018-10-17 04:24:28] Mean target coverages: 112X (tumor) 99X (normal).
WARN [2018-10-17 04:24:28] Allosome coverage missing, cannot determine sex.
WARN [2018-10-17 04:24:28] Allosome coverage missing, cannot determine sex.
INFO [2018-10-17 04:24:29] Removing 228 intervals with missing log.ratio.
INFO [2018-10-17 04:24:29] Removing 15 low/high GC targets.
INFO [2018-10-17 04:24:29] Removing 21 small (< 5bp) targets.
INFO [2018-10-17 04:24:29] Removing 19 targets with low total coverage in normal (< 150.00 reads).
WARN [2018-10-17 04:24:29] No normalDB provided. Provide one for better results.
INFO [2018-10-17 04:24:29] Removing 232 low coverage (< 15.0000X) targets.
INFO [2018-10-17 04:24:29] Using 9534 intervals (9534 on-target, 0 off-target).
INFO [2018-10-17 04:24:29] No off-target intervals. If this is hybrid-capture data, consider adding them.
INFO [2018-10-17 04:24:29] AT/GC dropout: 1.03 (tumor), 1.03 (normal).
INFO [2018-10-17 04:24:29] Loading VCF...
INFO [2018-10-17 04:24:29] Found 2331 variants in VCF file.
INFO [2018-10-17 04:24:29] 2282 (97.9%) variants annotated as likely germline (DB INFO flag).
INFO [2018-10-17 04:24:30] LIB-02240e4 is tumor in VCF file.
INFO [2018-10-17 04:24:30] No homozygous variants in VCF, provide unfiltered VCF.
INFO [2018-10-17 04:24:30] Removing 22 non heterozygous (in matched normal) germline SNPs.
INFO [2018-10-17 04:24:31] Initial testing for significant sample cross-contamination: unlikely
INFO [2018-10-17 04:24:31] Removing 48 variants with AF < 0.030 or AF >= 1.000 or less than 4 supporting reads or depth < 15.
INFO [2018-10-17 04:24:31] Removing 16 low quality variants with BQ < 25.
INFO [2018-10-17 04:24:31] Total size of targeted genomic region: 1.63Mb (2.58Mb with 50bp padding).
INFO [2018-10-17 04:24:31] 1.0% of targets contain variants.
INFO [2018-10-17 04:24:31] Removing 2122 variants outside intervals.
INFO [2018-10-17 04:24:31] Found SOMATIC annotation in VCF.
INFO [2018-10-17 04:24:31] Setting somatic prior probabilities for somatic variants to 0.999000 or to 0.000100 otherwise.
INFO [2018-10-17 04:24:32] Found SOMATIC annotation in VCF. Setting mapping bias to 0.973.
INFO [2018-10-17 04:24:32] Sample sex: ?
INFO [2018-10-17 04:24:32] Segmenting data...
INFO [2018-10-17 04:24:32] Loading pre-computed boundaries for DNAcopy...
INFO [2018-10-17 04:24:32] Setting undo.SD parameter to 1.000000.
INFO [2018-10-17 04:24:33] Setting prune.hclust.h parameter to 0.200000.
INFO [2018-10-17 04:24:33] Found 54 segments with median size of 24.88Mb.
INFO [2018-10-17 04:24:33] Using 123 variants.
INFO [2018-10-17 04:24:33] Mean standard deviation of log-ratios: 0.37
INFO [2018-10-17 04:24:33] 2D-grid search of purity and ploidy...
INFO [2018-10-17 04:24:34] Local optima: 0.63/1.9, 0.4/2.1
INFO [2018-10-17 04:24:34] Testing local optimum 1/2 at purity 0.63 and total ploidy 1.90...
INFO [2018-10-17 04:24:36] Fitting variants for purity 0.65, tumor ploidy 1.74 and contamination 0.01.
INFO [2018-10-17 04:24:37] Optimized purity: 0.65
INFO [2018-10-17 04:24:37] Testing local optimum 2/2 at purity 0.40 and total ploidy 2.10...
INFO [2018-10-17 04:24:37] Recalibrating log-ratios...
INFO [2018-10-17 04:24:37] Testing local optimum 2/2 at purity 0.40 and total ploidy 2.10...
INFO [2018-10-17 04:24:38] Recalibrating log-ratios...
INFO [2018-10-17 04:24:38] Testing local optimum 2/2 at purity 0.40 and total ploidy 2.10...
INFO [2018-10-17 04:24:38] Recalibrating log-ratios...
INFO [2018-10-17 04:24:38] Testing local optimum 2/2 at purity 0.40 and total ploidy 2.10...
INFO [2018-10-17 04:24:39] Done.
INFO [2018-10-17 04:24:39] ------------------------------------------------------------
INFO [2018-10-17 04:24:39] ------------------------------------------------------------
INFO [2018-10-17 04:24:39] PureCN 1.10.0
INFO [2018-10-17 04:24:39] ------------------------------------------------------------
INFO [2018-10-17 04:24:39] Loading coverage files...
INFO [2018-10-17 04:24:39] Mean target coverages: 112X (tumor) 99X (normal).
WARN [2018-10-17 04:24:39] Allosome coverage missing, cannot determine sex.
WARN [2018-10-17 04:24:39] Allosome coverage missing, cannot determine sex.
INFO [2018-10-17 04:24:40] Removing 228 intervals with missing log.ratio.
INFO [2018-10-17 04:24:40] Removing 15 low/high GC targets.
INFO [2018-10-17 04:24:40] Removing 21 small (< 5bp) targets.
INFO [2018-10-17 04:24:40] Removing 19 targets with low total coverage in normal (< 150.00 reads).
WARN [2018-10-17 04:24:40] No normalDB provided. Provide one for better results.
INFO [2018-10-17 04:24:40] Removing 232 low coverage (< 15.0000X) targets.
INFO [2018-10-17 04:24:40] Using 9534 intervals (9534 on-target, 0 off-target).
INFO [2018-10-17 04:24:40] No off-target intervals. If this is hybrid-capture data, consider adding them.
INFO [2018-10-17 04:24:40] AT/GC dropout: 1.03 (tumor), 1.03 (normal).
INFO [2018-10-17 04:24:40] Loading VCF...
INFO [2018-10-17 04:24:40] Found 2331 variants in VCF file.
INFO [2018-10-17 04:24:40] 2282 (97.9%) variants annotated as likely germline (DB INFO flag).
INFO [2018-10-17 04:24:41] LIB-02240e4 is tumor in VCF file.
INFO [2018-10-17 04:24:41] No homozygous variants in VCF, provide unfiltered VCF.
INFO [2018-10-17 04:24:41] Removing 22 non heterozygous (in matched normal) germline SNPs.
INFO [2018-10-17 04:24:42] Initial testing for significant sample cross-contamination: unlikely
INFO [2018-10-17 04:24:42] Removing 48 variants with AF < 0.030 or AF >= 1.000 or less than 4 supporting reads or depth < 15.
INFO [2018-10-17 04:24:42] Removing 16 low quality variants with BQ < 25.
INFO [2018-10-17 04:24:42] Total size of targeted genomic region: 1.63Mb (2.58Mb with 50bp padding).
INFO [2018-10-17 04:24:42] 1.0% of targets contain variants.
INFO [2018-10-17 04:24:42] Removing 2122 variants outside intervals.
INFO [2018-10-17 04:24:42] Found SOMATIC annotation in VCF.
INFO [2018-10-17 04:24:42] Setting somatic prior probabilities for somatic variants to 0.999000 or to 0.000100 otherwise.
INFO [2018-10-17 04:24:42] Found SOMATIC annotation in VCF. Setting mapping bias to 0.973.
INFO [2018-10-17 04:24:42] Sample sex: ?
INFO [2018-10-17 04:24:42] Segmenting data...
INFO [2018-10-17 04:24:43] Loading pre-computed boundaries for DNAcopy...
INFO [2018-10-17 04:24:43] Setting undo.SD parameter to 1.000000.
INFO [2018-10-17 04:24:44] Setting prune.hclust.h parameter to 0.200000.
INFO [2018-10-17 04:24:44] Found 54 segments with median size of 24.88Mb.
INFO [2018-10-17 04:24:44] Using 123 variants.
INFO [2018-10-17 04:24:44] Mean standard deviation of log-ratios: 0.37
INFO [2018-10-17 04:24:44] 2D-grid search of purity and ploidy...
INFO [2018-10-17 04:24:45] Local optima: 0.63/1.9, 0.4/2.1
INFO [2018-10-17 04:24:45] Testing local optimum 1/2 at purity 0.63 and total ploidy 1.90...
INFO [2018-10-17 04:24:47] Fitting variants for purity 0.65, tumor ploidy 1.74 and contamination 0.01.
INFO [2018-10-17 04:24:48] Optimized purity: 0.65
INFO [2018-10-17 04:24:48] Testing local optimum 2/2 at purity 0.40 and total ploidy 2.10...
INFO [2018-10-17 04:24:49] Recalibrating log-ratios...
INFO [2018-10-17 04:24:49] Testing local optimum 2/2 at purity 0.40 and total ploidy 2.10...
INFO [2018-10-17 04:24:50] Recalibrating log-ratios...
INFO [2018-10-17 04:24:50] Testing local optimum 2/2 at purity 0.40 and total ploidy 2.10...
INFO [2018-10-17 04:24:50] Recalibrating log-ratios...
INFO [2018-10-17 04:24:50] Testing local optimum 2/2 at purity 0.40 and total ploidy 2.10...
INFO [2018-10-17 04:24:51] Done.
INFO [2018-10-17 04:24:51] ------------------------------------------------------------
FATAL [2018-10-17 04:24:51] chr1 not valid chromosome name(s). Valid names are:
FATAL [2018-10-17 04:24:51] 1,2,3,4,5,6,7,8,9,10,11,12,13,14,15,16,17,18,19,20,21,22
FATAL [2018-10-17 04:24:51]
FATAL [2018-10-17 04:24:51] This is most likely a user error due to invalid input data or
FATAL [2018-10-17 04:24:51] parameters (PureCN 1.10.0).
INFO [2018-10-17 04:24:51] ------------------------------------------------------------
INFO [2018-10-17 04:24:51] PureCN 1.10.0
INFO [2018-10-17 04:24:51] ------------------------------------------------------------
INFO [2018-10-17 04:24:51] Loading coverage files...
WARN [2018-10-17 04:24:51] Provided sampleid (Sample2) does not match Sample1 found in segmentation.
WARN [2018-10-17 04:24:51] Allosome coverage missing, cannot determine sex.
WARN [2018-10-17 04:24:51] Allosome coverage missing, cannot determine sex.
INFO [2018-10-17 04:24:51] Removing 10 intervals with missing log.ratio.
INFO [2018-10-17 04:24:51] Using 10039 intervals (10039 on-target, 0 off-target).
INFO [2018-10-17 04:24:51] No off-target intervals. If this is hybrid-capture data, consider adding them.
INFO [2018-10-17 04:24:51] Loading VCF...
INFO [2018-10-17 04:24:51] Found 2331 variants in VCF file.
INFO [2018-10-17 04:24:51] 2282 (97.9%) variants annotated as likely germline (DB INFO flag).
INFO [2018-10-17 04:24:52] LIB-02240e4 is tumor in VCF file.
INFO [2018-10-17 04:24:52] No homozygous variants in VCF, provide unfiltered VCF.
INFO [2018-10-17 04:24:52] Removing 22 non heterozygous (in matched normal) germline SNPs.
INFO [2018-10-17 04:24:53] Initial testing for significant sample cross-contamination: unlikely
INFO [2018-10-17 04:24:53] Removing 48 variants with AF < 0.030 or AF >= 1.000 or less than 4 supporting reads or depth < 15.
INFO [2018-10-17 04:24:53] Removing 16 low quality variants with BQ < 25.
INFO [2018-10-17 04:24:53] Total size of targeted genomic region: 1.71Mb (2.70Mb with 50bp padding).
INFO [2018-10-17 04:24:53] 1.0% of targets contain variants.
INFO [2018-10-17 04:24:54] Removing 2120 variants outside intervals.
INFO [2018-10-17 04:24:54] Found SOMATIC annotation in VCF.
INFO [2018-10-17 04:24:54] Setting somatic prior probabilities for somatic variants to 0.999000 or to 0.000100 otherwise.
INFO [2018-10-17 04:24:54] Found SOMATIC annotation in VCF. Setting mapping bias to 0.976.
INFO [2018-10-17 04:24:54] Sample sex: ?
INFO [2018-10-17 04:24:54] Segmenting data...
WARN [2018-10-17 04:24:54] Provided sampleid (Sample2) does not match Sample1 found in segmentation.
INFO [2018-10-17 04:24:54] Loading pre-computed boundaries for DNAcopy...
INFO [2018-10-17 04:24:54] Setting undo.SD parameter to 0.000000.
Setting multi-figure configuration
INFO [2018-10-17 04:24:54] Setting prune.hclust.h parameter to 0.200000.
INFO [2018-10-17 04:24:55] Found 54 segments with median size of 24.88Mb.
INFO [2018-10-17 04:24:55] Using 125 variants.
INFO [2018-10-17 04:24:55] Mean standard deviation of log-ratios: 0.40
INFO [2018-10-17 04:24:55] 2D-grid search of purity and ploidy...
INFO [2018-10-17 04:24:56] Local optima: 0.6/1.9, 0.4/2.1
INFO [2018-10-17 04:24:56] Testing local optimum 1/2 at purity 0.60 and total ploidy 1.90...
INFO [2018-10-17 04:24:57] Fitting variants for purity 0.65, tumor ploidy 1.74 and contamination 0.01.
INFO [2018-10-17 04:24:58] Optimized purity: 0.65
INFO [2018-10-17 04:24:58] Testing local optimum 2/2 at purity 0.40 and total ploidy 2.10...
INFO [2018-10-17 04:24:59] Recalibrating log-ratios...
INFO [2018-10-17 04:24:59] Testing local optimum 2/2 at purity 0.40 and total ploidy 2.10...
INFO [2018-10-17 04:24:59] Recalibrating log-ratios...
INFO [2018-10-17 04:24:59] Testing local optimum 2/2 at purity 0.40 and total ploidy 2.10...
INFO [2018-10-17 04:25:00] Recalibrating log-ratios...
INFO [2018-10-17 04:25:00] Testing local optimum 2/2 at purity 0.40 and total ploidy 2.10...
INFO [2018-10-17 04:25:00] Done.
INFO [2018-10-17 04:25:00] ------------------------------------------------------------
WARN [2018-10-17 04:25:01] Allosome coverage missing, cannot determine sex.
WARN [2018-10-17 04:25:01] Allosome coverage missing, cannot determine sex.
Setting multi-figure configuration
FATAL [2018-10-17 04:25:10] Segmentation file expected with colnames ID, chrom, loc.start, loc.end,
FATAL [2018-10-17 04:25:10] num.mark, seg.mean
FATAL [2018-10-17 04:25:10]
FATAL [2018-10-17 04:25:10] This is most likely a user error due to invalid input data or
FATAL [2018-10-17 04:25:10] parameters (PureCN 1.10.0).
FATAL [2018-10-17 04:25:10] seg.file contains multiple samples and sampleid missing.
FATAL [2018-10-17 04:25:10]
FATAL [2018-10-17 04:25:10] This is most likely a user error due to invalid input data or
FATAL [2018-10-17 04:25:10] parameters (PureCN 1.10.0).
FATAL [2018-10-17 04:25:10] seg.file contains multiple samples and sampleid does not match any.
FATAL [2018-10-17 04:25:10]
FATAL [2018-10-17 04:25:10] This is most likely a user error due to invalid input data or
FATAL [2018-10-17 04:25:10] parameters (PureCN 1.10.0).
WARN [2018-10-17 04:25:11] Allosome coverage missing, cannot determine sex.
WARN [2018-10-17 04:25:11] Allosome coverage missing, cannot determine sex.
Setting multi-figure configuration
WARN [2018-10-17 04:25:21] Provided sampleid (example_tumor.txt) does not match Sample1 found in segmentation.
WARN [2018-10-17 04:25:21] Allosome coverage missing, cannot determine sex.
WARN [2018-10-17 04:25:21] Allosome coverage missing, cannot determine sex.
WARN [2018-10-17 04:25:21] Provided sampleid (example_tumor.txt) does not match Sample1 found in segmentation.
WARN [2018-10-17 04:25:27] Allosome coverage missing, cannot determine sex.
WARN [2018-10-17 04:25:27] Allosome coverage missing, cannot determine sex.
WARN [2018-10-17 04:25:29] No normalDB provided. Provide one for better results.
WARN [2018-10-17 04:25:41] Allosome coverage missing, cannot determine sex.
WARN [2018-10-17 04:25:41] Allosome coverage missing, cannot determine sex.
WARN [2018-10-17 04:25:46] Provided sampleid (Sample.1) does not match Sample1 found in segmentation.
Setting multi-figure configuration
WARN [2018-10-17 04:27:41] Allosome coverage missing, cannot determine sex.
WARN [2018-10-17 04:27:41] Allosome coverage missing, cannot determine sex.
WARN [2018-10-17 04:27:41] No normalDB provided. Provide one for better results.
WARN [2018-10-17 04:27:41] Sampleid looks like a normal in VCF, not like a tumor.
WARN [2018-10-17 04:27:53] Allosome coverage missing, cannot determine sex.
WARN [2018-10-17 04:27:53] Allosome coverage missing, cannot determine sex.
WARN [2018-10-17 04:27:53] No normalDB provided. Provide one for better results.
WARN [2018-10-17 04:28:01] Allosome coverage missing, cannot determine sex.
WARN [2018-10-17 04:28:01] Allosome coverage missing, cannot determine sex.
FATAL [2018-10-17 04:28:01] normalDB not a valid normalDB object. Use createNormalDatabase to
FATAL [2018-10-17 04:28:01] create one.
FATAL [2018-10-17 04:28:01]
FATAL [2018-10-17 04:28:01] This is most likely a user error due to invalid input data or
FATAL [2018-10-17 04:28:01] parameters (PureCN 1.10.0).
WARN [2018-10-17 04:28:01] You are likely not using the correct baits file!
WARN [2018-10-17 04:28:01] Allosome coverage missing, cannot determine sex.
WARN [2018-10-17 04:28:01] Allosome coverage missing, cannot determine sex.
WARN [2018-10-17 04:28:03] Allosome coverage missing, cannot determine sex.
WARN [2018-10-17 04:28:03] Allosome coverage missing, cannot determine sex.
FATAL [2018-10-17 04:28:04] normalDB appears to be empty.
FATAL [2018-10-17 04:28:04]
FATAL [2018-10-17 04:28:04] This is most likely a user error due to invalid input data or
FATAL [2018-10-17 04:28:04] parameters (PureCN 1.10.0).
WARN [2018-10-17 04:28:04] Allosome coverage missing, cannot determine sex.
WARN [2018-10-17 04:28:04] Allosome coverage missing, cannot determine sex.
WARN [2018-10-17 04:28:04] Intervals in coverage and interval.file have conflicting on/off-target annotation.
FATAL [2018-10-17 04:28:34] min.normals must be >=2.
FATAL [2018-10-17 04:28:34]
FATAL [2018-10-17 04:28:34] This is most likely a user error due to invalid input data or
FATAL [2018-10-17 04:28:34] parameters (PureCN 1.10.0).
FATAL [2018-10-17 04:28:46] The normal.panel.vcf.file contains only a single sample.
FATAL [2018-10-17 04:28:46]
FATAL [2018-10-17 04:28:46] This is most likely a user error due to invalid input data or
FATAL [2018-10-17 04:28:46] parameters (PureCN 1.10.0).
== testthat results ===========================================================
OK: 293 SKIPPED: 0 FAILED: 0
>
> proc.time()
user system elapsed
444.50 3.68 502.23
|
PureCN.Rcheck/tests_x64/testthat.Rout
R version 3.5.1 Patched (2018-07-24 r75005) -- "Feather Spray"
Copyright (C) 2018 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 [2018-10-17 04:29:12] Attempted adding gene symbols to intervals. Heuristics have been used to pick symbols for overlapping genes.
WARN [2018-10-17 04:29:12] Chromosome naming style of txdb file (UCSC) was different from interval file (NCBI).
WARN [2018-10-17 04:29:18] Attempted adding gene symbols to intervals. Heuristics have been used to pick symbols for overlapping genes.
FATAL [2018-10-17 04:29:20] tumor.coverage.file and interval.file do not align.
FATAL [2018-10-17 04:29:20]
FATAL [2018-10-17 04:29:20] This is most likely a user error due to invalid input data or
FATAL [2018-10-17 04:29:20] parameters (PureCN 1.10.0).
FATAL [2018-10-17 04:29:20] Need either f or purity and ploidy.
FATAL [2018-10-17 04:29:20]
FATAL [2018-10-17 04:29:20] This is most likely a user error due to invalid input data or
FATAL [2018-10-17 04:29:20] parameters (PureCN 1.10.0).
FATAL [2018-10-17 04:29:20] f not in expected range.
FATAL [2018-10-17 04:29:20]
FATAL [2018-10-17 04:29:20] This is most likely a user error due to invalid input data or
FATAL [2018-10-17 04:29:20] parameters (PureCN 1.10.0).
FATAL [2018-10-17 04:29:20] coverage not in expected range (>=2)
FATAL [2018-10-17 04:29:20]
FATAL [2018-10-17 04:29:20] This is most likely a user error due to invalid input data or
FATAL [2018-10-17 04:29:20] parameters (PureCN 1.10.0).
FATAL [2018-10-17 04:29:20] purity not in expected range.
FATAL [2018-10-17 04:29:20]
FATAL [2018-10-17 04:29:20] This is most likely a user error due to invalid input data or
FATAL [2018-10-17 04:29:20] parameters (PureCN 1.10.0).
FATAL [2018-10-17 04:29:20] ploidy not in expected range.
FATAL [2018-10-17 04:29:20]
FATAL [2018-10-17 04:29:20] This is most likely a user error due to invalid input data or
FATAL [2018-10-17 04:29:20] parameters (PureCN 1.10.0).
FATAL [2018-10-17 04:29:20] cell.fraction not in expected range.
FATAL [2018-10-17 04:29:20]
FATAL [2018-10-17 04:29:20] This is most likely a user error due to invalid input data or
FATAL [2018-10-17 04:29:20] parameters (PureCN 1.10.0).
INFO [2018-10-17 04:29:24] ------------------------------------------------------------
INFO [2018-10-17 04:29:24] PureCN 1.10.0
INFO [2018-10-17 04:29:24] ------------------------------------------------------------
INFO [2018-10-17 04:29:24] Loading coverage files...
INFO [2018-10-17 04:29:24] Mean target coverages: 112X (tumor) 99X (normal).
WARN [2018-10-17 04:29:24] Allosome coverage missing, cannot determine sex.
WARN [2018-10-17 04:29:24] Allosome coverage missing, cannot determine sex.
INFO [2018-10-17 04:29:24] Removing 228 intervals with missing log.ratio.
INFO [2018-10-17 04:29:24] Removing 22 small (< 5bp) targets.
INFO [2018-10-17 04:29:24] Removing 19 targets with low total coverage in normal (< 150.00 reads).
WARN [2018-10-17 04:29:24] No normalDB provided. Provide one for better results.
INFO [2018-10-17 04:29:24] Removing 233 low coverage (< 15.0000X) targets.
INFO [2018-10-17 04:29:24] Using 9547 intervals (9547 on-target, 0 off-target).
INFO [2018-10-17 04:29:24] No off-target intervals. If this is hybrid-capture data, consider adding them.
INFO [2018-10-17 04:29:24] No interval.file provided. Cannot check if data was GC-normalized. Was it?
INFO [2018-10-17 04:29:24] Loading VCF...
INFO [2018-10-17 04:29:24] Found 2331 variants in VCF file.
INFO [2018-10-17 04:29:25] 2282 (97.9%) variants annotated as likely germline (DB INFO flag).
INFO [2018-10-17 04:29:25] LIB-02240e4 is tumor in VCF file.
INFO [2018-10-17 04:29:25] No homozygous variants in VCF, provide unfiltered VCF.
INFO [2018-10-17 04:29:25] Removing 22 non heterozygous (in matched normal) germline SNPs.
INFO [2018-10-17 04:29:27] Initial testing for significant sample cross-contamination: unlikely
INFO [2018-10-17 04:29:27] Removing 48 variants with AF < 0.030 or AF >= 1.000 or less than 4 supporting reads or depth < 15.
INFO [2018-10-17 04:29:27] Removing 16 low quality variants with BQ < 25.
INFO [2018-10-17 04:29:27] Total size of targeted genomic region: 1.63Mb (2.58Mb with 50bp padding).
INFO [2018-10-17 04:29:28] 1.0% of targets contain variants.
INFO [2018-10-17 04:29:28] Removing 2122 variants outside intervals.
INFO [2018-10-17 04:29:28] Found SOMATIC annotation in VCF.
INFO [2018-10-17 04:29:28] Setting somatic prior probabilities for somatic variants to 0.999000 or to 0.000100 otherwise.
INFO [2018-10-17 04:29:28] Found SOMATIC annotation in VCF. Setting mapping bias to 0.973.
INFO [2018-10-17 04:29:28] Sample sex: ?
INFO [2018-10-17 04:29:28] Segmenting data...
INFO [2018-10-17 04:29:28] Loading pre-computed boundaries for DNAcopy...
INFO [2018-10-17 04:29:28] Setting undo.SD parameter to 1.000000.
INFO [2018-10-17 04:29:30] Setting prune.hclust.h parameter to 0.200000.
INFO [2018-10-17 04:29:31] Found 54 segments with median size of 24.88Mb.
INFO [2018-10-17 04:29:31] Using 123 variants.
INFO [2018-10-17 04:29:31] Mean standard deviation of log-ratios: 0.37
INFO [2018-10-17 04:29:31] 2D-grid search of purity and ploidy...
INFO [2018-10-17 04:29:34] Local optima: 0.63/1.9, 0.4/2.2, 0.5/2
INFO [2018-10-17 04:29:34] Testing local optimum 1/3 at purity 0.63 and total ploidy 1.90...
INFO [2018-10-17 04:29:36] Fitting variants for purity 0.65, tumor ploidy 1.74 and contamination 0.01.
INFO [2018-10-17 04:29:37] Optimized purity: 0.65
INFO [2018-10-17 04:29:37] Testing local optimum 2/3 at purity 0.40 and total ploidy 2.20...
INFO [2018-10-17 04:29:39] Fitting variants for purity 0.40, tumor ploidy 2.34 and contamination 0.01.
INFO [2018-10-17 04:29:41] Optimized purity: 0.40
INFO [2018-10-17 04:29:41] Testing local optimum 3/3 at purity 0.50 and total ploidy 2.00...
INFO [2018-10-17 04:29:42] Fitting variants for purity 0.65, tumor ploidy 1.74 and contamination 0.01.
INFO [2018-10-17 04:29:45] Optimized purity: 0.65
INFO [2018-10-17 04:29:45] Done.
INFO [2018-10-17 04:29:45] ------------------------------------------------------------
FATAL [2018-10-17 04:29:47] exclude not a GRanges object.
FATAL [2018-10-17 04:29:47]
FATAL [2018-10-17 04:29:47] This is most likely a user error due to invalid input data or
FATAL [2018-10-17 04:29:47] parameters (PureCN 1.10.0).
FATAL [2018-10-17 04:29:48] callable not a GRanges object.
FATAL [2018-10-17 04:29:48]
FATAL [2018-10-17 04:29:48] This is most likely a user error due to invalid input data or
FATAL [2018-10-17 04:29:48] parameters (PureCN 1.10.0).
WARN [2018-10-17 04:30:02] tumor.coverage.file and interval.file do not align.
FATAL [2018-10-17 04:30:04] No gc_bias column in interval.file.
FATAL [2018-10-17 04:30:04]
FATAL [2018-10-17 04:30:04] This is most likely a user error due to invalid input data or
FATAL [2018-10-17 04:30:04] parameters (PureCN 1.10.0).
FATAL [2018-10-17 04:30:05] Provided coverage is zero, most likely due to a corrupt BAM file.
FATAL [2018-10-17 04:30:05]
FATAL [2018-10-17 04:30:05] This is most likely a user error due to invalid input data or
FATAL [2018-10-17 04:30:05] parameters (PureCN 1.10.0).
INFO [2018-10-17 04:30:16] Reading C:\Users\biocbuild\bbs-3.7-bioc\tmpdir\RtmpwlNiM1\file271c6c8e817.rds...
INFO [2018-10-17 04:30:16] Reading C:\Users\biocbuild\bbs-3.7-bioc\tmpdir\RtmpwlNiM1\file271c6c8e817.rds...
INFO [2018-10-17 04:30:16] Reading C:\Users\biocbuild\bbs-3.7-bioc\tmpdir\RtmpwlNiM1\file271c6c8e817.rds...
INFO [2018-10-17 04:30:16] Reading C:\Users\biocbuild\bbs-3.7-bioc\tmpdir\RtmpwlNiM1\file271c6c8e817.rds...
INFO [2018-10-17 04:30:17] Reading C:\Users\biocbuild\bbs-3.7-bioc\tmpdir\RtmpwlNiM1\file271c6c8e817.rds...
INFO [2018-10-17 04:30:17] Reading C:\Users\biocbuild\bbs-3.7-bioc\tmpdir\RtmpwlNiM1\file271c6c8e817.rds...
FATAL [2018-10-17 04:30:17] Purity or Ploidy not numeric or in expected range.
FATAL [2018-10-17 04:30:17]
FATAL [2018-10-17 04:30:17] This is most likely a user error due to invalid input data or
FATAL [2018-10-17 04:30:17] parameters (PureCN 1.10.0).
INFO [2018-10-17 04:30:17] Reading C:\Users\biocbuild\bbs-3.7-bioc\tmpdir\RtmpwlNiM1\file271c6c8e817.rds...
INFO [2018-10-17 04:30:17] Reading C:\Users\biocbuild\bbs-3.7-bioc\tmpdir\RtmpwlNiM1\file271c6c8e817.rds...
FATAL [2018-10-17 04:30:18] 'Failed' column in
FATAL [2018-10-17 04:30:18] C:\Users\biocbuild\bbs-3.7-bioc\tmpdir\RtmpwlNiM1\file271c6c8e817.csv
FATAL [2018-10-17 04:30:18] not logical(1).
FATAL [2018-10-17 04:30:18]
FATAL [2018-10-17 04:30:18] This is most likely a user error due to invalid input data or
FATAL [2018-10-17 04:30:18] parameters (PureCN 1.10.0).
INFO [2018-10-17 04:30:19] 576 on-target bins with low coverage in all samples.
WARN [2018-10-17 04:30:19] You are likely not using the correct baits file!
WARN [2018-10-17 04:30:19] Allosome coverage missing, cannot determine sex.
WARN [2018-10-17 04:30:19] Allosome coverage missing, cannot determine sex.
INFO [2018-10-17 04:30:19] Processing on-target regions...
INFO [2018-10-17 04:30:20] Removing 930 targets with low coverage in normalDB.
INFO [2018-10-17 04:30:20] Removing 1 targets with zero coverage in more than 3% of normalDB.
INFO [2018-10-17 04:30:27] 576 on-target bins with low coverage in all samples.
WARN [2018-10-17 04:30:27] You are likely not using the correct baits file!
WARN [2018-10-17 04:30:27] Allosome coverage missing, cannot determine sex.
WARN [2018-10-17 04:30:27] Allosome coverage missing, cannot determine sex.
INFO [2018-10-17 04:30:27] Processing on-target regions...
INFO [2018-10-17 04:30:28] Removing 930 targets with low coverage in normalDB.
INFO [2018-10-17 04:30:28] Removing 1 targets with zero coverage in more than 3% of normalDB.
INFO [2018-10-17 04:30:30] 576 on-target bins with low coverage in all samples.
WARN [2018-10-17 04:30:30] You are likely not using the correct baits file!
WARN [2018-10-17 04:30:30] Allosome coverage missing, cannot determine sex.
WARN [2018-10-17 04:30:30] Allosome coverage missing, cannot determine sex.
INFO [2018-10-17 04:30:30] Processing on-target regions...
INFO [2018-10-17 04:30:31] Removing 930 targets with low coverage in normalDB.
INFO [2018-10-17 04:30:31] Removing 1 targets with zero coverage in more than 3% of normalDB.
INFO [2018-10-17 04:30:34] 576 on-target bins with low coverage in all samples.
WARN [2018-10-17 04:30:34] You are likely not using the correct baits file!
WARN [2018-10-17 04:30:34] Allosome coverage missing, cannot determine sex.
WARN [2018-10-17 04:30:34] Allosome coverage missing, cannot determine sex.
FATAL [2018-10-17 04:30:34] Length of normal.coverage.files and sex different
FATAL [2018-10-17 04:30:34]
FATAL [2018-10-17 04:30:34] This is most likely a user error due to invalid input data or
FATAL [2018-10-17 04:30:34] parameters (PureCN 1.10.0).
WARN [2018-10-17 04:30:40] Target intervals were not sorted.
INFO [2018-10-17 04:30:40] 560 on-target bins with low coverage in all samples.
WARN [2018-10-17 04:30:40] You are likely not using the correct baits file!
WARN [2018-10-17 04:30:40] Allosome coverage missing, cannot determine sex.
WARN [2018-10-17 04:30:40] Allosome coverage missing, cannot determine sex.
WARN [2018-10-17 04:30:40] Allosome coverage missing, cannot determine sex.
INFO [2018-10-17 04:30:40] Processing on-target regions...
INFO [2018-10-17 04:30:40] Removing 978 targets with low coverage in normalDB.
INFO [2018-10-17 04:30:40] Removing 11 targets with zero coverage in more than 3% of normalDB.
FATAL [2018-10-17 04:30:43] tumor.coverage.file and normalDB do not align.
FATAL [2018-10-17 04:30:43]
FATAL [2018-10-17 04:30:43] This is most likely a user error due to invalid input data or
FATAL [2018-10-17 04:30:43] parameters (PureCN 1.10.0).
INFO [2018-10-17 04:30:43] ------------------------------------------------------------
INFO [2018-10-17 04:30:43] PureCN 1.10.0
INFO [2018-10-17 04:30:43] ------------------------------------------------------------
INFO [2018-10-17 04:30:43] Loading coverage files...
INFO [2018-10-17 04:30:43] Mean target coverages: 112X (tumor) 99X (normal).
WARN [2018-10-17 04:30:43] Allosome coverage missing, cannot determine sex.
WARN [2018-10-17 04:30:43] Allosome coverage missing, cannot determine sex.
INFO [2018-10-17 04:30:43] Removing 228 intervals with missing log.ratio.
FATAL [2018-10-17 04:30:43] normalDB incompatible with this PureCN version. Please re-run
FATAL [2018-10-17 04:30:43] NormalDB.R.
FATAL [2018-10-17 04:30:43]
FATAL [2018-10-17 04:30:43] This is most likely a user error due to invalid input data or
FATAL [2018-10-17 04:30:43] parameters (PureCN 1.10.0).
FATAL [2018-10-17 04:30:43] normalDB incompatible with this PureCN version. Please re-run
FATAL [2018-10-17 04:30:43] NormalDB.R.
FATAL [2018-10-17 04:30:43]
FATAL [2018-10-17 04:30:43] This is most likely a user error due to invalid input data or
FATAL [2018-10-17 04:30:43] parameters (PureCN 1.10.0).
INFO [2018-10-17 04:30:44] Removing 22 non heterozygous (in matched normal) germline SNPs.
INFO [2018-10-17 04:30:45] Initial testing for significant sample cross-contamination: unlikely
INFO [2018-10-17 04:30:45] Removing 48 variants with AF < 0.030 or AF >= 1.000 or less than 4 supporting reads or depth < 15.
INFO [2018-10-17 04:30:45] Removing 16 low quality variants with BQ < 25.
INFO [2018-10-17 04:30:46] Removing 22 non heterozygous (in matched normal) germline SNPs.
INFO [2018-10-17 04:30:47] Initial testing for significant sample cross-contamination: unlikely
INFO [2018-10-17 04:30:47] Removing 48 variants with AF < 0.030 or AF >= 1.000 or less than 4 supporting reads or depth < 15.
INFO [2018-10-17 04:30:47] Removing 6 blacklisted variants.
INFO [2018-10-17 04:30:47] Removing 16 low quality variants with BQ < 25.
INFO [2018-10-17 04:30:48] Removing 22 non heterozygous (in matched normal) germline SNPs.
INFO [2018-10-17 04:30:49] Initial testing for significant sample cross-contamination: unlikely
INFO [2018-10-17 04:30:49] Removing 48 variants with AF < 0.030 or AF >= 1.000 or less than 4 supporting reads or depth < 15.
INFO [2018-10-17 04:30:49] Removing 16 low quality variants with BQ < 25.
WARN [2018-10-17 04:30:49] MuTect stats file lacks contig and position columns.
INFO [2018-10-17 04:30:50] Removing 22 non heterozygous (in matched normal) germline SNPs.
INFO [2018-10-17 04:30:51] Initial testing for significant sample cross-contamination: unlikely
INFO [2018-10-17 04:30:51] Removing 48 variants with AF < 0.030 or AF >= 1.000 or less than 4 supporting reads or depth < 15.
INFO [2018-10-17 04:30:51] Removing 16 low quality variants with BQ < 25.
WARN [2018-10-17 04:30:51] MuTect stats file and VCF file do not align perfectly. Will remove 2325 unmatched variants.
WARN [2018-10-17 04:30:51] MuTect stats file lacks failure_reasons column. Keeping all variants listed in stats file.
INFO [2018-10-17 04:30:52] Removing 0 non heterozygous (in matched normal) germline SNPs.
INFO [2018-10-17 04:30:53] Removing 0 variants with AF < 0.030 or AF >= 1.000 or less than 4 supporting reads or depth < 15.
INFO [2018-10-17 04:30:53] Removing 0 low quality variants with BQ < 25.
INFO [2018-10-17 04:30:53] Found 11 variants in VCF file.
WARN [2018-10-17 04:30:53] vcf.file has no DB info field for membership in germline databases. Guessing it based on available dbSNP ID.
INFO [2018-10-17 04:30:53] 1 (9.1%) variants annotated as likely germline (DB INFO flag).
INFO [2018-10-17 04:30:59] Found 11 variants in VCF file.
WARN [2018-10-17 04:30:59] vcf.file has no DB info field for membership in germline databases. Guessing it based on available dbSNP ID.
INFO [2018-10-17 04:30:59] 1 (9.1%) variants annotated as likely germline (DB INFO flag).
WARN [2018-10-17 04:31:00] Allosome coverage missing, cannot determine sex.
WARN [2018-10-17 04:31:00] Allosome coverage missing, cannot determine sex.
INFO [2018-10-17 04:31:00] Mean coverages: chrX: 108.58, chrY: 108.58, chr1-22: 104.87.
INFO [2018-10-17 04:31:00] Mean coverages: chrX: 108.58, chrY: 2.17, chr1-22: 104.87.
INFO [2018-10-17 04:31:01] Mean coverages: chrX: 108.58, chrY: 5.17, chr1-22: 104.87.
INFO [2018-10-17 04:31:01] No homozygous variants in VCF, provide unfiltered VCF.
INFO [2018-10-17 04:31:01] No germline variants in VCF.
FATAL [2018-10-17 04:31:01] No solution with id hello
FATAL [2018-10-17 04:31:01]
FATAL [2018-10-17 04:31:01] This is most likely a user error due to invalid input data or
FATAL [2018-10-17 04:31:01] parameters (PureCN 1.10.0).
FATAL [2018-10-17 04:31:01] No solution with id 10
FATAL [2018-10-17 04:31:01]
FATAL [2018-10-17 04:31:01] This is most likely a user error due to invalid input data or
FATAL [2018-10-17 04:31:01] parameters (PureCN 1.10.0).
FATAL [2018-10-17 04:31:03] all.data and w have different lengths.
FATAL [2018-10-17 04:31:03]
FATAL [2018-10-17 04:31:03] This is most likely a user error due to invalid input data or
FATAL [2018-10-17 04:31:03] parameters (PureCN 1.10.0).
WARN [2018-10-17 04:31:03] No mappability scores provided.
WARN [2018-10-17 04:31:03] No reptiming scores provided.
INFO [2018-10-17 04:31:03] Calculating GC-content...
WARN [2018-10-17 04:31:03] No mappability scores provided.
WARN [2018-10-17 04:31:03] No reptiming scores provided.
INFO [2018-10-17 04:31:03] Calculating GC-content...
INFO [2018-10-17 04:31:04] Removing 1 targets overlapping with exclude.
WARN [2018-10-17 04:31:04] No mappability scores provided.
WARN [2018-10-17 04:31:04] No reptiming scores provided.
INFO [2018-10-17 04:31:04] Calculating GC-content...
FATAL [2018-10-17 04:31:04] off.target.padding must be negative.
FATAL [2018-10-17 04:31:04]
FATAL [2018-10-17 04:31:04] This is most likely a user error due to invalid input data or
FATAL [2018-10-17 04:31:04] parameters (PureCN 1.10.0).
FATAL [2018-10-17 04:31:04] Interval coordinates should start at 1, not at 0
FATAL [2018-10-17 04:31:04]
FATAL [2018-10-17 04:31:04] This is most likely a user error due to invalid input data or
FATAL [2018-10-17 04:31:04] parameters (PureCN 1.10.0).
FATAL [2018-10-17 04:31:04] No off-target regions after filtering for mappability and
FATAL [2018-10-17 04:31:04] off.target.padding
FATAL [2018-10-17 04:31:04]
FATAL [2018-10-17 04:31:04] This is most likely a user error due to invalid input data or
FATAL [2018-10-17 04:31:04] parameters (PureCN 1.10.0).
WARN [2018-10-17 04:31:05] No mappability scores provided.
INFO [2018-10-17 04:31:05] Calculating GC-content...
INFO [2018-10-17 04:31:05] Splitting 1 large targets to an average width of 200.
WARN [2018-10-17 04:31:05] No mappability scores provided.
WARN [2018-10-17 04:31:05] No reptiming scores provided.
INFO [2018-10-17 04:31:05] Calculating GC-content...
WARN [2018-10-17 04:31:05] No mappability scores provided.
WARN [2018-10-17 04:31:05] No reptiming scores provided.
INFO [2018-10-17 04:31:05] Calculating GC-content...
INFO [2018-10-17 04:31:05] Tiling off-target regions to an average width of 200000.
WARN [2018-10-17 04:31:06] No mappability scores provided.
WARN [2018-10-17 04:31:06] No reptiming scores provided.
INFO [2018-10-17 04:31:06] Calculating GC-content...
WARN [2018-10-17 04:31:06] Intervals contain off-target regions. Will not change intervals.
WARN [2018-10-17 04:31:06] No mappability scores provided.
WARN [2018-10-17 04:31:06] No reptiming scores provided.
INFO [2018-10-17 04:31:06] Calculating GC-content...
WARN [2018-10-17 04:31:06] No reptiming scores provided.
INFO [2018-10-17 04:31:06] Calculating GC-content...
WARN [2018-10-17 04:31:06] 1 intervals without mapping score.
INFO [2018-10-17 04:31:06] Removing 1 targets with low mappability score (<0.50).
WARN [2018-10-17 04:31:06] No reptiming scores provided.
INFO [2018-10-17 04:31:06] Calculating GC-content...
INFO [2018-10-17 04:31:07] Tiling off-target regions to an average width of 200000.
WARN [2018-10-17 04:31:07] No reptiming scores provided.
INFO [2018-10-17 04:31:07] Calculating GC-content...
WARN [2018-10-17 04:31:08] No mappability scores provided.
WARN [2018-10-17 04:31:08] No reptiming scores provided.
INFO [2018-10-17 04:31:08] Calculating GC-content...
WARN [2018-10-17 04:31:08] Chromosome naming style of interval file (NCBI) was different from reference (UCSC).
WARN [2018-10-17 04:31:08] No mappability scores provided.
WARN [2018-10-17 04:31:08] No reptiming scores provided.
INFO [2018-10-17 04:31:08] Calculating GC-content...
FATAL [2018-10-17 04:31:08] Chromosome naming style of interval file unknown, should be UCSC.
FATAL [2018-10-17 04:31:08]
FATAL [2018-10-17 04:31:08] This is most likely a user error due to invalid input data or
FATAL [2018-10-17 04:31:08] parameters (PureCN 1.10.0).
WARN [2018-10-17 04:31:09] Chromosome naming style of interval file (NCBI) was different from reference (UCSC).
WARN [2018-10-17 04:31:09] Chromosome naming style of mappability file (NCBI) was different from reference (UCSC).
WARN [2018-10-17 04:31:09] No reptiming scores provided.
INFO [2018-10-17 04:31:09] Calculating GC-content...
INFO [2018-10-17 04:31:11] Loading coverage data...
INFO [2018-10-17 04:31:12] Mean target coverages: 71X (tumor) 99X (normal).
INFO [2018-10-17 04:31:12] Mean target coverages: 71X (tumor) 43X (normal).
INFO [2018-10-17 04:31:13] ------------------------------------------------------------
INFO [2018-10-17 04:31:13] PureCN 1.10.0
INFO [2018-10-17 04:31:13] ------------------------------------------------------------
INFO [2018-10-17 04:31:13] Using BiocParallel for parallel optimization.
INFO [2018-10-17 04:31:13] Loading coverage files...
INFO [2018-10-17 04:31:14] Mean target coverages: 112X (tumor) 99X (normal).
WARN [2018-10-17 04:31:14] Allosome coverage missing, cannot determine sex.
WARN [2018-10-17 04:31:14] Allosome coverage missing, cannot determine sex.
INFO [2018-10-17 04:31:14] Removing 228 intervals with missing log.ratio.
INFO [2018-10-17 04:31:14] Removing 22 small (< 5bp) targets.
INFO [2018-10-17 04:31:14] Removing 19 targets with low total coverage in normal (< 150.00 reads).
WARN [2018-10-17 04:31:14] No normalDB provided. Provide one for better results.
INFO [2018-10-17 04:31:14] Removing 233 low coverage (< 15.0000X) targets.
INFO [2018-10-17 04:31:14] Using 9547 intervals (9547 on-target, 0 off-target).
INFO [2018-10-17 04:31:14] No off-target intervals. If this is hybrid-capture data, consider adding them.
INFO [2018-10-17 04:31:14] No interval.file provided. Cannot check if data was GC-normalized. Was it?
INFO [2018-10-17 04:31:14] Sample sex: ?
INFO [2018-10-17 04:31:14] Segmenting data...
INFO [2018-10-17 04:31:14] Target weights found, will use weighted CBS.
INFO [2018-10-17 04:31:14] Loading pre-computed boundaries for DNAcopy...
INFO [2018-10-17 04:31:14] Setting undo.SD parameter to 1.000000.
Setting multi-figure configuration
INFO [2018-10-17 04:31:16] Found 54 segments with median size of 25.38Mb.
INFO [2018-10-17 04:31:16] Mean standard deviation of log-ratios: 0.37
INFO [2018-10-17 04:31:16] 2D-grid search of purity and ploidy...
INFO [2018-10-17 04:31:16] 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 [2018-10-17 04:32:14] Done.
INFO [2018-10-17 04:32:14] ------------------------------------------------------------
INFO [2018-10-17 04:32:14] Reading C:\Users\biocbuild\bbs-3.7-bioc\tmpdir\RtmpwlNiM1\file271c1b8e14e6.rds...
FATAL [2018-10-17 04:32:14] runAbsoluteCN was run without a VCF file.
FATAL [2018-10-17 04:32:14]
FATAL [2018-10-17 04:32:14] This is most likely a user error due to invalid input data or
FATAL [2018-10-17 04:32:14] parameters (PureCN 1.10.0).
FATAL [2018-10-17 04:32:14] runAbsoluteCN was run without a VCF file.
FATAL [2018-10-17 04:32:14]
FATAL [2018-10-17 04:32:14] This is most likely a user error due to invalid input data or
FATAL [2018-10-17 04:32:14] parameters (PureCN 1.10.0).
INFO [2018-10-17 04:32:14] ------------------------------------------------------------
INFO [2018-10-17 04:32:14] PureCN 1.10.0
INFO [2018-10-17 04:32:14] ------------------------------------------------------------
INFO [2018-10-17 04:32:14] Loading coverage files...
FATAL [2018-10-17 04:32:14] Need a normal coverage file if log.ratio and seg.file are not provided.
FATAL [2018-10-17 04:32:14]
FATAL [2018-10-17 04:32:14] This is most likely a user error due to invalid input data or
FATAL [2018-10-17 04:32:14] parameters (PureCN 1.10.0).
INFO [2018-10-17 04:32:14] ------------------------------------------------------------
INFO [2018-10-17 04:32:14] PureCN 1.10.0
INFO [2018-10-17 04:32:14] ------------------------------------------------------------
FATAL [2018-10-17 04:32:14] min.ploidy or max.ploidy not within expected range.
FATAL [2018-10-17 04:32:14]
FATAL [2018-10-17 04:32:14] This is most likely a user error due to invalid input data or
FATAL [2018-10-17 04:32:14] parameters (PureCN 1.10.0).
INFO [2018-10-17 04:32:14] ------------------------------------------------------------
INFO [2018-10-17 04:32:14] PureCN 1.10.0
INFO [2018-10-17 04:32:14] ------------------------------------------------------------
FATAL [2018-10-17 04:32:14] min.ploidy or max.ploidy not within expected range.
FATAL [2018-10-17 04:32:14]
FATAL [2018-10-17 04:32:14] This is most likely a user error due to invalid input data or
FATAL [2018-10-17 04:32:14] parameters (PureCN 1.10.0).
INFO [2018-10-17 04:32:14] ------------------------------------------------------------
INFO [2018-10-17 04:32:14] PureCN 1.10.0
INFO [2018-10-17 04:32:14] ------------------------------------------------------------
INFO [2018-10-17 04:32:14] ------------------------------------------------------------
INFO [2018-10-17 04:32:14] PureCN 1.10.0
INFO [2018-10-17 04:32:14] ------------------------------------------------------------
INFO [2018-10-17 04:32:14] ------------------------------------------------------------
INFO [2018-10-17 04:32:14] PureCN 1.10.0
INFO [2018-10-17 04:32:14] ------------------------------------------------------------
FATAL [2018-10-17 04:32:14] test.num.copy not within expected range.
FATAL [2018-10-17 04:32:14]
FATAL [2018-10-17 04:32:14] This is most likely a user error due to invalid input data or
FATAL [2018-10-17 04:32:14] parameters (PureCN 1.10.0).
INFO [2018-10-17 04:32:14] ------------------------------------------------------------
INFO [2018-10-17 04:32:14] PureCN 1.10.0
INFO [2018-10-17 04:32:14] ------------------------------------------------------------
WARN [2018-10-17 04:32:14] test.num.copy outside recommended range.
FATAL [2018-10-17 04:32:14] max.non.clonal not within expected range or format.
FATAL [2018-10-17 04:32:14]
FATAL [2018-10-17 04:32:14] This is most likely a user error due to invalid input data or
FATAL [2018-10-17 04:32:14] parameters (PureCN 1.10.0).
INFO [2018-10-17 04:32:14] ------------------------------------------------------------
INFO [2018-10-17 04:32:14] PureCN 1.10.0
INFO [2018-10-17 04:32:14] ------------------------------------------------------------
WARN [2018-10-17 04:32:14] test.num.copy outside recommended range.
FATAL [2018-10-17 04:32:14] max.non.clonal not within expected range or format.
FATAL [2018-10-17 04:32:14]
FATAL [2018-10-17 04:32:14] This is most likely a user error due to invalid input data or
FATAL [2018-10-17 04:32:14] parameters (PureCN 1.10.0).
INFO [2018-10-17 04:32:14] ------------------------------------------------------------
INFO [2018-10-17 04:32:14] PureCN 1.10.0
INFO [2018-10-17 04:32:14] ------------------------------------------------------------
FATAL [2018-10-17 04:32:14] test.purity not within expected range.
FATAL [2018-10-17 04:32:14]
FATAL [2018-10-17 04:32:14] This is most likely a user error due to invalid input data or
FATAL [2018-10-17 04:32:14] parameters (PureCN 1.10.0).
INFO [2018-10-17 04:32:14] ------------------------------------------------------------
INFO [2018-10-17 04:32:14] PureCN 1.10.0
INFO [2018-10-17 04:32:14] ------------------------------------------------------------
INFO [2018-10-17 04:32:14] Loading coverage files...
FATAL [2018-10-17 04:32:15] Tumor and normal are identical. This won't give any meaningful results
FATAL [2018-10-17 04:32:15] and I'm stopping here.
FATAL [2018-10-17 04:32:15]
FATAL [2018-10-17 04:32:15] This is most likely a user error due to invalid input data or
FATAL [2018-10-17 04:32:15] parameters (PureCN 1.10.0).
INFO [2018-10-17 04:32:15] ------------------------------------------------------------
INFO [2018-10-17 04:32:15] PureCN 1.10.0
INFO [2018-10-17 04:32:15] ------------------------------------------------------------
INFO [2018-10-17 04:32:15] Loading coverage files...
FATAL [2018-10-17 04:32:15] Length of log.ratio different from tumor coverage.
FATAL [2018-10-17 04:32:15]
FATAL [2018-10-17 04:32:15] This is most likely a user error due to invalid input data or
FATAL [2018-10-17 04:32:15] parameters (PureCN 1.10.0).
INFO [2018-10-17 04:32:15] ------------------------------------------------------------
INFO [2018-10-17 04:32:15] PureCN 1.10.0
INFO [2018-10-17 04:32:15] ------------------------------------------------------------
FATAL [2018-10-17 04:32:15] prior.purity must have the same length as test.purity.
FATAL [2018-10-17 04:32:15]
FATAL [2018-10-17 04:32:15] This is most likely a user error due to invalid input data or
FATAL [2018-10-17 04:32:15] parameters (PureCN 1.10.0).
INFO [2018-10-17 04:32:15] ------------------------------------------------------------
INFO [2018-10-17 04:32:15] PureCN 1.10.0
INFO [2018-10-17 04:32:15] ------------------------------------------------------------
FATAL [2018-10-17 04:32:15] min.gof not within expected range or format.
FATAL [2018-10-17 04:32:15]
FATAL [2018-10-17 04:32:15] This is most likely a user error due to invalid input data or
FATAL [2018-10-17 04:32:15] parameters (PureCN 1.10.0).
INFO [2018-10-17 04:32:15] ------------------------------------------------------------
INFO [2018-10-17 04:32:15] PureCN 1.10.0
INFO [2018-10-17 04:32:15] ------------------------------------------------------------
FATAL [2018-10-17 04:32:15] prior.purity not within expected range or format.
FATAL [2018-10-17 04:32:15]
FATAL [2018-10-17 04:32:15] This is most likely a user error due to invalid input data or
FATAL [2018-10-17 04:32:15] parameters (PureCN 1.10.0).
INFO [2018-10-17 04:32:15] ------------------------------------------------------------
INFO [2018-10-17 04:32:15] PureCN 1.10.0
INFO [2018-10-17 04:32:15] ------------------------------------------------------------
FATAL [2018-10-17 04:32:15] prior.purity must add to 1. Sum is 1.5
FATAL [2018-10-17 04:32:15]
FATAL [2018-10-17 04:32:15] This is most likely a user error due to invalid input data or
FATAL [2018-10-17 04:32:15] parameters (PureCN 1.10.0).
INFO [2018-10-17 04:32:15] ------------------------------------------------------------
INFO [2018-10-17 04:32:15] PureCN 1.10.0
INFO [2018-10-17 04:32:15] ------------------------------------------------------------
FATAL [2018-10-17 04:32:15] max.homozygous.loss not within expected range or format.
FATAL [2018-10-17 04:32:15]
FATAL [2018-10-17 04:32:15] This is most likely a user error due to invalid input data or
FATAL [2018-10-17 04:32:15] parameters (PureCN 1.10.0).
INFO [2018-10-17 04:32:15] ------------------------------------------------------------
INFO [2018-10-17 04:32:15] PureCN 1.10.0
INFO [2018-10-17 04:32:15] ------------------------------------------------------------
FATAL [2018-10-17 04:32:15] prior.K not within expected range or format.
FATAL [2018-10-17 04:32:15]
FATAL [2018-10-17 04:32:15] This is most likely a user error due to invalid input data or
FATAL [2018-10-17 04:32:15] parameters (PureCN 1.10.0).
INFO [2018-10-17 04:32:15] ------------------------------------------------------------
INFO [2018-10-17 04:32:15] PureCN 1.10.0
INFO [2018-10-17 04:32:15] ------------------------------------------------------------
FATAL [2018-10-17 04:32:15] prior.contamination not within expected range or format.
FATAL [2018-10-17 04:32:15]
FATAL [2018-10-17 04:32:15] This is most likely a user error due to invalid input data or
FATAL [2018-10-17 04:32:15] parameters (PureCN 1.10.0).
INFO [2018-10-17 04:32:15] ------------------------------------------------------------
INFO [2018-10-17 04:32:15] PureCN 1.10.0
INFO [2018-10-17 04:32:15] ------------------------------------------------------------
FATAL [2018-10-17 04:32:15] Iterations not in the expected range from 10 to 250.
FATAL [2018-10-17 04:32:15]
FATAL [2018-10-17 04:32:15] This is most likely a user error due to invalid input data or
FATAL [2018-10-17 04:32:15] parameters (PureCN 1.10.0).
INFO [2018-10-17 04:32:15] ------------------------------------------------------------
INFO [2018-10-17 04:32:15] PureCN 1.10.0
INFO [2018-10-17 04:32:15] ------------------------------------------------------------
FATAL [2018-10-17 04:32:15] Iterations not in the expected range from 10 to 250.
FATAL [2018-10-17 04:32:15]
FATAL [2018-10-17 04:32:15] This is most likely a user error due to invalid input data or
FATAL [2018-10-17 04:32:15] parameters (PureCN 1.10.0).
INFO [2018-10-17 04:32:15] ------------------------------------------------------------
INFO [2018-10-17 04:32:15] PureCN 1.10.0
INFO [2018-10-17 04:32:15] ------------------------------------------------------------
INFO [2018-10-17 04:32:16] ------------------------------------------------------------
INFO [2018-10-17 04:32:16] PureCN 1.10.0
INFO [2018-10-17 04:32:16] ------------------------------------------------------------
INFO [2018-10-17 04:32:16] Loading coverage files...
FATAL [2018-10-17 04:32:16] Interval files in normal and tumor different.
FATAL [2018-10-17 04:32:16]
FATAL [2018-10-17 04:32:16] This is most likely a user error due to invalid input data or
FATAL [2018-10-17 04:32:16] parameters (PureCN 1.10.0).
INFO [2018-10-17 04:32:16] ------------------------------------------------------------
INFO [2018-10-17 04:32:16] PureCN 1.10.0
INFO [2018-10-17 04:32:16] ------------------------------------------------------------
INFO [2018-10-17 04:32:16] Loading coverage files...
INFO [2018-10-17 04:32:17] Mean target coverages: 112X (tumor) 99X (normal).
WARN [2018-10-17 04:32:17] Allosome coverage missing, cannot determine sex.
WARN [2018-10-17 04:32:17] Allosome coverage missing, cannot determine sex.
INFO [2018-10-17 04:32:17] Removing 228 intervals with missing log.ratio.
INFO [2018-10-17 04:32:17] Removing 22 small (< 5bp) targets.
INFO [2018-10-17 04:32:17] Removing 19 targets with low total coverage in normal (< 150.00 reads).
WARN [2018-10-17 04:32:17] No normalDB provided. Provide one for better results.
INFO [2018-10-17 04:32:17] Removing 233 low coverage (< 15.0000X) targets.
INFO [2018-10-17 04:32:17] Using 9547 intervals (9547 on-target, 0 off-target).
INFO [2018-10-17 04:32:17] No off-target intervals. If this is hybrid-capture data, consider adding them.
INFO [2018-10-17 04:32:17] No interval.file provided. Cannot check if data was GC-normalized. Was it?
INFO [2018-10-17 04:32:17] Loading VCF...
INFO [2018-10-17 04:32:17] Found 2331 variants in VCF file.
INFO [2018-10-17 04:32:17] 2282 (97.9%) variants annotated as likely germline (DB INFO flag).
INFO [2018-10-17 04:32:17] LIB-02240e4 is tumor in VCF file.
INFO [2018-10-17 04:32:18] No homozygous variants in VCF, provide unfiltered VCF.
INFO [2018-10-17 04:32:18] Removing 22 non heterozygous (in matched normal) germline SNPs.
INFO [2018-10-17 04:32:19] Initial testing for significant sample cross-contamination: unlikely
INFO [2018-10-17 04:32:19] Removing 48 variants with AF < 0.030 or AF >= 1.000 or less than 4 supporting reads or depth < 15.
INFO [2018-10-17 04:32:19] Removing 16 low quality variants with BQ < 25.
INFO [2018-10-17 04:32:19] Total size of targeted genomic region: 1.63Mb (2.58Mb with 50bp padding).
INFO [2018-10-17 04:32:19] 1.0% of targets contain variants.
INFO [2018-10-17 04:32:19] Removing 2122 variants outside intervals.
INFO [2018-10-17 04:32:19] Found SOMATIC annotation in VCF.
INFO [2018-10-17 04:32:19] Setting somatic prior probabilities for somatic variants to 0.999000 or to 0.000100 otherwise.
INFO [2018-10-17 04:32:20] Found SOMATIC annotation in VCF. Setting mapping bias to 0.973.
INFO [2018-10-17 04:32:20] Sample sex: ?
INFO [2018-10-17 04:32:20] Segmenting data...
INFO [2018-10-17 04:32:20] Loading pre-computed boundaries for DNAcopy...
INFO [2018-10-17 04:32:20] Setting undo.SD parameter to 1.000000.
INFO [2018-10-17 04:32:21] Setting prune.hclust.h parameter to 0.200000.
INFO [2018-10-17 04:32:21] Found 54 segments with median size of 24.88Mb.
INFO [2018-10-17 04:32:21] Using 123 variants.
INFO [2018-10-17 04:32:21] Mean standard deviation of log-ratios: 0.37
INFO [2018-10-17 04:32:21] 2D-grid search of purity and ploidy...
INFO [2018-10-17 04:32:23] Local optima: 0.62/1.9, 0.3/1.9, 0.32/2
INFO [2018-10-17 04:32:23] Testing local optimum 1/3 at purity 0.62 and total ploidy 1.90...
INFO [2018-10-17 04:32:24] Fitting variants for purity 0.65, tumor ploidy 1.74 and contamination 0.01.
INFO [2018-10-17 04:32:25] Optimized purity: 0.65
INFO [2018-10-17 04:32:25] Testing local optimum 2/3 at purity 0.30 and total ploidy 1.90...
INFO [2018-10-17 04:32:25] Recalibrating log-ratios...
INFO [2018-10-17 04:32:25] Testing local optimum 2/3 at purity 0.30 and total ploidy 1.90...
INFO [2018-10-17 04:32:26] Recalibrating log-ratios...
INFO [2018-10-17 04:32:26] Testing local optimum 2/3 at purity 0.30 and total ploidy 1.90...
INFO [2018-10-17 04:32:26] Recalibrating log-ratios...
INFO [2018-10-17 04:32:26] Testing local optimum 2/3 at purity 0.30 and total ploidy 1.90...
INFO [2018-10-17 04:32:27] Testing local optimum 3/3 at purity 0.32 and total ploidy 2.00...
INFO [2018-10-17 04:32:28] Recalibrating log-ratios...
INFO [2018-10-17 04:32:28] Testing local optimum 3/3 at purity 0.32 and total ploidy 2.00...
INFO [2018-10-17 04:32:28] Recalibrating log-ratios...
INFO [2018-10-17 04:32:28] Testing local optimum 3/3 at purity 0.32 and total ploidy 2.00...
INFO [2018-10-17 04:32:29] Recalibrating log-ratios...
INFO [2018-10-17 04:32:29] Testing local optimum 3/3 at purity 0.32 and total ploidy 2.00...
INFO [2018-10-17 04:32:29] Done.
INFO [2018-10-17 04:32:29] ------------------------------------------------------------
INFO [2018-10-17 04:32:30] ------------------------------------------------------------
INFO [2018-10-17 04:32:30] PureCN 1.10.0
INFO [2018-10-17 04:32:30] ------------------------------------------------------------
INFO [2018-10-17 04:32:30] Loading coverage files...
INFO [2018-10-17 04:32:30] Mean target coverages: 112X (tumor) 99X (normal).
WARN [2018-10-17 04:32:30] Allosome coverage missing, cannot determine sex.
WARN [2018-10-17 04:32:30] Allosome coverage missing, cannot determine sex.
INFO [2018-10-17 04:32:30] Removing 228 intervals with missing log.ratio.
INFO [2018-10-17 04:32:30] Removing 22 small (< 5bp) targets.
INFO [2018-10-17 04:32:30] Removing 19 targets with low total coverage in normal (< 150.00 reads).
WARN [2018-10-17 04:32:30] No normalDB provided. Provide one for better results.
INFO [2018-10-17 04:32:30] Removing 233 low coverage (< 15.0000X) targets.
INFO [2018-10-17 04:32:30] Using 9547 intervals (9547 on-target, 0 off-target).
INFO [2018-10-17 04:32:30] No off-target intervals. If this is hybrid-capture data, consider adding them.
INFO [2018-10-17 04:32:30] No interval.file provided. Cannot check if data was GC-normalized. Was it?
INFO [2018-10-17 04:32:30] Loading VCF...
INFO [2018-10-17 04:32:30] Found 2331 variants in VCF file.
WARN [2018-10-17 04:32:30] DB INFO flag contains NAs
INFO [2018-10-17 04:32:31] 2277 (97.7%) variants annotated as likely germline (DB INFO flag).
WARN [2018-10-17 04:32:31] DP FORMAT field contains NAs. Removing 3 variants.
INFO [2018-10-17 04:32:31] LIB-02240e4 is tumor in VCF file.
INFO [2018-10-17 04:32:31] No homozygous variants in VCF, provide unfiltered VCF.
INFO [2018-10-17 04:32:33] Initial testing for significant sample cross-contamination: unlikely
INFO [2018-10-17 04:32:33] Removing 62 variants with AF < 0.030 or AF >= 0.970 or less than 4 supporting reads or depth < 15.
INFO [2018-10-17 04:32:33] Removing 16 low quality variants with BQ < 25.
INFO [2018-10-17 04:32:33] Total size of targeted genomic region: 1.63Mb (2.58Mb with 50bp padding).
INFO [2018-10-17 04:32:33] 1.0% of targets contain variants.
INFO [2018-10-17 04:32:33] Setting somatic prior probabilities for dbSNP hits to 0.000500 or to 0.500000 otherwise.
INFO [2018-10-17 04:32:33] Sample sex: ?
INFO [2018-10-17 04:32:33] Segmenting data...
INFO [2018-10-17 04:32:33] Loading pre-computed boundaries for DNAcopy...
INFO [2018-10-17 04:32:33] Setting undo.SD parameter to 1.000000.
Setting multi-figure configuration
INFO [2018-10-17 04:32:36] Setting prune.hclust.h parameter to 0.200000.
INFO [2018-10-17 04:32:36] Found 59 segments with median size of 17.67Mb.
INFO [2018-10-17 04:32:36] Removing 280 variants outside segments.
INFO [2018-10-17 04:32:36] Using 1970 variants.
INFO [2018-10-17 04:32:36] Mean standard deviation of log-ratios: 0.37
INFO [2018-10-17 04:32:37] 2D-grid search of purity and ploidy...
INFO [2018-10-17 04:32:39] Local optima: 0.62/1.9, 0.3/1.9, 0.32/2
INFO [2018-10-17 04:32:39] Testing local optimum 1/3 at purity 0.62 and total ploidy 1.90...
INFO [2018-10-17 04:32:40] Fitting variants for purity 0.65, tumor ploidy 1.74 and contamination 0.01.
INFO [2018-10-17 04:32:50] Optimized purity: 0.65
INFO [2018-10-17 04:32:50] Testing local optimum 2/3 at purity 0.30 and total ploidy 1.90...
INFO [2018-10-17 04:32:51] Recalibrating log-ratios...
INFO [2018-10-17 04:32:51] Testing local optimum 2/3 at purity 0.30 and total ploidy 1.90...
INFO [2018-10-17 04:32:51] Recalibrating log-ratios...
INFO [2018-10-17 04:32:51] Testing local optimum 2/3 at purity 0.30 and total ploidy 1.90...
INFO [2018-10-17 04:32:51] Recalibrating log-ratios...
INFO [2018-10-17 04:32:51] Testing local optimum 2/3 at purity 0.30 and total ploidy 1.90...
INFO [2018-10-17 04:32:52] Testing local optimum 3/3 at purity 0.32 and total ploidy 2.00...
INFO [2018-10-17 04:32:52] Recalibrating log-ratios...
INFO [2018-10-17 04:32:52] Testing local optimum 3/3 at purity 0.32 and total ploidy 2.00...
INFO [2018-10-17 04:32:53] Recalibrating log-ratios...
INFO [2018-10-17 04:32:53] Testing local optimum 3/3 at purity 0.32 and total ploidy 2.00...
INFO [2018-10-17 04:32:53] Recalibrating log-ratios...
INFO [2018-10-17 04:32:53] Testing local optimum 3/3 at purity 0.32 and total ploidy 2.00...
INFO [2018-10-17 04:32:53] Done.
INFO [2018-10-17 04:32:53] ------------------------------------------------------------
INFO [2018-10-17 04:32:54] ------------------------------------------------------------
INFO [2018-10-17 04:32:54] PureCN 1.10.0
INFO [2018-10-17 04:32:54] ------------------------------------------------------------
INFO [2018-10-17 04:32:54] Loading coverage files...
INFO [2018-10-17 04:32:54] Mean target coverages: 112X (tumor) 99X (normal).
WARN [2018-10-17 04:32:54] Allosome coverage missing, cannot determine sex.
WARN [2018-10-17 04:32:54] Allosome coverage missing, cannot determine sex.
INFO [2018-10-17 04:32:54] No Gene column in interval.file. You won't get gene-level calls.
INFO [2018-10-17 04:32:56] Removing 228 intervals with missing log.ratio.
INFO [2018-10-17 04:32:56] Removing 15 low/high GC targets.
INFO [2018-10-17 04:32:56] Removing 21 small (< 5bp) targets.
INFO [2018-10-17 04:32:56] Removing 19 targets with low total coverage in normal (< 150.00 reads).
WARN [2018-10-17 04:32:56] No normalDB provided. Provide one for better results.
INFO [2018-10-17 04:32:56] Removing 232 low coverage (< 15.0000X) targets.
INFO [2018-10-17 04:32:56] Using 9534 intervals (9534 on-target, 0 off-target).
INFO [2018-10-17 04:32:56] No off-target intervals. If this is hybrid-capture data, consider adding them.
INFO [2018-10-17 04:32:56] AT/GC dropout: 1.03 (tumor), 1.03 (normal).
INFO [2018-10-17 04:32:56] Loading VCF...
INFO [2018-10-17 04:32:56] Found 2331 variants in VCF file.
INFO [2018-10-17 04:32:56] 2282 (97.9%) variants annotated as likely germline (DB INFO flag).
INFO [2018-10-17 04:32:56] LIB-02240e4 is tumor in VCF file.
INFO [2018-10-17 04:32:56] No homozygous variants in VCF, provide unfiltered VCF.
INFO [2018-10-17 04:32:57] Removing 22 non heterozygous (in matched normal) germline SNPs.
INFO [2018-10-17 04:32:57] Initial testing for significant sample cross-contamination: unlikely
INFO [2018-10-17 04:32:57] Removing 48 variants with AF < 0.030 or AF >= Inf or less than 4 supporting reads or depth < 15.
INFO [2018-10-17 04:32:58] Removing 16 low quality variants with BQ < 25.
INFO [2018-10-17 04:32:58] Total size of targeted genomic region: 1.63Mb (2.58Mb with 50bp padding).
INFO [2018-10-17 04:32:58] 1.0% of targets contain variants.
INFO [2018-10-17 04:32:58] Removing 2122 variants outside intervals.
INFO [2018-10-17 04:32:58] Found SOMATIC annotation in VCF.
INFO [2018-10-17 04:32:58] Setting somatic prior probabilities for somatic variants to 0.999000 or to 0.000100 otherwise.
INFO [2018-10-17 04:32:58] Found SOMATIC annotation in VCF. Setting mapping bias to 0.973.
INFO [2018-10-17 04:32:58] Sample sex: ?
INFO [2018-10-17 04:32:58] Segmenting data...
INFO [2018-10-17 04:32:58] Loading pre-computed boundaries for DNAcopy...
INFO [2018-10-17 04:32:58] Setting undo.SD parameter to 1.000000.
Setting multi-figure configuration
INFO [2018-10-17 04:32:59] Setting prune.hclust.h parameter to 0.200000.
INFO [2018-10-17 04:33:00] Found 54 segments with median size of 24.88Mb.
INFO [2018-10-17 04:33:00] Using 123 variants.
INFO [2018-10-17 04:33:00] Mean standard deviation of log-ratios: 0.37
INFO [2018-10-17 04:33:00] 2D-grid search of purity and ploidy...
INFO [2018-10-17 04:33:00] 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 [2018-10-17 04:33:00] Testing local optimum 1/6 at purity 0.65 and total ploidy 1.90...
INFO [2018-10-17 04:33:01] Fitting variants for purity 0.65, tumor ploidy 1.74 and contamination 0.01.
INFO [2018-10-17 04:33:02] Optimized purity: 0.65
INFO [2018-10-17 04:33:02] Testing local optimum 2/6 at purity 0.33 and total ploidy 1.90...
INFO [2018-10-17 04:33:02] Recalibrating log-ratios...
INFO [2018-10-17 04:33:02] Testing local optimum 2/6 at purity 0.33 and total ploidy 1.90...
INFO [2018-10-17 04:33:03] Recalibrating log-ratios...
INFO [2018-10-17 04:33:03] Testing local optimum 2/6 at purity 0.33 and total ploidy 1.90...
INFO [2018-10-17 04:33:03] Recalibrating log-ratios...
INFO [2018-10-17 04:33:03] Testing local optimum 2/6 at purity 0.33 and total ploidy 1.90...
INFO [2018-10-17 04:33:04] Testing local optimum 3/6 at purity 0.38 and total ploidy 2.20...
INFO [2018-10-17 04:33:05] Fitting variants for purity 0.40, tumor ploidy 2.34 and contamination 0.01.
INFO [2018-10-17 04:33:06] Optimized purity: 0.40
INFO [2018-10-17 04:33:06] Testing local optimum 4/6 at purity 0.27 and total ploidy 2.10...
INFO [2018-10-17 04:33:06] Recalibrating log-ratios...
INFO [2018-10-17 04:33:06] Testing local optimum 4/6 at purity 0.27 and total ploidy 2.10...
INFO [2018-10-17 04:33:07] Recalibrating log-ratios...
INFO [2018-10-17 04:33:07] Testing local optimum 4/6 at purity 0.27 and total ploidy 2.10...
INFO [2018-10-17 04:33:07] Recalibrating log-ratios...
INFO [2018-10-17 04:33:07] Testing local optimum 4/6 at purity 0.27 and total ploidy 2.10...
INFO [2018-10-17 04:33:08] Testing local optimum 5/6 at purity 0.29 and total ploidy 2.00...
INFO [2018-10-17 04:33:08] Recalibrating log-ratios...
INFO [2018-10-17 04:33:08] Testing local optimum 5/6 at purity 0.29 and total ploidy 2.00...
INFO [2018-10-17 04:33:08] Recalibrating log-ratios...
INFO [2018-10-17 04:33:08] Testing local optimum 5/6 at purity 0.29 and total ploidy 2.00...
INFO [2018-10-17 04:33:09] Recalibrating log-ratios...
INFO [2018-10-17 04:33:09] Testing local optimum 5/6 at purity 0.29 and total ploidy 2.00...
INFO [2018-10-17 04:33:09] Testing local optimum 6/6 at purity 0.78 and total ploidy 2.00...
INFO [2018-10-17 04:33:10] Fitting variants for purity 0.65, tumor ploidy 1.74 and contamination 0.01.
INFO [2018-10-17 04:33:11] Optimized purity: 0.65
INFO [2018-10-17 04:33:11] Done.
INFO [2018-10-17 04:33:11] ------------------------------------------------------------
FATAL [2018-10-17 04:33:11] This function requires gene-level calls. Please add a column 'Gene'
FATAL [2018-10-17 04:33:11] containing gene symbols to the interval.file.
FATAL [2018-10-17 04:33:11]
FATAL [2018-10-17 04:33:11] This is most likely a user error due to invalid input data or
FATAL [2018-10-17 04:33:11] parameters (PureCN 1.10.0).
INFO [2018-10-17 04:33:13] ------------------------------------------------------------
INFO [2018-10-17 04:33:13] PureCN 1.10.0
INFO [2018-10-17 04:33:13] ------------------------------------------------------------
INFO [2018-10-17 04:33:13] Loading coverage files...
INFO [2018-10-17 04:33:13] Mean target coverages: 112X (tumor) 99X (normal).
WARN [2018-10-17 04:33:14] Allosome coverage missing, cannot determine sex.
WARN [2018-10-17 04:33:14] Allosome coverage missing, cannot determine sex.
INFO [2018-10-17 04:33:14] Removing 228 intervals with missing log.ratio.
INFO [2018-10-17 04:33:14] Removing 22 small (< 5bp) targets.
INFO [2018-10-17 04:33:14] Removing 19 targets with low total coverage in normal (< 150.00 reads).
WARN [2018-10-17 04:33:14] No normalDB provided. Provide one for better results.
INFO [2018-10-17 04:33:14] Removing 233 low coverage (< 15.0000X) targets.
INFO [2018-10-17 04:33:14] Using 9547 intervals (9547 on-target, 0 off-target).
INFO [2018-10-17 04:33:14] No off-target intervals. If this is hybrid-capture data, consider adding them.
INFO [2018-10-17 04:33:14] No interval.file provided. Cannot check if data was GC-normalized. Was it?
INFO [2018-10-17 04:33:14] Loading VCF...
INFO [2018-10-17 04:33:14] Found 2331 variants in VCF file.
INFO [2018-10-17 04:33:14] 2282 (97.9%) variants annotated as likely germline (DB INFO flag).
FATAL [2018-10-17 04:33:14] Different chromosome names in coverage and VCF.
FATAL [2018-10-17 04:33:14]
FATAL [2018-10-17 04:33:14] This is most likely a user error due to invalid input data or
FATAL [2018-10-17 04:33:14] parameters (PureCN 1.10.0).
INFO [2018-10-17 04:33:14] ------------------------------------------------------------
INFO [2018-10-17 04:33:14] PureCN 1.10.0
INFO [2018-10-17 04:33:14] ------------------------------------------------------------
INFO [2018-10-17 04:33:14] Loading coverage files...
INFO [2018-10-17 04:33:14] Mean target coverages: 112X (tumor) 99X (normal).
WARN [2018-10-17 04:33:14] Allosome coverage missing, cannot determine sex.
WARN [2018-10-17 04:33:14] Allosome coverage missing, cannot determine sex.
FATAL [2018-10-17 04:33:14] tumor.coverage.file and interval.file do not align.
FATAL [2018-10-17 04:33:14]
FATAL [2018-10-17 04:33:14] This is most likely a user error due to invalid input data or
FATAL [2018-10-17 04:33:14] parameters (PureCN 1.10.0).
INFO [2018-10-17 04:33:15] ------------------------------------------------------------
INFO [2018-10-17 04:33:15] PureCN 1.10.0
INFO [2018-10-17 04:33:15] ------------------------------------------------------------
INFO [2018-10-17 04:33:15] Loading coverage files...
INFO [2018-10-17 04:33:15] Mean target coverages: 112X (tumor) 99X (normal).
WARN [2018-10-17 04:33:15] Allosome coverage missing, cannot determine sex.
WARN [2018-10-17 04:33:15] Allosome coverage missing, cannot determine sex.
INFO [2018-10-17 04:33:16] Removing 228 intervals with missing log.ratio.
INFO [2018-10-17 04:33:16] Removing 15 low/high GC targets.
INFO [2018-10-17 04:33:16] Removing 21 small (< 5bp) targets.
INFO [2018-10-17 04:33:16] Removing 19 targets with low total coverage in normal (< 150.00 reads).
WARN [2018-10-17 04:33:16] No normalDB provided. Provide one for better results.
INFO [2018-10-17 04:33:16] Removing 232 low coverage (< 15.0000X) targets.
INFO [2018-10-17 04:33:16] Using 9534 intervals (9534 on-target, 0 off-target).
INFO [2018-10-17 04:33:16] No off-target intervals. If this is hybrid-capture data, consider adding them.
INFO [2018-10-17 04:33:17] AT/GC dropout: 1.03 (tumor), 1.03 (normal).
INFO [2018-10-17 04:33:17] Loading VCF...
INFO [2018-10-17 04:33:17] Found 2331 variants in VCF file.
INFO [2018-10-17 04:33:17] 2282 (97.9%) variants annotated as likely germline (DB INFO flag).
INFO [2018-10-17 04:33:17] LIB-02240e4 is tumor in VCF file.
INFO [2018-10-17 04:33:17] No homozygous variants in VCF, provide unfiltered VCF.
INFO [2018-10-17 04:33:17] Removing 22 non heterozygous (in matched normal) germline SNPs.
INFO [2018-10-17 04:33:18] Initial testing for significant sample cross-contamination: unlikely
INFO [2018-10-17 04:33:19] Removing 48 variants with AF < 0.030 or AF >= 1.000 or less than 4 supporting reads or depth < 15.
INFO [2018-10-17 04:33:19] Removing 16 low quality variants with BQ < 25.
INFO [2018-10-17 04:33:19] Total size of targeted genomic region: 1.63Mb (2.58Mb with 50bp padding).
INFO [2018-10-17 04:33:19] 1.0% of targets contain variants.
INFO [2018-10-17 04:33:19] Removing 2122 variants outside intervals.
INFO [2018-10-17 04:33:19] Found SOMATIC annotation in VCF.
INFO [2018-10-17 04:33:19] Setting somatic prior probabilities for somatic variants to 0.999000 or to 0.000100 otherwise.
INFO [2018-10-17 04:33:19] Found SOMATIC annotation in VCF. Setting mapping bias to 0.973.
INFO [2018-10-17 04:33:19] Sample sex: ?
INFO [2018-10-17 04:33:19] Segmenting data...
INFO [2018-10-17 04:33:20] Loading pre-computed boundaries for DNAcopy...
INFO [2018-10-17 04:33:20] Setting undo.SD parameter to 1.000000.
INFO [2018-10-17 04:33:21] Setting prune.hclust.h parameter to 0.200000.
INFO [2018-10-17 04:33:22] Found 54 segments with median size of 24.88Mb.
INFO [2018-10-17 04:33:22] Using 123 variants.
INFO [2018-10-17 04:33:22] Mean standard deviation of log-ratios: 0.37
INFO [2018-10-17 04:33:22] 2D-grid search of purity and ploidy...
INFO [2018-10-17 04:33:23] Local optima: 0.63/1.9, 0.4/2.1
INFO [2018-10-17 04:33:23] Testing local optimum 1/2 at purity 0.63 and total ploidy 1.90...
INFO [2018-10-17 04:33:24] Fitting variants for purity 0.65, tumor ploidy 1.74 and contamination 0.01.
INFO [2018-10-17 04:33:26] Optimized purity: 0.65
INFO [2018-10-17 04:33:26] Testing local optimum 2/2 at purity 0.40 and total ploidy 2.10...
INFO [2018-10-17 04:33:26] Recalibrating log-ratios...
INFO [2018-10-17 04:33:26] Testing local optimum 2/2 at purity 0.40 and total ploidy 2.10...
INFO [2018-10-17 04:33:27] Recalibrating log-ratios...
INFO [2018-10-17 04:33:27] Testing local optimum 2/2 at purity 0.40 and total ploidy 2.10...
INFO [2018-10-17 04:33:27] Recalibrating log-ratios...
INFO [2018-10-17 04:33:27] Testing local optimum 2/2 at purity 0.40 and total ploidy 2.10...
INFO [2018-10-17 04:33:28] Done.
INFO [2018-10-17 04:33:28] ------------------------------------------------------------
INFO [2018-10-17 04:33:28] ------------------------------------------------------------
INFO [2018-10-17 04:33:28] PureCN 1.10.0
INFO [2018-10-17 04:33:28] ------------------------------------------------------------
INFO [2018-10-17 04:33:28] Loading coverage files...
INFO [2018-10-17 04:33:28] Mean target coverages: 112X (tumor) 99X (normal).
WARN [2018-10-17 04:33:28] Allosome coverage missing, cannot determine sex.
WARN [2018-10-17 04:33:28] Allosome coverage missing, cannot determine sex.
INFO [2018-10-17 04:33:29] Removing 228 intervals with missing log.ratio.
INFO [2018-10-17 04:33:29] Removing 15 low/high GC targets.
INFO [2018-10-17 04:33:29] Removing 21 small (< 5bp) targets.
INFO [2018-10-17 04:33:29] Removing 19 targets with low total coverage in normal (< 150.00 reads).
WARN [2018-10-17 04:33:29] No normalDB provided. Provide one for better results.
INFO [2018-10-17 04:33:29] Removing 232 low coverage (< 15.0000X) targets.
INFO [2018-10-17 04:33:29] Using 9534 intervals (9534 on-target, 0 off-target).
INFO [2018-10-17 04:33:29] No off-target intervals. If this is hybrid-capture data, consider adding them.
INFO [2018-10-17 04:33:29] AT/GC dropout: 1.03 (tumor), 1.03 (normal).
INFO [2018-10-17 04:33:29] Loading VCF...
INFO [2018-10-17 04:33:29] Found 2331 variants in VCF file.
INFO [2018-10-17 04:33:29] 2282 (97.9%) variants annotated as likely germline (DB INFO flag).
INFO [2018-10-17 04:33:30] LIB-02240e4 is tumor in VCF file.
INFO [2018-10-17 04:33:30] No homozygous variants in VCF, provide unfiltered VCF.
INFO [2018-10-17 04:33:30] Removing 22 non heterozygous (in matched normal) germline SNPs.
INFO [2018-10-17 04:33:31] Initial testing for significant sample cross-contamination: unlikely
INFO [2018-10-17 04:33:31] Removing 48 variants with AF < 0.030 or AF >= 1.000 or less than 4 supporting reads or depth < 15.
INFO [2018-10-17 04:33:31] Removing 16 low quality variants with BQ < 25.
INFO [2018-10-17 04:33:31] Total size of targeted genomic region: 1.63Mb (2.58Mb with 50bp padding).
INFO [2018-10-17 04:33:32] 1.0% of targets contain variants.
INFO [2018-10-17 04:33:32] Removing 2122 variants outside intervals.
INFO [2018-10-17 04:33:32] Found SOMATIC annotation in VCF.
INFO [2018-10-17 04:33:32] Setting somatic prior probabilities for somatic variants to 0.999000 or to 0.000100 otherwise.
INFO [2018-10-17 04:33:32] Found SOMATIC annotation in VCF. Setting mapping bias to 0.973.
INFO [2018-10-17 04:33:32] Sample sex: ?
INFO [2018-10-17 04:33:32] Segmenting data...
INFO [2018-10-17 04:33:32] Loading pre-computed boundaries for DNAcopy...
INFO [2018-10-17 04:33:32] Setting undo.SD parameter to 1.000000.
INFO [2018-10-17 04:33:34] Setting prune.hclust.h parameter to 0.200000.
INFO [2018-10-17 04:33:34] Found 54 segments with median size of 24.88Mb.
INFO [2018-10-17 04:33:34] Using 123 variants.
INFO [2018-10-17 04:33:34] Mean standard deviation of log-ratios: 0.37
INFO [2018-10-17 04:33:34] 2D-grid search of purity and ploidy...
INFO [2018-10-17 04:33:35] Local optima: 0.63/1.9, 0.4/2.1
INFO [2018-10-17 04:33:35] Testing local optimum 1/2 at purity 0.63 and total ploidy 1.90...
INFO [2018-10-17 04:33:36] Fitting variants for purity 0.65, tumor ploidy 1.74 and contamination 0.01.
INFO [2018-10-17 04:33:38] Optimized purity: 0.65
INFO [2018-10-17 04:33:38] Testing local optimum 2/2 at purity 0.40 and total ploidy 2.10...
INFO [2018-10-17 04:33:39] Recalibrating log-ratios...
INFO [2018-10-17 04:33:39] Testing local optimum 2/2 at purity 0.40 and total ploidy 2.10...
INFO [2018-10-17 04:33:39] Recalibrating log-ratios...
INFO [2018-10-17 04:33:39] Testing local optimum 2/2 at purity 0.40 and total ploidy 2.10...
INFO [2018-10-17 04:33:39] Recalibrating log-ratios...
INFO [2018-10-17 04:33:39] Testing local optimum 2/2 at purity 0.40 and total ploidy 2.10...
INFO [2018-10-17 04:33:40] Done.
INFO [2018-10-17 04:33:40] ------------------------------------------------------------
FATAL [2018-10-17 04:33:40] chr1 not valid chromosome name(s). Valid names are:
FATAL [2018-10-17 04:33:40] 1,2,3,4,5,6,7,8,9,10,11,12,13,14,15,16,17,18,19,20,21,22
FATAL [2018-10-17 04:33:40]
FATAL [2018-10-17 04:33:40] This is most likely a user error due to invalid input data or
FATAL [2018-10-17 04:33:40] parameters (PureCN 1.10.0).
INFO [2018-10-17 04:33:40] ------------------------------------------------------------
INFO [2018-10-17 04:33:40] PureCN 1.10.0
INFO [2018-10-17 04:33:40] ------------------------------------------------------------
INFO [2018-10-17 04:33:40] Loading coverage files...
WARN [2018-10-17 04:33:41] Provided sampleid (Sample2) does not match Sample1 found in segmentation.
WARN [2018-10-17 04:33:41] Allosome coverage missing, cannot determine sex.
WARN [2018-10-17 04:33:41] Allosome coverage missing, cannot determine sex.
INFO [2018-10-17 04:33:41] Removing 10 intervals with missing log.ratio.
INFO [2018-10-17 04:33:41] Using 10039 intervals (10039 on-target, 0 off-target).
INFO [2018-10-17 04:33:41] No off-target intervals. If this is hybrid-capture data, consider adding them.
INFO [2018-10-17 04:33:41] Loading VCF...
INFO [2018-10-17 04:33:41] Found 2331 variants in VCF file.
INFO [2018-10-17 04:33:41] 2282 (97.9%) variants annotated as likely germline (DB INFO flag).
INFO [2018-10-17 04:33:41] LIB-02240e4 is tumor in VCF file.
INFO [2018-10-17 04:33:41] No homozygous variants in VCF, provide unfiltered VCF.
INFO [2018-10-17 04:33:42] Removing 22 non heterozygous (in matched normal) germline SNPs.
INFO [2018-10-17 04:33:43] Initial testing for significant sample cross-contamination: unlikely
INFO [2018-10-17 04:33:43] Removing 48 variants with AF < 0.030 or AF >= 1.000 or less than 4 supporting reads or depth < 15.
INFO [2018-10-17 04:33:43] Removing 16 low quality variants with BQ < 25.
INFO [2018-10-17 04:33:43] Total size of targeted genomic region: 1.71Mb (2.70Mb with 50bp padding).
INFO [2018-10-17 04:33:43] 1.0% of targets contain variants.
INFO [2018-10-17 04:33:43] Removing 2120 variants outside intervals.
INFO [2018-10-17 04:33:43] Found SOMATIC annotation in VCF.
INFO [2018-10-17 04:33:44] Setting somatic prior probabilities for somatic variants to 0.999000 or to 0.000100 otherwise.
INFO [2018-10-17 04:33:44] Found SOMATIC annotation in VCF. Setting mapping bias to 0.976.
INFO [2018-10-17 04:33:44] Sample sex: ?
INFO [2018-10-17 04:33:44] Segmenting data...
WARN [2018-10-17 04:33:44] Provided sampleid (Sample2) does not match Sample1 found in segmentation.
INFO [2018-10-17 04:33:44] Loading pre-computed boundaries for DNAcopy...
INFO [2018-10-17 04:33:44] Setting undo.SD parameter to 0.000000.
Setting multi-figure configuration
INFO [2018-10-17 04:33:44] Setting prune.hclust.h parameter to 0.200000.
INFO [2018-10-17 04:33:44] Found 54 segments with median size of 24.88Mb.
INFO [2018-10-17 04:33:45] Using 125 variants.
INFO [2018-10-17 04:33:45] Mean standard deviation of log-ratios: 0.40
INFO [2018-10-17 04:33:45] 2D-grid search of purity and ploidy...
INFO [2018-10-17 04:33:46] Local optima: 0.6/1.9, 0.4/2.1
INFO [2018-10-17 04:33:46] Testing local optimum 1/2 at purity 0.60 and total ploidy 1.90...
INFO [2018-10-17 04:33:47] Fitting variants for purity 0.65, tumor ploidy 1.74 and contamination 0.01.
INFO [2018-10-17 04:33:48] Optimized purity: 0.65
INFO [2018-10-17 04:33:48] Testing local optimum 2/2 at purity 0.40 and total ploidy 2.10...
INFO [2018-10-17 04:33:48] Recalibrating log-ratios...
INFO [2018-10-17 04:33:48] Testing local optimum 2/2 at purity 0.40 and total ploidy 2.10...
INFO [2018-10-17 04:33:49] Recalibrating log-ratios...
INFO [2018-10-17 04:33:49] Testing local optimum 2/2 at purity 0.40 and total ploidy 2.10...
INFO [2018-10-17 04:33:49] Recalibrating log-ratios...
INFO [2018-10-17 04:33:49] Testing local optimum 2/2 at purity 0.40 and total ploidy 2.10...
INFO [2018-10-17 04:33:50] Done.
INFO [2018-10-17 04:33:50] ------------------------------------------------------------
WARN [2018-10-17 04:33:50] Allosome coverage missing, cannot determine sex.
WARN [2018-10-17 04:33:50] Allosome coverage missing, cannot determine sex.
Setting multi-figure configuration
FATAL [2018-10-17 04:34:01] Segmentation file expected with colnames ID, chrom, loc.start, loc.end,
FATAL [2018-10-17 04:34:01] num.mark, seg.mean
FATAL [2018-10-17 04:34:01]
FATAL [2018-10-17 04:34:01] This is most likely a user error due to invalid input data or
FATAL [2018-10-17 04:34:01] parameters (PureCN 1.10.0).
FATAL [2018-10-17 04:34:02] seg.file contains multiple samples and sampleid missing.
FATAL [2018-10-17 04:34:02]
FATAL [2018-10-17 04:34:02] This is most likely a user error due to invalid input data or
FATAL [2018-10-17 04:34:02] parameters (PureCN 1.10.0).
FATAL [2018-10-17 04:34:02] seg.file contains multiple samples and sampleid does not match any.
FATAL [2018-10-17 04:34:02]
FATAL [2018-10-17 04:34:02] This is most likely a user error due to invalid input data or
FATAL [2018-10-17 04:34:02] parameters (PureCN 1.10.0).
WARN [2018-10-17 04:34:02] Allosome coverage missing, cannot determine sex.
WARN [2018-10-17 04:34:02] Allosome coverage missing, cannot determine sex.
Setting multi-figure configuration
WARN [2018-10-17 04:34:13] Provided sampleid (example_tumor.txt) does not match Sample1 found in segmentation.
WARN [2018-10-17 04:34:13] Allosome coverage missing, cannot determine sex.
WARN [2018-10-17 04:34:13] Allosome coverage missing, cannot determine sex.
WARN [2018-10-17 04:34:13] Provided sampleid (example_tumor.txt) does not match Sample1 found in segmentation.
WARN [2018-10-17 04:34:18] Allosome coverage missing, cannot determine sex.
WARN [2018-10-17 04:34:18] Allosome coverage missing, cannot determine sex.
WARN [2018-10-17 04:34:19] No normalDB provided. Provide one for better results.
WARN [2018-10-17 04:34:31] Allosome coverage missing, cannot determine sex.
WARN [2018-10-17 04:34:31] Allosome coverage missing, cannot determine sex.
WARN [2018-10-17 04:34:35] Provided sampleid (Sample.1) does not match Sample1 found in segmentation.
Setting multi-figure configuration
WARN [2018-10-17 04:37:10] Allosome coverage missing, cannot determine sex.
WARN [2018-10-17 04:37:10] Allosome coverage missing, cannot determine sex.
WARN [2018-10-17 04:37:10] No normalDB provided. Provide one for better results.
WARN [2018-10-17 04:37:11] Sampleid looks like a normal in VCF, not like a tumor.
WARN [2018-10-17 04:37:24] Allosome coverage missing, cannot determine sex.
WARN [2018-10-17 04:37:24] Allosome coverage missing, cannot determine sex.
WARN [2018-10-17 04:37:24] No normalDB provided. Provide one for better results.
WARN [2018-10-17 04:37:34] Allosome coverage missing, cannot determine sex.
WARN [2018-10-17 04:37:34] Allosome coverage missing, cannot determine sex.
FATAL [2018-10-17 04:37:34] normalDB not a valid normalDB object. Use createNormalDatabase to
FATAL [2018-10-17 04:37:34] create one.
FATAL [2018-10-17 04:37:34]
FATAL [2018-10-17 04:37:34] This is most likely a user error due to invalid input data or
FATAL [2018-10-17 04:37:34] parameters (PureCN 1.10.0).
WARN [2018-10-17 04:37:35] You are likely not using the correct baits file!
WARN [2018-10-17 04:37:35] Allosome coverage missing, cannot determine sex.
WARN [2018-10-17 04:37:35] Allosome coverage missing, cannot determine sex.
WARN [2018-10-17 04:37:38] Allosome coverage missing, cannot determine sex.
WARN [2018-10-17 04:37:38] Allosome coverage missing, cannot determine sex.
FATAL [2018-10-17 04:37:38] normalDB appears to be empty.
FATAL [2018-10-17 04:37:38]
FATAL [2018-10-17 04:37:38] This is most likely a user error due to invalid input data or
FATAL [2018-10-17 04:37:38] parameters (PureCN 1.10.0).
WARN [2018-10-17 04:37:39] Allosome coverage missing, cannot determine sex.
WARN [2018-10-17 04:37:39] Allosome coverage missing, cannot determine sex.
WARN [2018-10-17 04:37:39] Intervals in coverage and interval.file have conflicting on/off-target annotation.
FATAL [2018-10-17 04:38:11] min.normals must be >=2.
FATAL [2018-10-17 04:38:11]
FATAL [2018-10-17 04:38:11] This is most likely a user error due to invalid input data or
FATAL [2018-10-17 04:38:11] parameters (PureCN 1.10.0).
FATAL [2018-10-17 04:38:26] The normal.panel.vcf.file contains only a single sample.
FATAL [2018-10-17 04:38:26]
FATAL [2018-10-17 04:38:26] This is most likely a user error due to invalid input data or
FATAL [2018-10-17 04:38:26] parameters (PureCN 1.10.0).
== testthat results ===========================================================
OK: 293 SKIPPED: 0 FAILED: 0
>
> proc.time()
user system elapsed
519.06 2.31 579.23
|