联合检测SNV、CNV和FUSION变异的方法和装置与流程

本发明属于基因检测
技术领域：
，尤其涉及检测样本中目标基因的cnv的方法、系统和装置，以及联合检测样本中的snv、cnv和fusion变异的方法、系统和装置。
背景技术：
：细胞的dna通过凋亡、分泌或吞噬等多种机制进入血液循环系统，这种dna碎片称之为细胞游离dna(cellfreedna,cfdna)，大小通常为160-180bp。对于肿瘤患者而言，血浆中的cfdna除了来自正常的细胞外，还有部分来源于肿瘤细胞，这部分携带肿瘤细胞特异信息的dna被称为循环肿瘤dna(circulatingtumordna,ctdna)。ctdna在cfdna中的占比一般为0.1％-10％，并且随着病情阶段的不同差异很大。在肿瘤患者的ctdna中，常见的变异类型包括点突变(singlenucleotidevariant,snv)、插入缺失(insertionanddeletion,indel)、拷贝数变异(copynumbervariation,cnv)、基因融合(fusion)等。ctdna的相对含量与肿瘤的负荷和对治疗的反应是相关的，可以用于鉴定驱动基因、指导临床治疗、检测临床治疗效果及对癌症复发进行动态监控等，因此cfdna的液体活检越来越受到关注。相比传统影像学，cfdna检测可以更早的检测到癌症是否复发，此外cfdna样本收集相对容易，对于一些晚期患者很难取到组织样本。目前检测cnv的方法主要有荧光原位杂交(fish)、southern印迹杂交、数字pcr以及二代测序等。荧光原位杂交技术检测特异性高，但样本处理周期长，成本高(探针价格昂贵)，也无法做到高通量，结果判读专业性和主观性较强；数字pcr可以实现对扩增进行绝对定量，但对于样本的基因组要求较高，基因组紊乱情况下无法给出正确结果，甚至会误报；southern印记杂交技术也可以检测cnv，但是操作复杂繁琐，并且容易出现假阳性，临床推广比较困难；已有的一些二代检测技术无法在cfdna水平进行准确检测，主要表现在灵敏度不足，在肿瘤占比低的血液样本中假阴性率高。由上述可知，由于血浆中带有基因拷贝数变异的游离dna(cfdna)的浓度较低，目前检测cfdna中的cnv的方法存在灵敏度不高、特异性较低、准确度较低以及操作繁琐等缺点。因此，本领域对于有效检测低ctdna含量样本中的snv、cnv及fusion变异，尤其是低拷贝数扩增的cnv变异的方法存在持续需求。因此，本领域需要一种以高灵敏度和高特异性准确检测cnv的改进方法，从而有效提高cfdna中cnv的检出率；还需要一种联合检测cfdna样本中目标基因的snv、cnv和fusion变异的方法。技术实现要素：本发明提供了一种基于baf+logr信息，高灵敏度高特异性检测cnv的方法，该方法结合baf与测序深度的信息，可以有效提高cfdna中cnv的阳性检出率。在此基础上，本发明提供了一种联合检测snv、cnv和fusion变异的方法和装置，以有效检测ctdna占比极低的样本中的snv、cnv及fusion变异，尤其是低拷贝数扩增的cnv变异。本申请发明人发现，通过在目标基因区域及基因组上特定snp位点设计探针，并开发了利用baf和logr信息双重检测cnv的新技术，可以准确检出ctdna占比极低的样本(例如血液样本)中低频cnv扩增，克服现有技术存在的不足之处并具有极大的应用价值。在第一方面，本公开内容提供了一种检测样本中目标基因的cnv的方法，所述方法包括以下步骤：(aa)将样本原始测序数据比对到参考基因组上，生成比对结果；(a)计算每个snp位点的baf值；(b)对baf进行样本间矫正；(c)鉴定目标基因上的baf是否存在分离，以及各染色体上是否存在baf分离以评估全基因组的倍型完整程度；(d)统计原始测序深度并计算目标基因的logr值；(e)矫正logr值；(f)计算logr的背景噪音水平；和(g)按照如下规则判定目标基因是否存在cnv：当样本为整倍体时，此时如目标基因存在baf分离现象表明该基因存在cnv，此时如该基因的中位logr大于该基因的背景噪音水平或高于给定阈值t1则判定该基因存在扩增，如该基因的中位logr小于该基因的背景噪音水平的相反数或小于给定阈值t2则判定该基因存在缺失；或此时如果目标基因不存在baf分离现象，此时如该基因的logr大于该基因的背景噪音水平且高于给定阈值t3则判定该基因存在扩增，如该基因的logr小于该基因的背景噪音水平的相反数且小于给定阈值t4则判定该基因存在缺失；当样本为非整倍体时，此时如该基因的logr大于该基因的背景噪音水平且高于给定阈值t5则判定该基因存在扩增，如该基因的logr小于该基因的背景噪音水平的相反数且小于给定阈值t6则判定该基因存在缺失。在一个实施方案中，所述阈值t1为0.1至0.7范围内的任意值，t2为-0.7至-0.1范围内的任意值，t3为0.1至0.7范围内的任意值，t4为-0.7至-0.1范围内的任意值，t5为0.1至0.7范围内的任意值，和t6为-0.7至-0.1范围内的任意值；并且t1、t2、t3、t4、t5和t6可相同或不同。在一个实施方案中，所述阈值t1例如为0.1、0.2、0.3、0.4、0.5、0.6或0.7，t2例如为-0.7、-0.6、-0.5、-0.4、-0.3、-0.2或-0.1，t3例如为0.1、0.2、0.3、0.4、0.5、0.6或0.7，t4例如为-0.7、-0.6、-0.5、-0.4、-0.3、-0.2或-0.1，t5例如为0.1、0.2、0.3、0.4、0.5、0.6或0.7，和t6例如为-0.7、-0.6、-0.5、-0.4、-0.3、-0.2或-0.1。在一个实施方案中，所述阈值t1例如为0.2，t2例如为-0.2，t3例如为0.2，t4例如为-0.2，t5例如为0.3，和t6例如为-0.3。在一个实施方案中，在步骤(aa)之前，还包括选择snp位点并设计cnv探针的步骤。在一个实施方案中，步骤(b)包括：利用预先生成的baf基线文件对计算的baf值进行矫正，使得各杂合snp的baf在大量样本间的中值在0.5。在一个实施方案中，对于完整二倍体样本的杂合snp，其baf值与0.5的差值大于给定阈值则表明存在cnv，所述给定阈值例如0.05-0.2范围内的任意值，例如0.05、0.1、0.15、0.2等。在一个实施方案中，对于肿瘤样本，其baf偏离程度与目标基因的倍型和肿瘤纯度相关。在一个实施方案中，当肿瘤纯度为100%时，baf分布在0.33或0.67附近时表明该基因存在3个拷贝；而当肿瘤纯度为50%时，baf分布在0.33或0.67附近时表明该基因存在loh。在一个实施方案中，步骤(e)包括：利用预先生成的logr基线文件及广义相加模型(gam)对生成的logr值进行矫正，并对矫正后的logr的中位值平移到0值，输出矫正后的logr值。在一个实施方案中，步骤(f)包括：利用logr矫正模块输出的矫正后的logr值信息及cnv检测目标基因集作为输入评估每个目标基因的背景噪音水平，统计除去该目标基因所有探针后logr的sd值作为该目标基因的背景噪音水平并输出。在一个实施方案中，如果基因组上存在大量区域的baf分离现象，比如大于4条(例如大于4、5、6、7、8、9、10、11、12、13、14、15、16、17、18、19或20条)染色体上同时存在baf分离现象，则表明该样本的基因组为非整倍体。在一个实施方案中，如果基因组上不存在大量区域的baf分离现象，比如小于或等于4条(例如小于或等于4、5、6、7、8、9、10、11、12、13、14、15、16、17、18、19或20条)染色体上同时存在baf分离现象，则表明该样本的基因组为整倍体例如二倍体。在一个实施方案中，所述目标基因选自egfr、met、kras、nras、braf、erbb2、alk、kit、tp53、rb1、ret、ros1、fgfr2、ntrk1、ntrk2和ntrk3。在一个实施方案中，所述样本为血浆cfdna。在第一方面的一个具体实施方案中，提供一种检测样本，例如cfdna样本中目标基因的cnv的方法，所述方法包括以下步骤：(aa)将样本原始测序数据比对到参考基因组上，生成bam文件；(a)从bam文件中统计cnv相关snp位点的基因型信息，并计算每个snp位点的baf值；(b)利用预先生成的baf基线文件对计算的baf值进行矫正，使得各杂合snp的baf在大样本间的中值在0.5；(c)鉴定目标基因上的baf是否存在分离，以及各染色体上是否存在baf分离以评估全基因组的倍型完整程度；(d)从bam文件统计cnv相关探针的原始测序深度，并对测序深度进行以2为底的log转换，生成logr值；(e)利用预先生成的logr基线文件及广义相加模型(gam)对生成的logr值进行矫正，并对矫正后的logr的中位值平移到0值；(f)对于矫正后的logr值，在去除目标基因相关探针后计算余下所有cnv相关探针logr值的sd值作为该目标基因的背景噪音水平；(g)按照如下规则判定目标基因是否存在cnv：当样本为整倍体时，此时如目标基因存在baf分离现象表明该基因存在cnv，此时如该基因的中位logr大于该基因的背景噪音水平或高于给定阈值t1则判定该基因存在扩增，如该基因的中位logr小于该基因的背景噪音水平的相反数或小于给定阈值t2则判定该基因存在缺失；或此时如果目标基因不存在baf分离现象，此时如该基因的logr大于该基因的背景噪音水平且高于给定阈值t3则判定该基因存在扩增，如该基因的logr小于该基因的背景噪音水平的相反数且小于给定阈值t4则判定该基因存在缺失；当待检测样本存在多条染色体上大量snp位点的baf分离现象时表明该样本存在非整倍体扩增或缺失现象(即非整倍体时)，此时如该基因的logr大于该基因的背景噪音水平且高于给定阈值t5则判定该基因存在扩增，如该基因的logr小于该基因的背景噪音水平的相反数且小于给定阈值t6则判定该基因存在缺失，其中所述阈值t1例如为0.2，t2例如为-0.2，t3例如为0.2，t4例如为-0.2，t5例如为0.3，和t6例如为-0.3。在第二方面，本公开内容提供了一种联合检测样本中目标基因的snv、cnv和fusion的方法，所述方法包括snv检测、cnv检测和fusion变异检测，其中cnv检测根据第一方面所述的任一种方法进行。在第三方面，本公开内容提供了一种用于联合检测样本中的snv、cnv和fusion变异的装置，所述联合检测根据第二方面所述的任一种方法进行。在第三方面的一个实施方案中，所述用于联合检测样本中的snv、cnv和fusion变异的装置包括：测序数据读入模块，用于将原始测序数据比对到参考基因组并排序及去冗余；snv检测模块，用于检测样本中的所有snv；cnv检测模块，用于检测样本中目标基因的cnv；fusion变异检测模块，用于检测样本中的所有fusion变异；和结果输出模块，用于整合snv检测模块、cnv检测模块和fusion变异检测模块的结果并输出，其中cnv检测模块包括以下模块：baf计算模块，用于实施步骤(a)计算每个snp位点的baf值；baf矫正模块，用于实施步骤(b)对baf进行样本间矫正；baf分离鉴定模块，用于实施步骤(c)鉴定目标基因上的baf是否存在分离，以及各染色体上是否存在baf分离以评估全基因组的倍型完整程度；测序深度计算模块，用于实施步骤(d)统计原始测序深度并计算目标基因的logr值；logr矫正模块，用于实施步骤(e)矫正logr值；logr背景噪音计算模块，用于实施步骤(f)计算logr的背景噪音水平；和cnv判定模块，用于实施步骤(g)按照如下规则判定目标基因是否存在cnv：当样本为整倍体时，此时如目标基因存在baf分离现象表明该基因存在cnv，此时如该基因的中位logr大于该基因的背景噪音水平或高于给定阈值t1则判定该基因存在扩增，如该基因的中位logr小于该基因的背景噪音水平的相反数或小于给定阈值t2则判定该基因存在缺失；或此时如果目标基因不存在baf分离现象，此时如该基因的logr大于该基因的背景噪音水平且高于给定阈值t3则判定该基因存在扩增，如该基因的logr小于该基因的背景噪音水平的相反数且小于给定阈值t4则判定该基因存在缺失；当样本为非整倍体时，此时如该基因的logr大于该基因的背景噪音水平且高于给定阈值t5则判定该基因存在扩增，如该基因的logr小于该基因的背景噪音水平的相反数且小于给定阈值t6则判定该基因存在缺失。在一个实施方案中，所述阈值t1为0.1至0.7范围内的任意值，t2为-0.7至-0.1范围内的任意值，t3为0.1至0.7范围内的任意值，t4为-0.7至-0.1范围内的任意值，t5为0.1至0.7范围内的任意值，和t6为-0.7至-0.1范围内的任意值；并且t1、t2、t3、t4、t5和t6可相同或不同。在一个实施方案中，所述阈值t1例如为0.1、0.2、0.3、0.4、0.5、0.6或0.7，t2例如为-0.7、-0.6、-0.5、-0.4、-0.3、-0.2或-0.1，t3例如为0.1、0.2、0.3、0.4、0.5、0.6或0.7，t4例如为-0.7、-0.6、-0.5、-0.4、-0.3、-0.2或-0.1，t5例如为0.1、0.2、0.3、0.4、0.5、0.6或0.7，和t6例如为-0.7、-0.6、-0.5、-0.4、-0.3、-0.2或-0.1。在一个实施方案中，所述阈值t1例如为0.2，t2例如为-0.2，t3例如为0.2，t4例如为-0.2，t5例如为0.3，和t6例如为-0.3。在一个实施方案中，在将原始测序数据比对到参考基因组并排序及去冗余之前，还包括选择snp位点并设计cnv探针的步骤。在一个实施方案中，步骤(c)包括：利用预先生成的baf基线文件对计算的baf值进行矫正，使得各杂合snp的baf在大量样本间的中值在0.5。在一个实施方案中，对于完整二倍体样本的杂合snp，其baf值与0.5的差值大于给定阈值则表明存在cnv，所述给定阈值例如0.05-0.2范围内的任意值，例如0.05、0.1、0.15、0.2等。在一个实施方案中，对于肿瘤样本，其baf偏离程度与目标基因的倍型和肿瘤纯度相关。在一个实施方案中，当肿瘤纯度为100%时，baf分布在0.33或0.67附近时表明该基因存在3个拷贝；而当肿瘤纯度为50%时，baf分布在0.33或0.67附近时表明该基因存在loh。在一个实施方案中，步骤(e)包括：利用预先生成的logr基线文件及广义相加模型(gam)对生成的logr值进行矫正，并对矫正后的logr的中位值平移到0值，输出矫正后的logr值。在一个实施方案中，步骤(f)包括：利用logr矫正模块输出的矫正后的logr值信息及cnv检测目标基因集作为输入评估每个目标基因的背景噪音水平，统计除去该目标基因所有探针后logr的sd值作为该目标基因的背景噪音水平并输出。在一个实施方案中，如果基因组上存在大量区域的baf分离现象，比如大于4条(例如大于4、5、6、7、8、9、10、11、12、13、14、15、16、17、18、19或20条)染色体上同时存在baf分离现象，则表明该样本的基因组为非整倍体。在一个实施方案中，如果基因组上不存在大量区域的baf分离现象，比如小于或等于4条(例如小于或等于4、5、6、7、8、9、10、11、12、13、14、15、16、17、18、19或20条)染色体上同时存在baf分离现象，则表明该样本的基因组为整倍体例如二倍体。在一个实施方案中，所述目标基因选自egfr、met、kras、nras、braf、erbb2、alk、kit、tp53、rb1、ret、ros1、fgfr2、ntrk1、ntrk2和ntrk3。在一个实施方案中，所述样本为血浆cfdna。在第四方面，本公开内容提供了一种用于实施根据第一方面和第二方面所述的任一种方法的计算机系统，其包括：输入设备，用于输入样本核酸序列的测序信息；计算机存储器，用于存储计算机程序指令；计算机处理器，用于执行所述计算机程序指令，其中所述计算机程序指令实施单独的cnv检测或snv、cnv和fusion变异的联合检测，并将结果传输到输出设备；和输出设备，用于显示所述结果。在第五方面，本公开内容提供了一种计算机可读介质，其中所述计算机可读介质存储有计算机程序，其中所述计算机程序能被计算机处理器执行以实施根据第一方面和第二方面所述的任一种方法。下列描述和实施例详细阐述了本发明的实施方案。要理解的是，本发明不限于本文所述的具体实施方案并因此可改动。本领域技术人员将认识的是，存在本发明的许多变动和修改，所述变动和修改均包含在其范围之内。附图说明图1是根据本公开内容的一个实施方案的基于baf+logr的cnv检测方法；图2是根据本公开内容的一个实施方案的多基因突变联检的功能模块构成，用于联合检测snv、cnv和fusion变异。具体实施方式本申请发明人发现，通过在目标基因区域及基因组上特定snp位点设计探针，并开发了利用baf和logr信息双重检测cnv的新技术，可以准确检出ctdna占比极低的样本(例如血液样本)中低频cnv扩增，克服现有技术存在的不足之处并具有极大的应用价值。因此，本发明提供了一种基于baf+logr信息，高灵敏度高特异性检测cnv的方法，该方法结合baf与测序深度的信息，可以有效提高cfdna中cnv的阳性检出率。在此基础上，本发明提供了一种联合检测snv、cnv和fusion变异的方法和装置，以有效检测ctdna占比极低的样本中的snv、cnv及fusion变异，尤其是低拷贝数扩增的cnv变异。为了提高cnv检测的灵敏度，本发明在目标基因附近一定基因组范围内选择了若干在人群中高度杂合且gc含量在一定范围内的snp位点并设计cnv探针，同时在全基因组范围内选择了若干在人群中高度杂合且极少存在cnv且gc含量在一定范围内的snp位点并设计cnv探针，所选的snp位点用于统计目标基因及基因组范围内的baf值；同时针对目标基因的所有外显子区域设计cnv探针用于计算测序深度，经矫正后利用baf+logr算法检测目标基因的cnv。baf即b等位基因频率(ballelefrequency)，在人的基因组中，每条染色体都存在两个拷贝，因而对于每一个snp(单核苷酸多态性位点)，其在基因组上存在三种可能的形式，即aa、ab、bb型，对应的baf分别为0、0.5、1。如果该区域发生拷贝数变异，例如其中一条染色体的该区域发生扩增使得该区域变为3个拷贝，则存在aaa、aab、abb、bbb四种倍型，对应的baf分别为0、0.33、0.67、1。因而在等位基因snp杂合的情况下，偏离0.5附近的baf可以表征目标区域是否存在cnv及其可能的倍型，结合测序深度信息可以更为准确的检测cnv。在本公开内容的第一方面的一个具体实施方案中，提供了一种检测样本中目标基因的cnv的方法，所述方法包括以下步骤：(aa)将样本原始测序数据比对到参考基因组上，生成bam文件；(a)从bam文件中统计cnv相关snp位点的基因型信息，并计算每个snp位点的baf值；(c)利用预先生成的baf基线文件对计算的baf值进行矫正，使得各杂合snp的baf在大样本间的中值在0.5；(c)鉴定目标基因上的baf是否存在分离，以及各染色体上是否存在baf分离以评估全基因组的倍型完整程度；(d)从bam文件统计cnv相关探针的原始测序深度，并对测序深度进行以2为底的log转换，生成logr值；(e)利用预先生成的logr基线文件及广义相加模型(gam)对生成的logr值进行矫正，并对矫正后的logr的中位值平移到0值；(f)对于矫正后的logr值，在去除目标基因相关探针后计算余下所有cnv相关探针logr值的sd值作为该目标基因的背景噪音水平；(g)按照如下规则判定目标基因是否存在cnv：当样本为整倍体时，此时如目标基因存在baf分离现象表明该基因存在cnv，此时如该基因的中位logr大于该基因的背景噪音水平或高于给定阈值t1则判定该基因存在扩增，如该基因的中位logr小于该基因的背景噪音水平的相反数或小于给定阈值t2则判定该基因存在缺失；或此时如果目标基因不存在baf分离现象，此时如该基因的logr大于该基因的背景噪音水平且高于给定阈值t3则判定该基因存在扩增，如该基因的logr小于该基因的背景噪音水平的相反数且小于给定阈值t4则判定该基因存在缺失；当待检测样本存在多条染色体上大量snp位点的baf分离现象时表明该样本存在非整倍体扩增或缺失现象(即非整倍体时)，此时如该基因的logr大于该基因的背景噪音水平且高于给定阈值t5则判定该基因存在扩增，如该基因的logr小于该基因的背景噪音水平的相反数且小于给定阈值t6则判定该基因存在缺失。在一个实施方案中，所述阈值t1为0.1至0.7范围内的任意值，t2为-0.7至-0.1范围内的任意值，t3为0.1至0.7范围内的任意值，t4为-0.7至-0.1范围内的任意值，t5为0.1至0.7范围内的任意值，和t6为-0.7至-0.1范围内的任意值；并且t1、t2、t3、t4、t5和t6可相同或不同。在一个实施方案中，所述阈值t1例如为0.1、0.2、0.3、0.4、0.5、0.6或0.7，t2例如为-0.7、-0.6、-0.5、-0.4、-0.3、-0.2或-0.1，t3例如为0.1、0.2、0.3、0.4、0.5、0.6或0.7，t4例如为-0.7、-0.6、-0.5、-0.4、-0.3、-0.2或-0.1，t5例如为0.1、0.2、0.3、0.4、0.5、0.6或0.7，和t6例如为-0.7、-0.6、-0.5、-0.4、-0.3、-0.2或-0.1。参考图1，上述方法可以以包括若干功能模块的装置呈现，作为一个非限制性实例，这些功能模块包括：测序数据读入模块、baf计算模块、baf矫正模块、baf分离鉴定模块、测序深度计算模块、logr矫正模块、logr背景噪音计算模块、cnv判定模块、结果输出模块。以下详细描述这些功能模块。测序数据读入模块：用于读入本发明所提供的试剂盒探针捕获的文库经ngs测序后产生的原始测序数据，进而利用比对软件如bwa进行比对生成bam文件并进行排序、去冗余。baf计算模块：利用测序数据读入模块所产生的bam文件及cnv相关snp位点文件作为输入，计算每个位点的baf值并输出。baf矫正模块：利用baf计算模块输出的baf值及本发明所提供试剂盒在大规模样本中针对每个cnvsnp位点的baf所产生的基线文件作为输入，对样本中每个snp位点的baf进行样本间矫正，排除实验因素引入的误差，并输出矫正后的baf值。baf分离鉴定模块：利用baf矫正模块所产生的baf值作为输入，鉴定目标基因相关的snp位点是否存在baf分离的情况以及基因组倍型完整度。当目标基因上存在杂合snp时，杂合snp的baf如果不在0.5附近，则表明该目标基因存在拷贝数变异，其baf偏离0.5的程度与目标基因的倍型及肿瘤纯度相关，例如当肿瘤纯度为100%时，baf分布在0.33或0.67附近时表明该基因存在3个拷贝；当肿瘤纯度为50%时，baf分布在0.33或0.67附近时表明该基因存在loh。因而仅有baf的信息无法尚确定目标存在扩增或缺失。如果基因组上同时存在大量区域的baf分离现象，则表明该样本的基因组为非整倍体，即可能存在非整倍体扩增现象。测序深度(depth)计算模块：利用测序数据读入模块所产生的bam文件及cnv相关探针文件作为输入，计算每条探针的原始测序深度信息并输出。logr矫正模块：利用测序深度计算模块输出的测序深度信息及本发明所提供试剂盒在大规模正常人样本中所产生的测序深度信息作为输入对原始测序深度利用gam进行矫正，同时在样本内进行logr转换，并利用样本内所有cnv探针logr的中位值进行归零矫正，输出矫正后的logr值。logr背景噪音计算模块：利用logr矫正模块输出的矫正后的logr值信息及cnv检测目标基因集作为输入评估每个目标基因的背景噪音水平，统计除去该目标基因所有探针后logr的sd值作为该目标基因的背景噪音水平并输出。cnv判定模块：利用baf矫正模块、logr矫正模块及logr背景噪音计算模块所产生的baf、logr值及噪音水平作为输入判定目标基因是否存在cnv及相对拷贝数。当待检测样本不存在多条染色体上大量snp位点的baf分离现象时表明该样本为完整的二倍体样本，此时如目标基因存在baf分离现象表明该基因存在cnv，此时如该基因的中位logr大于该基因的背景噪音水平或高于给定阈值则判定该基因存在扩增，如该基因的中位logr小于该基因的背景噪音水平的相反数或小于给定阈值则判定该基因存在缺失。当待检测样本为完整的二倍体样本，此时如果目标基因不存在baf分离现象，此时如该基因的logr大于该基因的背景噪音水平且高于给定阈值则判定该基因存在扩增，如该基因的logr小于该基因的背景噪音水平的相反数且小于给定阈值则判定该基因存在缺失。当待检测样本存在多条染色体上大量snp位点的baf分离现象时表明该样本存在非整倍体扩增或缺失现象，此时如该基因的logr大于该基因的背景噪音水平且高于给定阈值则判定该基因存在扩增，如该基因的logr小于该基因的背景噪音水平的相反数且小于给定阈值则判定该基因存在缺失。结果输出模块：用于绘制cnv相关分析图及输出最终cnv结果。在一个具体实施方案中，所述判断目标基因baf是否存在分离的规则如下：如目标基因上存在大于给定阈值数目的杂合snp位点，比如5，则判断目标基因所有杂合snp中是否存在连续的给定阈值数目的杂合snp位点，使得这些连续的杂合snp位点的baf同时分布在0.5的两侧，且分布在0.5两侧的任意两个snp对应baf之间的差值大于给定阈值，如0.1；或使得这些连续的杂合snp位点的baf仅分布在0.5的单侧，且任意单个snp对应baf与0.5之间的差值大于给定阈值，如0.1；如满足上述条件，即判定目标基因存在baf分离，否则判定不存在baf分离。在一个具体实施方案中，所述判定基因组倍型是否完整的规则如下：对于任意一条染色体，如其包含的所有杂合snp中存在任意snp的baf与0.5之间的差值大于给定阈值，如0.1，则判定该染色体存在cnv；如果同时存在cnv的染色体数据大于给定阈值，如5，则判定该样本基因组倍型不完整；反之，则倍型完整。在一个具体实施方案中，所述判定目标基因是否存在cnv的规则如下：如果该样本基因组倍型完整，如果目标基因的baf存在分离，此时如果目标基因对应所有探针的中位logr值大于该基因对应背景噪音水平或者大于给定阈值，如0.2，则判定该基因存在cnv扩增；此时如果目标基因对应所有探针的中位logr值小于该基因对应背景噪音水平的相反数或者小于给定阈值，如-0.2，则判定该基因存在cnv缺失；或如果目标基因的baf不存在分离，此时如果目标基因对应所有探针的中位logr值大于该基因对应背景噪音水平且大于给定阈值，如0.2，则判定该基因存在cnv扩增；此时如果目标基因对应所有探针的中位logr值小于该基因对应背景噪音水平的相反数且小于给定阈值，如-0.2，则判定该基因存在cnv缺失；或如果该样本基因组倍型不完整，此时如果目标基因对应所有探针的中位logr值大于该基因对应背景噪音水平且大于给定阈值，如0.3，则判定该基因存在cnv扩增；此时如果目标基因对应所有探针的中位logr值小于该基因对应背景噪音水平的相反数且小于给定阈值，如-0.3，则判定该基因存在cnv缺失。在本公开内容的第二方面的一个具体实施方案中，提供了提供了一种联合检测样本中目标基因的snv、cnv和fusion的方法，所述方法包括snv检测、cnv检测和fusion变异检测，其中cnv变异检测根据第一方面所述的任一种方法进行。在本公开内容的第三方面的一个具体实施方案中，提供了一种用于联合检测样本中的snv、cnv和fusion变异的装置，所述联合检测根据第二方面所述的任一种方法进行。在第三方面的一个实施方案中，所述用于联合检测样本中的snv、cnv和fusion变异的装置包括：测序数据读入模块，用于将原始测序数据比对到参考基因组并排序及去冗余；snv检测模块，用于检测样本中的所有snv；cnv检测模块，用于检测样本中目标基因的cnv；fusion变异检测模块，用于检测样本中的所有fusion变异；和结果输出模块，用于整合snv检测模块、cnv检测模块和fusion变异检测模块的结果并输出，其中cnv检测模块包括以下模块：baf计算模块，用于实施步骤(a)计算每个snp位点的baf值；baf矫正模块，用于实施步骤(b)对baf进行样本间矫正；baf分离鉴定模块，用于实施步骤(c)鉴定目标基因上的baf是否存在分离，以及各染色体上是否存在baf分离以评估全基因组的倍型完整程度；测序深度计算模块，用于实施步骤(d)统计原始测序深度并计算目标基因的logr值；logr矫正模块，用于实施步骤(e)矫正logr值；logr背景噪音计算模块，用于实施步骤(f)计算logr的背景噪音水平；和cnv判定模块，用于实施步骤(g)按照如下规则判定目标基因是否存在cnv：当样本为整倍体时，此时如目标基因存在baf分离现象表明该基因存在cnv，此时如该基因的中位logr大于该基因的背景噪音水平或高于给定阈值t1则判定该基因存在扩增，如该基因的中位logr小于该基因的背景噪音水平的相反数或小于给定阈值t2则判定该基因存在缺失；或此时如果目标基因不存在baf分离现象，此时如该基因的logr大于该基因的背景噪音水平且高于给定阈值t3则判定该基因存在扩增，如该基因的logr小于该基因的背景噪音水平的相反数且小于给定阈值t4则判定该基因存在缺失；当样本为非整倍体时，此时如该基因的logr大于该基因的背景噪音水平且高于给定阈值t5则判定该基因存在扩增，如该基因的logr小于该基因的背景噪音水平的相反数且小于给定阈值t6则判定该基因存在缺失，其中所述阈值t1例如为0.2，t2例如为-0.2，t3例如为0.2，t4例如为-0.2，t5例如为0.3，和t6例如为-0.3。参考图2，该第三方面的装置可以包括若干功能模块，作为一个非限制性实例，这些功能模块包括：测序数据读入模块，用于将原始测序数据比对到参考基因组并排序及去冗余；snv检测模块，用于检测样本中的所有snv；cnv检测模块，用于检测样本中目标基因的cnv；fusion变异检测模块，用于检测样本中的所有fusion变异；和结果输出模块，用于整合snv检测模块、cnv检测模块和fusion变异检测模块的结果并输出。在第四方面，本公开内容提供了一种用于实施根据第一方面和第二方面所述的任一种方法的计算机系统，其包括：输入设备，用于输入样本核酸序列的测序信息；计算机存储器，用于存储计算机程序指令；计算机处理器，用于执行所述计算机程序指令，其中所述计算机程序指令实施单独的cnv检测或snv、cnv和fusion变异的联合检测，并将结果传输到输出设备；和输出设备，用于显示所述结果。在第五方面，本公开内容提供了一种计算机可读介质，其中所述计算机可读介质存储有计算机程序，其中所述计算机程序能被计算机处理器执行以实施根据第一方面和第二方面所述的任一种方法。本文所用的术语仅以描述具体的实施方案为目的而不意图限制本发明。除非上下文另有明确指示，否则本文所用的单数形式“一”、“一个”和“所述”也意图包括复数形式。此外，开放式的表述“包括”和“包含”解释为还可以含有没有述及的结构组成部分或方法步骤，但需要注意的是，该开放式的表述也涵盖仅由所述的组分和方法步骤组成的情形(即涵盖了封闭式表述“由……组成”的情形)。如全文所用，范围用作描述该范围内的每个数值和所有数值的简写形式。范围内的任何数值例如整数值、以十分之一递增的值(当范围的端值为小数点后一位时)或以百分之一递增的值(当范围的端值为小数点后二位时)都可选做该范围的终点。例如，范围0.1-10用作描述该范围内的所有数值，例如0.1、0.2、0.3、0.4、0.5、0.6、0.7、0.8……9.5、9.6、9.7、9.8、9.9和10(以十分之一递增的值)，并且包括所有子范围，例如0.1-1.0、2.0-3.0、4.0-5.0、6.0-7.0、8.0-9.0等。本说明书中提及的所有科学技术术语具有与本领域技术人员通常理解的含义相同的含义，在冲突的情况下以本说明书中的定义为准。为了使本发明的描述更易于理解，以下先对一些术语进行解释说明。本文所述的术语“建库”又称文库构建，是指，对于血液、体液或粪便等样本中存在的cfdna进行修复并连接到一段已知dna片段即adapter序列(也称为接头)上，从而可以用于在illumina设备上进行高通量dna测序的过程。本发明中所称“建库”是指用于高通量测序的建库。本文所述的术语“高通量测序”又可以称为下一代测序技术(nextgenerationsequencing,ngs)、大规模平行测序(massivelyparallelsequencing(mps))，是指采用“边合成边测序”的原理、对于几十万到几百万dna分子同时进行平行的测序反应，然后通过生物信息学分析所得到的原始图像数据或电化学信号、最终得到待测样品的核酸序列或拷贝数等信息的测序技术，又称为高通量测序、深度测序、二代测序等。高通量测序的基本程序是将待测dna随机打断成小片段，经末端修复、连接接头序列、pcr等步骤进行文库构建，最后使用illumina，iontorrent等测序仪进行测序。本文所述的术语“捕获测序”是指利用生物素标记的dna或rna探针，对dna样本中的目标片段进行捕获，并进行测序的技术。本文所述的术语“cfdna”，又称游离dna(cellfreedna)，是指存在血浆或血清、脑脊液等细胞外游离状态的核酸片段，约160-180bp，是细胞dna在生理或病理条件下的产物。cfdna可通过分泌或细胞死亡过程释放到循环中，例如细胞坏死或凋亡。某些cfdna为ctdna(见下文)。本文所述的术语“循环肿瘤dna(circulatingtumordna,ctdna)”是指起源于肿瘤的游离dna(cfdna)级分。本文所述的术语“突变”是指在生物或细胞的基因组中的遗传改变。例如，目的突变可以是相对于生物生殖系的变化，例如癌细胞特异性变化。突变可包括单核苷酸变异(snv)、拷贝数变异(cnv)、单核苷酸多态性(snp)、插入、缺失和重排(例如fusion(融合))。本文所述的术语“snv(singlenucleotidevariant，单核苷酸变异)”是一种体细胞单核苷酸突变。本文所述的术语“cnv(copynumbervariation，拷贝数变异)”是指由基因组发生重排而导致的，一般指长度为1kb以上的基因组大片段的拷贝数增加或者减少。本文所述的术语“fusion(融合)”是指两个或多个基因的编码区或非编码区断裂，重新首尾相连，置于同一套调控序列。本文所述的术语“snp”，即单核苷酸多态性(singlenucleotidepolymorphism)，是指由于单个核苷酸改变而导致的核酸序列多态性。一般来说，一个snp位点只有两种等位基因，因此又叫双等位基因。snp在人类基因组中的发生频率比较高，大约平均每1000个碱基中就有一个多态性位点。本文所述的术语“癌症”和“肿瘤”在本文中可交换使用以及可以指表现出自发的不受调控生长的细胞或组织，由此所述细胞或组织表现出以增加的细胞增殖(例如显著的细胞增殖)为特征的异常生长。用于本申请的检测、分析或处理的目的细胞可包括但不限于癌前细胞(例如良性细胞)、恶性细胞、转移前细胞、转移性细胞和非转移性细胞。事实上各组织的癌症均为本领域技术人员已知的，包括实体瘤，例如肉瘤、成胶质细胞瘤、黑素瘤、淋巴瘤、骨髓瘤等，以及循环癌症，例如白血病。癌症的实例包括但不限于卵巢癌、乳腺癌、结肠癌、肺癌、前列腺癌、肝细胞癌、胃癌、胰腺癌、子宫颈癌、卵巢癌、肝癌、膀胱癌、泌尿道癌、甲状腺癌、肾癌、癌瘤、黑素瘤、头颈癌和脑癌。虽然上文已描述了本发明的各种实施方案，但是应理解的是，其仅以实例的方式提供，而并非限制。对公开的实施方案的许多改变可依照本文的公开内容来进行，而不会背离本发明的精神或范围。因此，本发明的广度和范围不应受到任何上述的实施方案所限制。本文提及的所有文献都通过引用结合到本文中。本申请引用的所有出版物和专利文件都为所有目的而通过引用结合，引用程度如同单独地指出各个出版物或专利文件一样。实施例除非另外说明，否则本文实施例所用的材料均市购获得，用于进行实验的各种具体实验方法均为本领域常规的实验方法或者按照制造商所建议的步骤和条件，并能由本领域技术人员根据需要常规地确定。实施例1血液中多基因突变联检试剂盒的探针设计方案根据fda肿瘤治疗用药指南等公开信息报导选取了肺肠癌相关存在cnv的重要基因如met、erbb2，以met及erbb2基因范围为中心，在长度为1mb的基因组范围选择东亚人群频率在0.4-0.6之间、gc含量在0.3-0.7之间的snp位点；同时在全基因组范围内选择500个东亚人群频率在0.4-0.6之间、gc含量在0.4-0.6之间、且在人群间倍型高度稳定的snp位点。针对所选的这些snp位点及met、erbb2的外显子范围设计捕获探针，探针长度120nt，目标区域保持两层探针覆盖。根据fda肿瘤治疗用药指南等公开信息报导选取了肺肠癌相关存在snv的重要基因如egfr、met、kras、nras、braf、erbb2、alk、kit、tp53、rb1，针对其用药突变位点设计捕获探针，探针长度120nt，目标区域保持两层探针覆盖。根据fda肿瘤治疗用药指南等公开信息报导选取了肺肠癌相关存在fusion的重要基因如alk、ret、ros1、fgfr2、ntrk1、ntrk2、ntrk3，从cosmic等数据库中选择其常见融合断点所在内含子区域，针对这些内含子区域设计捕获探针，探针长度120nt，目标区域保持两层探针覆盖。上述探针合成后按相同浓度混合用于cfdna文库捕获，捕获的文库用于ngs测序。实施例2血液中多基因突变联检试剂盒的建库方案1.末端修复反应液配制先从冰箱-20℃取出下列试剂，解冻后震荡混匀，单个样本配制量参见表1。表1：体积cfdna25µl10x末端修复酶缓冲液3µl末端修复酶1.5µl总体积29.5µl2.末端修复反应向200ul离心管中加入4.5ul分装好的mix后按照表2程序进行反应。表2：步骤温度时间120℃30min24℃∞3.连接1反应mix配制先从冰箱-20℃取出下列试剂，解冻后震荡混匀，单个样本配制量参见表3。表3：体积连接缓冲液112.5µl连接接头11µl连接酶11.5µl总体积15µl配制后放置于冰盒上备用。4.纯化末端修复反应结束后使用1.8xampure磁珠进行纯化，使用15ul连接1反应mix回溶磁珠，并保证磁珠充分混匀。5.连接1反应将pcr管放置于pcr仪上，按照表4进行反应：表4：步骤温度时间120℃15min265℃15min34℃∞6.连接2反应mix配制先从冰箱-20℃取出下列试剂，解冻后震荡混匀，单个样本配制量参见表5。表5：体积连接缓冲液22.25µl连接接头22µl连接酶2a0.25µl连接酶2b0.5µl总体积5µl7.连接2反应于pcr管中加入5ul连接2反应mix后混匀保证磁珠充分重悬，并按照表6程序进行反应。表6：步骤温度时间165℃30min24℃∞8.纯化于pcr管中加入2.5xpeg/nacl溶液进行纯化，回溶于21ulnf水中并转移20ul至新的pcr管备用。9.pcr扩增反应液配制按照表7配制pcr反应液：表7：体积连接产物20uludiprimerpairs5µlhifihotstartreadymix25µl总体积50µl10.pcr反应按照表8的程序进行pcr反应：表8：11.纯化pcr反应结束后使用1.3xampure磁珠进行纯化，最后使用30ulnf水回溶。12.文库质检使用4150芯片和qubit进行文库片段大小和浓度的检验。13.杂交：13.1.杂交文库准备将待杂交文库放置于冰上融化后，取2ug置于1.5mlep管中，再加入5ulhumancot-1dna和1ulblockingoligomix后真空浓缩至干燥。13.2.杂交mix配制按照表9配制杂交反应mix。表9：体积2x杂交反应buffer8.5ul杂交反应增强液2.7ulprobe2ulnf水3.8ul总体积17ul13.3.杂交反应加入17ul杂交反应mix，吹吸混匀后离心将mix收集到管底，然后避光静置10min重溶。将重溶后的反应mix转移至0.2ml离心管中，按照表10的反应程序进行杂交：表10：程序时间95℃30s65℃16h65℃hold13.4.洗脱1)按照表11配制1x捕获buffer：表11：nf水buffer2xbeadswashbuffer16016010xwashbuffer12522810xwashbuffer21441610xwashbuffer31441610xstringentwashbuffer288322)将dynabeads®m270streptavidinc1从4℃冰箱取出，室温平衡30min；3)将m270涡旋均匀后取出50ul放置于1.5ml离心管，将离心管静置于磁力架上，beads完全分离后去掉上清；4)加入100ul1xbeadswashbuffer，涡旋均匀10sec后静置于磁力架上，beads完全分离后去掉上清，共清洗3次；5)按照表12配制磁珠重悬缓冲液。表12：nf水2x杂交反应buffer8.5ul杂交反应增强液2.7ulnf水5.8ul总体积17ul6)将17ul磁珠重悬缓冲液加入至磁珠中，涡旋均匀后分装到0.2mlpcr管中；7)杂交反应持续16h后，将17ul重悬后的m270磁珠加入到pcr仪上的杂交管中，涡旋混匀，然后放回pcr仪上65℃继续孵育45min，期间每12min短暂震荡约3s；8)孵育结束后加入100ul预热的1xwashbufferⅰ，短暂震荡约3s，然后转移到1.5mllobind离心管中，震荡大约3sec后短暂离心并立即静置于磁力架上，大约20s后beads完全分离，立即去掉上清；9)加入150ul预热的1xstringentwashbuffer，吹吸10次后短暂离心并立即静置于65℃金属浴上孵育5min，然后静置于磁力架上，大约20s后beads完全分离，立即去掉上清，重复该步骤一次；10)加入150ul1xwashbufferⅰ，振荡2min，后短暂轻离心，然后静置于磁力架上，beads完全分离后去掉上清；11)加入150ul1xwashbufferⅱ，振荡2min，后短暂轻离心，然后静置于磁力架上，beads完全分离后去掉上清；12)加入150ul1xwashbufferⅲ，振荡2min，后短暂轻离心，然后静置于磁力架上，beads完全分离后去掉上清；13)加入18ulnf水重悬磁珠备用。13.5.pcr扩增按照下表进行pcr反应mix的配制：表13：反应组分体积(ul)重悬的beads20kapahifihotstartreadymix25p5/p7primermix5总计50涡旋均匀后进行以下反应：表14：13.6.pcr产物回收加入1.5xampure磁珠进行产物回收，最后使用20ulnf水回溶。13.7.文库质检使用4150芯片和qubit进行文库片段大小和浓度的检验。14.上机测序构建的文库使用cn500进行测序。实施例3细胞系混合模拟cfdna样本中多基因突变检测将含有egfr19del突变的细胞系、含有krasg12d突变的细胞系、含有eml4-alk融合突变的细胞系以及含有met扩增的细胞系gdna打断，模拟cfdna，用阴性细胞系gdna稀释，使得egfr19del、krasg12d、eml4-alk的目标突变丰度分别在0.4%、0.2%及0.1%水平，使得met扩增的绝对拷贝数分别在10、6、3.5、2.5拷贝数，利用实施例1的试剂盒及实施例2的实验条件建库测序获得原始测序数据。原始测序数据经bwa比对到参考基因组hg19生成bam文件，利用sambamba排序去冗余。利用常规方法检测其中snv、fusion突变。利用本发明提供的cnv检测方法检测cnv。cnv检测具体步骤如下：1)从bam文件中统计cnv相关snp位点的基因型信息，并计算每个snp位点的baf值；2)利用预先生成的baf基线文件对计算的baf值进行矫正，使得各杂合snp的baf在大样本间的中值在0.5；3)鉴定目标基因上的baf是否存在分离，以及各染色体上是否存在baf分离以评估全基因组的倍型完整程度；4)从bam文件统计cnv相关探针的原始测序深度，并对测序深度进行以2为底的log转换，生成logr值；5)利用预先生成的logr基线文件及广义相加模型(gam)对生成的logr值进行矫正，并对矫正后的logr的中位值平移到0值；6)对于矫正后的logr值，在去除目标基因相关探针后计算余下所有cnv相关探针logr值的sd值作为该目标基因的背景噪音水平；7)按照如下规则判定目标基因是否存在cnv：当样本基因组倍型完整时，此时如目标基因存在baf分离现象表明该基因存在cnv，此时如该基因的中位logr大于该基因的背景噪音水平或高于给定阈值则判定该基因存在扩增，如该基因的中位logr小于该基因的背景噪音水平的相反数或小于给定阈值则判定该基因存在缺失；或此时如果目标基因不存在baf分离现象，此时如该基因的logr大于该基因的背景噪音水平且高于给定阈值则判定该基因存在扩增，如该基因的logr小于该基因的背景噪音水平的相反数且小于给定阈值则判定该基因存在缺失；当样本基因组倍型不完整时，此时如该基因的logr大于该基因的背景噪音水平且高于给定阈值则判定该基因存在扩增，如该基因的logr小于该基因的背景噪音水平的相反数且小于给定阈值则判定该基因存在缺失。所述预先生成的baf基线文件构建规则如下：利用本发明实施例1中提供的试剂盒对300例正常人样本建库测序，比对后统计各cnv相关snp位点的baf值，统计杂合位点baf分布的中值，该中值与0.5之间的差别作为对应snp位点baf的基线值。所述预先生成的logr基线文件构建规则如下：利用本发明实施例1中提供的试剂盒对30例正常人样本建库测序，比对后统计各cnv相关探针的logr值并进行样本内logr的中值归零矫正，统计每条探针logr的中值，作为对应探针logr的基线值。所述判断目标基因baf是否存在分离的规则如下：如目标基因上存在大于给定阈值数目的杂合snp位点，比如5，则判断目标基因所有杂合snp中是否存在连续的给定阈值数据的杂合snp位点，使得这些连续的杂合snp位点的baf同时分布在0.5的两侧，且分布在0.5两侧的任意两个snp对应baf之间的差值大于给定阈值，如0.1；或使得这些连续的杂合snp位点的baf仅分布在0.5的单侧，且任意单个snp对应baf与0.5之间的差值大于给定阈值，如0.1；如满足上述条件，即判定目标基因存在baf分离，否则判定不存在baf分离。所述判定基因组倍型是否完整的规则如下：对于任意一条染色体，如其包含的所有杂合snp中存在任意snp的baf与0.5之间的差值大于给定阈值，如0.1，则判定该染色体存在cnv；如果同时存在cnv的染色体数据大于给定阈值，如5，则判定该样本基因组倍型不完整；反之，则倍型完整。所述判定目标基因是否存在cnv的阈值如下：如果该样本基因组倍型完整，如果目标基因的baf存在分离，此时如果目标基因对应所有探针的中位logr值大于该基因对应背景噪音水平或者大于给定阈值，如0.2，则判定该基因存在cnv扩增；此时如果目标基因对应所有探针的中位logr值小于该基因对应背景噪音水平的相反数或者小于给定阈值，如-0.2，则判定该基因存在cnv缺失；或如果目标基因的baf不存在分离，此时如果目标基因对应所有探针的中位logr值大于该基因对应背景噪音水平且大于给定阈值，如0.2，则判定该基因存在cnv扩增；此时如果目标基因对应所有探针的中位logr值小于该基因对应背景噪音水平的相反数且小于给定阈值，如-0.2，则判定该基因存在cnv缺失；或如果该样本基因组倍型不完整，此时如果目标基因对应所有探针的中位logr值大于该基因对应背景噪音水平且大于给定阈值，如0.3，则判定该基因存在cnv扩增；此时如果目标基因对应所有探针的中位logr值小于该基因对应背景噪音水平的相反数且小于给定阈值，如-0.3，则判定该基因存在cnv缺失。表15显示了细胞系混样snv、fusion检测结果，其中阳性表示检出变异，阴性表示未检出。表16显示了细胞系混样cnv检测结果，其中扩增阳性表示检出met存在cnv扩增。表15：注：vaf为变异等位基因分数(variantallelefraction)。表16：突变类型ddpcr拷贝数ngs拷贝数检测结果met扩增1010.14扩增阳性met扩增1010.1扩增阳性met扩增5.86.21扩增阳性met扩增5.86.16扩增阳性met扩增3.53.7扩增阳性met扩增3.53.76扩增阳性met扩增2.52.78扩增阳性met扩增2.52.79扩增阳性实施例4临床血液样本中多基因突变联检将含有egfr19del突变的临床血液样本、含有egfrl858r突变的临床血液样本、含有egfrt790m突变的临床血液样本、含有krasg12d突变的临床血液样本、含有krasg13d突变的临床血液样本、含有brafv600e突变的临床血液样本、含有eml4-alk融合突变的临床血液样本以及含有met扩增的临床血液样本cfdna以及利用阴性细胞系提取的gdna打断模拟cfdna稀释，使用ddpcr对各突变进行定量，使用实施例1中提供的试剂盒及实施例2和3中提供的建库测序方法及检测方法检测。表17显示了临床样本snv、fusion检测结果，其中阳性表示检出变异，阴性表示未检出。表18显示了临床样本cnv检测结果，其中，扩增阳性表示检出met存在cnv扩增。这些实验结果表明，利用本发明的方法，成功实现了ctdna占比极低的样本中的snv、cnv及fusion变异，尤其是低拷贝数扩增的cnv变异的有效检测。表17：突变类型ddpcrvaf变异支持数ngsvaf检测结果brafv600e0.18%80.12%阳性krasg12d0.22%150.25%阳性krasg13d0.20%80.13%阳性egfrl858r0.20%140.20%阳性egfrl858r0.17%130.19%阳性egfrt790m0.19%80.13%阳性egfrt790m0.27%70.11%阳性egfrt790m0.26%80.12%阳性egfr19del0.18%220.38%阳性egfr19del0.18%120.20%阳性eml4-alk0.20%260.26%阳性eml4-alk0.40%490.49%阳性eml4-alk0.80%1101.14%阳性表18：突变类型ddpcr拷贝数ngs拷贝数检测结果met扩增2.993.13扩增阳性met扩增2.482.79扩增阳性met扩增5.224.55扩增阳性met扩增3.173.26扩增阳性met扩增2.492.88扩增阳性met扩增4.954.66扩增阳性met扩增2.973.22扩增阳性met扩增2.412.83扩增阳性虽然上文已描述了本发明的各种实施方案，但是应理解的是，其仅以实例的方式提供，而并非限制。在不脱离本发明精神和范围的前提下，本发明还会有各种变化和改进，这些变化和改进都将落入要求保护的本发明范围内。本发明要求保护的范围由所附的权利要求书及其等效物界定。当前第1页12