BayesCall: A model-based base-calling algorithm for high-throughput short-read sequencing

  1. Wei-Chun Kao1,
  2. Kristian Stevens2 and
  3. Yun S. Song1,3,4
  1. 1 Computer Science Division, University of California, Berkeley, Berkeley, California 94720, USA;
  2. 2 Department of Computer Science, University of California, Davis, Davis, California 95616, USA;
  3. 3 Department of Statistics, University of California, Berkeley, Berkeley, California 94720, USA

    Abstract

    Extracting sequence information from raw images of fluorescence is the foundation underlying several high-throughput sequencing platforms. Some of the main challenges associated with this technology include reducing the error rate, assigning accurate base-specific quality scores, and reducing the cost of sequencing by increasing the throughput per run. To demonstrate how computational advancement can help to meet these challenges, a novel model-based base-calling algorithm, BayesCall, is introduced for the Illumina sequencing platform. Being founded on the tools of statistical learning, BayesCall is flexible enough to incorporate various features of the sequencing process. In particular, it can easily incorporate time-dependent parameters and model residual effects. This new approach significantly improves the accuracy over Illumina's base-caller Bustard, particularly in the later cycles of a sequencing run. For 76-cycle data on a standard viral sample, phiX174, BayesCall improves Bustard's average per-base error rate by ∼51%. The probability of observing each base can be readily computed in BayesCall, and this probability can be transformed into a useful base-specific quality score with a high discrimination ability. A detailed study of BayesCall's performance is presented here.

    Footnotes

    | Table of Contents

    Preprint Server