NZ796095A - Phasing correction - Google Patents

Abstract

Memory efficient methods determine corrected color values from image data acquired by a nucleic acid sequencer during a base calling cycle. Such methods may: (a) obtain an image of a substrate (e.g., a portion of a flow cell) including a plurality of sites where nucleic acid bases are read; (b) measure color values of the plurality of sites from the image of the substrate; (c) store the color values in a processor buffer of the sequencer’s one or more processors; (d) retrieve partially phase-corrected color values of the plurality of sites, where the partially phase-corrected color values were stored in the sequencer’s memory during an immediately preceding base calling cycle; (e) determine a prephasing correction; and (f) determine the corrected color values. In various implementations, these operations are all performed during a single base calling cycle. In certain embodiments, the methods additionally include using the corrected color values to make base calls for the plurality of sites. Sequencers may be designed or configured to implement such methods.

Description

Memory efficient methods determine corrected color values from image data acquired by a nucleic acid sequencer during a base calling cycle. Such methods may: (a) obtain an image of a substrate (e.g., a portion of a flow cell) including a plurality of sites where nucleic acid bases are read; (b) measure color values of the plurality of sites from the image of the substrate; (c) store the color values in a processor buffer of the sequencer’s one or more processors; (d) retrieve lly phase-corrected color values of the plurality of sites, where the lly phase-corrected color values were stored in the sequencer’s memory during an immediately ing base calling cycle; (e) determine a prephasing tion; and (f) determine the corrected color values. In various implementations, these operations are all performed during a single base g cycle. In certain embodiments, the methods additionally include using the corrected color values to make base calls for the plurality of sites. Sequencers may be designed or configured to implement such methods.

NZ 796095 PHASING CORRECTION CROSS-REFERENCE TO RELATED APPLICATIONS This ation is a onal application of New d Patent Application No. 754912, a National Phase Entry of International Patent Application No. .

International Patent Application No. claims the benefits of U.S.

Provisional Patent Application No. 62/443,294, filed January 06, 2017, and entitled "PHASING CORRECTION," which is hereby incorporated herein by reference in its entirety and for all purposes. The content of New Zealand Patent ation No. 754912 is hereby incorporated by reference herein in its entirety and for all es.

BACKGROUND The disclosure relates to sequencing nucleic acids. More specifically, the disclosure relates to systems and methods for real time sequencing with phasing corrections.

At a particular site on a flow cell or other substrate, multiple copies of a nucleic acid molecule, all having the same sequence (possibly with limited ions unintentionally introduced by sample processing), are ed together. Enough copies are used to ensure that sufficient signal is produced to permit reliable base calling. The collection of nucleic acid molecules at a site is called a cluster.

Phasing represents an unintended ct that arises from sequencing multiple nucleic acid molecules within a cluster. Phasing is the rate at which signals such as fluorescence from single molecules within a cluster lose sync with each other. Often the term phasing is reserved for contaminating signal from some molecules that fall , and the term pre-phasing is used for contaminating signal from other molecules that go ahead. Together phasing and pre-phasing describe how well the sequencing apparatus and chemistry is performing.

SUMMARY n aspects of this disclosure pertain to methods of determining corrected color values from image data acquired by a nucleic acid sequencer during a base calling cycle, where the sequencer includes an image acquisition system, one or more processors, and memory. Such methods may be terized by the following operations:(a) obtaining an image of a substrate (e.g., a portion of a flow cell) including a plurality of sites where nucleic acid bases are read; (b) ing color values of the plurality of sites from the image of the substrate;(c) storing the color values in a processor buffer of the sequencer’s one or more processors; (d) retrieving partially phase-corrected color values of the plurality of sites, where the partially phase-corrected color values were stored in the sequencer’s memory during an immediately preceding base calling cycle; (e) determining a prephasing correction; and (f) determining the ted color values. In various implementations, these operations are all performed during a single base calling cycle. In certain embodiments, the methods onally include using the corrected color values to make base calls for the plurality of sites.

During sequencing, the sites exhibit colors enting nucleic acid base types.

The measured and stored color values may be intensity or other magnitude values at a particular wavelength or range of wavelengths. In some implementations, the color values are determined from only two channels of the sequencer. In some implementations, the color values are obtained from four channels of the sequencer. While this sure s on phasing correction of color signals, the concepts apply to other types of signals generated during sequencing clusters of c acids having identical sequences. Examples of such other signals include radiation outside the visible spectrum, ion concentration, etc.

In certain ments, determining the corrected color values in (f) uses (i) the color values in the processor buffer, (ii) the partially phase corrected values stored during the immediately preceding cycle, and (iii) the prephasing correction. In certain embodiments, determining the pre-phasing correction in (e) uses (i) the partially phase-corrected color values stored during the immediately ing base calling cycle, and (ii) the color values stored in the processor buffer.

In certain embodiments, the prephasing tion includes a weight. In such embodiments, the operation of determining the corrected color values may e multiplying the weight by the color values of the plurality of sites measured from the image of the substrate.

In certain implementations, the methods additionally include determining a g correction for the immediately succeeding base calling cycle. As an example, determining the phasing correction for the immediately succeeding base calling cycle es ing (i) the partially phase-corrected color values stored in the sequencer’s memory, and (ii) the color values stored in the sor buffer. In certain embodiments including determining a phasing correction for the immediately succeeding base calling cycle, the methods additionally include (i) producing partially phase-corrected color values for the immediately succeeding base calling cycle by applying the phasing correction to color values of the plurality of sites stored in the sequencer’s ; and (ii) storing the partially phase- corrected color values for the immediately succeeding base calling cycle in the sequencer’s memory. In certain embodiments, producing the partially phase-corrected color values for the immediately succeeding base calling cycle additionally includes summing (i) the phasing corrected color values of the plurality of sites, and (ii) the color values of the plurality of sites from the image of the ate measured in (b). In some implementations, storing the partially phase-corrected color values for the immediately succeeding base calling cycle stores the partially-corrected color values in tile buffers of the sequencer’s .

In certain embodiments, the methods are performed in real time during acquisition of sequence reads by the nucleic acid sequencer. In n embodiments, the nucleic acid sequencer sequences by synthesizing nucleic acids at the plurality of sites. In certain embodiments where the substrate includes a flow cell, the flow cell is logically d into tiles, and each tile ents a region of the flow cell comprising a subset of sites, which subset is captured in a single image from the image acquisition .

In some embodiments ing such systems, in ion (d) (retrieving partially corrected color values of the plurality of sites), the partially phase-corrected color values were previously stored in tile buffers of the sequencer’s memory, where the tile buffers are ated for storing data representing images of individual tiles on the ate.

In certain embodiments, the memory has a storage capacity of about 512 Gigabytes or less, or about 256 Gigabytes or less. In certain embodiments, for example, the memory has a storage ty of less than twice the capacity ed to store the data contained in the total number of tiles on two flow cells. In some embodiments, the sing described herein saves at least about 50 Gigabytes; in some embodiments it saves at least about 100 Gigabytes.

In some implementations, prior to operation (a) (obtaining an image of a substrate), the methods additionally include providing reagents to the flow cell and allowing the reagents to interact with sites to exhibit the colors representing nucleic acid base types during the base calling cycle. In such implementations, the method may additionally include, after operation (f) (determining the corrected color values): (i) providing fresh reagents to the flow cell and allowing the fresh reagents to interact with the sites to exhibit colors representing nucleic acid base types for a next base calling cycle; and (ii) repeating ions (a)-(e) for the next base calling cycle. Such methods may additionally include creating a first sor thread for performing operations (a)-(f) for the base calling cycle, and creating a second processor thread for performing operations ) for the next base calling cycle. In certain embodiments, the methods additionally include allocating the sor buffer and a second processor buffer, where the second sor buffer is used to determine the corrected color values in (f).

Certain other aspects of the disclosure pertain to nucleic acid sequencers which may be characterized by the following elements: an image acquisition system; ; and one or more processors designed or configured to: (a) obtain data representing an image of a substrate ing a plurality of sites where nucleic acid bases are read (the sites exhibit, e.g., colors representing nucleic acid base types); (b) obtain color values of the plurality of sites from the image of the substrate; (c) store the color values in a processor buffer; (d) retrieve partially phase-corrected color values of the plurality of sites for a base calling cycle (the partially phase-corrected color values were stored in the sequencer’s memory during an immediately preceding base calling cycle); (e) determine a prephasing correction; and ((f) determine corrected color values from, e.g., (i) the color values in the sor buffer, (ii) the partially phase corrected values stored during the immediately preceding cycle, and (iii) the prephasing correction.

The ctions or other uration for determining a prephasing correction may include configuration for determining the sing correction from (i) the partially phasecorrected color values stored during the immediately ing base calling cycle, and (ii) the color values stored in the processor buffer.

In certain embodiments, the memory is divided into a plurality of tile buffers, each designated for storing data enting a single image of a tile on the substrate. In certain embodiments, the memory has a storage capacity of less than about 550 Gigabytes (in some examples, this is less than twice the capacity required to store the data contained in the total number of tiles on two flow cells).

The processors may be configured to perform the recited operations in various ways such as receiving executable machine readable instructions. In some cases, the processors are programmed with firmware or custom processing cores such as digital signal processing cores. In various embodiments, the processor(s) are designed or configured to m (and/or control) any or more of the method operations described above.

In some implementations, phasing correction features disclosed herein substantially reduce the cost of a sequencing instrument by more efficiently utilizing memory (e.g., random access memory (RAM)). Some embodiments employ these phasing correction features in the context of real time analysis (RTA) on sequencing platforms.

] According to one aspect of the t disclosure, there is provided a nucleic acid sequencer comprising: an image acquisition system; memory; and one or more processors designed or configured to perform a first base calling cycle by: (a) obtaining data representing an image of a substrate comprising a plurality of sites where nucleic acid bases are read by the image acquisition system, n the sites exhibit colors representing nucleic acid base types; (b) obtaining first color values of the plurality of sites from the image of the substrate; (c) storing the first color values in a processor buffer; (d) retrieving second color values of the ity of sites, wherein the second color values were stored in the memory during a second base calling cycle immediately preceding the first base g cycle; (e) retrieving third color values of the plurality of sites, wherein the third color values were stored in the memory during a third base calling cycle immediately preceding the second base calling cycle; and (f) determining corrected color values for the first base calling cycle from the first color values in the processor buffer, the second color values, and the third color values. [0017B] According to one aspect of the present disclosure, there is ed A method of determining corrected color values from image data acquired, during a base calling cycle, by a nucleic acid cer comprising an image acquisition system, one or more processors, and memory, the method comprising: (a) obtaining data representing an image of a substrate comprising a plurality of sites where nucleic acid bases are read by the image acquisition system, n the sites exhibit colors enting c acid base types; (b) obtaining first color values of the plurality of sites from the image of the substrate; (c) storing the first color values in a processor buffer; (d) ving second color values of the plurality of sites, wherein the second color values were stored in the memory during a second base calling cycle immediately preceding the first base calling cycle; (e) retrieving third color values of the plurality of sites, wherein the third color values were stored in the memory during a third base g cycle immediately preceding the second base calling cycle; and (f) determining corrected color values for the first base calling cycle from the first color values in the processor buffer, the second color values, and the third color values.

These and other features of the disclosure will be presented in greater detail below, with reference to the ated drawings.

BRIEF DESCRIPTION OF THE DRAWINGS Figure 1 is a block m of a cer with hardware for real time analysis of image data taken from nucleic acid clusters.

Figure 2 is an illustration of two channel sequencing data used to illustrate the concepts of phasing and pre-phasing.

Figure 3 depicts a flow cell architecture including a plurality of tiles, each containing many clusters.

Figure 4 depicts a data array containing magnitude data for clusters in a tile or other imaged portion of a flow cell; the magnitude data may be light intensity values for each of two or more color channels.

Figure 5 schematically depicts a first processing configuration and methodology for conducting g correction in real time.

Figure 6 presents a flowchart of a base calling process that may employ the sor and memory configuration depicted in Figure 5.

Figure 7 tically depicts a second processing configuration and methodology for conducting phasing correction in real time. This configuration reduces the requirements on system memory.

Figure 8 schematically depicts a third processing configuration and ology for conducting phasing correction in real time. This configuration r reduces the requirements on system memory.

Figure 9 presents a high-level flowchart of the first few processing cycles that may be employed with the processor and memory configuration of Figure 8 and, in some implementations, Figure 7.

Figure 10 presents a flow chart of processing cycles that conduct fully phasing corrected base calling. Such cycle may be performed in the third and subsequent processing cycles when sequencing clusters of a tile.

Figure 11 presents comparative data for phasing tion methods, one using a reduced main memory thm.

DETAILED DESCRIPTION DEFINITIONS Numeric ranges are inclusive of the numbers defining the range. It is intended that every maximum numerical tion given throughout this specification includes every lower numerical limitation, as if such lower cal limitations were sly written herein. Every m numerical limitation given throughout this specification will include every higher numerical limitation, as if such higher numerical limitations were expressly written herein. Every numerical range given throughout this specification will include every narrower numerical range that falls within such broader numerical range, as if such narrower cal ranges were all expressly written herein.

The headings provided herein are not intended to limit the disclosure.

Unless defined otherwise herein, all technical and scientific terms used herein have the same meaning as commonly understood by one of ordinary skill in the art. Various scientific naries that include the terms included herein are well known and available to those in the art. gh any methods and materials similar or equivalent to those described herein find use in the practice or testing of the embodiments disclosed herein, some methods and materials are described.

The terms defined immediately below are more fully described by nce to the specification as a whole. It is to be understood that this sure is not limited to the particular methodology, protocols, and reagents described, as these may vary, ing upon the context they are used by those of skill in the art.

As used herein, the singular terms “a,” “an,” and “the” include the plural nce unless the context y indicates otherwise. The term “plurality” refers to more than one element. For example, the term is used herein in reference to a number of reads to produce phased island using the methods disclosed herein.

The term “portion” is used herein in reference to the amount of sequence information of genome, chromosome, or haplotype in a ical sample that in sum amount to less than the sequence information of one complete genome, one complete chromosome, or one complete haplotype, as apparent from context.

The term “sample” herein refers to a sample, typically d from a ical fluid, cell, tissue, organ, or organism containing a nucleic acid or a mixture of nucleic acids containing at least one nucleic acid sequence that is to be sequenced. Such samples include, but are not limited to sputum/oral fluid, ic fluid, cerebrospinal fluid, blood, a blood fraction (e.g., serum or plasma), fine needle biopsy samples (e.g., surgical biopsy, fine needle biopsy, etc.), urine, saliva, semen, sweat, tears, peritoneal fluid, l fluid, lavage fluid tissue explant, organ culture and any other tissue or cell preparation, or fraction or derivative thereof or isolated therefrom.

Although the sample is often taken from a human subject (e.g., patient), samples can be taken from any organism having chromosomes, including, but not limited to dogs, cats, horses, goats, sheep, , pigs, etc. The sample may be used directly as obtained from the biological source or following a pretreatment to modify the character of the sample. For example, such pretreatment may include preparing plasma from blood, diluting viscous fluids and so forth. Methods of pretreatment may also involve, but are not limited to, filtration, precipitation, dilution, distillation, mixing, fugation, freezing, lyophilization, concentration, amplification, nucleic acid fragmentation, inactivation of interfering components, the addition of ts, lysing, etc. If such s of pretreatment are employed with respect to the sample, such pretreatment methods are typically such that the nucleic acid(s) of interest remain in the test sample, sometimes at a concentration proportional to that in an untreated test sample (e.g., namely, a sample that is not subjected to any such pretreatment method(s)). Such ed” or “processed” samples are still ered to be biological “test” samples with respect to the methods described herein.

The terms “polynucleotide,” ic acid” and “nucleic acid molecules” are used interchangeably and refer to a covalently linked sequence of nucleotides (i.e., ribonucleotides for RNA and deoxyribonucleotides for DNA) in which the 3’ position of the pentose of one nucleotide is joined by a phosphodiester group to the 5’ position of the pentose of the next.

The nucleotides include sequences of any form of nucleic acid, including, but not limited to RNA and DNA molecules. The term “polynucleotide” es, without limitation, single - and double-stranded polynucleotide.

Single ed polynucleotide molecules can have originated in single-stranded form, as DNA or RNA or have ated in double-stranded DNA (dsDNA) form (e.g. genomic DNA segments, PCR and amplification products and the like). Thus a single stranded polynucleotide may be the sense or antisense strand of a polynucleotide duplex. s of preparation of single stranded polynucleotide molecules suitable for use in the described methods using rd techniques are well known in the art. The precise ce of the primary cleotide molecules is generally not material to the disclosed embodiments and may be known or unknown. The single stranded polynucleotide molecules can ent genomic DNA molecules (e.g., human genomic DNA) including both intron and exon sequences (coding sequence), as well as ding regulatory sequences such as promoter and enhancer sequences.

The nucleic acid described herein can be of any length suitable for use in the provided methods. For example, the target nucleic acids can be at least 10, at least 20, at least , at least 40, at least 50, at least 75, at least 100, at least 150, at least 200, at least 250, at least 500, or at least 1000 kb in length or .

In the context of a flow cell or other substrate for sequencing, the term “site” refers to small region where sequencing takes place. In many embodiments, a site contains multiple, typically numerous, copies of a single nucleic acid sequence from which sequencing data is obtained. The sequence data obtained from a site may be a “read.” The term “polymorphism” or ic polymorphism” is used herein in reference to the occurrence in the same population of two or more alleles at one genetic locus. Various forms of polymorphism include single nucleotide polymorphisms, tandem repeats, microdeletions , ions, indels, and other polymorphisms.

A “base call” is an assigned base (nucleotide type) to sequence data for a particular on in a polynucleotide sequence. A base call may be output by a sequencer for each position in nucleic acid being sequenced. A y of the call is sometimes ascribed to a base call.

The term “read” refers to a sequence read from a portion of a nucleic acid sample.

Typically, though not necessarily, a read represents a short sequence of contiguous base pairs in the . The read may be represented ically by the base pair sequence (in ATCG) of the sample portion. It may be stored in a memory device and processed as appropriate to determine whether it s a reference sequence or meets other criteria. A read may be obtained directly from a cing apparatus or indirectly from stored sequence information concerning the sample. In some cases, a read is a DNA sequence of sufficient length (e.g., at least about 25 bp) that can be used to identify a larger ce or region, e.g., that can be aligned and specifically assigned to a chromosome or genomic region or gene.

The term “Next Generation Sequencing (NGS)” herein refers to sequencing methods that allow for massively parallel sequencing of clonally amplified molecules and of single nucleic acid molecules. Non-limiting examples of NGS include sequencing-bysynthesis using reversible dye ators, and sequencing-by-ligation.

The term “parameter” herein refers to a numerical value that characterizes a physical property or a representation of that property. In some situations, a parameter numerically terizes a quantitative data set and/or a numerical relationship between tative data sets. For example, the mean and variance of a standard distribution fit to a histogram are parameters.

The terms “threshold” herein refer to any number that is used as a cutoff to characterize a sample, a c acid, or portion thereof (e.g., a read). The threshold may be compared to a measured or calculated value to determine whether the source giving rise to such value suggests should be classified in a particular manner. Threshold values can be identified empirically or analytically. The choice of a threshold is dependent on the level of confidence that the user wishes to have to make the classification. Sometimes they are chosen for a particular purpose (e.g., to balance sensitivity and selectivity).

Real time analysis refers to a process and system in which sing and data analysis are performed in the background of data acquisition during a DNA sequencing run.

An example of a real time analysis system is bed in US Patent No. 8,965,076, which is incorporated herein by reference in its entirety.

CONTEXT FOR PHASING Sequencing apparatus Figure 1 shows a block diagram of some features of a typical nucleic acid sequencer 100 or a system including such sequencer. Notably, the system 100 includes a flow cell 101, and image acquisition system 103, one or more processors 105 with one or more s 107, and system memory (sometimes referred to as main memory) 109 including a plurality of tile buffers 111. Typically, system memory 109 is provided on device that is not part of an integrated circuit ning any of the one or more processor(s) 105. In certain embodiments, the system memory is volatile memory such as Random Access Memory or RAM, e.g., DRAM, a solid state hard drive, or a hard disk drive.

The flow cell and image ition system contain components designed or ured in ance with principles tood in the field of c acid sequencing, and they will not be described in detail herein. le image analysis systems and associated flow cells are employed in nucleic acid sequencers such as the MiSeq and HiSeq series of sequencers available from Illumina, Inc. of San Diego, California. For additional information, see US patent number 573, US patent number 9,193,996, and US patent number 8,951,781, each of which is incorporated herein by reference in its entirety.

In general, nucleic acid sequences le for use with the disclosed methods provide rapid and efficient detection of a plurality of target nucleic acid in parallel. They can include fluidic components capable of delivering amplification reagents and/or sequencing reagents to one or more immobilized DNA fragments, the system including components such as pumps, valves, reservoirs, fluidic lines and the like. A flow cell can be configured and/or used in an integrated system for detection of target nucleic acids. Exemplary flow cells are described, for example, in US 2010/0111768 A1 and US Ser. No. 13/273,666, each of which is incorporated herein by nce in its entirety. As exemplified for flow cells, one or more of the c components of an integrated system can be used for both an amplification method and for a detection method. For example, one or more of the fluidic components of an integrated system can be used for an amplification method and for the ry of sequencing reagents in a sequencing method. Alternatively, an integrated system can include separate fluidic systems to carry out amplification methods and to carry out detection methods.

For purposes of this disclosure, it is sufficient to understand that the flow cell first receives and immobilizes or otherwise es a nucleic acid sample which is to be sequenced and then exposed to various reagents associated with the cing process. In certain embodiments, the sequencing process is a sequence by synthesis process, although other sequencing technologies may be ed.

The image acquisition system 103 includes l components such as scence tion components (e.g., a laser and associated mirrors and lenses) for illuminating sites on the flow cell where sequencing is taking place and image capture components for capturing images of fluorescence on portions of the flow cell having multiple sites. The data captured by the image acquisition system contains information le for determining which nucleotide is being read on any given site at any given sequencing cycle.

To allow for real-time analysis, the sequencer 100 typically includes onboard processors and memory that interpret and store image data from the image ition system 103. Examples of suitable processors for the sequencer include s Xeon E5 class.

Typically, the processor 105 includes multiple buffers 107 that temporarily store image data taken during a single image acquisition cycle. In the ed embodiment, the processor buffers are allocated in the system . A given processor buffer may be associated with a particular processor thread created to analyze image data of a region of the flow cell during real time analysis. In certain embodiments, the image data analyzed by a thread is that of a single tile (described below), captured during a single image acquisition cycle. In certain embodiments, the buffer can store about 400 Gigabytes of data. As used herein, a thread is an ordered sequence of instructions that tells the processor what operations to e. The instructions configure the processor using executable machine code selected from a specific machine language instruction set, or “native instructions,” designed into the hardware processor.

The machine ge instruction set, or native instruction set, is known to, and essentially built into, the hardware processor(s), or CPUs. This is the “language” by which the system and application software icates with the hardware processors. Each native instruction is a discrete code that is recognized by the processing architecture and that can specify particular registers for arithmetic, addressing, or control functions; particular memory locations or offsets; and particular addressing modes used to interpret ds. More complex operations are built up by combining these simple native instructions, which are executed sequentially, or as otherwise directed by control flow instructions.

System memory 109 includes le tile buffers 111, each configured to store a portion of the image data acquired from the flow cell during a single image acquisition cycle.

Tile buffers in this example are referred to as such e they are configured to hold a single tile's worth of image data. As explained more fully below, a tile is a region of a flow cell that can be captured in a single image taken during a single image acquisition cycle. Tile buffers 111 are intended to store image data over a longer period of time than processor buffers 107. In n embodiments, tile buffers 111 store image data for at least two image acquisition cycles. While this application describes buffers that buffer data from a tile of a flow cell, the disclosed embodiments are not limited to buffers storing this amount of data.

Unless otherwise stated or clear from context, references to “tile buffers” are understood to include any type of buffer that stores image data from a portion of a flow cell, which image data is processed as a unit as described herein.

To make base calls, the one or more processors 105 acts on data provided from system memory 109 and data stored in processor s 107. Typically, a single base call is made for a single site during a single image acquisition cycle.

As shown, the one or more processors 105 and the main memory 109 share data ctionally.

Additionally, the one or more processors 105 e image data from image ition system 103. In certain embodiments, image acquisition system 103 obtains data from flow cell 101 by exciting the sequencing sites on flow cell 101 and receiving optical s from those sites. In certain embodiments, the signal received by image acquisition system 103 is a fluorescence signal created when system 103 illuminates flow cell 101 with light at appropriate wavelengths. In such embodiments, the scence signal is provided as intensity values for a plurality of colors.

The concept of a cycle is used throughout this sure. A single sequencing cycle involves reading a single nucleotide from each of one or more sites captured on an image.

The reading is referred to as making a base call. In various ments described herein, a single computational cycle--from the perspective of the processor(s) and memory--performs both base calling and image capture but for different nucleotides, with the base calling g image capture in the sequence of nucleotides being read or called. For example, in a single computational cycle, the one or more processors t base calling for a nucleotide in sequencing cycle n and concurrently conduct image capture for nucleotide in sequencing cycle n +1. Thus, in a single computational cycle, the sequencer (a) stores and processes unmodified image data for nucleotides in sequencing cycle n +1 and (b) makes a base call for nucleotides in sequencing cycle n. The use of the processor buffers and tile buffers in this cycle-by-cycle processing will be described in more detail below.

Phasing Generally At a particular site on a flow cell or other substrate, multiple copies of a c acid molecule, all having the same sequence (possibly with d variations unintentionally introduced by sample processing), are analyzed together. Enough copies are used to ensure that sufficient signal is produced to permit reliable base calling. The collection of nucleic acid molecules at a site is called a cluster. In some cases, an unsequenced cluster contains only single stranded nucleic acid molecules.

Phasing represents an unintended artifact that arises from sequencing le nucleic acid molecules within a cluster. g is the rate at which signals such as fluorescence from single molecules within a r lose sync with each other. Often the term g is reserved for contaminating signal from some molecules that fall behind, and the term pre-phasing is used for contaminating signal from other molecules that go ahead.

Together phasing and pre-phasing be how well the sequencing apparatus and chemistry is performing.

Low numbers are better. Values of .10 mean 0.10% of the molecules in a cluster are both falling behind and 0.10% are running ahead at each base calling cycle. In other words 0.20% of the true signal is lost each cycle and will therefore contribute to noise. r example, 0.20/0.20 means that 0.4% of the true signal is lost per cycle, in which case after 250 cycles (without correction) the noise would be equal to the signal.

A real time analysis component of a sequencer may determine g and prephasing in order to apply the correct level of phasing correction as sequencing proceeds. This works by artificially pushing signal in or out of each sequencer channel based on base calls before or after the t cycle. usly, phasing and pre-phasing were estimated over a defined number of cycles (e.g., the first 12 cycles of each read) and then applied to all subsequent cycles. Some recent sequencers employ an algorithm called cal phasing correction to optimize the phasing correction at every cycle by trying a range of corrections and selecting the one which results in the highest chastity (signal purity). While empirical phasing correction provides improved performance, it requires greater computational resources.

In conventional sequencers, each base has a unique fluorescent dye color; e.g., green to thymine, red for ne, blue for guanine, and yellow for e. To capture information for base calling, a four channel sequencer takes four images of a tile or other portion of a flow cell. Some sequencers now have only two channels, and therefore take only two images of the same portion of the flow cell. A two-channel cer uses a mix of dyes for each base and uses red and green filters for the two images. In an example of a two channel sequencer, clusters seen in red or green images are interpreted as C and T bases, respectively. Clusters observed in in both red and green images are d as A bases, while unlabeled clusters are fied as G bases.

Figure 2 illustrates phasing during sequencing of a nucleic r having the sequence . . . ACGTAAG . . . . As rated, the during the base calling cycle for the first G, 98.4% of the fluorescence signal originates from sequences currently generating signal for G, while 1.5% of the fluorescence signal originates from ces currently producing signal for the prior base C, and 1.1% of the fluorescence signal ates from sequences currently producing signal for the next base T. The signal contribution for the prior base C is from phasing and the signal contribution from the next base T is from pre-phasing.

Phasing correction for this G base call is reflected in the graph on the right side of Figure 2. As shown for a two-channel sequencer, the fluorescence signal can be represented on a two-dimensional plot, with maximal intensity signal on a “green axis” representing T, maximal intensity on a “red axis” representing C, maximal intensity mid-way between the axes representing A, and minimal intensity on both axes representing G. Without phasing error, the signal for G should have zero intensity on both the red and green axes. Instead, with the phasing error discussed, the fluorescence signal has some intensity contribution on both the green and red axes. In this example, pre-phasing correction reduces the signal intensity to zero on the green axis and phasing correction reduces the signal intensity to zero on the red axis. Similar corrections may be made on base calls for the bases T, C, and A.

Tiles and Flow Cells As explained, a flow cell contains multiple sites where cing information is collected. In certain embodiments, each site of a flow cell ns a cluster of singlestranded nucleic acids sharing the same ce. A single image used in real time sequencing may contain ns of such clusters. A typical flow cell is so large that it requires hundreds or even thousands of separate images to cover its entire area. In certain embodiments, the processor and associated memory employed for real-time analysis processes all these images currently to make base calls for a single cycle. In some implementations, the sor and memory concurrently process all images ed over two or more flow cells during a single base calling cycle. Figure 3 schematically depicts a flow cell architecture used in some sequencers from na, Inc. In the depicted example, the sequencer makes concurrent base calls on two flow cells, Flow Cell 1 and Flow Cell 2. In certain embodiments, each flow cell has sequencing sites on each of two surfaces, a top surface in the bottom e. In such cases, the sequencer images both the top and bottom surfaces during each base calling cycle. As depicted in Figure 3, each flow cell surface includes four lanes, L1, L2, L3, and L4; of course other numbers are possible. Each lane of each surface may have multiple subdivisions referred to as swaths. Each swath is in turn divided into multiple tiles. For example, there may be approximately 120 tiles per swath.

Considering two flow cells, each having two surfaces, with each surface having four lanes, each lane having six , and each swath having 120 tiles, several thousand tiles of data need to be analyzed per cycle. In s embodiments, each tile image (or other image from a portion of a flow cell) is acted on by a single processor thread. In certain embodiments, a sequencer employing a flow cell having the architecture depicted in Figure 3 processes 8000 or more tiles of data in each base calling cycle. In such cases, the real time processing logic would employ 8000 or more processor threads in each base calling cycle The data from a single tile captured during a single cycle can be stored in the memory as an array, with each entry in the array representing a color value for each channel of a single r in the tile. An array for a two-channel arrangement is depicted in Figure 4.

As an e, a color intensity or can generate signal counts between about 400 and 1500 for each channel. A tile buffer in the system memory is configured to store all the information in the array, in other words the color values of all clusters on a tile at a single base calling cycle. A sor buffer may be similarly configured to store all the information in the array.

Phasing Process A significant memory burden of real time analysis of sequence data stems from the requirement in phasing correction that two or three cycles of cluster intensities must be saved for every tile for the full length of the run. On an Illumina HiSeqX with a 700nm flowcell, this takes up 73 Gigabytes of memory. This burden is sufficiently large that most of the data (on this platform) is cached to a solid state hard drive.

As explained, phasing correction adjusts the ity values of an image to address out of phase sequencing of some nucleic acid stands in a r. Phasing correction accomplishes this by starting with the measured cluster color intensity values (or other signals ed by with the sequencing ) for a current base calling cycle and adding or cting a correction value using measured intensity values from the previous base calling cycle and/or using measured intensity values from the subsequent base calling cycle.

In various implementations, a phasing corrected intensity value for making a base call applies an expression as shown in the bottom of Figure 5. As shown there, phasing corrected intensity values for a current base calling cycle in an image equal the measured intensity values for the current base calling cycle minus the product of a first coefficient and the measured intensity values at the immediately previous base calling cycle and minus the product of a second coefficient and measured intensity values at the immediately successive base calling cycle: Corrected Intensity = -a.In-1 + In - b.In+1 where In-1, In, and In+1 are the intensity values of clusters in a tile at the immediately preceding base calling cycle, at the current base calling cycle, and the immediately ding base calling cycle respectively. The coefficients a and b are the g and prephasing coefficients (sometimes called weights), respectively. These may be calculated anew for each base calling cycle of a tile.

Returning to Figure 2, the measured intensity value for the third base in the depicted sequence (for a single cluster in an image) is shown as dot in the graph on the right side of Figure 2. The pre -phasing correction to this measured intensity value is reflected by the vertical arrow from the measured intensity value down to the horizontal axis. In the expression for phasing corrected intensity values, this pre-phasing correction is represented by the t of the cient b and the intensity value measured for the next successive base calling cycle. In addition, the measured intensity value is corrected by a phasing correction represented by the horizontal arrow on the graph. This phasing correction is ented by subtracting from the measured intensity value, the t of a coefficient a and the measured intensity value for the immediately preceding base calling cycle. The coefficients a and b may be determined by numerous methods, but in many implementations, they are calculated fresh for each base calling cycle. A description of methods for determining the coefficients to be used in phasing correction is described in International Patent Application having Publication Number WO2015/084985 by Belitz et al. and published on June 11, 2015, which is incorporated herein by nce in its entirety.

In n embodiments, the g algorithm determines g coefficients empirically by maximizing the cumulative chastity (or similar metric) of the cluster ity data during a base g cycle. One implementation of the thm iterates over all or many phasing coefficients and determines which ones give the best results. For example, the phasing algorithm may optimize a and b at every cycle using a pattern search employing a cost function that counts the number of clusters that fail a chastity . Thus, a and b are ed to maximize the data quality.

In some embodiments, phasing coefficients are determined as an on-going analysis throughout a sequencing run (e.g., during generate of a read). As a result of this approach, an inaccurate phasing estimation made during early cycles will not adversely affect later cycles.

Some methods determine chastity of a cluster intensity value as a function of relative distances to Gaussian centroids for the other cluster intensity values determined for the same base calling cycle. The centroids ideally align with expected locations of the A, T, C, and G intensities for two channels (see Figure 2), assuming that a two-channel system is used. In certain ments, chastity can be calculated using the expression: chastity= 1 - Dl/(Dl + D2), where D1 is the distance to the nearest Gaussian centroid, and D2 is the distance to the next nearest centroid. Utilizing this approach, when the mean chastity (quality) of intensity values are maximized, the correct values of a and b are chosen,. Once these values are identified, then a correction can be applied to all cluster values and base calling can occur directly. s of g Gaussian distributions to a two-channel data set are described in International Patent Application having ation Number /084985, previously incorporated by reference.

In some embodiments, a phasing correction is calculated at nearly every cycle during a sequencing run. In some embodiments, a phasing correction is calculated at every cycle during a sequencing run. In some embodiments, a separate phasing correction is calculated for different locations of an imaged surface at the same cycle. For example, in some embodiments, a separate phasing correction is calculated for every individual lane of an imaged surface, such as an individual flow cell lane. In some embodiments a te phasing correction is calculated for every subset of a lane, such as an imaging swath within a flow cell lane. In some embodiments, a separate phasing correction is calculated for each individual image, such as, for e, every tile. In certain embodiments, a separate phasing correction is ated for every tile at every cycle.

As reads get , higher order terms can become more important in g correction. Thus, in particular ments, to correct for this, a second order empirical phasing correction can be calculated. For example, in some embodiments, the method comprises a second order phasing correction as d by the following: I( cycle)= -a*I(cycle-2)- A *I( I)+ I(cycle)- B*I(cycle+ 1 )-b*I(cycle+2) where I represents intensity and a, A, B, and b represent the first and second order terms to the phasing correction. In particular embodiments, the calculation is optimized over a, A, B, and b.

Figure 5 schematically depicts a processing configuration and methodology for conducting phasing correction in real time. In the depicted embodiment, a processor 502 creates a new processing thread 503 when the processor is called upon to make base calls from clusters in an image, e.g., an image of a tile. A new thread may be generated for each base calling cycle for each tile. In the depicted embodiment, the processor 502 makes available a single processor buffer 505 for each base calling cycle of a tile (and the designated processing thread). The processor buffer temporarily stores intensity values that are computationally manipulated by the processor to conduct phasing correction for a current base calling cycle n. In the depicted embodiment, the sor interfaces with a system memory 507 containing three buffers, one each for g image data captured for a ular base calling cycle. In the case of the flow cell architecture depicted in Figure 3, each buffer stores image data for the clusters of a single tile; hence the buffers are ed to as tile s. Of course, for other flow cell architectures and/or image acquisition systems, the buffers may store more or less cluster data. For convenience, the specification will refer to tile buffers. Each tile buffer stores data for a single tile (or other portion of a flow cell) captured during a single base calling cycle. The image data may be provided as an array of data such as shown in Figure 4.

As ed, system memory 507 includes a tile buffer 509 which temporarily stores intensity values for the immediately previous base calling cycle (in comparison to the current base calling cycle handled by the processor), a tile buffer 511 which stores intensity values measured for the current base calling cycle, and a tile buffer 513 which stores intensity values for the immediately succeeding base calling cycle. Again, each of the tile buffers 509, 511, and 513 contain measured data of a single tile for a single base calling cycle n.

As shown, thread 503 makes use of the intensity values in each of the tile buffers 509, 511, and 513 during a single base calling cycle. The intensity values are sively loaded into processor buffer 505 and manipulated to implement the phasing correction expression presented at the bottom of Figure 5. After the base calling process is completed as depicted in the processor and memory configuration of Figure 5, the processor buffer holds adjusted ity values used to make a phasing corrected base call.

Figure 6 presents a flowchart of a base calling process that may employ the processor and memory configuration depicted in Figure 5. As shown in Figure 6, a process 601 initiates a new base calling cycle by creating a sor thread and allocating a processor buffer to that thread. See process block 603. Thereafter, the processor extracts intensity data from an image of a flow cell tile (or other appropriate portion of the flow cell) taken rently with the current processing cycle. In the ed implementation, the ed image and associated intensity values are the primary ity values for the next successive base calling cycle, not the current base calling cycle (the current processing iteration). In other words, the current processing cycle performs a base call for image data collected in an immediately preceding processing cycle. Thus, as depicted in a process block 605 of process 601, the ted intensity values are given the reference In +1, where n ents the current base calling cycle. Stated another way, a processing cycle both (i) calls bases for base calling cycle n, and (ii) captures image data for base calling cycle n+1.

The newly extracted intensity data, which may be provided in the form of an array as depicted in Figure 4, is stored in an available tile buffer on the system memory (e.g., tile buffer 513). In certain embodiments, this tile buffer is one that stored intensity data that was previously used but is no longer necessary for base calling.

In the current processing cycle, process 601 also retrieves intensity data stored during a computational cycle previous to the current computational cycle. See process block 607. The retrieved intensity data is for the current base calling cycle and is given reference In.

The retrieved intensity data is ed from an appropriate tile buffer such as tile buffer 511 of the system memory as shown in Figure 5.

In on, process 601 retrieves intensity data that was stored two cycles previous to the current base calling cycle. See process block 609. As an example, with reference to Figure 5, such ity data may be ed from a tile buffer 509 of the system memory.

The array of intensity values retrieved in operation 609 is identified by In -1.

While operations 605, 607, and 609 are shown as ing sequentially, this order of operations is flexible and the process can be implemented such that any order is acceptable, so long as it is consistent with base calling that incorporates phasing correction.

Upon retrieving the intensity values for the t base calling cycle (process block 607) and the intensity values for the immediately preceding base calling cycle (processing block 609), the processor has available all intensity values it needs to perform a phasing correction. It does this by first determining the pre-phasing correction weight b and the phasing correction weight a for the current base calling cycle. See process block 611, which illustrates that this may be accomplished using the extracted intensity values for the next succeeding base calling cycle along with the intensity values for the current and ately preceding base calling cycles. Then, using the phasing and asing tion weights, the processor ates phasing corrected intensity values for the current base calling cycle as depicted in process block 613. The corrected values are for the clusters in the tile under consideration. The calculation may employ the expression depicted in block 613. Using the phasing corrected intensity values, the sor makes calls for the current base calling cycle as depicted in process block 615.

At this point, the processing for the current base calling cycle is complete and the next iteration of base calling may be executed. The decision of whether to conduct another base g cycle is depicted in a block 617 which determines whether there are any further nucleotides to be sequenced in the clusters of the tile under consideration. If there are none, the s is completed as depicted at block 619. If there are, process control is handed to a process block 621 where the processor increments a cycle count. This effectively indexes the intensity values for the current base calling cycle In to intensity values for the immediately ing base calling cycle In -1. At the same time, the intensity values for the immediately next base calling cycle (In+1) become the intensity values for the new current base calling cycle (In). These increments are made with respect to the indexes applied to the intensity data stored in the tile buffers.

Phasing Process (Reduced Main ) The approach of Figures 5 and 6 can work fine so long as the sequencer and its associated real-time analysis system is not memory constrained. r, given the amount of data that must be processed in certain modern sequencers, such as those employed to m whole genome sequencing, insufficient memory may be ble, particularly at a commercially viable cost. Therefore, storing three times the amount of data required to fully image the flow cell (or flow cells) during a base calling cycle can present a serious neck.

A phasing algorithm such as represented in Figures 5 and 6 is an important contribution to real time analysis, in that it significantly improves sequencing results, particularly on non-standard samples, e.g. low diversity samples. However, the imposed memory burden becomes greater as the throughput of next tion sequencing s grows. The following embodiments reduce memory burden by using phasing weights learned from data that was already partially phasing corrected. The phasing and pre-phasing weights can be learned independently and still e high quality sequencing results. In some es, the main memory requirement is less than twice the capacity required to store the data contained in the total number of tiles on two flow cells.

In certain embodiments, the processor and memory configuration for phasing corrected base calling is adjusted to reduce the requirements on system memory. One example of how this works is ed in Figure 7. Intensity values are corrected as described above, e.g., phasing and pre-phasing weights are calculated and applied to the immediately preceding and immediately ding . However, in the example of Figure 7, system memory 707 employs only two tile buffers for phasing correction: tile buffer 709 and tile buffer 711. In this example, a sor 702 employs a processing thread 703 which, contrary to the example of Figure 5, has two associated sor buffers: a processor buffer 705 for storing and operating on the intensity values retrieved from memory 707 and a processor buffer 706 for storing and using the newly captured image intensity values In+1. In the depicted example, the processor buffers are ted in main memory, but this is not always required. In some embodiments, the sor buffers are allocated in a different physical memory or even on the processor chip.

Replacing tile buffers with processor buffers effectively reduces the total memory ements. By using multiple sors and/or multithreaded processing, a few processors handle many tiles. As an example, the number of tiles in a system may be on the order of 1000-2000, while the number of processors handling all these tiles is about twenty.

In , such system can realize a memory reduction on the order of 50x. In some implementations, the reduction is on the order of 20x.

In this implementation, the intensity values captured from tile images in the current processing cycle (In+1) are stored locally on the processor and used to calculate the phasing and pre-phasing weights and subsequently make a base call. In some implementations, only after this process is complete are the most recently captured intensity values (In+1) stored in a tile buffer on system memory 707.

In some embodiments, a processor and system memory are configured as depicted in Figure 8. As with the process or/memory configuration in Figure 7, a processor 802 employs sing threads 803, each associated with two processor s: a processor buffer 805 for temporarily storing intensity values from a system memory 807 (tile buffer 811), and a processor buffer 806 for temporarily storing intensity values captured during the current processing cycle (In+1). In order to allow this configuration to work efficiently and effectively, the intensity values stored in tile buffer 811 must be partially phasing corrected.

Examples of mechanisms for accomplishing this are described below. Processor buffer 705 in Figure 7 and processor buffer 805 in Figure 8 load intensities from main memory and then manipulate those intensities to generate the corrected intensities which are ed for base calling. In the depicted example, the processor buffers are allocated in main memory, but this is not always required. In some embodiments, the processor buffers are allocated in a different al memory or even on the processor chip.

Figure 9 presents a high-level view of a s 901 that may be employed with the processor and memory configuration of Figure 8 and, in some implementations, Figure 7. As illustrated in Figure 9, the first and second processing cycles employ insufficient information to t full phasing correction on clusters imaged in a tile. However, phasing is not a significant problem in the very first cycles.

To t full phasing correction, the cer requires three consecutive cycles of image data. In the first processing cycle, the sequencer does not make a base call; it merely stores intensity data for the next processing, i.e., the cycle in which the first base call is made.

As depicted, the process 901 begins at a process block 903 where a thread is d for the first processing cycle. The instructions in this thread direct extraction of intensity data from an image of the clusters during the first sequencing cycle (I1), i.e. the cycle during which the first nucleotides of the clusters are read. See process block 905. The image data is stored in a tile buffer in system memory. At this point, the first sing cycle is effectively complete.

The process continues at a process block 907 where a new thread is d in preparation for the second processing cycle. In this process, first and second processor buffers are allocated for the second processing cycle. See block 907. Collectively, process blocks 907, 909, 911, 913, 915, 917, 919, 921, and 923 are performed during the second processing cycle, which executes using the thread and processor buffers generated at process block 907.

As depicted, the processor extracts intensity data from the image for the next base calling cycle (I2) and stores that data in a first processor buffer. See process block 909. Next, during the second processing cycle, the processor retrieves the intensity data stored in the tile buffer during the first processing cycle, which intensity data is for the current base calling cycle (I1). See block 911. Using the intensity data collected during the first and second processing cycles, the processor can calculate a pre-phasing weight b for the current base calling cycle (i.e., the first base calls in the reads). See process block 913. With the intensity values for the first two cycles and the asing weight, the processor calculates corrected intensity data values for the second base calling cycle (I2). The corrected intensity data values may be stored in the second processor buffer. See process block 915. Next, the processor makes the base calls for the second base calling cycle using the corrected intensity data values obtained in block 915. See process block 917.

At this point, the sequencing process is ready to begin preparing for the next base calling cycle. It starts at a process block 919 by determining a phasing tion weight a using the next (or second) base calling cycle intensity data (I2) and the current base calling cycle data (I1), which was stored in the tile buffer. Using the phasing correction weight a, the processor next calculates phasing corrected (but not pre-phasing corrected) intensity data values from the currently ected intensity data (I2) ted during this second processing cycle and the intensity data values for the first processing cycle (I1) according to the expression presented in s block 921. This results in a partially corrected intensity value array (I2(partially corrected)) for the second base calling cycle. The cer will have to await the next processing cycle before conducting pre-phasing correction. However, at this point much of the calculation is completed and the array data for a single image can be stored in a tile buffer for use in the next base calling cycle. To this end, the processor stores the phasing corrected (but not pre-phasing corrected) ity data in the tile buffer (such that I2(partially corrected) es I1 in the tile buffer). See process block 923.

At this point, the first and second processing cycles are completed and base calls are made for the first base calling cycle, which is the second processing cycle. Subsequent base calling cycles may be performed with full phasing correction as described in Figure 10. See process block 925.

Figure 10 s a sequence of operations it may be performed during a processing cycle that conducts fully g corrected base calling. Such cycle may be performed in the third and subsequent processing cycles when sequencing clusters of a tile. In certain embodiments, the sequence of operations depicted in Figure 10 ponds to process block 925 of Figure 9.

As depicted, the process begins by allocating a thread and associated first and second processor buffers. See process block 1003. Next, the processor extracts intensity data values from an image for the next base calling cycle (In+1) and stores t hose values in a first processor buffer. See process block 1005. Concurrently, the processor retrieves the partially corrected intensity data values that were stored during the previous base calling cycle (as a non-limiting example, I2(partially corrected) in the embodiment of or In – a(In-1)). These values now represent the intensity values for the current base calling cycle (In). They were previously stored in the system memory’s tile buffer and are now retrieved therefrom. See process block 1007. With the lly corrected intensity data values for the current base calling cycle, which were phasing corrected, the processor need only conduct pre-phasing correction to complete the correction of the intensity data and make the necessary base calls for the current base calling cycle. To this end, the processor determines the pre-phasing correction weight b for the current base calling cycle. It does this using extracted intensity data that it just retrieved from the image data, for the next cycle , along with the previously partially corrected intensity data for the current base calling cycle. Recall that this partially corrected data that was just retrieved from the tile buffer. The lly corrected intensity data may be represented by the expression In – a(In-1). See s block 1009.

With the pre-phasing correction weight b calculated for the current base calling cycle, the processor has all it needs to calculate a fully phasing corrected intensity data array for the t base calling cycle (In). The calculation is conducted as depicted in process block 1009. The resulting fully corrected intensity data values are stored in the second sor buffer. See process block 1011. fter, the processor makes the base calls for the current base g cycle using the corrected intensity data values stored in the second processor buffer. See process block 1013.

The current processing cycle can begin ing for the next base calling cycle which will be executed during the next processing cycle. In the depicted embodiment, the processor determines the phasing correction weight a for the next base calling cycle using ity data available for the current base calling cycle. See process block 1015. Recall that the next base g cycle intensity data was extracted and stored in the first sor buffer at s operation 1005. Partially corrected intensity values for the current base g cycle were ved from the tile buffer for purposes of making the current base calls. The same partially corrected intensity values are now used to calculate the phasing correction weight a for the next base calling cycle. With the g correction weight for the next base calling cycle now calculated, the processor calculates phasing corrected (but not pre-phasing corrected) intensity data values as depicted in process block 1017. Processor then stores these phasing corrected ity data values for the next base calling cycle in the tile buffer. See process block 1019.

Before this ion, it was assumed that that base calling accuracy would suffer by learning prephasing weights from phasing corrected intensities. However, the results herein show that little or no inaccuracy results. In some implementations, the image data is compressed (e.g., lossy compression) and even the lly phase corrected data is compressed. In both cases, it has been demonstrated that the compression could be performed t loss of cy. As an example, without compression, an implementation uses two float buffers for each tile (a float buffer is 4 bytes in size). With compression, an implementation uses a single byte buffer, thus realizing 4x less memory.

At this point the current processing cycle is effectively complete, so the sor determines whether there are any more cycles that need be conducted in sequencing the clusters of the t tile. See decision block 1021. If no further bases need be read from the clusters, the process is complete and no further processing cycles are ted. However, if one or more additional sequencing cycles are required, process control is directed to a process block 1023 where the processor increments the t cycle at which point the partially corrected intensity data values stored in the tile buffer become current; i.e., they become the values for the new base calling cycle. Process control then returns to process block 1003 where the next processing cycle begins.

EXAMPLE As explained, certain embodiments reduce memory burden by using phasing s learned from data that was already partially g corrected. However, it was not clear that the phasing and pre-phasing weights can be learned independently and still provide high quality sequencing results. The example presented in Figure 11 establishes that they As shown, two comparisons were made, each using a baseline process (e.g., a process of Figures 5 and 6) and a new process that was optimized to reduce main memory requirements (e.g., a process of s 8 and 10). In each comparison, the same sequencer and sample were employed. Specifically, an na HiSeqX ment was converted to use 2 dye try. The sequencer’s output images were saved and the two phasing algorithms were both tested on the same sequencing images, providing a completely controlled test. The “Clusters PF” indicates the throughput delivered by the sequencer; the %Aligned indicates the number of clusters that successfully aligned to the nce , and the r Rate” indicates the mean error rate of the sequences called by the software compared to the reference genome.

The sequencing results trate that the memory-efficient phasing algorithm is comparable to the baseline algorithm. In this example, the memory efficient process produced an approximately 3% se in error rate, which is offset by a reduction in main memory (estimated be from 420 Gigabytes to 340 Gigabytes in some implementations).

SEQUENCING METHODS As indicated above, the disclosure pertains to sequencing nucleic acid samples.

Any of a number of sequencing technologies using one or more channels of information for base calling, particularly l channels, may be used. Particularly applicable techniques are those where nucleic acids are attached at fixed locations in an array (e.g., as cluster) and where the array is repeatedly imaged. Embodiments in which images are obtained in different color channels, for example, coinciding with ent labels used to distinguish one nucleotide base type from another are particularly applicable. In some embodiments, the process to determine the tide sequence of a target nucleic acid can be an automated process. Certain embodiments e sequencing-by-synthesis ("SBS") techniques. While sequencing by synthesis techniques are emphasized here, other sequencing technologies may be employed.

In many implementations, SBS techniques involve the enzymatic extension of a nascent nucleic acid strand through the iterative addition of nucleotides against a template strand. In traditional s of SBS, a single nucleotide monomer may be provided to a target nucleotide in the presence of a polymerase in each ry. However, in the methods described herein, more than one type of nucleotide monomer can be provided to a target nucleic acid in the presence of a polymerase in a delivery.

SBS can utilize tide monomers that have a terminator moiety or those that lack any terminator moieties. Methods utilizing nucleotide monomers lacking terminators include, for example, pyrosequencing and sequencing using γ-phosphate-labeled nucleotides.

In methods using nucleotide monomers lacking terminators, the number of nucleotides added in each cycle is generally variable and dependent upon the template sequence and the mode of nucleotide delivery. For SBS techniques that utilize nucleotide monomers having a terminator moiety, the terminator can be effectively irreversible under the sequencing conditions used as is the case for traditional Sanger sequencing which utilizes dideoxynucleotides, or the terminator can be reversible as is the case for sequencing methods ped by Solexa (now Illumina, Inc.).

SBS techniques can utilize nucleotide monomers that have a label moiety or those that lack a label moiety. Accordingly, incorporation events can be detected based on a characteristic of the label, such as fluorescence of the label; a characteristic of the nucleotide monomer such as molecular weight or charge; a byproduct of incorporation of the nucleotide, such as e of pyrophosphate; or the like. In embodiments, where two or more different nucleotides are present in a cing reagent, the different nucleotides can be distinguishable from each other, or atively, the two or more ent labels can be the indistinguishable under the detection techniques being used. For example, the different nucleotides t in a sequencing reagent can have different labels and they can be distinguished using appropriate optics as exemplified by the sequencing s developed by Solexa (now Illumina, Inc.).

Some ments include pyrosequencing techniques. Pyrosequencing detects the release of inorganic pyrophosphate (PPi) as particular nucleotides are incorporated into the nascent strand (Ronaghi, M., Karamohamed, S., Pettersson, B., Uhlen, M. and Nyren, P. (1996) "Real-time DNA sequencing using ion of osphate release." ical Biochemistry 242(1), 84-9; Ronaghi, M. (2001) "Pyrosequencing sheds light on DNA sequencing." Genome Res. 11(1), 3-11; Ronaghi, M., Uhlen, M. and Nyren, P. (1998) "A sequencing method based on real-time pyrophosphate." Science 281(5375), 363; U.S. Pat.

No. 6,210,891; U.S. Pat. No. 568 and U.S. Pat. No. 6,274,320, the disclosures of which are incorporated herein by reference in their entireties). In pyrosequencing, released PPi can be detected by being immediately converted to adenosine triphosphate (ATP) by ATP sulfurylase, and the level of ATP generated is ed via luciferase-produced photons. The nucleic acids to be sequenced can be attached to features in an array and the array can be imaged to capture the chemiluminescent signals that are produced due to incorporation of tides at the features of the array. An image can be obtained after the array is treated with a ular nucleotide type (e.g. A, T, C or G). Images obtained after on of each nucleotide type will differ with regard to which features in the array are detected. These differences in the image t the different sequence content of the features on the array.

However, the relative locations of each feature will remain unchanged in the images. The images can be stored, processed and analyzed using the methods set forth herein. For example, images obtained after treatment of the array with each different nucleotide type can be handled in the same way as exemplified herein for images obtained from different detection channels for ible terminator-based sequencing methods.] In another exemplary type of SBS, cycle sequencing is accomplished by stepwise addition of reversible terminator nucleotides containing, for example, a cleavable or photobleachable dye label as bed, for e, in WO 04/018497 and U.S. Pat. No. 7,057,026, the disclosures of which are incorporated herein by reference. This approach is being commercialized by Solexa (now na Inc.), and is also described in WO 78 and WO 07/123,744, each of which is orated herein by reference. The availability of fluorescently-labeled ators in which both the ation can be ed and the fluorescent label cleaved facilitates efficient cyclic reversible termination (CRT) sequencing.

Polymerases can also be co-engineered to efficiently incorporate and extend from these modified nucleotides.

In reversible terminator-based sequencing embodiments, the labels may not substantially inhibit extension under SBS reaction conditions. However, the detection labels can be removable, for example, by cleavage or degradation. Images can be captured following incorporation of labels into d nucleic acid features. In particular embodiments, each cycle involves simultaneous delivery of four different nucleotide types to the array and each nucleotide type has a spectrally distinct label. Four images can then be obtained, each using a detection channel that is selective for one of the four different labels.

Alternatively, different nucleotide types can be added sequentially and an image of the array can be obtained between each addition step. In such embodiments each image will show nucleic acid features that have incorporated nucleotides of a particular type. Different features will be present or absent in the different images due the different sequence content of each e. However, the relative position of the features will remain unchanged in the images. Images ed from such reversible terminator-SBS methods can be stored, processed and analyzed as set forth herein. Following the image capture step, labels can be removed and reversible terminator moieties can be removed for subsequent cycles of nucleotide addition and detection. Removal of the labels after they have been detected in a particular cycle and prior to a subsequent cycle can e the advantage of ng background signal and crosstalk between cycles.

In particular embodiments some or all of the nucleotide monomers can include reversible terminators. In such embodiments, ible terminators/cleavable fluors can include fluor linked to the ribose moiety via a 3' ester linkage (Metzker, Genome Res. :1767-1776 (2005), which is incorporated herein by reference). Other approaches have ted the terminator chemistry from the cleavage of the fluorescence label (Ruparel et al., Proc Natl Acad Sci USA 102: 5932-7 (2005), which is incorporated herein by reference in its entirety). Ruparel et al described the development of reversible terminators that used a small 3' allyl gr oup to block extension, but could easily be deblocked by a short treatment with a palladium catalyst. The fluorophore was attached to the base via a photocleavable linker that could easily be cleaved by a 30 second exposure to long wavelength UV light. Thus, either disulfide reduction or photocleavage can be used as a cleavable linker. Another approach to reversible termination is the use of natural termination that ensues after ent of a bulky dye on a dNTP. The ce of a charged bulky dye on the dNTP can act as an effective terminator through steric and/or electrostatic nce. The presence of one incorporation event prevents further incorporations unless the dye is removed. Cleavage of the dye removes the fluor and effectively reverses the termination. Examples of modified nucleotides are also described in U.S. Pat. No. 7,427,673, and U.S. Pat. No. 7,057,026, the disclosures of which are incorporated herein by reference in their entireties.

Additional exemplary SBS systems and methods which can be utilized with the methods and systems described herein are described in U.S. Patent Application Publication No. 166705, U.S. Patent Application ation No. 2006/0188901, U.S. Pat. No. 7,057,026, U.S. Patent Application Publication No. 2006/0240439, U.S. Patent Application Publication No. 2006/0281109, PCT Publication No. WO 05/065814, U.S. Patent ation Publication No. 2005/0100900, PCT Publication No. WO 06/064199, PCT Publication No.

WO 07/010,251, U.S. Patent ation Publication No. 2012/0270305 and U.S. Patent Application Publication No. 2013/0260372, the disclosures of which are incorporated herein by reference in their entireties.

Some embodiments can utilize detection of four different nucleotides using fewer than four different labels. For example, SBS can be performed utilizing methods and systems described in the incorporated materials of U.S. Patent Application Publication No. 2013/0079232. As a first example, a pair of nucleotide types can be detected at the same wavelength, but distinguished based on a difference in intensity for one member of the pair compared to the other, or based on a change to one member of the pair (e.g. via chemical modification, photochemical cation or physical modification) that causes apparent signal to appear or disappear compared to the signal detected for the other member of the pair. As a second example, three of four different nucleotide types can be detected under particular conditions while a fourth nucleotide type lacks a label that is detectable under those conditions, or is minimally detected under those conditions (e.g., minimal ion due to background scence, etc). Incorporation of the first three nucleotide types into a c acid can be determined based on presence of their respective signals and incorporation of the fourth nucleotide type into the nucleic acid can be determined based on absence or minimal detection of any signal. As a third example, one nucleotide type can include label(s) that are ed in two different channels, s other nucleotide types are detected in no more than one of the ls. The aforementioned three exemplary configurations are not considered mutually exclusive and can be used in various combinations. An exemplary embodiment that es all three examples, is a fluorescent-based SBS method that uses a first nucleotide type that is detected in a first channel (e.g. dATP having a label that is detected in the first channel when excited by a first tion wavelength), a second nucleotide type that is detected in a second channel (e.g. dCTP having a label that is ed in the second channel when excited by a second excitation wavelength), a third tide type that is detected in both the first and the second channel (e.g. dTTP having at least one label that is ed in both channels when excited by the first and/or second excitation wavelength) and a fourth nucleotide type that lacks a label that is not, or minimally, detected in either channel (e.g. dGTP having no .

Further, as described in the incorporated materials of U.S. Patent Application Publication No. 2013/0079232, sequencing data can be obtained using a single channel. In such so-called e sequencing approaches, the first nucleotide type is labeled but the label is removed after the first image is generated, and the second nucleotide type is labeled only after a first image is generated. The third tide type retains its label in both the first and second images, and the fourth nucleotide type remains unlabeled in both images.

Some embodiments can utilize sequencing by ligation techniques. Such techniques utilize DNA ligase to incorporate oligonucleotides and identify the incorporation of such oligonucleotides. The oligonucleotides typically have different labels that are correlated with the identity of a particular nucleotide in a sequence to which the oligonucleotides ize.

As with other SBS methods, images can be obtained following treatment of an array of nucleic acid es with the labeled sequencing reagents. Each image will show nucleic acid es that have incorporated labels of a particular type. Different features will be present or absent in the ent images due the different sequence content of each feature, but the relative position of the es will remain unchanged in the images. Images obtained from ligation-based cing methods can be , processed and analyzed as set forth herein.

Exemplary SBS systems and methods which can be utilized with the methods and systems described herein are described in U.S. Pat. No. 6,969,488, U.S. Pat. No. 6,172,218, and U.S.

Pat. No. 6,306,597, the disclosures of which are incorporated herein by reference in their entireties.

Some embodiments can utilize nanopore sequencing (Deamer, D. W. & Akeson, M.

"Nanopores and nucleic acids: prospects for apid sequencing." Trends Biotechnol. 18, 147-151 (2000); Deamer, D. and D. Branton, "Characterization of nucleic acids by nanopore analysis". Acc. Chem. Res. 35:817-825 (2002); Li, J., M. Gershow, D. Stein, E. Brandin, and J. A. Golovchenko, "DNA molecules and configurations in a solid-state nanopore microscope" Nat. Mater. 2:611-615 (2003), the sures of which are incorporated herein by reference in their entireties). In such embodiments, the target nucleic acid passes through a nanopore. The nanopore can be a synthetic pore or biological membrane n, such as α- sin. As the target nucleic acid passes through the nanopore, each base-pair can be identified by measuring fluctuations in the electrical conductance of the pore. (U.S. Pat. No. 7,001,792; Soni, G. V. & Meller, "A. Progress toward ast DNA cing using solidstate nanopores." Clin. Chem. 53, 1996-2001 (2007); Healy, K. "Nanopore-based singlemolecule DNA analysis." Nanomed. 2, 459-481 (2007); Cockroft, S. L., Chu, J., Amorin, M. & Ghadiri, M. R. "A single-molecule re device detects DNA polymerase activity with -nucleotide resolution." J. Am. Chem. Soc. 130, 818-820 (2008), the disclosures of which are incorporated herein by reference in their entireties). Data obtained from re sequencing can be , processed and analyzed as set forth . In particular, the data can be d as an image in accordance with the exemplary treatment of optical images and other images that is set forth herein.

Some embodiments can utilize methods involving the real-time monitoring of DNA polymerase ty. Nucleotide orations can be detected through fluorescence resonance energy transfer (FRET) interactions between a fluorophore-bearing polymerase and γ-phosphate-labeled nucleotides as described, for example, in U.S. Pat. No. 7,329,492 and U.S. Pat. No. 7,211,414 (each of which is orated herein by reference) or nucleotide incorporations can be ed with zero-mode waveguides as described, for example, in U.S.

Pat. No. 7,315,019 (which is incorporated herein by reference) and using fluorescent nucleotide s and engineered polymerases as described, for example, in U.S. Pat. No. 7,405,281 and U.S. Patent Application Publication No. 2008/0108082 (each of which is incorporated herein by reference). The illumination can be restricted to a zeptoliter-scale volume around a surface-tethered polymerase such that incorporation of fluorescently labeled nucleotides can be observed with low background (Levene, M. J. et al. "Zero-mode waveguides for single-molecule analysis at high concentrations." Science 299, 682-686 (2003); Lundquist, P. M. et al. "Parallel confocal detection of single molecules in real time." Opt. Lett. 33, 1026-1028 (2008); Korlach, J. et al. "Selective um passivation for targeted immobilization of single DNA polymerase molecules in zero-mode waveguide nano structures." Proc. Natl. Acad. Sci. USA 105, 1176-1181 (2008), the disclosures of which are incorporated herein by reference in their entireties). Images obtained from such s can be stored, processed and analyzed as set forth .

Some SBS embodiments include detection of a proton released upon incorporation of a nucleotide into an extension product. For example, sequencing based on detection of released protons can use an electrical detector and associated techniques that are commercially available from Ion t (Guilford, CT, a Life Technologies subsidiary) or sequencing methods and systems described in US 2009/0026082 A1; US 2009/0127589 A1; US 2010/0137143 A1; or US 2010/0282617 A1, each of which is orated herein by reference. Methods set forth herein for amplifying target nucleic acids using kinetic ion can be readily applied to ates used for detecting protons. More specifically, methods set forth herein can be used to produce clonal populations of amplicons that are used to detect protons.

The above SBS methods can be advantageously carried out in multiplex formats such that multiple different target nucleic acids are manipulated simultaneously. In particular embodiments, different target nucleic acids can be treated in a common reaction vessel or on a surface of a particular substrate. This allows ient ry of sequencing reagents, removal of unreacted reagents and detection of oration events in a multiplex manner. In ments using surface-bound target nucleic acids, the target nucleic acids can be in an array format. In an array format, the target nucleic acids can be typically bound to a surface in a spatially distinguishable manner. The target nucleic acids can be bound by direct nt attachment, attachment to a bead or other particle or binding to a polymerase or other molecule that is attached to the surface. The array can include a single copy of a target nucleic acid at each site (also referred to as a feature) or le copies having the same sequence can be present at each site or feature. Multiple copies can be produced by amplification methods such as, bridge amplification or emulsion PCR.

The methods set forth herein can use arrays having features at any of a variety of densities including, for e, at least about 10 features/cm2, 100 features/cm2, 500 features/cm2, 1,000 features/cm2, 5,000 features/cm2, 10,000 features/cm2, 50,000 features/cm2, 100,000 features/cm2, 1,000,000 features/cm2, 5,000,000 features/cm2, or higher.

The methods set forth herein can provide for rapid and ent detection of a plurality of target nucleic acid in parallel. ingly the present disclosure provides integrated systems capable of preparing and detecting nucleic acids using techniques known in the art such as those exemplified above. Thus, an integrated system of the present disclosure can include fluidic components capable of delivering amplification reagents and/or sequencing reagents to one or more immobilized DNA nts, the system comprising ents such as pumps, valves, reservoirs, fluidic lines and the like. A flow cell can be configured and/or used in an integrated system for detection of target nucleic acids. Exemplary flow cells are described, for example, in US 2010/0111768 A1 and US Ser.

No. 13/273,666, each of which is incorporated herein by reference. As exemplified for flow cells, one or more of the c components of an integrated system can be used for an amplification method and for a detection method. Taking a c acid sequencing embodiment as an example, one or more of the fluidic components of an integrated system can be used for an amplification method set forth herein and for the delivery of sequencing ts in a sequencing method such as those exemplified above. Alternatively, an integrated system can include separate fluidic systems to carry out amplification s and to carry out detection methods. Examples of integrated sequencing systems that are capable of creating amplified nucleic acids and also determining the sequence of the nucleic acids include, without limitation, the MiSeqTM platform (Illumina, Inc., San Diego, CA) and devices described in US Ser. No. 13/273,666, which is incorporated herein by reference.

In some embodiments of the methods bed herein, the mapped sequence tags se sequence reads of about 20bp, about 25bp, about 30bp, about 35bp, about 40bp, about 45bp, about 50bp, about 55bp, about 60bp, about 65bp, about 70bp, about 75bp, about 80bp, about 85bp, about90bp, about 95bp, about 100bp, about 110bp, about 120bp, about 130, about 140bp, about 150bp, about 200bp, about 250bp, about 300bp, about 350bp, about 400bp, about 450bp, or about 500bp. In some cases, single-end reads of greater than 500bp are employed for reads of greater than about 1000bp when paired end reads are ted.

Mapping of the sequence tags is ed by comparing the ce of the tag with the sequence of the reference to determine the somal origin of the sequenced nucleic acid molecule, and specific genetic sequence information is not needed. A small degree of ch (0-2 mismatches per sequence tag) may be allowed to account for minor polymorphisms that may exist between the reference genome and the genomes in the mixed sample.

SYSTEMS AND APPARATUS FOR REAL TIME ANALYSIS OF SEQUENCING DATA Analysis of the sequencing data is typically performed using various computer executed algorithms and programs. Therefore, certain embodiments employ processes involving data stored in or transferred through one or more er s or other processing systems. Embodiments disclosed herein also relate to apparatus for performing these operations. This tus may be specially constructed for the required purposes, or it may be a general-purpose computer (or a group of ers) selectively activated or reconfigured by a computer program and/or data structure stored in the computer. In some embodiments, a group of processors ms some or all of the recited analytical operations collaboratively (e.g., via a network or cloud computing) and/or in parallel. A processor or group of processors for performing the methods described herein may be of various types including microcontrollers and microprocessors such as programmable devices (e.g., CPLDs and FPGAs) and non-programmable devices such as gate array ASICs or general purpose microprocessors.

In addition, certain embodiments relate to tangible and/or non-transitory computer le media or computer program ts that include program instructions and/or data (including data structures) for performing various computer-implemented operations.

Examples of computer-readable media include, but are not limited to, semiconductor memory devices, magnetic media such as disk drives, magnetic tape, optical media such as CDs, magneto-optical media, and hardware s that are specially configured to store and perform program instructions, such as nly memory devices (ROM) and random access memory (RAM). The er readable media may be directly controlled by an end user or the media may be indirectly controlled by the end user. Examples of directly controlled media include the media located at a user facility and/or media that are not shared with other entities. Examples of indirectly controlled media include media that is indirectly accessible to the user via an external network and/or via a service providing shared resources such as the “cloud.” es of program instructions include both machine code, such as produced by a compiler, and files containing higher level code that may be executed by the er using an interpreter.

In various ments, the data or information employed in the disclosed s and apparatus is ed in an electronic format. Such data or information may include reads derived from a nucleic acid sample, counts or densities of such tags that align with particular regions of a reference sequence (e.g., that align to a chromosome or chromosome t), separation distances between adjacent reads or fragments, distributions of such separation distances, ses, and the like. As used herein, data or other information provided in electronic format is available for storage on a machine and transmission between machines. Conventionally, data in electronic format is provided digitally and may be stored as bits and/or bytes in various data structures, lists, databases, etc. The data may be embodied electronically, optically, etc.

One ment provides a computer program product for determining phasing and pre-phasing coefficients, as well as phasing corrected ude values and associated base calls. The computer product may n instructions for performing any one or more of the above-described methods for phasing and base calling. As ned, the computer product may include a non-transitory and/or le computer readable medium having a computer executable or compilable logic (e.g., ctions) recorded thereon for enabling a processor to align reads, fy fragments and/or islands from aligned reads, identify alleles, including indel alleles, of heterozygous polymorphisms, phase portions of chromosomes, and haplotype chromosomes and genomes. In one example, the computer product includes (1) a computer readable medium having a computer executable or compilable logic (e.g., instructions) stored thereon for ng a processor conduct phasing correction on magnitude data (e.g., color intensity data from two or more channels) on nucleic acid samples; (2) computer assisted logic for making base calls of the nucleic acid samples; and (3) an output procedure for generating an output characterizing the nucleic acid samples.

It should be understood that it is not cal, or even possible in most cases, for an unaided human being to perform the computational operations of the methods disclosed herein. For e, generating phasing coefficients for even a single tile during a single base calling cycle might require years of effort without the assistance of a computational apparatus. Of course, the problem is compounded because le NGS sequencing generally require g correction and base calling for at least nds or even millions of reads.

The methods disclosed herein can be performed using a system for sequencing nucleic acid samples. The system may include: (a) a sequencer for receiving c acids from the test sample providing nucleic acid sequence information from the sample; (b) a processor; and (c) one or more computer-readable e media having stored n instructions for execution on the processor to evaluate data from the sequencer. The computer-readable storage media may also store partially phasing corrected magnitude data from the clusters on a flow cell.

In some embodiments, the methods are instructed by a computer-readable medium having stored thereon computer-readable ctions for carrying out a method for determining the phase of a sequence. Thus one embodiment provides a computer program product include one or more computer-readable non-transitory storage media having stored thereon computer-executable instructions that, when executed by one or more processors of a computer system, cause the er system to ent a method for sequencing a DNA sample. The method includes: (a) obtain data representing an image (e.g., the image itself) of a substrate comprising a plurality of sites where nucleic acid bases are read; (b) obtain color values (or other values representing dual bases/nucleotides) of the plurality of sites from the image of the substrate; (c) store the color values in a processor buffer; (d) retrieve partially phase-corrected color values of the plurality of sites for a base g cycle, where the partially phase-corrected color values were stored in the sequencer’s memory during an immediately preceding base calling cycle; (e) determine a prephasing correction from (i) the lly corrected color values stored during the immediately preceding base calling cycle, and (ii) the color values stored in the processor buffer; and (f) determine corrected color values from (i) the color values in the processor buffer, (ii) the partially phase corrected values stored during the immediately preceding cycle, and (iii) the prephasing tion.

Sequence or other data can be input into a computer or stored on a computer readable medium either directly or indirectly. In various embodiments, a computer system is on board or directly coupled to a sequencing device that reads and/or analyzes sequences of c acids from samples. ces or other information from such tools are provided to the computer system (or simply on board processing hardware) via a data transmission interface. In addition, the memory device may store reads, base calling quality information, phasing coefficients information, etc. The memory may also store various routines and/or programs for analyzing and presenting the sequence data. Such programs/routines may include programs for performing statistical analyses, etc.

In one example, a user provides a sample into a sequencing apparatus. Data is collected and/or analyzed by the sequencing tus which is connected to a computer.

Software on the computer allows for data collection and/or is. Data can be stored, displayed (via a monitor or other r device), and/or sent to another location. The computer may be connected to the internet which is used to transmit data to a handheld device utilized by a remote user (e.g., a physician, scientist or analyst). It is understood that the data can be stored and/or analyzed prior to ittal. In some ments, raw data is ted and sent to a remote user or apparatus that will analyze and/or store the data. For e, reads may be transmitted as they are generated, or soon thereafter, and aligned and other analyzed remotely. Transmittal can occur via the internet, but can also occur via satellite or other connection. Alternately, data can be stored on a computer-readable medium and the medium can be shipped to an end user (e.g., via mail). The remote user can be in the same or a different geographical location including, but not limited to a building, city, state, country or continent.

In some embodiments, the methods also include ting data regarding a plurality of polynucleotide sequences (e.g., reads) and sending the data to a computer or other computational system. For example, the computer can be connected to tory equipment, e.g., a sample collection apparatus, a polynucleotide amplification apparatus, or a nucleotide sequencing apparatus. The data collected or stored can be transmitted from the computer to a remote location, e.g., via a local network or a wide area network such as the internet. At the remote location various operations can be performed on the transmitted data.

In some ments of any of the systems provided herein, the sequencer is configured to perform next generation sequencing (NGS). In some embodiments, the sequencer is configured to perform massively el sequencing using sequencing-bysynthesis with reversible dye terminators. In other embodiments, the sequencer is configured to perform single le cing.

CONCLUSION The present disclosure may be embodied in other specific forms without departing from its spirit or essential characteristics. The described embodiments are to be considered in all respects only as illustrative and not restrictive. The scope of the sure is, therefore, indicated by the ed claims rather than by the foregoing description. All changes which come within the meaning and range of equivalency of the claims are to be embraced within their scope.

Claims (28)

Claims:

1. A nucleic acid sequencer comprising: an image acquisition system; ; and one or more processors ed or configured to perform a first base calling cycle by: (a) obtaining data representing an image of a substrate comprising a plurality of sites where nucleic acid bases are read by the image acquisition system, n the sites t colors representing nucleic acid base types; (b) obtaining first color values of the plurality of sites from the image of the substrate; (c) storing the first color values in a sor buffer; (d) retrieving second color values of the plurality of sites, wherein the second color values were stored in the memory during a second base calling cycle immediately preceding the first base g cycle; (e) retrieving third color values of the plurality of sites, wherein the third color values were stored in the memory during a third base calling cycle immediately preceding the second base calling cycle; and (f) determining ted color values for the first base calling cycle from the first color values in the processor buffer, the second color values, and the third color values.

2. The nucleic acid sequencer of claim 1, wherein the one or more processors are further designed or configured to perform the first base calling cycle by using the corrected color values to make base calls for the plurality of sites.

3. The nucleic acid sequencer of claim 1, wherein the one or more processors are further designed or configured to determine a pre-phasing correction from: the first color values stored in the processor buffer, and the second color values stored during the second base calling cycled immediately preceding the first base calling cycle.

4. The nucleic acid sequencer of claim 1, wherein the one or more processors are further designed or configured to determine a phasing correction from: the second color values stored during the second base g cycled ately preceding the first base calling cycle, and the third color stored during the third base g cycled ately preceding the second base calling cycle.

5. The nucleic acid sequencer of claim 1, wherein the one or more processors are further designed or configured to store the first color values in memory.

6. The nucleic acid cer of claim 1, wherein the one or more processors are further designed or configured to overwrite the third color values stored in memory with the first color values after determining the corrected color values.

7. The nucleic acid sequencer of claim 1, wherein the memory has a storage capacity of about 512 Gigabytes or less.

8. The nucleic acid sequencer of claim 1, wherein the memory is divided into a plurality of tile s, each designated for storing data enting a single image of a tile on the substrate.

9. The nucleic acid sequencer of claim 1, wherein the one or more processors are further designed or configured to m (a) – (e) in real time during base calling.

10. The nucleic acid sequencer of claim 1, wherein the nucleic acid sequencer synthesizes nucleic acids at the plurality of sites.

11. The nucleic acid sequencer of claim 1, wherein the color values are ined from only two channels of the sequencer.

12. The nucleic acid sequencer of claim 1, wherein the color values are obtained from four channels of the sequencer.

13. The nucleic acid sequencer of claim 1, wherein the substrate comprises a flow cell, wherein the flow cell is logically divided into tiles, and wherein each tile represents a region of the flow cell comprising a subset of sites, which subset is captured in a single image from the image acquisition system.

14. The nucleic acid sequencer of claim 1, further comprising, prior to operation (a), providing reagents to the flow cell and allowing the reagents to ct with sites to exhibit the colors representing nucleic acid base types during the base calling cycle.

15. A method of ining corrected color values from image data acquired, during a base g cycle, by a nucleic acid sequencer comprising an image acquisition system, one or more processors, and memory, the method comprising: (a) obtaining data representing an image of a substrate comprising a plurality of sites where nucleic acid bases are read by the image acquisition system, wherein the sites exhibit colors representing nucleic acid base types; (b) obtaining first color values of the plurality of sites from the image of the substrate; (c) storing the first color values in a processor buffer; (d) ving second color values of the plurality of sites, wherein the second color values were stored in the memory during a second base calling cycle immediately preceding the first base calling cycle; (e) ving third color values of the plurality of sites, wherein the third color values were stored in the memory during a third base calling cycle immediately preceding the second base g cycle; and (f) determining ted color values for the first base g cycle from the first color values in the processor buffer, the second color values, and the third color values.

16. The method of claim 15, further comprising performing the first base calling cycle by using the corrected color values to make base calls for the plurality of sites.

17. The method of claim 15, further comprising determining a pre-phasing correction from: the first color values stored in the processor buffer, and the second color values stored during the second base calling cycled immediately preceding the first base calling cycle.

18. The method of claim 15, further comprising determining a phasing correction from: the second color values stored during the second base calling cycled immediately preceding the first base calling cycle, and the third color stored during the third base calling cycled immediately preceding the second base calling cycle.

19. The method of claim 15, further comprising g the first color values in memory.

20. The method of claim 15, further comprising overwriting the third color values stored in memory with the first color values after determining the corrected color values.

21. The method of claim 15, wherein the memory has a storage ty of about 512 Gigabytes or less.

22. The method of claim 15, wherein the memory is divided into a plurality of tile buffers, each ated for storing data representing a single image of a tile on the substrate.

23. The method of claim 15, further comprising performing (a) – (e) in real time during base calling.

24. The method of claim 15, wherein the nucleic acid sequencer synthesizes nucleic acids at the plurality of sites.

25. The method of claim 15, wherein the color values are determined from only two channels of the sequencer.

26. The method of claim 15, wherein the color values are obtained from four channels of the sequencer.

27. The method of claim 15, n the substrate comprises a flow cell, wherein the flow cell is logically divided into tiles, and wherein each tile represents a region of the flow cell sing a subset of sites, which subset is captured in a single image from the image ition system.

28. The method of claim 15, further comprising, prior to operation (a), providing reagents to the flow cell and allowing the reagents to interact with sites to exhibit the colors representing c acid base types during the base calling cycle. Illumina, Inc. Patent Attorneys for the Applicant/Nominated Person SPRUSON & FERGUSON khcgmfi Exam “x mE. uoi “masm Mg. Emwgm Emma ................................... m0? mmrcwmmuai £me Emwgm Emu W/ comflmmsvmgx a 30E mwmmmmmm fimEmm on; mmmEm 1.5% (ohaoing contribution) 98u4% (fluorescence g) 11% (prowphasiog contribution) Green Base 3 image ooior WProphasing correction axis; Phaing correction C; Homo 2 Fan Cafii ”E Fan Caii 2 Emaga Data Par Cyaia Praaaaaaa in Raai Time 2 Raw Caiia 2 Sarfaaaaiflriaw Caii a LaaaaISarfaaa 6 lLaaa 12G TiEaaiSwath 8009+ Tiflaa af Data Par flyaia Hgara 3 as”; mguv m3m> E335. 535:2 Emou Lomou Mmmmccmxuv aQEwE Vic c gm mmmmmcmfim mwcm wmmwmcmfim Exam Awmwmmgmwcm Ewmmufim wmfimtamv mmmfih 3%me mcmwmwufinm mmmmv Start e thread and eitoeete a processor 6‘33 butter to the thread t intensity date from image for 605 the next oyoie (Emmi) and store in en ~ eveitebte tiie butter ousiy storing oidest intensity date, it any) increment eyeie, at Retrieve intensity dete (stored during which time in previous eyeie) for the oorrent eyeie (in) heeernes ) from appropriate tiie butter 3: Emu) becomes in in the tits butters 609 ve intensity data (stored two eyeies previous) for the previeus oyoie (in_1)t’rorn appropriate tiie hotter .. ...Btt Determine prephesing correction eight (h) and phasing oorreotion weight (a) for the current eyeie using extracted intensity dete tor the next oyoie (th) eiong with the current end previous intensity date List retrieved from the enroariete tiie butters Cetouiete tuiiy corrected ity date for ourrent eyeie using previous, current, and next intensity date with eppiied phasing and prephesing weights Erim-arrested) : in(frorntiie buffer) '“ b