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