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
相关推荐
恒点虚拟仿真13 小时前
XR数字融合工作站赋能新能源汽车专业建设的创新路径
人工智能·汽车·xr·虚拟现实·虚拟仿真·新能源汽车·ai+虚拟仿真
Hody914 天前
【XR硬件系列】Vivo Vision 与 Apple VisionPro 深度技术对比:MR 时代的轻量化革命与生态霸权
xr·mr
平行云4 天前
Paraverse平行云实时云渲染助力第82届威尼斯电影节XR沉浸式体验
unity·云原生·ue5·xr·实时云渲染
爱看科技4 天前
苹果Vision Air蓝图或定档2027,三星/微美全息加速XR+AI核心生态布局卡位
人工智能·xr
北京阿法龙科技有限公司4 天前
AR眼镜在智能制造的应用方向和场景用例|阿法龙XR云平台
ar·xr·制造
Hody911 个月前
梳理一下实现3D显示的途径有哪些?
xr
安卓开发者1 个月前
Jetpack Compose for XR:构建下一代空间UI的完整指南
ui·xr
沧海归城1 个月前
Unity_XR控制手部动画
unity·游戏引擎·xr
道可云2 个月前
道可云人工智能每日资讯|天津市人工智能(AI+信创)创新生态联盟成立
人工智能·百度·ar·xr·deepseek
X-mj2 个月前
Unity URP + XR 自定义 Skybox 在真机变黑问题全解析与解决方案(支持 Pico、Quest 等一体机)
unity·游戏引擎·xr