isovic / racon

Ultrafast consensus module for raw de novo genome assembly of long uncorrected reads. http://genome.cshlp.org/content/early/2017/01/18/gr.214270.116 Note: This was the original repository which will no longer be officially maintained. Please use the new official repository here:
https://github.com/lbcb-sci/racon
MIT License
257 stars 48 forks source link

error: unable to find subsampled sequences! #211

Open ziczhang opened 2 years ago

ziczhang commented 2 years ago

Hi,

I've found a small bug on racon_wrapper.py.

When I used the racon_wrapper as below: racon_wrapper --subsample 3000000000 20 ${NAME}.fq.gz mapped.paf racon_0.fa > racon_1.fa

it caused

[RaconWrapper::run] preparing data with rampler [RaconWrapper::run] error: unable to find subsampled sequences!

This was because of my fastq file's name, which was "fq.gz" but not "fastq.gz" and on my racon's work_directory, "${NAME}_20x.fq" was generated, but not "${NAME}_20x.fastq"

Specifically, this bug was mainly caused by lines 72-76 in racon_wrapper.py.

Best,

Zicong.