This is a test version of Biostars. For the public version, visit https://www.biostars.org.
How to create scaffolds from contigs?

Hi all!

I was using velveth and velvetg to create contigs.fa file. It includes all sequences with kmer=31, min contig length 100 and expected coverage 5.

Now I would like to create scaffolds, which tool do you recommend for this purpose?

Thanks in advance,

Best,

Agata

de novo
• 6,597 views • link
written by agata88 agata88

Spades and IDBA (which I know) but why don´t you play with kmer lenght? With that programs you can make assembly using different kmer lenghts.

*That programs makes scaffolds after contigs assembly but using read`s information.

• 0 views • link
written by Buffo Buffo

Thanks but velvet does not support higher kmer length than 31. I see that Spades and IDBA need fastq files as input not contig file. Do you know any tool that can concatenate contings into scaffolds with the use of contig fasta file?

• 0 views • link
written by agata88 agata88

Sorry I dont know, Reads are used to infer distance between contigs and then make scaffolds, so If you dont have them, it is complicated issue, but; If you have a reference genome (or genome from related specie) you can map to it these contigs ( with Bowtie2 or BWA) and then make scaffolds ussing coordinates if it is posible after mapping.

• 0 views • link
written by Buffo Buffo

Even you can extract consensus sequence from IGV after mapping contigs to reference genome.

• 0 views • link
written by Buffo Buffo

Thanks I'll try that, although I would like to do it as much automatic as possible :) extracting manually consensus sequence is not comfortable in this case.

• 0 views • link
written by agata88 agata88

You can extract consensus sequence automatically from IGV.

• 0 views • link
written by Buffo Buffo

Thanks, SPADES is doing really good job :) I used fastq raw files and quit velvet. Best, Agata

• 0 views • link
written by agata88 agata88

0 answers

No answers yet.

Log in to answer this question.

More posts like this