Spectrum Alignment: Efficient Resequencing by Hybridization
Proceedings of the Eighth International Conference on Intelligent Systems for Molecular Biology
Handling long targets and errors in sequencing by hybridization
Proceedings of the sixth annual international conference on Computational biology
Restricting SBH Ambiguity via Restriction Enzymes
WABI '02 Proceedings of the Second International Workshop on Algorithms in Bioinformatics
De novo repeat classification and fragment assembly
RECOMB '04 Proceedings of the eighth annual international conference on Resaerch in computational molecular biology
Restricting SBH ambiguity via restriction enzymes
Discrete Applied Mathematics
Hi-index | 0.01 |
Sequencing by hybridization is a method for reconstructing a DNA sequence based on its k-mer content. This content, called the spectrum of the sequence, can be obtained from hybridization with a universal DNA chip. However, even with a sequencing chip containing all 49 9-mers and assuming no hybridization errors, only about 400 bases-long sequences can be reconstructed unambiguously.Drmanac et al. suggested sequencing long DNA targets by obtaining spectra of many short overlapping fragments of the target, inferring their relative positions along the target and computing spectra of subfragments that are short enough to be uniquely recoverable. Drmanac et al. do not treat the realistic case of errors in the hybridization process. In this paper we study the effect of such errors. We show that the probability of ambiguous reconstruction in the presence of (false negative) errors is close to the probability in the errorless case. More precisely, the ratio between these probabilities is 1 + &Ogr;(p/(1 - p)4 · 1/d) where d is the average distance between neighboring subfragments, and p is the probability of a false negative.We also obtain lower and upper bounds for the probability of unambiguous reconstruction based on errorless spectrum. For realistic chip sizes, these bounds are tighter than those given by Arratia et al. Finally, we report results on simulations with real DNA sequences, showing that even in the presence of 50% false negative errors, a target of cosmid length can be recovered with less than 0.1% miscalled bases.