SPAdes contigs

[vpi]

Hi,

I am new to SPAdes and I have some questions. I wonder about reads, like how many reads make a contig. I don’t see any information in the SPAdes.log file. Is there a way to see the number of reads for each contig?

[ctseto]

Velvet used to store read information differently; which then allowed extraction of that information via a perl script.

It may be worth mapping your original reads back to your contigs anyways-to check for changes in coverage, or to see if there are spurious gaps. Programs like Opera and REAPr rely on this mapping of the starting reads to the finished product as part of their quality control.

[Mario Moreno]

Hello all!
I am new to SPAdes too,
I have the contigs number but I dont have info about the number of read used for make the contigs. there any way to know this information without using external programs?
Thanks

[Dmitry.Antipov]

Hi.
SPAdes do not output the exact number of reads mapped to the contig. However you can estimate it with help of contig length and average contig coverage (NB, it is kmer coverage for the last k used for assembly) that can be found in contig's fasta header.