计算机应用研究

作者：胡爽,陈长波,

单位：1.中国科学院重庆绿色智能技术研究院生物计算安全重庆市重点实验院,重庆400714；2.重庆邮电大学计算机科学与技术学院,重庆400065；3.中国科学院大学重庆学院,重庆400714；

关键词：序列比对,罗伯斯-惠勒变换,第二代测序,BWA,古代DNA,

基金：国家重点研发计划资助项目(2020YFA0712300)；重庆英才计划青年拔尖项目(2021000263)；重庆市院士牵头科技创新引导专项资助项目(cstc2021yszx-jcyjX0004,2022YSZX-JCX0011CSTB,CSTB2023YSZX-JCX0008)；；

序列比对是基因组数据分析的关键一环，提高其比对效率对推动测序技术在医学、古生物学等领域的应用具有重要意义。针对BWA算法两大步骤之一的SAMSE存在冗余读取索引导致效率欠佳的问题，提出了新的算法BWA*。其通过运用流程优化消除了SAMSE中参考序列及其索引的冗余读取，在此基础上叠加运用关键参数值调整和多线程优化等技术，进一步提高了计算效率。公开数据库中的真实序列测试表明，BWA*的SAMSE步骤的性能是BWA中SAMSE性能的7.11~8.61倍，平均为7.84倍，BWA*的整体性能是BWA的1.25~1.70倍，平均1.47倍。针对实际应用中的古代DNA序列比对，实验表明和另一常用工具BWA-MEM相比，优化后的BWA*在继承原有BWA高精度特性的同时，实现了对BWA-MEM速度的超越。

来源：2024年第12期

《计算机应用研究》期刊编辑部

查看计算机应用研究杂志2024年第12期