lmrodriguezr / nonpareil

Estimate metagenomic coverage and sequence diversity
http://enve-omics.ce.gatech.edu/nonpareil/
Other
42 stars 11 forks source link

nonpareil failing with too few reads #60

Closed fgvieira closed 3 weeks ago

fgvieira commented 1 year ago

I am running nonpareil as:

nonpareil -T kmer -s input.fastq -f fastq -F -o out.npo -a out.npa -C out.npc

through a pipeline and, sometimes, there are files with not that many reads. For those, I am getting the error:

Nonpareil v3.401
Fatal error:
Reduce the number of query reads (-X) to ≤ 10% of total reads
 [      0.0] Fatal error: Reduce the number of query reads (-X) to ≤ 10% of total reads

Looking at the code, it seems that nonpareil requires -X to be less than 10% of the total reads. I am a bit confused with this options since, either: