I have applied two de novo genome assembly tools to data from the run SRR12707453, corresponding to a phage (I downloaded the data and have no relation with the authors of the study).

Using rnaSPAdes I get many long contigs (>50.000 bp). When analyzing the contigs using Quast against a reference genome (MN813694.1) most long contigs are discarded and the seventh-longest one is the only one selected. If I blast the other contigs I find correspondence with other phages, which makes me conclude that many phages were present in the same probe.

This seventh contig has 47686 bp length, while the reference has 47713 bp length. But they are not equivalently arranged. Apparently, the genome was circularized and then it was broken around the position 32400 in the reference. Also, the sequence does not match between positions 45000 and 45951 (in the reference), and a little misassembly seems to occur around position 46300.

See here some screen captures for the Icarus visualization obtained with Quast:

enter image description here enter image description here

In contrast, using Quast with the contigs obtained from the other assembly tool, I get some differences. There is also a contig with almost the same size as the reference genome (47573 bp) that has exactly the same structure than the contig obtained with rnaSPAdes: circularized, broken around position 32400, does not match between position 45000 and 45951, and small misassemblie aronund position 46300. But there are also some contigs in which there is no break around position 32400.

enter image description here

May I interpret that these contigs indicate that the phage genome is sometimes circularized and sometimes broken at position ~32400 and that's why some contigs don't find a sequence break at that position?

Am I interpreting anything wrong? Any other conclusion can be obtained from the data?

You can (temporaly) download here the Quast results.


1 Answer 1


some additional information to your request.


it contain assembled by reference graph using ncbi tool saute(skesa) which will be public later.

install bandage app and load this graph https://rrwick.github.io/Bandage/

picture looks intriguing,I would say , run contain many errors located in one position. Reason? I don’t really know.

  • $\begingroup$ Yep... I don't know how circularization works. I guess it could be that binding the sequence generates every time some garbage in the union and so every RNA chain has a slightly different sequence at this point. $\endgroup$
    – juanjo75es
    Commented Nov 4, 2020 at 21:14

Your Answer

By clicking “Post Your Answer”, you agree to our terms of service and acknowledge you have read our privacy policy.

Not the answer you're looking for? Browse other questions tagged or ask your own question.