PureCN: Tool for estimation of tumor purity and ploidy
1
0
Entering edit mode
3.2 years ago
sruthi ▴ 30

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.

software error next-gen • 1.1k views
ADD COMMENT
1
Entering edit mode
3.2 years ago

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 COMMENT
0
Entering edit mode

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 REPLY
0
Entering edit mode

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

ADD REPLY
0
Entering edit mode

Thanks Markus. I will open an GitHub issue.

ADD REPLY

Login before adding your answer.

Traffic: 1935 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6