Figbird: a probabilistic method for filling gaps in genome assemblies

Abstract

Motivation Advances in sequencing technologies have led to the sequencing of genomes of a multitude of organisms. However, draft genomes of many of these organisms contain a large number of gaps due to the repeats in genomes, low sequencing coverage and limitations in sequencing technologies. Although there exist several tools for filling gaps, many of these do not utilize all information relevant to gap filling. Results Here, we present a probabilistic method for filling gaps in draft genome assemblies using second generation reads based on a generative model for sequencing that takes into account information on insert sizes and sequencing errors. Our method is based on the expectation-maximization (EM) algorithm unlike the graph based methods adopted in the literature. Experiments on real biological datasets show that this novel approach can fill up large portions of gaps with small number of errors and misassemblies compared to other state of the art gap filling tools. Availability and Implementation The method is implemented using C ++ in a software named “Filling Gaps by Iterative Read Distribution (Figbird)”, which is available at: https://github.com/SumitTarafder/Figbird. Supplementary information Supplementary data are available at Bioinformatics online.