Problem: kmergenie has a minimum number of input of reads?
Hi all, I'm the following problem with kmergenie:
Warning! using max number of read files (2000)
error opening file: #
fitting model to histograms to estimate best k
could not predict a best k value
Execution of decide failed (return code 0)
How can I increase the read number of input??
Other question is how much of memory I need to run kmergenie?
Thanks,
Leandro
• 5,845 views
•
link
updated
by
Ram
•
written
by
Leandro de Mattos
•
written
by
Leandro de Mattos
0 answers
No answers yet.
Log in to answer this question.
More posts like this
-
not best k value found
written by shaileshdesai76 •
-
kmergenie show nothing
written by zhangmengjie •
-
kmergenie does not predicted a best k value
written by StudentBio •
-
No best K value found
written by kspata
-
Question about k-mer genie output
written by stacy734 •
-
Unable to get kmergenie to work, but make check seems fine
written by jeremycfd •
-
kmergenie - Error 1 - It is doesnt work when I typed 'make'
written by margxenscienculo
-
The problem to run kmergenie
written by ziyingliu1996 •
-
Execution of decide failed (return code 0)
written by Anaïs Vittu •
-
No could not fit >15 [Assemblied by SOAPdenovo]
written by margxenscienculo
Hi!
What command line did you use? Also, what operating system?
Dear Rayan, I used :
./kmergenie *.csfastaOperating systems: Biolinux (last version)
I have installed the R and python.
Thanks. How many
*.csfastafiles do you have? Kmergenie indeed has a limit on the number of input files (2000), as mentioned in the error. You could try merging them, like this:then run
Ah also: Kmergenie doesn't (yet) work with an input like
*.fasta. (it might in the future; right now is version 1.6950)If you have a list of fasta, please do the following:
Hi, I use this command it doesn't work. Instead, it shows:
here is my list file:
This looks like a working directory problem. The ~/data/list file does not seem to contain absolute paths, thus you need to run kmergenie inside the ~/data/ folder.
Thanks a lot :)
Hi Rayal,
I have a side question, would you please clarify it for me:
if I run kmergenie for a pair-end read set (contain Read 1 and Read 2 fastq files), do I need to find a way to translate Read 2 into its compliment sequence before combining with Read 1 for kmergenie run? (because all sequence infomation in Read 2 is compliment to Read 1). If NOT, would it double the number of distinct kmer in statistical calculation of kmergenie? Overall, what we want to know is only 01 single strain of DNA only, isn't it?
Sorry, I am very new to this field. Thank you very much in advance!
Phuong
Hi Phuong,
No need. Kmergenie does not care if a read is in forward or reverse orientation, also does not care about reads are paired-end or single-end or mate-pairs. Just input all the fastq files that you would give to an assembler, in any order.
It won't double the number of kmers, as, kmergenie considers that a kmer and its reverse complement are the same object.
Thank you very much, this really enlightens me, especially the fact that kmergenie considers a kmer and its reverse complement are the same object.
Hi Hian, but I used only one input, I have one file.
I see.. Can you please paste the output of the following commands?
(By the way, Biostars encourages that you respond in a reply, not in a separate response, which is reserved for when an answer to the original problem is found)