Skip to main content

S1.2 Contig graph reads alignment visualization

S1.2 节 “Contig graph reads alignment visualization” 讲的是：GFA 2 允许用 fragment lines 记录 reads 比对到 contigs 的情况，而 GfaViz 能够把这些 reads 到 contigs 的比对关系，放在组装图结构的背景下进行可视化。这样可以帮助用户理解 contig 的哪些区域有足够的 reads 支持，哪些区域支持度较低。

具体解释如下：

1. GFA 2 的 fragment lines 是什么？​

GFA 2 引入了一种叫 fragment line 的记录方式，用来描述 reads 如何比对到 contigs 上。这些比对信息反映了在组装过程中，哪些 reads 贡献了该 contig 的哪些区域。换句话说，它记录了 contig 的“覆盖度来源”。

2. 为什么要可视化这些比对？​

通过观察 reads 到 contigs 的比对，可以判断 contig 的某些区域是否被足够多的 reads 支持。如果某个区域只有很少的 reads 比对，说明该区域的组装可能不可靠。因此，可视化这些比对有助于评估组装质量。

3. 传统方法与 GfaViz 的区别​

  • 传统 pileup 可视化：通常是把 reads 比对到一个线性参考序列上，每个 contig 单独展示。也就是说，你只能看到 reads 相对于某一条 contig 的排列情况，看不到 contig 之间的图结构关系。
  • GfaViz 的可视化：GfaViz 允许在组装图结构的背景下展示比对区域。也就是说，你不仅能看到 reads 比对到哪个 contig，还能看到这些 contig 在图中如何相互连接。通过交互模式，用户还可以定位并识别比对到图中特定区域的 reads。

4. 视觉表示方式​

在 GfaViz 中，reads 到 contigs 的比对用以下方式表示：

  • 绿色线：表示来自**正向链（forward strand）**的 read。
  • 红色线：表示来自**反向链（reverse strand）**的 read。
  • 蓝色高亮：表示 contig 上被 read 比对覆盖的区域。
  • 透明度：与 S1.1 类似，如果多个 reads 比对到同一区域，蓝色会叠加得更深，颜色饱和度更高；如果只有一个 read 比对，颜色就较浅。这样可以直观看出多重覆盖。

5. Fig. S2 的例子​

图 S2 展示的是 GAGE S. aureus 数据集中，Abyss2 contig 图的一个连通分量。该分量包含 3 个 contigs：214、398 和 127。

  • 这 3 个 contigs 作为图中的节点。
  • 比对到这些 contigs 的 reads 用 GFA 2 的 fragment lines 表示。
  • 在图中，这些 reads 被画成绿色线（正向链）或红色线（反向链）。
  • contig 上被比对覆盖的区域用蓝色高亮显示。
  • 通过颜色饱和度可以看出哪些区域有多个 reads 比对。

6. 这一节想说明什么？​

S1.2 主要说明 GfaViz 能够：

  • 利用 GFA 2 的 fragment lines 记录 reads 到 contigs 的比对；
  • 在组装图结构的背景下可视化这些比对，而不仅仅是在线性参考上；
  • 通过绿/红线区分 reads 的方向；
  • 通过蓝色高亮和透明度展示比对区域和多重覆盖；
  • 支持交互式探索，帮助用户识别特定区域的 reads 支持情况，从而评估 contig 组装的可靠性。

简言之，S1.2 展示的是：GfaViz 可以把 reads 到 contigs 的比对画在组装图里，用绿/红线和蓝色高亮直观地显示哪些 reads 以什么方向支持了 contig 的哪些区域。