Question: SPAdes error "Missing @SQ header"
gravatar for Angelo Victoria
4.1 years ago by
Philippines/University of the Philippines Diliman
Angelo Victoria30 wrote:


I am trying to do a de-novo assembly of a prokaryotic genome using SPAdes 3.5.0. I have been getting the following error

[samopen] no @SQ lines in the header.

[sam_read1] missing header? Abort!

I tried doing a subset of the reads and it worked. Thinking that it was a memory issue, I increased the memory allocation from 128GB to 350GB to give it ample memory but, unfortunately, it still didn't work with the full data set.

I think it has something to do with the mismatch correction (--careful) because the assembly was finished but the correction was aborted.

It also created a file named "@SQ SN:NODE_87368_length_231_cov_2.38312_ID_174735 LN:231" in the home directory (outside the spades folder), if that would be of any help.


ADD COMMENTlink modified 4.1 years ago by Philipp Bayer6.6k • written 4.1 years ago by Angelo Victoria30
gravatar for Philipp Bayer
4.1 years ago by
Philipp Bayer6.6k
Philipp Bayer6.6k wrote:

In --careful, Spades uses BWA to align the reads back to the assembly in order to correct minor indels. This step generates the BAM files about which it complains.

There's a good chance that this breaks for your for some reason - could you post the whole STDOUT to or something like that? I have a hunch there's a tiny BWA error hidden somewhere, so that the BAM file isn't generated correctly, maybe there's a bogus read in the whole dataset, etc.

ADD COMMENTlink written 4.1 years ago by Philipp Bayer6.6k
Please log in to add an answer.


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