CN106355045B - 一种基于扩增子二代测序小片段插入缺失检测的方法及装置 - Google Patents
一种基于扩增子二代测序小片段插入缺失检测的方法及装置 Download PDFInfo
- Publication number
- CN106355045B CN106355045B CN201610769361.1A CN201610769361A CN106355045B CN 106355045 B CN106355045 B CN 106355045B CN 201610769361 A CN201610769361 A CN 201610769361A CN 106355045 B CN106355045 B CN 106355045B
- Authority
- CN
- China
- Prior art keywords
- sample
- deletion
- base
- tested
- small fragment
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
Classifications
-
- G—PHYSICS
- G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
- G16B—BIOINFORMATICS, i.e. INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR GENETIC OR PROTEIN-RELATED DATA PROCESSING IN COMPUTATIONAL MOLECULAR BIOLOGY
- G16B25/00—ICT specially adapted for hybridisation; ICT specially adapted for gene or protein expression
Abstract
本发明公开了一种基于扩增子二代测序小片段插入缺失检测的方法及装置。该方法包括以下步骤:S1,分别提取待测样本和阴性对照样本的DNA,多重扩增子扩增与目标性状相关突变的目标区域;S2,通过二代测序得到目标区域的序列;S3,将目标区域的序列与参考基因组相比对,与参考基因组未比对上的碱基形成错配;S4,将由阴性对照样本得到的目标性状相关热点短片段插入缺失的背景噪音,以二项分布作为模型,在每个碱基上区分待测样本中热点短片段插入缺失与背景噪音,如果待测样本的错配碱基所占参考基因型比例与背景噪音有显著区别,那么将错配碱基作为小片段插入缺失阳性检出。应用本发明,提高了小片段插入缺失检测的准确性。
Description
技术领域
本发明涉及生物学领域,具体而言,涉及一种基于扩增子二代测序小片段插入缺失检测的方法及装置。
背景技术
多重扩增子二代测序是将感兴趣的基因组区域定制成特异性扩增引物与基因组DNA进行特异性扩增,将目标基因组区域的DNA片段进行富集后再利用第二代测序技术进行测序的研究策略。多重扩增子二代测序是目前基因组学研究中的一个热点技术,主要原因是该技术消耗少量的成本和时间。在相同成本下,研究者可以研究到更多的样本数量和测到更深的深度。作为一个强大、有效的技术,它在新一代高通量测序中发挥独特之处,应用领域越来越广泛。
赛默飞公司推出的Ion AmpliSeqTM是扩增子二代测序典型产品。通过扩增富集目标区域序列,然后使用proton进行二代测序得到序列具体情况,然后使用统计学检验对样本中低频率的热点短片段插入缺失与测序错误进行区分。但是这种方法检出效果的准确性还待进一步地提高。
发明内容
本发明旨在提供一种基于扩增子二代测序小片段插入缺失检测的方法及装置,以提高基于扩增子二代测序小片段插入缺失检测的准确性。
为了实现上述目的,根据本发明的一个方面,提供了一种基于扩增子二代测序小片段插入缺失检测的方法。该方法包括以下步骤:S1,分别提取待测样本和阴性对照样本的DNA,多重扩增子扩增与目标性状相关突变的目标区域,其中,阴性对照样本为目标性状相关基因未突变的样本;S2,通过二代测序得到目标区域的序列;S3,将目标区域的序列与参考基因组相比对,与参考基因组未比对上的碱基形成错配;S4,将由阴性对照样本得到的目标性状相关热点短片段插入缺失的背景噪音,以二项分布作为模型,在每个碱基上区分待测样本中热点短片段插入缺失与背景噪音,如果待测样本的错配碱基所占参考基因型比例与背景噪音有显著区别,那么将错配碱基作为小片段插入缺失阳性检出。
进一步的,S4中,如果待测样本的错配碱基所占参考基因型比例与背景噪音相比p-value<1e-6,那么将错配碱基作为小片段插入缺失阳性检出。
进一步的,S3中数据信息处理包括:S31,使用TMAP比对软件将目标区域的序列与参考基因组相比对,与参考基因组未比对上的碱基形成错配;S32,使用samtools软件建立比对文件的索引。
根据本发明的另一方面,提供了一种基于扩增子二代测序小片段插入缺失检测的装置。该装置包括:样本处理装置,用于分别提取待测样本和阴性对照样本的DNA,并多重扩增子扩增与目标性状相关突变的目标区域,其中,阴性对照样本为目标性状相关基因未突变的样本;测序装置,用于通过二代测序得到目标区域的序列;错配碱基获取装置,用于将目标区域的序列与参考基因组相比对,与参考基因组未比对上的碱基形成错配;小片段插入缺失检出装置,用于将由阴性对照样本得到的目标性状相关热点短片段插入缺失的背景噪音,以二项分布作为模型,在每个碱基上区分待测样本中热点短片段插入缺失与背景噪音,如果待测样本的错配碱基所占参考基因型比例与背景噪音有显著区别,那么将错配碱基作为小片段插入缺失阳性检出。
进一步的,小片段插入缺失检出装置中,如果待测样本的错配碱基所占参考基因型比例与背景噪音相比p-value<1e-6,那么将错配碱基作为小片段插入缺失阳性检出。
进一步的,错配碱基获取装置中包括:TMAP比对软件,用于将目标区域的序列与参考基因组相比对,与参考基因组未比对上的碱基形成错配;samtools软件,用于建立比对文件的索引。
应用本发明的技术方案,将由阴性对照样本得到的目标性状相关热点短片段插入缺失的背景噪音,以二项分布作为模型,在每个碱基上区分待测样本中热点短片段插入缺失与背景噪音,如果待测样本的错配碱基所占参考基因型比例与背景噪音有显著区别,那么将错配碱基作为小片段插入缺失阳性检出。在某些短片段插入缺失有特定位置和突变型,其中,较长的短片段插入缺失由于其在检测时由于测序或扩增原因造成假阳性检出的概率较低(<0.01%),所以对于不同位置不同长度的短片段插入缺失需要有不同检出阈值,阈值通过阴性对照样本在每个靶向位点的检出情况建立模型得到,对每一个突变以阴性样本中的发生率拟合二项分布作为模型,p-value为1e-6作为阈值,检测其余样本中的突变,这样可以最大限度地检测待测样本中的低频率(0.1%)热点短片段插入缺失,而不产生假阳性检出。
附图说明
构成本申请的一部分的说明书附图用来提供对本发明的进一步理解,本发明的示意性实施例及其说明用于解释本发明,并不构成对本发明的不当限定。在附图中:
图1示出了根据本发明一具体实施方式的基于扩增子二代测序小片段插入缺失检测的方法流程示意图。
具体实施方式
需要说明的是,在不冲突的情况下,本申请中的实施例及实施例中的特征可以相互组合。下面将参考附图并结合实施例来详细说明本发明。
在二代测序中,一次测序可以对几十个以上的目标性状相关热点短片段插入缺失进行检测。发明人发现,现有技术对所有位点采取同样阈值,而明显由于基因组上序列的差异,在扩增和测序的过程中不同长短的短片段插入缺失的背景噪音会有一定区别,采取同样阈值不能最大限区分序列原因造成的背景噪音与真实突变。
本发明使用阴性对照样本测序结果得到所有热点短片段插入缺失上的背景噪音分布,以区分在扩增和测序过程中产生的噪音和样本中真实突变,使待测样本中低频率突变的灵敏度和特异度得到提升。
根据本发明一种典型的实施方式,提供一种基于扩增子二代测序小片段插入缺失检测的方法。该方法包括以下步骤:S1,分别提取待测样本和阴性对照样本的DNA,多重扩增子扩增与目标性状相关突变的目标区域,其中,阴性对照样本为目标性状相关基因未突变的样本;S2,通过二代测序得到目标区域的序列;S3,将目标区域的序列与参考基因组相比对,与参考基因组未比对上的碱基形成错配;S4,将由阴性对照样本得到的目标性状相关热点短片段插入缺失的背景噪音,以二项分布作为模型,在每个碱基上区分待测样本中热点短片段插入缺失与背景噪音,如果待测样本的错配碱基所占参考基因型比例与背景噪音有显著区别,那么将错配碱基作为小片段插入缺失阳性检出。
优选的,S4中,如果待测样本的错配碱基所占参考基因型比例与背景噪音相比p-value<1e-6,那么将错配碱基作为小片段插入缺失阳性检出。
根据本发明一种典型的实施方式,S2中数据信息处理包括:S31,使用TMAP比对软件将目标区域的序列与参考基因组相比对,与参考基因组未比对上的碱基形成错配;S32,使用samtools软件建立比对文件的索引。
根据本发明一种典型的实施方式,提供一种基于扩增子二代测序小片段插入缺失检测的装置。该装置包括:样本处理装置,用于分别提取待测样本和阴性对照样本的DNA,并多重扩增子扩增与目标性状相关突变的目标区域,其中,阴性对照样本为目标性状相关基因未突变的样本;测序装置,用于通过二代测序得到目标区域的序列;错配碱基获取装置,用于将目标区域的序列与参考基因组相比对,与参考基因组未比对上的碱基形成错配;小片段插入缺失检出装置,用于将由阴性对照样本得到的目标性状相关热点短片段插入缺失的背景噪音,以二项分布作为模型,在每个碱基上区分待测样本中热点短片段插入缺失与背景噪音,如果待测样本的错配碱基所占参考基因型比例与背景噪音有显著区别,那么将错配碱基作为小片段插入缺失阳性检出。
优选的,小片段插入缺失检出装置中,如果待测样本的错配碱基所占参考基因型比例与背景噪音相比p-value<1e-6,那么将错配碱基作为小片段插入缺失阳性检出。
根据本发明一种典型的实施方式,错配碱基获取装置中包括:TMAP比对软件,用于将目标区域的序列与参考基因组相比对,与参考基因组未比对上的碱基形成错配;samtools软件,用于建立比对文件的索引。
根据本发明一种典型的实施方式,基于扩增子二代测序小片段插入缺失检测的方法如图1所示,主要包括小片段插入缺失检测程序外完成和小片段插入缺失检测程序内完成两大部分,其中,前一部分包括样本采集、样本处理测序及下机数据;后者包括阴性对照样本数据处理得到背景噪音,待测样本数据与背景噪音进行二相分布检测得到检测结果。
根据本发明一种典型的实施方式,基于扩增子二代测序拷贝数变异检测的方法包括以下步骤:1)样本处理:目标区域扩增、引物消化、末端修复、加接头、磁珠纯化、模板富集、测序;2)数据处理:使用TMAP比对软件将二代测序序列比对到参考基因组上,与参考基因组未比对上的碱基形成错配,以及使用samtools软件建立比对文件的索引;3)输入与目标性状相关的热点短片段插入缺失文件和使用阴性对照样本相同实验方法得到的每个碱基上背景噪音文件,与未比对上参考基因组碱基的错配进行统计学检验,如果背景噪音与错配碱基所占参考基因型比例有显著区别(p-value<1e-6),那么输出阳性检出,否则输出阴性结果。
下面将结合实施例进一步说明本发明的有益效果。本实施例中没有详细描述的步骤可采用本领域的常规技术手段实现。
实施例1
本实施例中,待测样本是NCI-H1650细胞系(美国模式培养物集存库购买),为一例细胞系样本,经过3d-pcr测试带有0.21%频率的EGFR:p.E746_A750delELREA阳性突变。
本实施例的具体步骤如下(所有试剂都是从赛默飞公司购买):
1)提取待测样本DNA,利用荧光定量计(Qubit)进行定量,进行浓缩或加无核酸酶的水稀释,使其浓度为5ng/ul,体积为6ul。
2)利用多重PCR技术对目标区域进行扩增,PCR反应体系如表1所示。
表1
试剂 | 体积 |
高保真多重扩增反应混和液 | 4μl |
引物混合液 | 10μl |
血浆游离DNA | 6μl |
其中,引物具体序列见表2。
表2
多重PCR反应条件如表3所示。
表3
3)利用尿嘧啶-DNA糖基化酶对带有尿嘧啶的引物进行消化,利用末端修复酶将DNA片段的末端补平。
消化及补平反应体系如表4所示。
表4
引物消化和末端修复条件如表5所示:
表5
4)将引物消化后的目标区域DNA用连接酶和测序接头进行连接,反应体系如表6:
表6
试剂 | 体积 |
1号接头溶液 | 0.5μl |
2号接头溶液 | 0.5μl |
连接酶 | 2μl |
连接反应混合液 | 4μl |
无核酸酶的水 | 1μl |
上步反应得到的DNA | 22ul |
其中,1号,2号接头序列见表7。
表7
*代表T碱基硫代修饰(phosphorothioate bond)。
连接条件如表8:
表8
5)将上一步得到的连接产物进行磁珠纯化,然后利用qPCR定量得到文库的摩尔浓度。
6)根据每个文库的浓度和目标数据量计算上机需要的体积,按计算结果吸取文库混合,然后采用乳液PCR的方法使模板在测序微珠上进行单克隆扩增,将微珠纯化后加入测序芯片。
7)测序,在proton基因测序仪上完成测序,测序平台将得到的电信号转化为碱基序列下机数据为bam文件存储所有测序结果。
8)使用TMAP将测序结果比对上参考基因组,平均测序深度为67043X,质控信息如表9所示。
表9
其中,目标区域序列为比对上目标区域的序列除以总序列,为有效数据率。
9)遍历每一个热点位点,得到支持突变的序列数和支持参考序列的序列数。与阴性样本对照集进行对比,如果有显著性差异(二项检验p<1e-6)将其提取出来。通常,本领域技术人员可以根据实际需要找到所关注的突变。本实施例中所关注的突变检出结果如表10所示。
表10
对比:本实施例的样本使用赛默飞公司的Torrent-Variant-Caller软件使用默认参数进行检测没有检出结果,这是因为赛默飞公司的软件对所有位点采用同一阈值,就需要使用最保守阈值以对噪音最大的位点进行过滤。
目前,多重扩增子测序数据使用TVC软件进行变异检测,可以检测数十个热点短片段插入缺失,对所有位点有相同阈值,并没有考虑到热点附近序列和插入缺失长度对测序造成的影响,这样组织样本中高频率(>1%)的热点短片段插入缺失并没有问题,但是对于血液中的接近背景噪音低频率(~0.1%)的热点短片段插入缺失会造成对于易错碱基的假阳性或不易测错的碱基的假阴性。另外短片段插入缺失的背景噪音错误率与长度也有一定关系,理论上越长的插入缺失可以与之混淆的背景噪音就越低(0.01%>)。通过这个方法也可以利用这部分信息进而检测低频突变。使用本发明的技术方案,利用阴性对照样本在所有热点短片段插入缺失的背景噪音信息,判断待测样本中的每个热点短片段插入缺失是否能与背景噪音区分,最大化检测血浆中低频突变的能力,同时不会产生假阳性结果。
以上所述仅为本发明的优选实施例而已,并不用于限制本发明,对于本领域的技术人员来说,本发明可以有各种更改和变化。凡在本发明的精神和原则之内,所作的任何修改、等同替换、改进等,均应包含在本发明的保护范围之内。
Claims (6)
1.一种基于扩增子二代测序小片段插入缺失检测的方法,其特征在于,包括以下步骤:
S1,分别提取待测样本和阴性对照样本的DNA,多重扩增子扩增与目标性状相关突变的目标区域,其中,阴性对照样本为目标性状相关基因未突变的样本;
S2,通过二代测序得到所述目标区域的序列;
S3,将所述目标区域的序列与参考基因组相比对,与所述参考基因组未比对上的碱基形成错配;
S4,将由所述阴性对照样本得到的目标性状相关热点短片段插入缺失的背景噪音,以二项分布作为模型,在每个碱基上区分所述待测样本中热点短片段插入缺失与所述背景噪音,如果所述待测样本的错配碱基所占参考基因型比例与所述背景噪音有显著区别,那么将所述错配碱基作为小片段插入缺失阳性检出。
2.根据权利要求1所述的方法,其特征在于,所述S4中,如果所述待测样本的错配碱基所占参考基因型比例与所述背景噪音相比p-value<1e-6,那么将所述错配碱基作为小片段插入缺失阳性检出。
3.根据权利要求1所述的方法,其特征在于,所述S3中数据信息处理包括:
S31,使用TMAP比对软件将所述目标区域的序列与所述参考基因组相比对,与所述参考基因组未比对上的碱基形成错配;
S32,使用samtools软件建立比对文件的索引。
4.一种基于扩增子二代测序小片段插入缺失检测的装置,其特征在于,包括:
样本处理装置,用于分别提取待测样本和阴性对照样本的DNA,并多重扩增子扩增与目标性状相关突变的目标区域,其中,阴性对照样本为目标性状相关基因未突变的样本;
测序装置,用于通过二代测序得到所述目标区域的序列;
错配碱基获取装置,用于将所述目标区域的序列与参考基因组相比对,与所述参考基因组未比对上的碱基形成错配;
小片段插入缺失检出装置,用于将由所述阴性对照样本得到的目标性状相关热点短片段插入缺失的背景噪音,以二项分布作为模型,在每个碱基上区分所述待测样本中热点短片段插入缺失与所述背景噪音,如果所述待测样本的错配碱基所占参考基因型比例与所述背景噪音有显著区别,那么将所述错配碱基作为小片段插入缺失阳性检出。
5.根据权利要求4所述的装置,其特征在于,所述小片段插入缺失检出装置中,如果所述待测样本的错配碱基所占参考基因型比例与所述背景噪音相比p-value<1e-6,那么将所述错配碱基作为小片段插入缺失阳性检出。
6.根据权利要求4所述的装置,其特征在于,所述错配碱基获取装置中包括:
TMAP比对软件,用于将所述目标区域的序列与所述参考基因组相比对,与所述参考基因组未比对上的碱基形成错配;
samtools软件,用于建立比对文件的索引。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610769361.1A CN106355045B (zh) | 2016-08-30 | 2016-08-30 | 一种基于扩增子二代测序小片段插入缺失检测的方法及装置 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610769361.1A CN106355045B (zh) | 2016-08-30 | 2016-08-30 | 一种基于扩增子二代测序小片段插入缺失检测的方法及装置 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN106355045A CN106355045A (zh) | 2017-01-25 |
CN106355045B true CN106355045B (zh) | 2019-03-15 |
Family
ID=57858287
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201610769361.1A Active CN106355045B (zh) | 2016-08-30 | 2016-08-30 | 一种基于扩增子二代测序小片段插入缺失检测的方法及装置 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN106355045B (zh) |
Families Citing this family (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108690871B (zh) * | 2018-03-29 | 2022-05-20 | 深圳裕策生物科技有限公司 | 基于二代测序的插入缺失突变检测方法、装置和存储介质 |
CN108470114B (zh) * | 2018-04-27 | 2020-02-28 | 元码基因科技(北京)股份有限公司 | 基于单样本的二代测序数据分析肿瘤突变负荷的方法 |
CN109266729B (zh) * | 2018-09-29 | 2020-11-27 | 中国科学院遗传与发育生物学研究所 | 一种基于基因组二代测序的大片段缺失的检测方法 |
CN109754843B (zh) * | 2018-12-04 | 2021-02-19 | 志诺维思(北京)基因科技有限公司 | 一种探测基因组小片段插入缺失的方法及装置 |
CN109767810B (zh) * | 2019-01-10 | 2021-04-27 | 上海思路迪生物医学科技有限公司 | 高通量测序数据分析方法及装置 |
CN111508559B (zh) * | 2020-04-21 | 2021-08-13 | 北京橡鑫生物科技有限公司 | 检测目标区域cnv的方法及装置 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104662168A (zh) * | 2012-06-21 | 2015-05-27 | 香港中文大学 | 用于癌症检测的血浆dna突变分析 |
CN104762402A (zh) * | 2015-04-21 | 2015-07-08 | 广州定康信息科技有限公司 | 超快速检测人类基因组单碱基突变和微插入缺失的方法 |
CN105518151A (zh) * | 2013-03-15 | 2016-04-20 | 莱兰斯坦福初级大学评议会 | 循环核酸肿瘤标志物的鉴别和用途 |
CN103874767B (zh) * | 2011-10-14 | 2016-08-17 | 深圳华大基因股份有限公司 | 对核酸样本中预定区域进行基因分型的方法和系统 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CA2507189C (en) * | 2002-11-27 | 2018-06-12 | Sequenom, Inc. | Fragmentation-based methods and systems for sequence variation detection and discovery |
-
2016
- 2016-08-30 CN CN201610769361.1A patent/CN106355045B/zh active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103874767B (zh) * | 2011-10-14 | 2016-08-17 | 深圳华大基因股份有限公司 | 对核酸样本中预定区域进行基因分型的方法和系统 |
CN104662168A (zh) * | 2012-06-21 | 2015-05-27 | 香港中文大学 | 用于癌症检测的血浆dna突变分析 |
CN105518151A (zh) * | 2013-03-15 | 2016-04-20 | 莱兰斯坦福初级大学评议会 | 循环核酸肿瘤标志物的鉴别和用途 |
CN104762402A (zh) * | 2015-04-21 | 2015-07-08 | 广州定康信息科技有限公司 | 超快速检测人类基因组单碱基突变和微插入缺失的方法 |
Also Published As
Publication number | Publication date |
---|---|
CN106355045A (zh) | 2017-01-25 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN106355045B (zh) | 一种基于扩增子二代测序小片段插入缺失检测的方法及装置 | |
CN106282356B (zh) | 一种基于扩增子二代测序点突变检测的方法及装置 | |
CN106372459B (zh) | 一种基于扩增子二代测序拷贝数变异检测的方法及装置 | |
JP6817259B2 (ja) | 癌の検出のための血漿dna中のサイズ及び数異常の使用 | |
CN107403074B (zh) | 一种突变蛋白的检测方法及装置 | |
EP2926288B1 (en) | Accurate and fast mapping of targeted sequencing reads | |
CN106462670A (zh) | 超深度测序中的罕见变体召集 | |
CN106497920A (zh) | 一种用于非小细胞肺癌基因突变检测的文库构建方法及试剂盒 | |
CN108588194A (zh) | 利用高通量测序数据检测肿瘤突变负荷的方法及装置 | |
CN105861710A (zh) | 测序接头、其制备方法及其在超低频变异检测中的应用 | |
CN112397151B (zh) | 基于靶向捕获测序的甲基化标志物筛选与评价方法及装置 | |
Babarinde et al. | Computational methods for mapping, assembly and quantification for coding and non-coding transcripts | |
CN113035273B (zh) | 一种快速、超高灵敏度的dna融合基因检测方法 | |
CN107058551A (zh) | 检测微卫星位点不稳定性的方法及装置 | |
CN108660200B (zh) | 一种检测短串联重复序列扩张的方法 | |
CN104846089A (zh) | 一种孕妇外周血中胎儿游离dna比例的定量方法 | |
CN111052249A (zh) | 确定预定染色体保守区域的方法、确定样本基因组中是否存在拷贝数变异的方法、系统和计算机可读介质 | |
CN113674803A (zh) | 一种拷贝数变异的检测方法及其应用 | |
CN115803447A (zh) | 染色体邻近实验中的结构变异检测 | |
CN109461473B (zh) | 胎儿游离dna浓度获取方法和装置 | |
CN108866173A (zh) | 一种标准序列的验证方法、装置及其应用 | |
CN109680054A (zh) | 一种低频dna突变的检测方法 | |
CN107885972B (zh) | 一种基于单端测序的融合基因检测方法及其应用 | |
CN112837748A (zh) | 一种区分不同解剖学起源肿瘤的系统及其方法 | |
WO2023124779A1 (zh) | 基于三代测序数据检测点突变的分析方法和装置 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |