bismark OT CTOT OB CTOB 以及mapping后的bam文件中的XG,XR列的含义

首先,OT,OB,CTOT,CTOB都是描述测序reads的,而不是描述参考基因组的。

bisul-fate建库会将DNA双链文库中非甲基化的C转化成U。转化结束后,被转化的U和互补链的G并不配对。此时正链(+,OT,original top strand)和反链(-,OB,original bottom strand)中,均为C to T转换的原始reads。

上述的reads经过PCR扩增后,正反单链均产生完全互补链。OT的互补链为CTOT(Complementary original top strand),OB的互补链为CTOB(Complementary bottom strand)。可知CTOT和CTOB均为GA转换。

图片引用自:

https://zhuanlan.zhihu.com/p/163495878

  • 正链(+) :是指 FASTA 文件中提供的原始序列本身,即参考基因组中记录的那条链。

  • 负链(-) :是指 FASTA 文件中记录的序列的反义链,也就是将参考序列取反向互补得到的链。

关于bismark比对,可以参考这篇:

https://www.zxzyl.com/archives/759/

比对生成的bam文件中,XR字段如果是CT,表示该reads是经过CT变换后匹配到了基因组中,即该reads属于OB或者OT;此时若XG字段为GA,表示它mapping到了参考基因组正链中的GA变换也就是反链中的CT变换,即表示它属于反链,即OB。如上图,即(2)对应的情况(reads上的C全部转换成T,然后mapping到了基因组正链GA转换。)

列个表供参考:

XR XG 说明 Strand
CT CT read 是 C→T(OB或者OT),基因组是 C→T(OT或CTOT) OT
CT GA read 是 C→T(OB或者OT),基因组是 G→A(OB或CTOB) OB
GA CT read 是 G→A(CTOT或者CTOB),基因组是 C→T(OT或CTOT) CTOT
GA GA read 是 G→A(CTOT或者CTOB),基因组是 G→A(OB或CTOB) CTOB
相关推荐
众智鸿图2 天前
XR驱动管网数据智能体检!众智鸿图精彩亮相2026燃气运营与安全研讨会
大数据·安全·xr
chaoyuanl3 天前
超元力玻璃剧场:裸眼3D、8K高清与悬空栈桥的成熟文旅解决方案
大数据·人工智能·3d·xr·vr
wasee0015 天前
亲子XR大空间为什么受欢迎?研学项目的客群与运营模式分析
xr
chaoyuanl21 天前
从平面到包裹:超元力飞行影院如何构建“人在景中、身随景动”的感知闭环
人工智能·科技·平面·xr·vr·娱乐
chaoyuanl2 个月前
超元力主题乐园飞行影院:把景区文化做成可持续更新的沉浸体验
大数据·3d·xr·娱乐·mr
XR技术研习社2 个月前
关于 PICO 串流测试中报错 IndexOutOfRangeException: renderPassIndex 的两个排障方法
unity·ar·xr·vr
Hody912 个月前
【XR硬件介绍】PICO Space Pro(Project Swan)深度解析:字节 2026 旗舰 MR 头显延期背后
xr·mr
Hody912 个月前
【XR硬件介绍】苹果N50智能眼镜技术解读:Vision Pro让位之后,“无屏AI眼镜”凭什么接棒
人工智能·xr
XR技术研习社2 个月前
从 PICO 文档看未来短期内的 Unity 版本选择
unity·游戏引擎·xr·vr
心前阳光3 个月前
Unity之XR Interaction Toolkit | IPointerDownHandler问题
unity·游戏引擎·xr