**50**wrote:

Hi folks. I am trying to use your software KmerGenie. I have some problems. Because it does not draw the histograms at >15 value. It is hilarious when i put without --diploid command. It could tell me the best k but without doing at >15. I tried from 15 to 55 because at more k (to121K) it is more slower (1h30min) at the result is the same "could not fit".

`sudo ./kmergenie /home/soba/out_WGS_k25_quitado.fasta --diploid -l 15 -k 55 -o prueba running histogram estimation Linear estimation: ~15 M distinct 38-mers are in the reads K-mer sampling: 1/2 |`

`processing | [going to estimate histograms for values of k: 55 45 35 25 15 ---------------------------------------------------------------------------------Total time Wallclock 34.736 s fitting model to histograms to estimate best k could not fit prueba-k15.histo could not fit prueba-k35.histo could not fit prueba-k25.histo could not fit prueba-k55.histo could not fit prueba-k45.histo could not predict a best k value Execution of decide failed (return code 0)`

```
sudo ./kmergenie /home/soba/Documentos/out_WGS_k25_quitado.fasta -l 15 -k 55 -o prueba
running histogram estimation
Linear estimation: ~15 M distinct 38-mers are in the reads
K-mer sampling: 1/2
| processing |
[going to estimate histograms for values of k: 55 45 35 25 15
---------------------------------------------------------------------------------Total time Wallclock 35.6302 s
fitting model to histograms to estimate best k
could not fit prueba-k25.histo
could not fit prueba-k35.histo
could not fit prueba-k45.histo
could not fit prueba-k55.histo
estimation of the best k so far: 15
refining estimation around [15; 21], with a step of 2
running histogram estimation
Linear estimation: ~19 M distinct 21-mers are in the reads
K-mer sampling: 1/3
| processing |
[going to estimate histograms for values of k: 21 19 17 15
---------------------------------------------------------------------------------Total time Wallclock 22.5893 s
fitting model to histograms to estimate best k
could not fit prueba-k19.histo
could not fit prueba-k17.histo
could not fit prueba-k21.histo
could not fit prueba-k25.histo
could not fit prueba-k45.histo
could not fit prueba-k35.histo
could not fit prueba-k55.histo
table of predicted num. of genomic k-mers: prueba.dat
recommended coverage cut-off for best k: 18
best k: 15
```

Just a note: try not to run programs as root (e.g. do not use sudo).

27kOk. Now it showed the line was failed.

26k• written 5.8 years ago by margxenscienculo •50this one is due to a failure to overwrite the (possibly root) file prueba_report.html

1.4kIt is me, the same. So, May it be a problem with R program? What are the correct installed packages?

26k• written 5.8 years ago by margxenscienculo •50Thanks for reporting this issue. Could you please send the

`prueba_report.html`

file? I would like to see what the histograms that didn't fit look like. It seems to be a tiny dataset.26k• written 5.8 years ago by Rayan Chikhi •1.4kYes, that is. It is a small dataset 25,770,300pb. I have sent email. Regards.

50