This is a test version of Biostars. For the public version, visit https://www.biostars.org.
BBnorm not outputting txt files when trying to determine ploidy on pb fastq.gz reads

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?

kmer khist bbnorm bbmap ploidy
• 1,397 views • link
updated by GenoMax GenoMax • written by kristina.mahan kristina.mahan

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.

• 0 views • link
written by GenoMax GenoMax

0 answers

No answers yet.

Log in to answer this question.

More posts like this