[发明专利]一种基于参考序列的基因压缩方法

权利要求书

1.一种基于参考序列的基因压缩方法，其特征在于，包括以下步骤：

步骤1)、遍历参考基因文件；将其中小写字符写入内存，此小写字符信息包括碱基、N字符和其他字符在内的所有小写字符信息；遍历过程中遇到大写的碱基字符时，将其写入内存；遇到小写的碱基字符时，将其转成大写字符，然后写入内存，并且使用开始位置与长度二元组记录参考基因文件小写字符的位置信息，得到参考基因文件小写字符二元组；

步骤2)、参考基因文件碱基信息Hash值计算；使用k-mer方法将基因分段，对碱基ACGT分别编码，并计算Hash值，使得每一个Hash值代表一个基因片段；

使用k-mer方法将基因分段，对碱基ACGT分别编码，并计算Hash值，使得每一个Hash值代表一个基因片段的方法：

首先将A、C、G、T信息编码，分别使用0、1、2、3表示；其次从首位开始，选取连续的11位字符，将第一位的数字与4的10次方相乘，将第二位的数字与4的9次方相乘，以此类推，直至最后一位数字与4的0次方相乘，然后将各个结果相加得到最终的结果即为Hash值，并存入Hash表中；

Hash值计算方法：

对于一段字符x₁₀x₉x₈x₇x₆x₅x₄x₃x₂x₁x₀，其中x_i对应的编码为code(x_i)，则计算Hash值的方式为

其中，HashValue表示Hash值；

步骤3)、遍历待压缩基因文件；将其中小写字符写入内存，此小写字符信息包括碱基、N字符和其他字符在内的所有小写字符信息；遍历过程中遇到大写的碱基字符、N字符和其他字符时，将其写入内存；遇到小写的碱基字符、N字符和其他字符时，将其转成大写字符，然后写入内存；同时将换行信息也写入内存；并且使用开始位置与长度二元组记录压缩基因文件小写字符的位置信息，得到压缩基因文件小写字符二元组；

步骤4)、小写字符写入压缩文件；对比已经生成的小写字符二元组，将待压缩基因文件的小写字符二元组对比参考基因文件的小写字符二元组，进一步将相同的部分用二元组表示，并写入压缩文件；不同的部分直接写入压缩文件；

步骤5)、换行信息、N字符和其他字符写入压缩文件；将带压缩文件中换行信息、N字符和其他字符以二元组形式表示并将二元组写入压缩文件；

步骤6)、碱基信息写入压缩文件；使用k-mer方法对带压缩基因文件每一段求取Hash值，并根据Hash值匹配得到与参考基因文件相同的k-mer段，记录该k-mer段的开始位置与连续重复最大长度为一个新二元组，直至遍历待压缩基因文件小写字符二元组完成，生成新二元组写入压缩文件。