Question: Velvet can't handle k-mers as long as X! We'll stick to 32 if you don't mind.
0
gravatar for emmaielle
4.3 years ago by
emmaielle0
Uruguay
emmaielle0 wrote:

I am trying to create the k-mer hash table with velveth, with this command line:

velveth directory 15,49 -fasta -shortPaired interleavedInputFile.fasta

My issue here is that I get "Velvet can't handle k-mers as long as 49! We'll stick to 32 if you don't mind". I searched in other posts and figured I need to change the Makefile from Velvet directory, where 'MAXKMERLENGTH=31', to the number I wished.

After changing this, there was no change.

I also tried by writing make ’MAXKMERLENGTH=57’ at the directory when compiling, but the terminal just doesn't respond.

What am I doing wrong?

Thanks!

kmer length velvet assembly • 1.9k views
ADD COMMENTlink modified 4.2 years ago by blanca10 • written 4.3 years ago by emmaielle0

The make MAXKMERLENGTH=57 step should fix it. What do you mean by "terminal won't respond"? 

ADD REPLYlink written 4.3 years ago by RamRS21k

I meant that it didn't compile and just showed me a ">" as in that it kept thinking indefinitely.

I just tried make MAXKMERLENGTH=57 without quotes and it compiled. I had taken it literally from Velvet Can'T Handle K-Mers As Long As ....

The problem is that it didn't do a thing, because the velvet directory still doesn't have velveth and velvetg, and Makefile stil has MAXKMERLENGTH?=31

ADD REPLYlink modified 4.3 years ago • written 4.3 years ago by emmaielle0

Could you check if the output of which velveth points to the directory with the velvet binaries you just compiled? Maybe the velvet you compiled is not the velvet being used.

ADD REPLYlink written 4.3 years ago by RamRS21k

Yep, it does point to the directory I just compiled.

ADD REPLYlink written 4.3 years ago by emmaielle0

Really weird.

ADD REPLYlink written 4.3 years ago by RamRS21k

If you say it's possible to compile by make MAXKMERLENGTH=57, then it is weird. I'll try again later just in case, though.

ADD REPLYlink modified 4.3 years ago • written 4.3 years ago by emmaielle0
1

Try make clean, then make MAXKMERLENGTH=57

ADD REPLYlink written 4.3 years ago by RamRS21k

Well, that worked perfectly. I don't know what had happened before. Thanks!

ADD REPLYlink written 4.3 years ago by emmaielle0
1

make decides which files to compile based on if the prereq files changed. Because none of them did, it did not recompile. make clean removed all previously compiled files and forced a re-compile from scratch, and took your MAXKMERLENGTH into account when it compiled this time.

This is a bug. make with added parameters should make clean automatically.

ADD REPLYlink written 4.3 years ago by RamRS21k
0
gravatar for emmaielle
4.3 years ago by
emmaielle0
Uruguay
emmaielle0 wrote:

Okay, so I figured what was I doing wrong:

I changed the MakeFile into MAXKMERLENGTH=57 before compiling, and then just used make.

Now it works perfectly. I didn't know the make basics very well.

ADD COMMENTlink written 4.3 years ago by emmaielle0

That is weird. make MAXKERLENGTH=57 should work. Hmmm.

ADD REPLYlink written 4.3 years ago by RamRS21k
0
gravatar for blanca
4.2 years ago by
blanca10
Spain
blanca10 wrote:

I am sorry I still have not clear how to compile Velvet with a longer K-mer. Either changing the MAXKMERLENGTH option in Makefile and compiling with "make" afterwards or changing that option with make MAXKMERLENGTH=57 (with or without quotes) don't seem to work (still showing the error "Velvet can't handle k-mers as long as ....").

I have make clean and even removed the directory each time so that Velvet does not get confused (because of that mentioned bug).

 

Does someone have any advice? any other option I have still not tried?

Thank you!

 

 

EDITED: Just find out. In case anyone have this same problem, I fixed it like mentioned in here: 

Question: Can't change MAXKMERLENGTH on Velvet

ADD COMMENTlink modified 4.2 years ago • written 4.2 years ago by blanca10

What's the output of this command:

which velveth

Is this the directory you're running make MAXKMERLENGTH=57 in?

ADD REPLYlink written 4.2 years ago by RamRS21k
Please log in to add an answer.

Help
Access

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 1315 users visited in the last hour