[发明专利]一种二代序列的短序列纠错的方法和装置有效
申请号: | 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序列。
该专利技术资料仅供研究查看技术是否侵权等信息,商用须获得专利权人授权。该专利全部权利属于深圳华大基因科技服务有限公司,未经深圳华大基因科技服务有限公司许可,擅自商用是侵权行为。如果您想购买此专利、获得商业授权和技术合作,请联系【客服】
本文链接:http://www.vipzhuanli.com/pat/books/201611138690.2/1.html,转载请声明来源钻瓜专利网。