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
相关推荐
2501_946490386 天前
媒体演播场景xR虚拟拍摄技术落地实操——Hecoos服务器全链路技术拆解
服务器·xr·媒体·hirender·hecoos
UhMITKJvJkro12 天前
含UPFC电力系统潮流计算程序复现之路
xr
2501_9464903817 天前
xR+AR双虚拟技术在晚会舞台的协同落地实践——以浙江卫视苏宁易购超级秀为例
ar·xr·hirender·hecoos
2501_9464903817 天前
xR虚拟制片技术在游戏IP展演中的高密度拍摄落地实践——基于hecoos xR+UE4的技术实现
tcp/ip·游戏·xr
爱看科技18 天前
OpenAI或入局AI眼镜赛道,Meta/微美全息XR+AI战略凸显引领产业新风向
人工智能·xr
2501_9464903820 天前
xR 虚拟演播室赋能时政传播新范式 —— 基于 hecoos xR/VP 的全国两会报道技术实践
xr·媒体服务器·虚拟制作·hirender
2501_9464903821 天前
xR虚拟制作在文化类节目中的高海拔落地实践——以《念奴娇・赤壁怀古》《书香曲水》为例
xr·多媒体播放器·虚拟制作·hirender
2501_9464903823 天前
xR虚拟拍摄在媒体访谈节目中的技术落地实践——以上海报业集团xR创作室为例
xr
蛙色3DVR1 个月前
告别“一次性体验”:如何用社交与竞技撬动XR的持续吸引力?
xr·vr·vr大空间
北京阿法龙科技有限公司1 个月前
AI视觉赋能+AR眼镜,破解制造业产线组装错漏难题|阿法龙XR云平台
人工智能·ar·xr