Question: PureCN: Tool for estimation of tumor purity and ploidy
0
gravatar for sruthi
9 months ago by
sruthi20
sruthi20 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 • 400 views
ADD COMMENTlink modified 9 months ago by markus.riester480 • written 9 months ago by sruthi20
1
gravatar for markus.riester
9 months ago by
markus.riester480 wrote:

Hi,

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 https://github.com/lima1/PureCN/issues, ideally including a minimal BED file to reproduce the issue.

Thanks a lot, Markus

ADD COMMENTlink written 9 months ago by markus.riester480

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 9 months ago by sruthi20

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 9 months ago by markus.riester480

Thanks Markus. I will open an GitHub issue.

ADD REPLYlink written 8 months ago by sruthi20
Please log in to add an answer.

Help
Access

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