Safe and Complete Contig Assembly Through Omnitigs

Alexandru I. Tomescu, Paul Medvedev

Research output: Contribution to journalArticlepeer-review

2 Scopus citations

Abstract

Contig assembly is the first stage that most assemblers solve when reconstructing a genome from a set of reads. Its output consists of contigs - a set of strings that are promised to appear in any genome that could have generated the reads. From the introduction of contigs 20 years ago, assemblers have tried to obtain longer and longer contigs, but the following question remains: given a genome graph G (e.g., a de Bruijn, or a string graph), what are all the strings that can be safely reported from G as contigs? In this article, we answer this question using a model in which the genome is a circular covering walk. We also give a polynomial-time algorithm to find such strings, which we call omnitigs. Our experiments show that omnitigs are 66%-82% longer on average than the popular unitigs, and 29% of dbSNP locations have more neighbors in omnitigs than in unitigs.

Original languageEnglish (US)
Pages (from-to)590-602
Number of pages13
JournalJournal of Computational Biology
Volume24
Issue number6
DOIs
StatePublished - Jun 2017

All Science Journal Classification (ASJC) codes

  • Modeling and Simulation
  • Molecular Biology
  • Genetics
  • Computational Mathematics
  • Computational Theory and Mathematics

Fingerprint Dive into the research topics of 'Safe and Complete Contig Assembly Through Omnitigs'. Together they form a unique fingerprint.

Cite this