【示例】CRISPR-Cas9 是怎么“剪”基因的

示例内容：这是一篇用于演示排版与功能的占位内容，可直接删除或替换为真实内容。

三个角色

• gRNA：约 20 nt 的向导序列，负责“找位置”；
• PAM：紧邻靶点的 NGG 序列，Cas9 只在它旁边动刀；
• Cas9：核酸酶，在 PAM 上游约 3 bp 处造成双链断裂。

一段找 PAM 的小代码

import re
seq = "ATGCTGACCGGTTAGCAGGCTTACGG"
for m in re.finditer(r"(?=([ACGT]{20})[ACGT]GG)", seq):
    print(m.start(), m.group(1))

编辑效率常用 $\text{indel\%}=\frac{\text{含插入/缺失的 reads}}{\text{总 reads}}\times 100\%$ 来衡量。