Dear scientists,
I have performed de novo assembly of bacterial WGS data using spades and abyss. However I am not sure which assembly approach is good to go because N50 value is not half or more than half of total assembly length.
All reads are QC passed (no adapter or primer sequences are found and base quality is above 30).
And also please suggest that how many contigs should I keep for further analysis. Length of reference genome is 5,240,075 bp.
Here i am attaching assembly reports please input your valuable suggestions and guidance.
Assembly with abyss
Assembly with spades
Assembly with spades (using --trusted-contigs flag where reference genome was used to guide assembly)
Assembly with unicycler
Your SPAdes assembly looks fine. N50 is not supposed to be half a genome. I would take a >5000 bp as a contig cutoff and that should be similar to your reference genome. This can be confirmed by mapping. You can lower the threshold maybe wo 2000 bp, but that will hardly be different from 5 kb.
Thank you @Mensur Dlakic for your valuable suggestion!
Please have a look on stats again I mistakenly pasted stats of spades assembly with --trusted-contig flag twice. Adding --trusted-contigs parameter brought difference in total contigs number and N50 value as well.
Please suggest which spades approach should I follow?
Not sure why you need me to answer this question, as the --trusted-contigs assembly is clearly the best. It has fewest contigs, largest N50 and the largest contig overall.
Thank you for your answer. Best,
Log in to answer this question.
Looks like you pasted the stats for plain
SPAdesassembly twice or vice versa. Since the last two stats are identical.Thanks for informing GenoMax. I mistakenly pasted stats of spades (with --trusted-contigs option) twice. I have edited the post.
N50 is the length of the size-ordered contig that puts you over that 50% threshold
http://jermdemo.blogspot.com/2008/11/calculating-n50-from-velvet-output.html
Thank you so much Jeremy for sharing useful link.