I've been using the latest version of Velvet (1.2.10) and the latest version of Metavelvet (1.2.02) to assemble metagenomic datasets of pair-end reads of about 120 bp length.
I use Velvetoptimiser to choose the appropriate k-mer then run vevleth again (as for some reason velvetoptimiser doesn't produce the Roadmaps file) and then try to run metavelvet. But when I do so I get a segmentation fault error like this:
[126.515876] Starting pebble resolution...
Segmentation fault (core dumped)
It's not a memory problem as I am running into a 512GB RAM machine (64 cpus) and the process is quite quick. (I've also checked with free -m)
I've compiled both velvet and metavelvet at 'MAKMERLENGTH=101', 'OPENMP=1'.
My server is running on CentOS 7.1.1503 with g++ version 4.8.5.
Any ideas anyone?
Oh and Merry Christmas!