voutcn / megahit

Ultra-fast and memory-efficient (meta-)genome assembler
http://www.ncbi.nlm.nih.gov/pubmed/25609793
GNU General Public License v3.0
591 stars 135 forks source link

Getting started #1

Closed accopeland closed 10 years ago

accopeland commented 10 years ago

Hi,

I cloned the repo and built megahit without problems, but am having some trouble getting a test run to complete (start) without errors. It's not quite clear from the docs if you require or support fasta, fastq, and/or fastq.gz so I've tried both fastq and fasta as follows:

python ../megahit -r uazx.sample.fa -m 60000000000 -l 150

Which produces the following error

Number of CPU threads 32 [Tue Sep 30 09:03:47 2014]: Extracting solid (k+1)-mers for k = 21 sdbg_builder_cpu: cx1_functions.cpp:192: void phase1::ReadInputFile(global_data_t&): Assertion `packed_reads != __null' failed. Error occurs when running "builder count"

Can you advise me on a correct command line and any specific input requirements e.g. paired, interleaved, fasta, etc. ? Thanks in advance.

voutcn commented 10 years ago

Hi Alex,

You command line is correct. MEGAHIT accept fasta or fastq file as input. The program crashed because of insufficient memory. In the command line, the -m option was set to be 60000000000 (~56GB). Make sure you have enough memory. You may try a smaller number for the -m option.

Regards, Dinghua

On 1 Oct, 2014, at 12:06 am, Alex Copeland notifications@github.com wrote:

Hi,

I cloned the repo and built megahit without problems, but am having some trouble getting a test run to complete (start) without errors. It's not quite clear from the docs if you require or support fasta, fastq, and/or fastq.gz so I've tried both fastq and fasta as follows:

python ../megahit -r uazx.sample.fa -m 60000000000 -l 150

Which produces the following error

Number of CPU threads 32 [Tue Sep 30 09:03:47 2014]: Extracting solid (k+1)-mers for k = 21 sdbg_builder_cpu: cx1_functions.cpp:192: void phase1::ReadInputFile(global_data_t&): Assertion `packed_reads != __null' failed. Error occurs when running "builder count"

Can you advise me on a correct command line and any specific input requirements e.g. paired, interleaved, fasta, etc. ? Thanks in advance.

— Reply to this email directly or view it on GitHub.

voutcn commented 10 years ago

Hi Alex,

Did MEGAHIT complete a test run successfully? If yes I can close this issue. BTW you may want to try the latest release of MEGAHIT.

Dinghua

accopeland commented 10 years ago

Yes, it did complete correctly thanks. Go ahead and close. I will try the latest version as well.

On Sun, Oct 12, 2014 at 1:04 AM, Dinghua Li notifications@github.com wrote:

Hi Alex,

Did MEGAHIT complete a test run successfully? If yes I can close this issue. BTW you may want to try the latest release of MEGAHIT.

Dinghua

— Reply to this email directly or view it on GitHub https://github.com/voutcn/megahit/issues/1#issuecomment-58776886.

Alex Copeland Genome Assembly Group Lead Joint Genome Institute, Walnut Creek, CA 94598