Entering edit mode
9.5 years ago
abujamel.t
•
0
I am using kmergenie version 1.6473 to estimate the optimal kmer for my assembly. The program runs perfect, but my problem is that it only generates .histo files but no histograms.
Here is the command I used:
kmergenie unligned_R1.fastq unligned_R2.fastq -o kmer -s 2 -l 13 -t 32
Linear estimation: ~51 M distinct 71-mers are in the reads
K-mer sampling: 1/15
| processing |
[going to estimate histograms for values of k: 121 119 117 115 113 111 109 107 105 103 101 99 97 95 93 91 89 87 85 83 81 79 77 75 73 71 69 67 65 63 61 59 57 55 53 51 49 47 45 43 41 39 37 35 33 31 29 27 25 23 21 19 17 15 13
-----------------------------------------------------------------------------------------------------------------------Total time Wallclock 628.354 s
The result was only a list of histo files from kmer-k13.histo to kmer-k121.histo BUT no graphs!!!
Could anyone help me solve this problem?
Cheers,
TJ
Hi, can you post the list of files that it generated, as well as the first few lines of one of the .histo files?
Also, did it output any error in stderr?
Hi Rayan,
the files list is as following:
kmer-k121.histoshows:there was no error
TJ
Thanks. Is this an actual sequencing dataset? There shouldn't be so few kmers that appear once, twice, 3 times, etc.. (that's the meaning of the hist files)
Or is this a simulated dataset? Kmergenie only works with real data, or simulated according to realistic sequencing scenario
it is an actual sequencing dataset of metagenomic sequencing generated from Illumina HiSeq2500 sequencing
TJ
Thanks. Histogram generation should work on that data. (Note, however, that Kmergenie is not designed for metagenomic data).
Let's see if the fault comes from your system or from that specific dataset. Could you try a simple (small, for quick execution) genomic dataset, using the same parameters
-o kmer -s 2 -l 13 -t 32, and see if kmergenie completes successfully on your system?Thanks Rayan for your reply,
I tried sample paird end fastq files with 750K reads and 36 reads length. Same problem.
TJ
Alright, so something is wrong with either the system or the software. I checked that command line, it works fine on my computer on a sample dataset.
What's the output of the "make check" command in kmergenie folder, does it print an error?
What's your system? (desktop, cluster?)
Could you try running kmergenie without any option except
-o kmeron that sample data?This is the output of
make checkMy system is a Dell workstation running Biolinux 8 (Ubuntu 14.04)
I tried running the command as you specified, same thing happened but it tested only two kmers: 21 and 31 (only histo files no graphs)
TJ
Thanks for this information. Your install looks normal, I'm getting puzzled.
One more debug I can think of, what's the output of the following command, executed in the folder where your .histo files reside, (replace XXX by the path to kmergenie folder):
(kmer is the string given to the
-oparameter)This should run the second phase of Kmergenie manually and print some debug information in stdout.
Hi Rayan,
it worked! it generated the graphs from the histo files.
So what do you think the issue is?
TJ
For some reason the "decide" script doesn't seen to be getting executed. That is very strange, as the histogram creation program is correctly executed; both programs are called the same way in the kmergenie main program. To be honest, I do not know why this happens, and havn't seen that with any other user. Can you please paste the full output of the decide command, maybe I'll find something unusual?
Was this ever resolved? I have the same issue, thanks.
I haven't heard back from TJ. Could you please execute the following command, and paste the full output, so that I may gain some insight into the problem?
where
[prefix]is the string given to the-oparameter.I get a list of histo files but no pdfs when I disown kmergenie and log out (also does not proceed into the second round). When I disown only and keep the terminal alive, it runs fine. Don't know if this might provide any insight.
Cheers, John