I have some 300 fasta genomes and 6 nanopore fastq files and I was trying to run mashtree_jackknife
. However I keep getting the following error
mashtree_jackknife.pl --reps 100 --numcpus 4 *.fastq *.fasta -- --min-depth 0 > mashtree.jackknife.dnd
Thread 1 terminated abnormally: mashtree_jackknife.pl: Mashtree::mashDist: Can't take log of 0
Stopped at /home/subudhak/miniconda3/envs/lib/perl5/site_perl/Mashtree.pm line 351.
mashtree_jackknife.pl: main::main: ERROR: one or more threads did not return an array of jack knife tree files as expected.
Stopped at /home/subudhak/miniconda3/envs/bin/mashtree_jackknife.pl line 153.
Perl exited with active threads:
3 running and unjoined
0 finished and unjoined
0 running and detached
However, mashtree runs successfully.
Besides, I wish to understand the use of --genomeSize
option. These genomes are on an average 8.3Kb long so should I set it to 8300 ?