This is a test version of Biostars. For the public version, visit https://www.biostars.org.
The --contig option for SPAdes

Hello everyone,

I am using the SPAdes with --trusted-contigs and --untrusted-contigs.

I want to know how these contigs exactly help with graph construction and path extension. Is any document mentioned the related algorithms?

assembly spades contig
• 6,001 views • link
updated by dark.lord • dark.lord • written by 934963534 • 934963534

Taken from the SPAdes manual:

Additional contigs

In case you have contigs of the same genome generated by other assembler(s) and you wish to merge them into SPAdes assembly, you can specify additional contigs using --trusted-contigs or --untrusted-contigs. First option is used when high quality contigs are available. These contigs will be used for graph construction, gap closure and repeat resolution. Second option is used for less reliable contigs that may have more errors or contigs of unknown quality. These contigs will be used only for gap closure and repeat resolution. The number of additional contigs is unlimited.

Note, that SPAdes does not perform assembly using genomes of closely-related species. Only contigs of the same genome should be specified.

• 0 views • link
written by Sej Modha Sej Modha

I have read it before and wonder how the contigs help with graph construction, gap closure and repeat resolution, since they should not be considered just as single read.

• 1 views • link
written by 934963534 • 934963534

1 answer

My understanding of this is that the contigs are used as a backbone for assembly. Using these scaffolds/contigs as a backbone, the assembly becomes more of a mapping - the raw (trimmed) reads stretch the contigs, and that's how it helps with gap closure. The same thing could happen with the mapping + re-assembly, but this would require many more steps because no graph is constructed, whereas in SPAdes it is.

In numbers - this also helped me to better understand - say we have a 1000 bp long contig. A 250-bp long read overlaps the terminal 100 bp of the contig. This means that there is a 150 bp long stretch over the tip of the 1000 bp contig, and on this stretch, other reads can overlap - something that doesn't happen when you map.

Hope it helps.

Cheers

Stefano

• 1 views • link
written by dark.lord • dark.lord

Log in to answer this question.

More posts like this