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