[发明专利]一种二代序列的短序列纠错的方法和装置有效

专利信息
申请号: 201611138690.2 申请日: 2016-12-12
公开(公告)号: CN108614954B 公开(公告)日: 2020-07-28
发明(设计)人: 邓天全;高强;盛琴;陈世璇;杨林峰 申请(专利权)人: 深圳华大基因科技服务有限公司
主分类号: G16B40/20 分类号: G16B40/20;G16B30/00
代理公司: 深圳鼎合诚知识产权代理有限公司 44281 代理人: 孙银行;彭家恩
地址: 518083 广东省深圳市*** 国省代码: 广东;44
权利要求书: 查看更多 说明书: 查看更多
摘要:
搜索关键词: 一种 二代 序列 纠错 方法 装置
【权利要求书】:

1.一种二代序列的短序列纠错的方法,其特征在于,所述方法包括:

获取二代序列的K-mer序列;

对所述K-mer序列进行分类,得到多个分类文件;

对每个所述分类文件,统计所述K-mer序列的频数;

筛选并获得低于预定频数的K-mer序列;

用所述预定频数的K-mer序列对所述二代序列进行纠错;

其中,所述纠错包括:将所述预定频数的K-mer序列比对回所述二代序列,若比对上,则从该K-mer头尾两端对应在所述二代序列上的位置将所述二代序列截断;保留截断后的长度大于预定长度的序列,舍弃截断后的长度小于所述预定长度的序列。

2.根据权利要求1所述的二代序列的短序列纠错的方法,其特征在于,所述方法在所述分类和所述统计步骤之间,还包括:压缩所述K-mer序列。

3.根据权利要求2所述的二代序列的短序列纠错的方法,其特征在于,所述压缩所述K-mer序列包括:对所述K-mer序列中的每三个碱基用1位ASCII字符替代。

4.根据权利要求1-3任一项所述的二代序列的短序列纠错的方法,其特征在于,所述K-mer序列的长度范围是17至75个碱基长度。

5.根据权利要求4所述的二代序列的短序列纠错的方法,其特征在于,所述K-mer序列的长度是3的整数倍。

6.根据权利要求5所述的二代序列的短序列纠错的方法,其特征在于,所述K-mer序列的长度是39个碱基长度。

7.根据权利要求1-3任一项所述的二代序列的短序列纠错的方法,其特征在于,所述对所述K-mer序列进行分类包括:截取所述K-mer序列的连续n位碱基序列,其中含N的K-mer序列去掉,将所述K-mer序列分类到4的n次方份文件中。

8.根据权利要求1-3任一项所述的二代序列的短序列纠错的方法,其特征在于,所述预定频数是4。

9.根据权利要求1-3任一项所述的二代序列的短序列纠错的方法,其特征在于,所述预定长度是80-120个碱基长度。

10.根据权利要求9所述的二代序列的短序列纠错的方法,其特征在于,所述预定长度是100个碱基长度。

11.一种二代序列的短序列纠错的装置,其特征在于,所述装置包括:

获取单元,用于获取二代序列的K-mer序列;

分类单元,用于对所述K-mer序列进行分类,得到多个分类文件;

统计单元,用于对每个所述分类文件,统计所述K-mer序列的频数;

筛选单元,用于筛选并获得低于预定频数的K-mer序列;

纠错单元,用于用所述预定频数的K-mer序列对所述二代序列进行纠错;

其中,所述纠错包括:将所述预定频数的K-mer序列比对回所述二代序列,若比对上,则从该K-mer头尾两端对应在所述二代序列上的位置将所述二代序列截断;保留截断后的长度大于预定长度的序列,舍弃截断后的长度小于所述预定长度的序列。

12.根据权利要求11所述的二代序列的短序列纠错的装置,其特征在于,所述装置还包括:压缩单元,用于压缩分类过的K-mer序列。

下载完整专利技术内容需要扣除积分,VIP会员可以免费下载。

该专利技术资料仅供研究查看技术是否侵权等信息,商用须获得专利权人授权。该专利全部权利属于深圳华大基因科技服务有限公司,未经深圳华大基因科技服务有限公司许可,擅自商用是侵权行为。如果您想购买此专利、获得商业授权和技术合作,请联系【客服

本文链接:http://www.vipzhuanli.com/pat/books/201611138690.2/1.html,转载请声明来源钻瓜专利网。

同类专利
专利分类
×

专利文献下载

说明:

1、专利原文基于中国国家知识产权局专利说明书;

2、支持发明专利 、实用新型专利、外观设计专利(升级中);

3、专利数据每周两次同步更新,支持Adobe PDF格式;

4、内容包括专利技术的结构示意图流程工艺图技术构造图

5、已全新升级为极速版,下载速度显著提升!欢迎使用!

请您登陆后,进行下载,点击【登陆】 【注册】

关于我们 寻求报道 投稿须知 广告合作 版权声明 网站地图 友情链接 企业标识 联系我们

钻瓜专利网在线咨询

周一至周五 9:00-18:00

咨询在线客服咨询在线客服
tel code back_top