BBnorm not outputting txt files when trying to determine ploidy on pb fastq.gz reads

Question

[kristina.mahan]

I'm trying to determine ploidy of my genome using BBnorm. I am not getting any error but no txt files are produced.

khist.sh in=reads.fq khist=khist.txt peaks=peaks.txt

or equivalently

bbnorm.sh in=reads.fq khist=khist.txt peaks=peaks.txt passes=1 prefilter minprob=0 minqual=0 mindepth=0

This is my input:

khist.sh in=A.subreads.fastq.gz khist=khist.txt peaks=peaks.txt

Can it use fastq.gz files? Can you use illumina reads?

[GenoMax]

All BBTools can use compressed files. They can use Illumina data for sure. Title of your post seems to suggest that you may have PacBio data? That may not work since with exception of mapPacBio.sh most BBTools programs are for short read data. That said can you try the following?

For khist.sh usage indicates :

 khist.sh in=<input> hist=<histogram output>

For bbnorm.sh:

 bbnorm.sh in=<input> out=<reads to keep> outt=<reads to toss> hist=<histogram output>

I assume you are using commands from BBNorm Guide page.