JavaScript可用于DNA序列分析，支持合法性检查、碱基统计、互补链获取及转录操作；通过findORFs函数识别开放阅读框；实现序列比对与相似度计算。
JavaScript 可以用于前端或 Node.js 环境下的轻量级 DNA 序列分析，尤其适合在网页应用中快速处理和可视化生物信息数据。虽然主流生物信息学多使用 Python 或 R，但通过 JavaScript 也能实现基础的 DNA 序列操作与分析。
DNA序列基本操作
常见的 DNA 碱基为 A（腺嘌呤）、T（胸腺嘧啶）、C（胞嘧啶）、G（鸟嘌呤）。JavaScript 可以对字符串形式的 DNA 序列进行处理：
• 检查序列合法性：确保只包含 ATCG 字符• 计算碱基组成：统计各碱基出现频率
• 获取互补链：A↔T，C↔G，反向互补
• 转录为 RNA：将 T 替换为 U
示例代码：
<font>
function validateDNA(seq) {
return /^[ATCG]+$/i.test(seq);
}
<p>function countBases(seq) {
const counts = { A: 0, T: 0, C: 0, G: 0 };
for (let base of seq.toUpperCase()) {
if (counts.hasOwnProperty(base)) counts[base]++;
}
return counts;
}</p><p>function complementStrand(seq) {
const compMap = { 'A': 'T', 'T': 'A', 'C': 'G', 'G': 'C' };
return seq.toUpperCase().split('').reverse().map(b => compMap[b]).join('');
}
</font>开放阅读框（ORF）查找
在 DNA 序列中识别可能编码蛋白质的区域，通常从起始密码子 ATG 开始，到终止密码子（TAA、TAG、TGA）结束。
立即学习“Java免费学习笔记（深入）”；
• 将 DNA 转为 RNA（T → U）• 在三个正向读码框中搜索起始与终止密码子
• 支持反向互补链分析
简化版 ORF 检测逻辑：
<font>
function findORFs(dnaSeq) {
const startCodon = 'AUG';
const stopCodons = ['UAA', 'UAG', 'UGA'];
const rna = dnaSeq.replace(/T/g, 'U');
const orfs = [];
<p>for (let frame = 0; frame < 3; frame++) {
for (let i = frame; i < rna.length - 2; i += 3) {
const codon = rna.slice(i, i + 3);
if (codon === startCodon) {
for (let j = i + 3; j < rna.length - 2; j += 3) {
const stop = rna.slice(j, j + 3);
if (stopCodons.includes(stop)) {
orfs.push(rna.slice(i, j + 3));
break;
}
}
}
}
}
return orfs;
}
</font>序列比对简易实现
对于短序列，可使用 JavaScript 实现简单的 Needleman-Wunsch 或滑动比对。
• 计算两序列的匹配率• 支持模糊匹配与错配计数
示例：计算相似度
<font>
function sequenceIdentity(seq1, seq2) {
if (seq1.length !== seq2.length) {
const minLen = Math.min(seq1.length, seq2.length);
seq1 = seq1.slice(0, minLen);
seq2 = seq2.slice(0, minLen);
}
let matches = 0;
for (let i = 0; i < seq1.length; i++) {
if (seq1[i].toUpperCase() === seq2[i].toUpperCase()) matches++;
}
return (matches / seq1.length) * 100;
}
</font>前端可视化集成
结合 HTML5 和图表库（如 Chart.js 或 D3.js），可在网页中展示碱基分布、GC 含量、ORF 位置等。
• 用柱状图显示碱基频率• 用进度条样式展示 GC 含量
• 高亮显示 ORF 区域
例如，计算 GC 含量：
<font>
function gcContent(seq) {
const gcCount = (seq.match(/[GC]/gi) || []).length;
return (gcCount / seq.length) * 100;
}
</font>基本上就这些。JavaScript 适合做交互式 DNA 分析工具原型，特别是在浏览器中运行的小型项目。复杂分析建议结合 Node.js 和专用库（如 biojs），或调用后端服务处理。不复杂但容易忽略的是大小写处理和边界检查。
