TransWikia.com

Viral genome finishing

Bioinformatics Asked on July 17, 2021

I have assembled poxvirus genome using Ray. The assembly is good. Out of several thousand contigs I got, I was able to get one scaffold using Contiguator tool, which is about 90% of my genome. I have well annotated reference sequence. If there any good tool for filling gaps and checking assembly quality?

Things I’ve tried so far:

  1. IMAGE from PAGIT. It asks for ace file which I guess is obtained from Newbler. So, I am not able to use this.

  2. gapfill_py script from broad: There are some gaps still missing.

I came to know that Consed is good. But, I have windows which is not compatible for Consed.

I checked my BAM file with Tablet. There are reads that cover almost entire gene. Are there any good tools or manual ways to finish my assembly so that I can get the full genome?

One Answer

Assuming that gap filling is a domain-independent and works the same for viruses as for bacteria or eukarytes.

There are two real options I am aware of: Sealer a successor of GapFiller from ABySS package and GapCloser from SOAPdenovo2 package. According to Sealer paper it's better, but I am not sure if any indepented benchmarking have confirmed their conclusion. On my insect data GapCloser worked better than GapFiller (back then I have not tried to compare it to Sealer).

Answered by Kamil S Jaron on July 17, 2021

Add your own answers!

Ask a Question

Get help from others!

© 2024 TransWikia.com. All rights reserved. Sites we Love: PCI Database, UKBizDB, Menu Kuliner, Sharing RPP