Base identification method and device, computer equipment and storage medium

Abstract

The invention relates to a base identification method and device, computer equipment and a storage medium, and relates to the technical field of signal processing. The method comprises the following steps: acquiring target channel signal values of a plurality of base channels of a nucleotide sequence cluster to be detected in a plurality of cycles of nucleic acid sequencing; determining a library type as an unbalanced library according to each target channel signal value, and correcting a target quantile value of the target channel signal value of each base channel; and normalizing the target channel signal value of each base channel according to the target quantile value of each base channel, and determining the base type of the nucleotide sequence cluster to be detected according to the normalized target channel signal value. According to the method and the device, smooth signal value correction processing on the target quantile of the base channel is ensured, the accuracy of the target quantile value corresponding to the target quantile is improved, and the problem that the target quantile value is not accurate due to deviation of the signal value corresponding to the quantile of the base channel in the base image is solved. And the accuracy of the gene base sequence obtained by analysis is relatively low.

Description

Technical Field

[0001] This application relates to the field of signal processing technology, and in particular to a base identification method, apparatus, computer device, and storage medium. Background Technology

[0002] Second-generation sequencing technology can use high-resolution microscopic imaging to obtain gene base images after being treated with fluorescent dyes. Then, the base channels of the base images are analyzed by GPU (graphics processing unit) four-color sequencing data processing algorithms to obtain the gene base sequence.

[0003] However, when analyzing base channels in a base image, the accuracy of the obtained gene base sequence is low due to deviations in the signal values ​​corresponding to the quantile sites in the base channels. Therefore, a solution is urgently needed to correct the signal values ​​of the quantile sites to ensure the accuracy of the gene base sequence. Summary of the Invention

[0004] Therefore, it is necessary to provide a base identification method, apparatus, computer equipment, and storage medium that can correct the signal values ​​of quantile sites to address the above-mentioned technical problems.

[0005] Firstly, this application provides a base recognition method. The method includes:

[0006] Obtain the target channel signal values ​​of multiple base channels in multiple cycles of nucleic acid sequencing for the nucleotide sequence cluster to be tested;

[0007] Based on the signal values ​​of each target channel, the library type is determined to be an unbalanced library, and the target quantile values ​​of the target channel signal values ​​of each base channel are corrected.

[0008] The target channel signal values ​​of each base channel are normalized based on the target quantile values ​​of each base channel, and the base type of the nucleotide sequence cluster to be tested is determined based on the normalized target channel signal values.

[0009] In one embodiment, the target quantile value for correcting the target channel signal value of each base channel includes:

[0010] For each base channel, determine the target channel signal value corresponding to the base channel in the target image;

[0011] Determine the relationship between the target channel signal value and the signal threshold;

[0012] Based on the size relationship and the reference quantile value of the target quantile, the signal value correction processing is performed on the target quantile of the base channel to obtain the target quantile value corresponding to the target quantile.

[0013] Determine the target channel signal value corresponding to the base channel in the target image;

[0014] Determine the relationship between the target channel signal value and the signal threshold;

[0015] Based on the size relationship and the reference quantile value of the target quantile, the signal value correction processing is performed on the target quantile of the base channel to obtain the target quantile value corresponding to the target quantile.

[0016] In one embodiment, the reference quantile value includes an initial signal value, a maximum signal value, and a minimum signal value. Based on the magnitude relationship and the reference quantile value of the target quantile, signal value correction processing is performed on the target quantile of the base channel to obtain the target quantile value corresponding to the target quantile, including:

[0017] If the target channel signal value is greater than or equal to the signal threshold, then the target quantile value corresponding to the target quantile is determined based on the initial signal value of the target quantile.

[0018] If the target channel signal value is less than the signal threshold, then the target quantile of the base channel is corrected based on the maximum and minimum signal values ​​to obtain the target quantile value corresponding to the target quantile.

[0019] In one embodiment, the target quantile includes an upper quantile and a lower quantile. Based on the maximum and minimum signal values, signal value correction processing is performed on the target quantile of the base channel to obtain the target quantile value corresponding to the target quantile, including:

[0020] Based on the minimum signal value, the upper quantile of the base channel is subjected to signal value correction processing to obtain the first target quantile value corresponding to the upper quantile contained in the target quantile value;

[0021] Based on the maximum signal value, the lower quantile of the base channel is subjected to signal value correction processing to obtain the second target quantile value corresponding to the lower quantile contained in the target quantile value.

[0022] In one embodiment, determining the target quantile value corresponding to the target quantile based on the initial signal value of the target quantile includes:

[0023] Use the initial signal value of the target quantile as the target quantile value; or,

[0024] If the initial signal value is inconsistent with the standard signal value of the target quantile, the standard signal value of the target quantile shall be used as the target quantile value corresponding to the target quantile.

[0025] In one embodiment, the number of base channels is four, and determining the target channel signal value corresponding to the base channels in the target image includes:

[0026] For each base channel, the initial channel signal values ​​contained in that base channel are sorted in ascending order to obtain the signal value sorting result;

[0027] The signal value corresponding to the 99th percentile in the signal value sorting results is taken as the candidate channel signal value for that base channel;

[0028] The smallest signal value among all candidate channel signal values ​​is taken as the target channel signal value.

[0029] In one embodiment, if the library type is a balanced library, the base type of the nucleotide sequence cluster to be tested is determined based on the target channel signal value of each base channel.

[0030] Secondly, this application also provides a quantile correction device. The device includes:

[0031] The acquisition module is used to acquire the target channel signal values ​​of multiple base channels in multiple cycles of nucleic acid sequencing for the nucleotide sequence cluster to be tested.

[0032] The first determining module is used to determine the library type as an unbalanced library based on the signal values ​​of each target channel, and to correct the target quantile values ​​of the target channel signal values ​​of each base channel.

[0033] The second determining module is used to normalize the target channel signal value of each base channel according to the target quantile value of each base channel, and determine the base type of the nucleotide sequence cluster to be tested based on the normalized target channel signal value.

[0034] Thirdly, this application also provides a computer device. The computer device includes a memory and a processor, the memory storing a computer program, and the processor executing the computer program to perform the following steps:

[0035] Determine the target channel signal value corresponding to the base channel in the target image;

[0036] Determine the relationship between the target channel signal value and the signal threshold;

[0037] Based on the size relationship and the reference quantile value of the target quantile, the signal value correction processing is performed on the target quantile of the base channel to obtain the target quantile value corresponding to the target quantile.

[0038] Fourthly, this application also provides a computer-readable storage medium. The computer-readable storage medium stores a computer program thereon, which, when executed by a processor, performs the following steps:

[0039] Determine the target channel signal value corresponding to the base channel in the target image;

[0040] Determine the relationship between the target channel signal value and the signal threshold;

[0041] Based on the size relationship and the reference quantile value of the target quantile, the signal value correction processing is performed on the target quantile of the base channel to obtain the target quantile value corresponding to the target quantile.

[0042] Fifthly, this application also provides a computer program product. The computer program product includes a computer program that, when executed by a processor, performs the following steps:

[0043] Determine the target channel signal value corresponding to the base channel in the target image;

[0044] Determine the relationship between the target channel signal value and the signal threshold;

[0045] Based on the size relationship and the reference quantile value of the target quantile, the signal value correction processing is performed on the target quantile of the base channel to obtain the target quantile value corresponding to the target quantile.

[0046] The aforementioned base identification method, apparatus, computer equipment, and storage medium acquire target channel signal values ​​for multiple base channels in multiple cycles of nucleic acid sequencing of the nucleotide sequence cluster to be tested; determine the library type as an unbalanced library based on each target channel signal value; correct the target quantile values ​​of the target channel signal values ​​for each base channel; normalize the target channel signal values ​​of each base channel based on the target quantile values ​​of each base channel; and determine the base type of the nucleotide sequence cluster to be tested based on the normalized target channel signal values. As can be seen from the above, after determining that the library type is an unbalanced library based on the signal values ​​of each target channel, this application corrects the target quantile values ​​of the target channel signals of each base channel. This normalizes the target channel signal values ​​of each base channel based on their target quantile values, and determines the base type of the nucleotide sequence cluster to be tested based on the normalized target channel signal values. This ensures the smooth progress of signal value correction processing for the target quantiles of the base channels, improves the accuracy of the target quantile values ​​corresponding to the target quantiles, and prevents the occurrence of low accuracy of the analyzed gene base sequences due to deviations in the signal values ​​corresponding to the quantiles of the base channels in the base image. Attached Figure Description

[0047] Figure 1 An application environment diagram of a base recognition method provided in the embodiments of this application;

[0048] Figure 2A flowchart illustrating a base recognition method provided in this application embodiment;

[0049] Figure 3 A flowchart of another base recognition method provided in the embodiments of this application;

[0050] Figure 4 This is a flowchart illustrating the process of determining the target quantile value corresponding to the target quantile, provided in an embodiment of this application.

[0051] Figure 5 This is a flowchart illustrating the process of determining the target channel signal value provided in an embodiment of this application.

[0052] Figure 6 This is a schematic diagram of the first quantile point provided in the embodiments of this application;

[0053] Figure 7 This is a schematic diagram of a second quantile location provided in an embodiment of this application;

[0054] Figure 8 A flowchart illustrating yet another base recognition method provided in the embodiments of this application;

[0055] Figure 9 This is a structural block diagram of the first quantile correction device provided in the embodiments of this application;

[0056] Figure 10 This is a structural block diagram of the second quantile correction device provided in the embodiments of this application;

[0057] Figure 11 A structural block diagram of the third quantile correction device provided in the embodiments of this application;

[0058] Figure 12 This is a structural block diagram of the fourth quantile correction device provided in the embodiments of this application;

[0059] Figure 13 This is an internal structural diagram of a computer device in one embodiment. Detailed Implementation

[0060] To make the objectives, technical solutions, and advantages of this application clearer, the following detailed description is provided in conjunction with the accompanying drawings and embodiments. It should be understood that the specific embodiments described herein are merely illustrative and not intended to limit the scope of this application.

[0061] It should be understood that the specific embodiments described herein are merely illustrative of this application and are not intended to limit this application. In the description of this application, references to terms such as "one embodiment," "some embodiments," "example," "specific example," or "some examples," etc., indicate that a specific feature, structure, material, or characteristic described in connection with that embodiment or example is included in at least one embodiment or example of this application. In this specification, the illustrative expressions of the above terms do not necessarily refer to the same embodiment or example. Furthermore, the specific features, structures, materials, or characteristics described may be combined in a suitable manner in any one or more embodiments or examples. Moreover, without contradiction, those skilled in the art can combine and integrate the different embodiments or examples described in this specification, as well as the features of different embodiments or examples.

[0062] Based on the above, the base recognition method provided in this application embodiment can be applied to, for example... Figure 1 In the application environment shown, in one embodiment, a computer device is provided, which may be a server, and its internal structure diagram may be as follows. Figure 1 As shown, the computer device includes a processor, memory, and a network interface connected via a system bus. The processor provides computational and control capabilities. The memory includes non-volatile storage media and internal memory. The non-volatile storage media stores the operating system, computer programs, and a database. The internal memory provides an environment for the operation of the operating system and computer programs in the non-volatile storage media. The database stores data acquired using a base recognition method. The network interface communicates with external terminals via a network connection. When executed by the processor, the computer program implements a base recognition method.

[0063] This application discloses a base identification method, apparatus, computer device, and storage medium, which may specifically include the following: determining the target channel signal value corresponding to the base channel, and then, based on the relationship between the target channel signal value and the signal threshold, performing signal value correction processing on the target quantile of the base channel to obtain the target quantile value corresponding to the target quantile.

[0064] In one embodiment, such as Figure 2 As shown, Figure 2 This is a flowchart illustrating a base identification method provided in an embodiment of this application. Figure 1 The base recognition method performed by the computer device in the process may include the following steps:

[0065] S201, obtain the target channel signal values ​​of multiple base channels of the nucleotide sequence cluster to be tested in multiple cycles of nucleic acid sequencing.

[0066] In this context, a nucleotide sequence cluster is a group of similar or identical molecules, nucleotide sequences, or DNA strands; for example, a cluster can be an amplified nucleotide, an oligonucleotide, or a polynucleotide having the same or similar sequences. In an embodiment, during a base sequencing cycle, the nucleotide sequence cluster is immobilized to the reaction site and / or the reaction chamber.

[0067] Furthermore, the reaction carrier for loading DNA is a microarray. The DNA reacts on the microarray, and the base sequence at that position is determined by measuring the fluorescent groups linked in the reaction. The fluorescence signal of the DNA on the microarray can be acquired using an optical system and a scientific camera and converted into a digital signal.

[0068] The core principle of next-generation sequencing is sequencing-by-synthesis. Each biochemical reaction synthesizes one base, and the sequencer needs to collect a fluorescence signal from the chip; one such cycle is called a sequence. As a result, in each cycle, one base of the DNA sequence on the chip array can be detected. Through several cycles, the DNA sequence on the chip array can be obtained.

[0069] Nucleic acid sequencing can include multiple cycles for a cluster of nucleotide sequences to be tested. For example, a cluster of nucleotide sequences to be tested may include 100 cycles in nucleic acid sequencing.

[0070] For a given FOV acquired within a loop, signal samples are taken from the four base channels of each DNB, resulting in multiple signal values ​​for each DNB across these four base channels. These signal values ​​can be sorted in ascending order to obtain four ascending sequence of signal values ​​for each of the four base channels. There can be four base signal channels. Candidate channel signal values ​​can be selected from the signal value sequences of each base channel in the target image, using preset selectable quantiles as candidate channel signal values. The smallest candidate channel signal value is then selected as the target channel signal value.

[0071] In another embodiment, there can be two base signal channels. The signal values ​​corresponding to 99% (or 98% or 90% or other optional values) of the initial channel signal values ​​contained in each base channel in the target image can be used as candidate channel signal values ​​for each base channel, and the smallest signal value among the candidate channel signal values ​​of each base channel can be used as the target channel signal value.

[0072] S202, determine the library type as an unbalanced library based on the signal values ​​of each target channel, and correct the target quantile values ​​of the target channel signal values ​​of each base channel.

[0073] Determining whether a library is unbalanced or balanced depends primarily on the distribution of its components (such as samples, data, and literature). In fields such as biology, genetics, and sequencing technology, especially in high-throughput sequencing, library balance typically refers to whether the relative proportions of the library's components (such as different genes or different sequence fragments) are close to or reach the expected proportions.

[0074] Specifically, a balanced library refers to a library in which the proportions of its various components (such as genes and sequence fragments) are relatively equal, without significant bias or enrichment. An unbalanced library, on the other hand, refers to a library in which the proportions of certain components are significantly higher or lower than those of other components, exhibiting significant bias or enrichment.

[0075] In one embodiment of this application, when it is necessary to determine whether a library is an unbalanced library, data analysis can be performed on the signal values ​​of each target channel. By analyzing the signal values ​​of the target channels, the proportion of each component in the library can be understood. If the proportion of some components is significantly higher or lower than that of other components, it may be an unbalanced library.

[0076] S203, normalize the target channel signal value of each base channel according to the target quantile value of each base channel, and determine the base type of the nucleotide sequence cluster to be tested based on the normalized target channel signal value.

[0077] It should be noted that when it is necessary to normalize the target channel signal values ​​of each base channel, the following may be included: Calculating the normalization factor using the statistical characteristics (such as mean, median, minimum, maximum, etc.) of the target quantile values ​​and target channel signal values. The specific calculation method for the normalization factor may vary depending on the experimental design and data analysis requirements; dividing the target channel signal value of each base channel by the corresponding normalization factor yields the normalized target channel signal value.

[0078] To further explain, when it is necessary to determine the base type of a nucleotide sequence cluster to be tested, the following steps may be included: Setting a threshold: Based on the normalized target channel signal value distribution and experimental requirements, a threshold is set to distinguish different base types. Comparing signal values: The normalized target channel signal value is compared with the set threshold to determine whether the signal of each nucleotide channel reaches or exceeds the threshold. Determining the base type: Based on the comparison results, each position in the nucleotide sequence cluster to be tested is assigned to the base channel with the highest signal value, thereby determining the base type at that position.

[0079] When performing normalization and determining base types, the complexity of the experimental design, the quality and reliability of the data, and the biological background knowledge should be fully considered.

[0080] To further explain, if the library type is a balanced library, the base type of the nucleotide sequence cluster to be tested is determined based on the target channel signal value of each base channel.

[0081] The aforementioned base identification method involves obtaining target channel signal values ​​for multiple base channels in multiple cycles of nucleic acid sequencing for the nucleotide sequence cluster to be tested; determining the library type as an unbalanced library based on the target channel signal values; correcting the target quantile values ​​of the target channel signal values ​​for each base channel; normalizing the target channel signal values ​​of each base channel based on the target quantile values ​​of each base channel; and determining the base type of the nucleotide sequence cluster to be tested based on the normalized target channel signal values. As can be seen from the above, after determining that the library type is an unbalanced library based on the signal values ​​of each target channel, this application corrects the target quantile values ​​of the target channel signals of each base channel. This normalizes the target channel signal values ​​of each base channel based on their target quantile values, and determines the base type of the nucleotide sequence cluster to be tested based on the normalized target channel signal values. This ensures the smooth progress of signal value correction processing for the target quantiles of the base channels, improves the accuracy of the target quantile values ​​corresponding to the target quantiles, and prevents the occurrence of low accuracy of the analyzed gene base sequences due to deviations in the signal values ​​corresponding to the quantiles of the base channels in the base image.

[0082] In one embodiment, such as Figure 3 As shown, when it is necessary to correct the target quantile value of the target channel signal value of each base channel, the following can be included:

[0083] Step 301: Determine the target channel signal value corresponding to the base channel in the target image.

[0084] It should be noted that when it is necessary to determine the target channel signal value, the following can be included: the signal values ​​corresponding to the 99%, 98%, or 90% percentiles of the initial channel signal values ​​contained in each base channel in the target image can be used as the candidate channel signal values ​​for each base channel, and the smallest signal value among the candidate channel signal values ​​of each base channel can be used as the target channel signal value.

[0085] To further clarify, the target image contains four base channels: A (adenine), G (guanine), T (thymine), and C (cytosine).

[0086] Step 302: Determine the relationship between the target channel signal value and the signal threshold.

[0087] The signal threshold can be set or modified according to the actual situation and the staff's historical work experience. The value of the signal threshold is not limited here.

[0088] Specifically, when it is necessary to determine the signal threshold, the intermediate segment value can be taken as the signal threshold based on the data characteristics of the balanced and unbalanced libraries. For example, the average value 'a' of the minimum V99 value of multiple cycles in the balanced library and the average value 'b' of the minimum V99 value of multiple cycles in the unbalanced library can be taken as the signal threshold.

[0089] In one embodiment of this application, when it is necessary to determine the relationship between the target channel signal value and the signal threshold, the target channel signal value and the signal threshold can be compared. If the target channel signal value is greater than or equal to the signal threshold, the relationship between the target channel signal value and the signal threshold is determined as: the target channel signal value is greater than or equal to the signal threshold; if the target channel signal value is less than the signal threshold, the relationship between the target channel signal value and the signal threshold is determined as: the target channel signal value is less than the signal threshold.

[0090] Step 303: Based on the size relationship and the reference quantile value of the target quantile, perform signal value correction processing on the target quantile of the base channel to obtain the target quantile value corresponding to the target quantile.

[0091] It should be noted that the target quantile includes upper quantile and lower quantile. Therefore, when performing signal value correction processing on the target quantile of the base channel, it can be divided into performing signal value correction processing on the upper quantile of the base channel and performing signal value correction processing on the lower quantile of the base channel.

[0092] The reference quantile values ​​include the initial signal value, the maximum signal value, and the minimum signal value. Specifically, the initial signal value refers to the initial signal value corresponding to the upper quantile in the target quantile, and the initial signal value corresponding to the lower quantile in the target quantile; the maximum signal value refers to the maximum signal value corresponding to the lower quantile in multiple channels; and the minimum signal value refers to the minimum signal value corresponding to the upper quantile in the base channel.

[0093] In one embodiment of this application, when it is necessary to perform signal value correction processing on the target quantile of the base channel, the specific process may include the following: if it is determined that the target channel signal value is less than the signal threshold, then the target quantile of the base channel is corrected according to the maximum and minimum signal values ​​contained in the reference quantile values ​​to obtain the target quantile value corresponding to the target quantile.

[0094] Specifically, when performing signal value correction processing on the target quantile of the base channel based on the maximum and minimum signal values ​​contained in the reference quantile values, the following steps are included: performing signal value correction processing on the upper and lower quantiles of the base channel based on the maximum and minimum signal values ​​contained in the reference quantile values ​​to obtain the target quantile values ​​corresponding to the upper and lower quantiles.

[0095] In one embodiment of this application, when it is necessary to perform signal value correction processing on the target quantile of the base channel, the following may be included: if it is determined that the relationship between the target channel signal value and the signal threshold is that the target channel signal value is greater than or equal to the signal threshold, then the target quantile of the base channel is subjected to signal value correction processing according to the initial signal value contained in the reference quantile value to obtain the target quantile value corresponding to the target quantile.

[0096] Specifically, when performing signal value correction processing on the target quantile of the base channel based on the initial signal value contained in the reference quantile value, the following steps are included: performing signal value correction processing on the upper and lower quantiles of the base channel based on the initial signal value corresponding to the upper quantile and the initial signal value corresponding to the lower quantile in the target quantile, to obtain the target quantile values ​​corresponding to the upper and lower quantiles.

[0097] The aforementioned base identification method determines the target channel signal value corresponding to the base channel, and then performs signal value correction processing on the target quantile of the base channel based on the relationship between the target channel signal value and the signal threshold, thereby obtaining the target quantile value corresponding to the target quantile. As can be seen from the above, this application first determines the target channel signal value corresponding to the base channel, and then performs signal value correction processing on the target quantile of the base channel based on the relationship between the target channel signal value and the signal threshold. Therefore, when performing signal value correction on the target quantile, this application considers not only the target channel signal value corresponding to the base channel, but also the relationship between the target channel signal value and the signal threshold. This ensures the smooth progress of signal value correction processing on the target quantile of the base channel, improves the accuracy of the target quantile value, and prevents the low accuracy of the analyzed gene base sequence due to deviations in the signal values ​​corresponding to the quantiles of the base channels in the base image.

[0098] In one embodiment, the accuracy of the obtained gene base sequence is low due to deviations in the signal values ​​corresponding to the cleavage sites of base channels in the base image. Therefore, a solution is urgently needed to correct the signal values ​​of the cleavage sites to ensure the accuracy of the gene base sequence. To solve the above problem, a solution such as... Figure 4As shown, based on the size relationship and the reference quantile value of the target quantile, signal value correction processing is performed on the target quantile of the base channel to obtain the target quantile value corresponding to the target quantile. Specifically, this may include the following:

[0099] Step 401: If the target channel signal value is greater than or equal to the signal threshold, then determine the target quantile value corresponding to the target quantile based on the initial signal value of the target quantile.

[0100] It should be noted that when determining the target quantile value based on the initial signal value of the target quantile, the following may be included: using the initial signal value of the target quantile as the target quantile value; or, if the initial signal value is inconsistent with the standard signal value of the target quantile, using the standard signal value of the target quantile as the target quantile value.

[0101] In one embodiment of this application, if the target channel signal value is greater than or equal to the signal threshold, the initial signal value of the target quantile is not modified. Specifically, for each base channel, the initial signal value of the upper quantile in the target quantile of the base channel is used as the first target quantile value corresponding to the upper quantile in the target quantile value; the initial signal value of the lower quantile in the target quantile of the base channel is used as the second target quantile value corresponding to the lower quantile in the target quantile value.

[0102] Specifically, when it is necessary to determine the signal threshold, the intermediate segment value can be taken as the signal threshold based on the data characteristics of the balanced and unbalanced libraries. For example, the average value 'a' of the minimum V99 value of multiple cycles in the balanced library and the average value 'b' of the minimum V99 value of multiple cycles in the unbalanced library can be taken as the signal threshold.

[0103] In another embodiment of this application, if the target channel signal value is greater than or equal to the signal threshold, the target quantile value corresponding to the target quantile is determined according to the standard signal value of the predetermined target quantile. Specifically, the standard signal value includes the standard signal value corresponding to the upper quantile in each base channel and the standard signal value corresponding to the lower quantile in each base channel. Therefore, for each base channel, it is determined whether the standard signal value corresponding to the upper quantile in the base channel is the same as the initial signal value of the upper quantile. If they are the same, no processing is performed; if they are different, the standard signal value corresponding to the upper quantile is used as the target quantile value of the upper quantile. It is also determined whether the standard signal value corresponding to the lower quantile in the base channel is the same as the initial signal value of the lower quantile. If they are the same, no processing is performed; if they are different, the standard signal value corresponding to the lower quantile is used as the target quantile value of the lower quantile.

[0104] Step 402: If the size relationship is that the target channel signal value is less than the signal threshold, perform signal value correction processing on the target quantile of the base channel according to the maximum signal value and the minimum signal value, and obtain the target quantile value corresponding to the target quantile.

[0105] It should be noted that when performing signal value correction processing on the target quantile of the base channel according to the maximum signal value and the minimum signal value, the following specific contents may be included: perform signal value correction processing on the upper quantile of the base channel according to the minimum signal value, and obtain the first target quantile value corresponding to the upper quantile included in the target quantile value; perform signal value correction processing on the lower quantile of the base channel according to the maximum signal value, and obtain the second target quantile value corresponding to the lower quantile included in the target quantile value.

[0106] Furthermore, the signal value correction processing of the upper quantile and the lower quantile of the base channel can be implemented based on the dichotomy method.

[0107] In an embodiment of the present application, for each base channel, the minimum signal value can be expressed as: V U =MIN(S[C][K*P U ), where C ∈ {A base channel, G base channel, T base channel, C base channel}; the maximum signal value can be expressed as: V D =MAX(S[C][K*P D ). Furthermore, the signal value correction processing of the upper quantile and the lower quantile of the base channel can be implemented based on the dichotomy method. Specifically, the first target quantile value corresponding to the upper quantile can be determined according to the calculation formula (1), and the second target quantile value corresponding to the lower quantile can be determined according to the calculation formula (2).

[0108] The calculation formula (1) is as follows:

[0109] P U o[C]= t / K, when S[C][t - 1]< V U < S[C][t]) (1)

[0110] Where 0 < t < K, C ∈ {A base channel, G base channel, T base channel, C base channel}, P U refers to the upper quantile, and P U o refers to the first target quantile value corresponding to the upper quantile.

[0111] The calculation formula (2) is as follows:

[0112] P Do[C]= t / K, when S[C][t - 1]< V D < S[C][t] (2)

[0113] Where 0 < t < K, C ∈ {A base channel, G base channel, T base channel, C base channel}, P D refers to the lower quantile, P D o refers to the second target quantile value corresponding to the lower quantile.

[0114] In the above base recognition method, based on the magnitude relationship between the target channel signal value and the signal threshold, through the initial signal value of the target quantile, or the maximum signal value and the minimum signal value of the target quantile, signal value correction processing is performed on the target quantile of the base channel to obtain the target quantile value corresponding to the target quantile.

[0115] In one embodiment, as Figure 5 shown, the number of base channels is 4. Therefore, when it is necessary to determine the target channel signal value corresponding to the base channel in the target image, the following specific contents may be included:

[0116] Step 501, determine the candidate channel signal values corresponding to each base channel in the target image.

[0117] It should be noted that when it is necessary to determine the candidate channel signal values corresponding to each base channel, the following specific contents may be included: For each base channel, sort the initial channel signal values included in this base channel in ascending order to obtain a signal value sorting result; use the signal value corresponding to the 99% quantile in the signal value sorting result as the candidate channel signal value corresponding to this base channel.

[0118] In one embodiment of the present application, the signal value corresponding to the qu99% quantile can be expressed as: V99[C]= S[C][K * 0.99], C ∈ {A base channel, G base channel, T base channel, C base channel}. Where k is the number of sampling points.

[0119] Among them, the 99% quantile refers to the value of the probability distribution function at the 99% position in a certain probability distribution. Specifically, if the probability distribution function is cut at the 99% quantile, then only 1% of the random variable values will be higher than this point, and the remaining 99% of the values will be lower than this point.

[0120] Step 502, use the smallest signal value among the candidate channel signal values as the target channel signal value.

[0121] In one embodiment of the present application, the target channel signal value can be represented by the calculation formula (3), where the calculation formula (3) is as follows:

[0122] V99 min =MIN(V99[C])(3)

[0123] Where C ∈ {A base channels, G base channels, T base channels, C base channels}, V99 min V99 refers to the target channel signal value, while V99 refers to the candidate channel signal value.

[0124] As an example, the 100% quantile, 99% quantile, 30% quantile, and 0% quantile of each channel in the balanced library data are as follows: Figure 6 As shown; the 100%, 99%, 30%, and 0% quantiles for each channel of the imbalanced library data are as follows: Figure 7 As shown. Among them, Figure 6 and Figure 7 In the middle, the upper and lower ends of the thin line correspond to the 0% and 100% quantiles; the upper and lower ends of the thick line correspond to PuPd. Furthermore, the V99min corresponding to the unbalanced library is smaller, while the V99min corresponding to the balanced library is larger.

[0125] The aforementioned base identification method determines the candidate channel signal value corresponding to each base channel, and then determines the target channel signal value based on the smallest signal value among the candidate channel signal values. By determining the target channel signal value, a data basis is provided for subsequent signal value correction processing of the target quantile of the base channel.

[0126] In one embodiment, such as Figure 8 As shown, when signal value correction processing is required for target quantile sites of base channels, the following can be included:

[0127] Step 801: For each base channel, sort the initial channel signal values ​​contained in that base channel in ascending order to obtain the signal value sorting result.

[0128] Step 802: The signal value corresponding to the 99th percentile in the signal value sorting result is taken as the candidate channel signal value corresponding to the base channel.

[0129] Step 803: Select the smallest signal value among the candidate channel signal values ​​as the target channel signal value.

[0130] Step 804: Determine the relationship between the target channel signal value and the signal threshold. If the target channel signal value is greater than or equal to the signal threshold, proceed to step 805; if the target channel signal value is less than the signal threshold, proceed to step 806.

[0131] Step 805: Use the initial signal value of the target quantile as the target quantile value corresponding to the target quantile; or, in the case where the initial signal value is inconsistent with the standard signal value of the target quantile, use the standard signal value of the target quantile as the target quantile value corresponding to the target quantile.

[0132] Step 806: According to the minimum signal value, perform signal value correction processing on the upper quantiles of the base channels to obtain the first target quantile values corresponding to the upper quantiles included in the target quantile values; according to the maximum signal value, perform signal value correction processing on the lower quantiles of the base channels to obtain the second target quantile values corresponding to the lower quantiles included in the target quantile values.

[0133] In an embodiment of the present application, the 4-channel signal sampling data D[C][K] (C ∈ {A, C, G, T}) are sorted in ascending order respectively to obtain S[C][K] (C ∈ {A, C, G, T}); take the signal values corresponding to the 99% quantiles for each channel data: V99[C] = S[C][K*0.99] (C ∈ {A, C, G, T}); if V99min is less than the comparison threshold VC, then take PU and PD as the upper and lower quantile values respectively (0 < PD < PU < 1), and find the maximum and minimum signal values corresponding to the 4 channels: among them, PU corresponds to the minimum signal: VU = MIN(S[C][K*PU]) (C ∈ {A, C, G, T}); PD corresponds to the maximum signal: VD = MAX(S[C][K*PD]) (C ∈ {A, C, G, T}); based on the bisection method, approximate and search for the actual corresponding quantile values in the 4-channel data according to the signal values VU and VD: PUo[C] = t / K, when S[C][t - 1] < VU < S[C][t] (0 < t < K, C ∈ {A, C, G, T}); PDo[C] = t / K, when S[C][t - 1] < VD < S[C][t] (0 < t < K, C ∈ {A, C, G, T}); output the actual upper and lower quantiles PUo[C] and PDo[C] (C ∈ {A, C, G, T}) of each channel for subsequent use; if V99min is not less than the comparison threshold VC, then the upper and lower quantiles of each channel take fixed values PUo[C] = 0.99, PDo[C] = 0.3, (C ∈ {A, C, G, T}) for subsequent use.

[0134] The aforementioned base identification method determines the target channel signal value corresponding to the base channel, and then performs signal value correction processing on the target quantile of the base channel based on the relationship between the target channel signal value and the signal threshold, thereby obtaining the target quantile value corresponding to the target quantile. As can be seen from the above, this application first determines the target channel signal value corresponding to the base channel, and then performs signal value correction processing on the target quantile of the base channel based on the relationship between the target channel signal value and the signal threshold. Therefore, when performing signal value correction on the target quantile, this application considers not only the target channel signal value corresponding to the base channel, but also the relationship between the target channel signal value and the signal threshold. This ensures the smooth progress of signal value correction processing on the target quantile of the base channel, improves the accuracy of the target quantile value, and prevents the low accuracy of the analyzed gene base sequence due to deviations in the signal values ​​corresponding to the quantiles of the base channels in the base image.

[0135] It should be understood that although the steps in the flowcharts of the above embodiments are shown sequentially according to the arrows, these steps are not necessarily executed in the order indicated by the arrows. Unless explicitly stated herein, there is no strict order restriction on the execution of these steps, and they can be executed in other orders. Moreover, at least some steps in the flowcharts of the above embodiments may include multiple steps or multiple stages. These steps or stages are not necessarily completed at the same time, but can be executed at different times. The execution order of these steps or stages is not necessarily sequential, but can be performed alternately or in turn with other steps or at least some of the steps or stages of other steps.

[0136] Based on the same inventive concept, this application also provides a quantile correction device for implementing the base identification method described above. The solution provided by this device is similar to the solution described in the above method; therefore, the specific limitations in one or more quantile correction device embodiments provided below can be found in the limitations of the base identification method described above, and will not be repeated here.

[0137] In one embodiment, such as Figure 9 As shown, a quantile correction device is provided, comprising: a first acquisition module 10, a first determination module 20, and a second determination module 30, wherein:

[0138] The acquisition module 10 is used to acquire the target channel signal values ​​of multiple base channels of the nucleotide sequence cluster to be tested in multiple cycles of nucleic acid sequencing.

[0139] The first determining module 20 is used to determine the library type as an unbalanced library based on the signal values ​​of each target channel, and to correct the target quantile values ​​of the target channel signal values ​​of each base channel.

[0140] The second determining module 30 is used to normalize the target channel signal value of each base channel according to the target quantile value of each base channel, and determine the base type of the nucleotide sequence cluster to be tested according to the normalized target channel signal value.

[0141] The aforementioned quantile correction device acquires the target channel signal values ​​of multiple base channels in multiple cycles of nucleic acid sequencing for the nucleotide sequence cluster to be tested; determines the library type as an unbalanced library based on the target channel signal values; corrects the target quantile values ​​of the target channel signal values ​​for each base channel; normalizes the target channel signal values ​​of each base channel based on the target quantile values ​​of each base channel; and determines the base type of the nucleotide sequence cluster to be tested based on the normalized target channel signal values. As can be seen from the above, after determining that the library type is an unbalanced library based on the signal values ​​of each target channel, this application corrects the target quantile values ​​of the target channel signals of each base channel. This normalizes the target channel signal values ​​of each base channel based on their target quantile values, and determines the base type of the nucleotide sequence cluster to be tested based on the normalized target channel signal values. This ensures the smooth progress of signal value correction processing for the target quantiles of the base channels, improves the accuracy of the target quantile values ​​corresponding to the target quantiles, and prevents the occurrence of low accuracy of the analyzed gene base sequences due to deviations in the signal values ​​corresponding to the quantiles of the base channels in the base image.

[0142] In one embodiment, such as Figure 10 As shown, a quantile correction device is provided. The first determining module 20 of the quantile correction device includes: a first determining unit 21, a second determining unit 22, and a processing unit 23, wherein:

[0143] The first determining unit 21 is used to determine the target channel signal value corresponding to the base channel in the target image.

[0144] The second determining unit 22 is used to determine the relationship between the target channel signal value and the signal threshold.

[0145] The processing unit 23 is used to perform signal value correction processing on the target quantile of the base channel according to the size relationship and the reference quantile value of the target quantile, so as to obtain the target quantile value corresponding to the target quantile.

[0146] In one embodiment, such as Figure 11As shown, a quantile correction device is provided. The processing unit 23 in this quantile correction device includes: a first determining subunit 231 and a second determining subunit 232, wherein:

[0147] The first determining subunit 231 is used to determine the target quantile value corresponding to the target quantile based on the initial signal value of the target quantile if the target channel signal value is greater than or equal to the signal threshold.

[0148] The second determining subunit 232 is used to perform signal value correction processing on the target quantile of the base channel based on the maximum and minimum signal values ​​if the target channel signal value is less than the signal threshold, so as to obtain the target quantile value corresponding to the target quantile.

[0149] The second determining subunit is specifically used to perform signal value correction processing on the upper quantile of the base channel based on the minimum signal value, to obtain the first target quantile value corresponding to the upper quantile contained in the target quantile value. Based on the maximum signal value, it performs signal value correction processing on the lower quantile of the base channel, to obtain the second target quantile value corresponding to the lower quantile contained in the target quantile value.

[0150] The second determining subunit is further specifically used to take the initial signal value of the target quantile as the target quantile value corresponding to the target quantile. If the initial signal value is inconsistent with the standard signal value of the target quantile, the standard signal value of the target quantile is taken as the target quantile value corresponding to the target quantile.

[0151] In one embodiment, such as Figure 12 As shown, a quantile correction device is provided. The first determining unit 21 in this device includes a third determining subunit 211 and a fourth determining subunit 212, wherein:

[0152] The third determining subunit 211 is used to determine the candidate channel signal values ​​corresponding to each base channel in the target image.

[0153] The third determining subunit is specifically used to sort the initial channel signal values ​​contained in each base channel in ascending order to obtain the signal value sorting result; and to take the signal value corresponding to the 99th percentile in the signal value sorting result as the candidate channel signal value corresponding to the base channel.

[0154] The fourth determining subunit 212 is used to take the smallest signal value among the candidate channel signal values ​​as the target channel signal value.

[0155] Each module in the aforementioned quantile correction device can be implemented entirely or partially through software, hardware, or a combination thereof. These modules can be embedded in the processor of a computer device in hardware form or independent of it, or stored in the memory of a computer device in software form, so that the processor can call and execute the operations corresponding to each module.

[0156] In one embodiment, a computer device is provided, which may be a terminal, and its internal structure diagram may be as follows: Figure 13 As shown, the computer device includes a processor, memory, input / output interface, communication interface, display unit, and input device. The processor, memory, and input / output interface are connected via a system bus, and the communication interface, display unit, and input device are also connected to the system bus via the input / output interface. The processor provides computing and control capabilities. The memory includes non-volatile storage media and internal memory. The non-volatile storage media stores the operating system and computer programs. The internal memory provides an environment for the operation of the operating system and computer programs in the non-volatile storage media. The input / output interface is used for exchanging information between the processor and external devices. The communication interface is used for wired or wireless communication with external terminals; wireless communication can be achieved through Wi-Fi, mobile cellular networks, NFC (Near Field Communication), or other technologies. When executed by the processor, the computer program implements a base recognition method. The display unit is used to form a visually visible image and can be a display screen, projection device, or virtual reality imaging device. The display screen can be an LCD screen or an e-ink screen. The input device of the computer device can be a touch layer covering the display screen, or buttons, trackballs, or touchpads set on the casing of the computer device, or external keyboards, touchpads, or mice, etc.

[0157] Those skilled in the art will understand that Figure 13 The structure shown is merely a block diagram of a portion of the structure related to the present application and does not constitute a limitation on the computer device to which the present application is applied. Specific computer devices may include more or fewer components than those shown in the figure, or combine certain components, or have different component arrangements.

[0158] In one embodiment, a computer device is provided, including a memory and a processor, wherein the memory stores a computer program, and the processor executes the computer program to perform the following steps:

[0159] Obtain the target channel signal values ​​of multiple base channels in multiple cycles of nucleic acid sequencing for the nucleotide sequence cluster to be tested;

[0160] Based on the signal values ​​of each target channel, the library type is determined to be an unbalanced library, and the target quantile values ​​of the target channel signal values ​​of each base channel are corrected.

[0161] The target channel signal values ​​of each base channel are normalized based on the target quantile values ​​of each base channel, and the base type of the nucleotide sequence cluster to be tested is determined based on the normalized target channel signal values.

[0162] In one embodiment, the target channel signal value corresponding to the base channel in the target image is determined;

[0163] Determine the relationship between the target channel signal value and the signal threshold;

[0164] Based on the size relationship and the reference quantile value of the target quantile, the signal value correction processing is performed on the target quantile of the base channel to obtain the target quantile value corresponding to the target quantile.

[0165] In one embodiment, when the processor executes the computer program, it further implements the following steps: The reference quantile value includes an initial signal value, a maximum signal value, and a minimum signal value; based on the magnitude relationship and the reference quantile value of the target quantile, signal value correction processing is performed on the target quantile of the base channel to obtain the target quantile value corresponding to the target quantile, including:

[0166] If the target channel signal value is greater than or equal to the signal threshold, then the target quantile value corresponding to the target quantile is determined based on the initial signal value of the target quantile.

[0167] If the target channel signal value is less than the signal threshold, then the target quantile of the base channel is corrected based on the maximum and minimum signal values ​​to obtain the target quantile value corresponding to the target quantile.

[0168] In one embodiment, when the processor executes the computer program, it further implements the following steps: the target quantile includes an upper quantile and a lower quantile; based on the maximum and minimum signal values, signal value correction processing is performed on the target quantile of the base channel to obtain the target quantile value corresponding to the target quantile, including:

[0169] Based on the minimum signal value, the upper quantile of the base channel is subjected to signal value correction processing to obtain the first target quantile value corresponding to the upper quantile contained in the target quantile value;

[0170] Based on the maximum signal value, the lower quantile of the base channel is subjected to signal value correction processing to obtain the second target quantile value corresponding to the lower quantile contained in the target quantile value.

[0171] In one embodiment, when the processor executes the computer program, it further performs the following steps: determining the target quantile value corresponding to the target quantile based on the initial signal value of the target quantile, including:

[0172] Use the initial signal value of the target quantile as the target quantile value; or,

[0173] If the initial signal value is inconsistent with the standard signal value of the target quantile, the standard signal value of the target quantile shall be used as the target quantile value corresponding to the target quantile.

[0174] In one embodiment, when the processor executes the computer program, it further performs the following steps: If the number of base channels is four, determine the target channel signal value corresponding to the base channels in the target image, including:

[0175] Determine the candidate channel signal values ​​corresponding to each base channel in the target image;

[0176] The smallest signal value among all candidate channel signal values ​​is taken as the target channel signal value.

[0177] In one embodiment, when the processor executes the computer program, it further performs the following steps: determining the candidate channel signal values ​​corresponding to each base channel in the target image, including:

[0178] For each base channel, the initial channel signal values ​​contained in that base channel are sorted in ascending order to obtain the signal value sorting result;

[0179] The signal value corresponding to the 99th percentile in the signal value sorting results is taken as the candidate channel signal value for that base channel.

[0180] In one embodiment, a computer-readable storage medium is provided having a computer program stored thereon, the computer program performing the following steps when executed by a processor:

[0181] Obtain the target channel signal values ​​of multiple base channels in multiple cycles of nucleic acid sequencing for the nucleotide sequence cluster to be tested;

[0182] Based on the signal values ​​of each target channel, the library type is determined to be an unbalanced library, and the target quantile values ​​of the target channel signal values ​​of each base channel are corrected.

[0183] The target channel signal values ​​of each base channel are normalized based on the target quantile values ​​of each base channel, and the base type of the nucleotide sequence cluster to be tested is determined based on the normalized target channel signal values.

[0184] In one embodiment, the target channel signal value corresponding to the base channel in the target image is determined;

[0185] Determine the relationship between the target channel signal value and the signal threshold;

[0186] Based on the size relationship and the reference quantile value of the target quantile, the signal value correction processing is performed on the target quantile of the base channel to obtain the target quantile value corresponding to the target quantile.

[0187] In one embodiment, when the computer program is executed by the processor, it further implements the following steps: The reference quantile value includes an initial signal value, a maximum signal value, and a minimum signal value; based on the magnitude relationship and the reference quantile value of the target quantile, signal value correction processing is performed on the target quantile of the base channel to obtain the target quantile value corresponding to the target quantile, including:

[0188] If the target channel signal value is greater than or equal to the signal threshold, then the target quantile value corresponding to the target quantile is determined based on the initial signal value of the target quantile.

[0189] If the target channel signal value is less than the signal threshold, then the target quantile of the base channel is corrected based on the maximum and minimum signal values ​​to obtain the target quantile value corresponding to the target quantile.

[0190] In one embodiment, when the computer program is executed by the processor, it further implements the following steps: the target quantile includes an upper quantile and a lower quantile; based on the maximum and minimum signal values, signal value correction processing is performed on the target quantile of the base channel to obtain the target quantile value corresponding to the target quantile, including:

[0191] Based on the minimum signal value, the upper quantile of the base channel is subjected to signal value correction processing to obtain the first target quantile value corresponding to the upper quantile contained in the target quantile value;

[0192] Based on the maximum signal value, the lower quantile of the base channel is subjected to signal value correction processing to obtain the second target quantile value corresponding to the lower quantile contained in the target quantile value.

[0193] In one embodiment, when the computer program is executed by a processor, it further performs the following steps: determining the target quantile value corresponding to the target quantile based on the initial signal value of the target quantile, including:

[0194] Use the initial signal value of the target quantile as the target quantile value; or,

[0195] If the initial signal value is inconsistent with the standard signal value of the target quantile, the standard signal value of the target quantile shall be used as the target quantile value corresponding to the target quantile.

[0196] In one embodiment, when the computer program is executed by a processor, it further performs the following steps: If the number of base channels is four, determine the target channel signal value corresponding to the base channels in the target image, including:

[0197] Determine the candidate channel signal values ​​corresponding to each base channel in the target image;

[0198] The smallest signal value among all candidate channel signal values ​​is taken as the target channel signal value.

[0199] In one embodiment, when the computer program is executed by a processor, it further performs the following steps: determining candidate channel signal values ​​corresponding to each base channel in the target image, including:

[0200] For each base channel, the initial channel signal values ​​contained in that base channel are sorted in ascending order to obtain the signal value sorting result;

[0201] The signal value corresponding to the 99th percentile in the signal value sorting results is taken as the candidate channel signal value for that base channel.

[0202] In one embodiment, a computer program product is provided, including a computer program that, when executed by a processor, performs the following steps:

[0203] Obtain the target channel signal values ​​of multiple base channels in multiple cycles of nucleic acid sequencing for the nucleotide sequence cluster to be tested;

[0204] Based on the signal values ​​of each target channel, the library type is determined to be an unbalanced library, and the target quantile values ​​of the target channel signal values ​​of each base channel are corrected.

[0205] The target channel signal values ​​of each base channel are normalized based on the target quantile values ​​of each base channel, and the base type of the nucleotide sequence cluster to be tested is determined based on the normalized target channel signal values.

[0206] In one embodiment, the target channel signal value corresponding to the base channel in the target image is determined;

[0207] Determine the relationship between the target channel signal value and the signal threshold;

[0208] Based on the size relationship and the reference quantile value of the target quantile, the signal value correction processing is performed on the target quantile of the base channel to obtain the target quantile value corresponding to the target quantile.

[0209] In one embodiment, when the computer program is executed by the processor, it further implements the following steps: The reference quantile value includes an initial signal value, a maximum signal value, and a minimum signal value; based on the magnitude relationship and the reference quantile value of the target quantile, signal value correction processing is performed on the target quantile of the base channel to obtain the target quantile value corresponding to the target quantile, including:

[0210] If the target channel signal value is greater than or equal to the signal threshold, then the target quantile value corresponding to the target quantile is determined based on the initial signal value of the target quantile.

[0211] If the target channel signal value is less than the signal threshold, then the target quantile of the base channel is corrected based on the maximum and minimum signal values ​​to obtain the target quantile value corresponding to the target quantile.

[0212] In one embodiment, when the computer program is executed by the processor, it further implements the following steps: the target quantile includes an upper quantile and a lower quantile; based on the maximum and minimum signal values, signal value correction processing is performed on the target quantile of the base channel to obtain the target quantile value corresponding to the target quantile, including:

[0213] Based on the minimum signal value, the upper quantile of the base channel is subjected to signal value correction processing to obtain the first target quantile value corresponding to the upper quantile contained in the target quantile value;

[0214] Based on the maximum signal value, the lower quantile of the base channel is subjected to signal value correction processing to obtain the second target quantile value corresponding to the lower quantile contained in the target quantile value.

[0215] In one embodiment, when the computer program is executed by a processor, it further performs the following steps: determining the target quantile value corresponding to the target quantile based on the initial signal value of the target quantile, including:

[0216] Use the initial signal value of the target quantile as the target quantile value; or,

[0217] If the initial signal value is inconsistent with the standard signal value of the target quantile, the standard signal value of the target quantile shall be used as the target quantile value corresponding to the target quantile.

[0218] In one embodiment, when the computer program is executed by a processor, it further performs the following steps: If the number of base channels is four, determine the target channel signal value corresponding to the base channels in the target image, including:

[0219] Determine the candidate channel signal values ​​corresponding to each base channel in the target image;

[0220] The smallest signal value among all candidate channel signal values ​​is taken as the target channel signal value.

[0221] In one embodiment, when the computer program is executed by a processor, it further performs the following steps: determining candidate channel signal values ​​corresponding to each base channel in the target image, including:

[0222] For each base channel, the initial channel signal values ​​contained in that base channel are sorted in ascending order to obtain the signal value sorting result;

[0223] The signal value corresponding to the 99th percentile in the signal value sorting results is taken as the candidate channel signal value for that base channel.

[0224] It should be noted that the user information (including but not limited to user device information, user personal information, etc.) and data (including but not limited to data used for analysis, data stored, data displayed, etc.) involved in this application are all information and data authorized by the user or fully authorized by all parties, and the collection, use and processing of related data must comply with the relevant laws, regulations and standards of the relevant countries and regions.

[0225] Those skilled in the art will understand that all or part of the processes in the methods of the above embodiments can be implemented by a computer program instructing related hardware. The computer program can be stored in a non-volatile computer-readable storage medium. When executed, the computer program can include the processes of the embodiments of the above methods. Any references to memory, databases, or other media used in the embodiments provided in this application can include at least one of non-volatile and volatile memory. Non-volatile memory can include read-only memory (ROM), magnetic tape, floppy disk, flash memory, optical memory, high-density embedded non-volatile memory, resistive random access memory (ReRAM), magnetic random access memory (MRAM), ferroelectric random access memory (FRAM), phase change memory (PCM), graphene memory, etc. Volatile memory can include random access memory (RAM) or external cache memory, etc. By way of illustration and not limitation, RAM can take many forms, such as Static Random Access Memory (SRAM) or Dynamic Random Access Memory (DRAM). The databases involved in the embodiments provided in this application may include at least one type of relational database and non-relational database. Non-relational databases may include, but are not limited to, blockchain-based distributed databases. The processors involved in the embodiments provided in this application may be general-purpose processors, central processing units, graphics processing units, digital signal processors, programmable logic devices, quantum computing-based data processing logic devices, etc., and are not limited to these.

[0226] The technical features of the above embodiments can be combined in any way. For the sake of brevity, not all possible combinations of the technical features in the above embodiments are described. However, as long as there is no contradiction in the combination of these technical features, they should be considered to be within the scope of this specification.

[0227] The above embodiments are merely illustrative of several implementation methods of this application, and their descriptions are relatively specific and detailed. However, they should not be construed as limiting the scope of this application. It should be noted that those skilled in the art can make various modifications and improvements without departing from the concept of this application, and these all fall within the protection scope of this application. Therefore, the protection scope of this application should be determined by the appended claims.

Claims

1. A base recognition method, characterized in that, The method includes: Obtain the target channel signal values ​​of multiple base channels in multiple cycles of nucleic acid sequencing for the nucleotide sequence cluster to be tested; Based on the signal values ​​of each target channel, the library type is determined to be an unbalanced library, and the target quantile values ​​of the target channel signal values ​​of each base channel are corrected. The target channel signal values ​​of each base channel are normalized according to the target quantile values ​​of each base channel, and the base type of the nucleotide sequence cluster to be tested is determined based on the normalized target channel signal values.

2. The method according to claim 1, characterized in that, The target quantile values ​​for correcting the target channel signal values ​​of each base channel include: For each base channel, determine the target channel signal value corresponding to the base channel in the target image; Determine the relationship between the target channel signal value and the signal threshold; Based on the size relationship and the reference quantile value of the target quantile, the target quantile of the base channel is subjected to signal value correction processing to obtain the target quantile value corresponding to the target quantile.

3. The method according to claim 2, characterized in that, The reference quantile value includes an initial signal value, a maximum signal value, and a minimum signal value. The step of performing signal value correction processing on the target quantile of the base channel based on the magnitude relationship and the reference quantile value of the target quantile to obtain the target quantile value corresponding to the target quantile includes: If the magnitude relationship is such that the target channel signal value is greater than or equal to the signal threshold, then the target quantile value corresponding to the target quantile is determined based on the initial signal value of the target quantile. If the target channel signal value is less than the signal threshold, then the target quantile of the base channel is subjected to signal value correction processing based on the maximum signal value and the minimum signal value to obtain the target quantile value corresponding to the target quantile.

4. The method according to claim 3, characterized in that, The target quantile includes an upper quantile and a lower quantile. The step of performing signal value correction processing on the target quantile of the base channel based on the maximum and minimum signal values ​​to obtain the target quantile value corresponding to the target quantile includes: Based on the minimum signal value, the upper quantile of the base channel is subjected to signal value correction processing to obtain the first target quantile value corresponding to the upper quantile contained in the target quantile value; Based on the maximum signal value, signal value correction processing is performed on the lower quantile of the base channel to obtain the second target quantile value corresponding to the lower quantile contained in the target quantile value.

5. The method according to claim 4, characterized in that, The step of determining the target quantile value corresponding to the target quantile based on the initial signal value of the target quantile includes: The initial signal value of the target quantile is used as the target quantile value corresponding to the target quantile; or, If the initial signal value is inconsistent with the standard signal value of the target quantile, the standard signal value of the target quantile shall be used as the target quantile value corresponding to the target quantile.

6. The method according to claim 2, characterized in that, The number of base channels is 4, and determining the target channel signal value corresponding to the base channels in the target image includes: For each base channel, the initial channel signal values ​​contained in that base channel are sorted in ascending order to obtain the signal value sorting result; The signal value corresponding to the 99th percentile in the signal value sorting results is taken as the candidate channel signal value for that base channel; The smallest signal value among all candidate channel signal values ​​is taken as the target channel signal value.

7. The method according to claim 1, characterized in that, The method further includes: If the library type is a balanced library, the base type of the nucleotide sequence cluster to be tested is determined based on the target channel signal value of each base channel.

8. A quantile correction device, characterized in that, The device includes: The acquisition module is used to acquire the target channel signal values ​​of multiple base channels in multiple cycles of nucleic acid sequencing for the nucleotide sequence cluster to be tested. The first determining module is used to determine the library type as an unbalanced library based on the signal values ​​of each target channel, and to correct the target quantile values ​​of the target channel signal values ​​of each base channel. The second determining module is used to normalize the target channel signal value of each base channel according to the target quantile value of each base channel, and determine the base type of the nucleotide sequence cluster to be tested according to the normalized target channel signal value.

9. A computer device comprising a memory and a processor, wherein the memory stores a computer program, characterized in that, When the processor executes the computer program, it implements the steps of the method according to any one of claims 1 to 7.

10. A computer-readable storage medium having a computer program stored thereon, characterized in that, When the computer program is executed by a processor, it implements the steps of the method according to any one of claims 1 to 7.

11. A computer program product, comprising a computer program, characterized in that, When the computer program is executed by a processor, it implements the steps of the method according to any one of claims 1 to 7.

Images