[发明专利]高通量测序数据后期处理方法

摘要

本发明涉及一种高通量测序数据后期处理方法，属于分子生物学技术领域。该方法根据用户需要，对高通量测序数据后期进行高效的数据平滑、数据缩放、组间标准化、特定数据行提取、组间数据量平衡等操作，并支持任何测序平台生成的实验数据，特别适用于存放有多个样本组数据的矩阵文件，从而能够从而减少数据分析人员的工作负担，降低数据处理的难度，且本发明的高通量测序数据后期处理方法的应用范围也较为广泛。

主权项

一种高通量测序数据后期处理方法，其特征在于，所述的方法包括以下步骤：(1)系统获得高通量测序数据的矩阵文件，通过用户设置的矩阵最大值和最小值，生成经尺寸变换处理的新的矩阵文件；(2)系统对所述的新的矩阵文件的列或行进行分位数规范化处理，形成统一的列向量或行向量之间的值域范围；(3)系统对所述的新的矩阵文件进行压缩处理，保持所述的新的矩阵文件特征不变，并降低新的矩阵文件的容量大小；(4)系统根据用户操作获取关键字段文件、模版文件及模版文件中关键字段所属列，通过对模版文件中特定关键字段的信息提取，筛选出符合条件的模版文件子集并输出；(5)系统根据用户操作高通量数据序列比对结果文件，并根据设定的随机重抽样文件的个数和单个重抽样文件的序列比对结果记录条数，生成预设数量的模拟格式高通量数据序列比对结果文件；(6)系统根据所述的新的矩阵文件、模版文件子集和模拟格式高通量数据序列比对结果文件生成前景数据处理结果文件和背景数据处理文件，通过逐一比较每个信号值在前景数据处理结果文件和背景数据处理结果文件中的数量分布，生成经验性假阳性率分布，并对每个前景信号的假阳性概率进行估计，生成前景信号假阳性概率；(7)系统根据特定格式的高通量测序数据的矩阵文件及预设的预期聚合酶链反应冗余片段在总样本中所占比例，生成聚合酶链反应冗余片段阈值并在所述的高通量测序数据的矩阵文件中标记所述的聚合酶链反应冗余片段位置；所述的特定格式的高通量测序数据的矩阵文件中其第一列为染色体序号，第二列为位置信息，第三列为位点正链上的序列个数，第四列为位点负链上的序列个数。