基于用户输入切换DNA序列比对函数中的Wobble变量

问题背景 #

需要实现含Wobble碱基的DNA序列比对并报告错配：

• DNA序列由A/T/C/G四种碱基组成
• Wobble字符（如R代表A/G、Y代表C/T等）可匹配多种碱基
• 当前代码需手动修改!N为!R/!Y等数组引用，无法通过用户输入动态切换；曾尝试合并所有数组用flat()处理，但因Wobble碱基重叠导致判断逻辑失效

当前代码的局限 #

现有代码硬编码了特定Wobble数组（如N），无法动态适配用户选择的Wobble类型，每次调整检测逻辑都要修改条件中的数组引用，灵活性差。

优化方案 #

使用Wobble映射对象存储字符与对应碱基数组的关联，通过用户输入的目标Wobble字符动态获取对应数组，实现逻辑的动态切换：

1. 定义映射对象，将每个Wobble字符映射到其对应的碱基集合
2. 接收用户输入的目标Wobble字符（如'N'/'Y'/'R'）
3. 在比对逻辑中动态调用对应数组，判断待测碱基是否属于该Wobble的允许范围

优化后代码示例 #

// 定义Wobble字符到对应碱基的映射
const wobbleMap = {
  R: ['A', 'G'],
  Y: ['C', 'T'],
  K: ['G', 'T'],
  N: ['A', 'G', 'C', 'T'],
  M: ['A', 'C'],
  S: ['C', 'G'],
  W: ['A', 'T']
};
const allWobbles = Object.keys(wobbleMap);

// 用户输入目标Wobble（示例：用户输入'N'或'Y'等）
const targetWobble = 'N'; // 可替换为用户输入值
const targetBases = wobbleMap[targetWobble];

let mismatchPositions = [];
const referenceLength = referenceSeq.length;

for (let i = 0; i < referenceLength; i++) {
  const refBase = referenceSeq[i];
  const testBase = seq.sequence[i];

  if (allWobbles.includes(refBase)) {
    // 参考碱基是Wobble：判断待测碱基是否不在当前目标Wobble的允许集合中
    if (!targetBases.includes(testBase) && refBase !== testBase) {
      mismatchPositions.push(`(${i+1} ${refBase} > ${testBase} Mismatch)`);
    }
  } else {
    // 参考碱基是普通碱基：直接比对
    if (refBase !== testBase) {
      mismatchPositions.push(`(${i+1} ${refBase} > ${testBase} Mismatch)`);
    }
  }
}

// 输出结果示例
if (mismatchPositions.length > 0) {
  console.log(`${seq.sequence} - ${mismatchPositions.join(' ')}`);
} else {
  console.log(`${seq.sequence} - No mismatches`);
}

测试效果 #

以参考序列GTAAGTTGNACG为例：

• 当targetWobble = 'N'时，第9位的N对应的错配不会被检测到，输出如：

GTAAGTTGAATG -(11 C > T Mismatch)
GTAAGTTAAACG -(8 G > A Mismatch)

• 当targetWobble = 'Y'时，若参考序列第9位为Y，待测碱基不在Y的允许集合（C/T）中时，会报告该位置错配。

内容的提问来源于stack exchange，提问作者Dmishra