This is a test version of Biostars. For the public version, visit https://www.biostars.org.
Alignment of mutli-contig fasta files

Hi,

I'm new to this so apologies if this question has been answered elsewhere:

I've recently sequenced 70 similar bacterial isolates using a MiSeq and assembled them using Spades. Each resulting fasta file has about 500 contigs.

My aim is to make each of these files into a single contig/scaffold, so that I can align with MUSCLE and then use FASTTREE to build the phylogeny.

Is there any way to do this?

assembly alignment next-gen
• 3,012 views • link
updated by vmicrobio vmicrobio • written by L Ron Hubbard • L Ron Hubbard

My aim is to make each of these files into a single contig/scaffold

If SPAdes was not able to create a single contig perhaps your sequence data is not complete. How can you fill-in data that is missing? Are you hoping to borrow from across samples i.e. essentially make single hybrid assembly?

• 0 views • link
written by GenoMax GenoMax

With bacteria, assemblies virtually always have multiple contigs. I think due to repeat elements. Can anyone confirm this?

• 0 views • link
written by L Ron Hubbard • L Ron Hubbard

1 answer

you can use scaffold builder (source forge or online) to map your contigs against a reference in order to get a scaffold

• 0 views • link
written by vmicrobio vmicrobio

Hi Haro, I tried that but the scaffold still produces the same number of contigs. Any ideas?

• 0 views • link
written by L Ron Hubbard • L Ron Hubbard

Hi, you may try quast to evaluate the quality of the assembly. Maybe the coverage is not homogenous or there is a lot of repeated sequences which may affect the assembly

• 0 views • link
written by vmicrobio vmicrobio

Log in to answer this question.

More posts like this