[发明专利]对目标序列进行模拟的方法、装置和电子设备

权利要求书

1.一种对目标序列进行模拟的方法，其特征在于，包括：

获取参比序列和待拟合序列，且所述待拟合序列包含所述目标序列；

基于所述目标序列的碱基对应的所述参比序列的碱基的基准参数和所述待拟合序列的碱基的基准参数，获得所述目标序列的碱基的比较参数，其中，所述待拟合序列的基准参数基于至少一个所述待拟合序列的碱基的比较参数获得；

基于所述目标序列的比较参数确定所述目标序列的碱基是否合格以获得所述目标序列的合格碱基和不合格碱基；

对所述不合格碱基的所述比较参数进行迭代消减以获得目标参数，所述迭代消减的迭代次数基于所述比较参数，且每次迭代消减与至少一个碱基的至少一个比较参数相关联；以及

基于所述合格碱基和所述不合格碱基的所述目标参数确定所述目标序列的模拟后的序列。

2.根据权利要求1所述的方法，其特征在于，所述基准参数为碱基深度。

3.根据权利要求2所述的方法，其特征在于，对所述不合格碱基的所述比较参数进行迭代消减以获得目标参数包括：

对所述目标序列的所述碱基进行评分以获得碱基分数；

基于所述碱基分数对所述目标序列的所述碱基中的各个碱基进行排序以获得排序后的不合格碱基；以及

对所述排序后的不合格碱基的所述比较参数进行迭代消减以获得目标参数。

4.根据权利要求3所述的方法，其特征在于，基于所述碱基分数对所述目标序列的所述不合格碱基中的各个碱基进行排序以获得排序后的不合格碱基包括：

按照所述目标序列的所述碱基分数的由低到高的顺序对所述各个碱基进行排序以获得排序后的不合格碱基。

5.根据权利要求4所述的方法，其特征在于，对所述目标序列的所述排序后的所述不合格碱基的所述比较参数进行迭代消减以获得目标参数包括：

将所述不合格碱基的各读段进行评分以获得读段分数；

按照所述读段分数的由高到低的顺序对所述不合格碱基的各读段依次进行迭代删除，所述迭代删除在满足下列条件之一的情况下终止：

消减后的所述比较参数小于参数阈值；以及

所述读段覆盖的合格碱基数大于碱基阈值；基于所述不合格碱基的所述比较参数减去删除的读段的数目以获得所述目标参数。

6.根据权利要求5所述的方法，其特征在于，所述目标序列的所述不合格碱基的碱基分数为所述不合格碱基的所述比较参数与所述基准参数的商与一之和。

7.根据权利要求6所述的方法，其特征在于，所述读段分数为所述读段覆盖的所有碱基的碱基分数的平均值；以及

所述不合格碱基的各所述读段按照所述读段的读段分数的由高到低的顺序进行所述迭代删除。

8.根据权利要求3所述的方法，其特征在于，对所述目标序列的所述不合格碱基的所述比较参数进行迭代消减以获得目标参数包括：

在每次对所述排序后的不合格碱基的所述比较参数进行迭代消减之后，基于所述消减后的所述比较参数对所述不合格碱基重新计算碱基分数并排序，所述重新计算的碱基分数是计算前的碱基分数与所述迭代次数的差值。

9.根据权利要求8所述的方法，其特征在于，在所述消减后的所述比较参数小于用于确定所述碱基是否合格的偏差值的情况下，所述重新计算的碱基分数是计算前的碱基分数与所述迭代次数与修订系数之积的差值。

10.根据权利要求5所述的方法，其特征在于，所述碱基阈值是大于等于三且小于等于五的自然数。

11.根据权利要求1所述的方法，其特征在于，所述待拟合序列通过在所述目标序列的至少一侧延长预定单位而获得。