Question: PureCN: Tool for estimation of tumor purity and ploidy
gravatar for sruthi
18 days ago by
sruthi0 wrote:

Hi, I'm trying to calculate tumor purity and ploidy using a tool called PureCN written in Rscript. The initial step is to generate an interval file based on GC content from hg38.bed and hg38.fa files. The command used for the process is given below: Rscript $PURECN/IntervalFile.R --infile /path/to/hg19.bed --fasta /path/to/hg19.fa --outfile /path/to/hg19_intervals.txt

The following error shows up:

Error in value[3L] : record 4423505 (chr1_KN196472v1_fix:111333-111727) failed file: /path/to/hg38.fa Calls: preprocessIntervals ... tryCatch -> tryCatchList -> tryCatchOne -> <anonymous>

It would be of great help if someone could help me out in troubleshooting this error. Thanks in advance.

next-gen software error • 111 views
ADD COMMENTlink modified 18 days ago by markus.riester320 • written 18 days ago by sruthi0
gravatar for markus.riester
18 days ago by
markus.riester320 wrote:


the infile needs to be a BED file with the capture baits. This file should only contain regular chromosomes. If it does, please open a GitHub issue at, ideally including a minimal BED file to reproduce the issue.

Thanks a lot, Markus

ADD COMMENTlink written 18 days ago by markus.riester320

Thanks Markus. I uploaded a BED file with capture baits, but the same error pops up. Is there something that I must look in the FASTA file? Also, which build is supported by the tool? hg19 or hg38.

ADD REPLYlink written 17 days ago by sruthi0

Both builds are supported. Please try to create a reproducible example, like a minimal BED file and open an GitHub issue. Thank you.

ADD REPLYlink written 16 days ago by markus.riester320

Thanks Markus. I will open an GitHub issue.

ADD REPLYlink written 9 days ago by sruthi0
Please log in to add an answer.


Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 1211 users visited in the last hour