This is a test version of Biostars. For the public version, visit https://www.biostars.org.
bacterial genome assembly output from canu

Hi,

I have done de novo genome assembly of a bacterial strain using canu. I want to find structural variants comparing with the reference genome. Since, the bacterial genome is circular, It is hard to find the origin of replication to align with the reference genome. Does anyone have suggestions how to find structural variants in the bacterial genome comparing with the reference genome?

Thanks

bacteria genome assembly canu

Yes, I tried to repeat the assembly after removing read lengths shorter than 2000 bp. The assembly result is similar. I think the issue is due to incorrect identification of origin of replication during genome assembly.

Where is the dnaA gene located in your assembly and where is it in the reference? Paper1 and Paper2.

2 answers

There are many ways to do this. I suggest to do the following:
1) Align your genome to the reference genome using pairwise megablast
2) From the alignment you'll be able to find the position of your genome that corresponds to the first position in the reference. Then change your genome sequence, so now it starts from this position.
3) Align your genome to the reference genome again and look at the dotplot. You'll be able to see structural differences on the dotplot.

I think Quast is a good tool for that. It already aligns the assembly to the reference independently of any issue with circularity. I think it's also useful to get two different assemblies with two different software. Sometimes it's just the assembler that fails.

Here you have a likely real rearrangement verified by three different algorithms (SPAdes, rnaSPAdes and Contignant s-aligner):

Alignment of reads obtained with SPAdes, rnaSPAdes and Contignant s-aligner

Here you have a false rearrangement detected with SPAdes

Alignment of reads obtained with SPAdes

Log in to answer this question.