当前位置: 首页 > news >正文

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转换。)

列个表供参考:

XRXG说明Strand
CTCTread 是 C→T(OB或者OT),基因组是 C→T(OT或CTOT)OT
CTGAread 是 C→T(OB或者OT),基因组是 G→A(OB或CTOB)OB
GACTread 是 G→A(CTOT或者CTOB),基因组是 C→T(OT或CTOT)CTOT
GAGAread 是 G→A(CTOT或者CTOB),基因组是 G→A(OB或CTOB)CTOB

相关文章:

  • Python基础:人生重开模拟器(小游戏)
  • 2024年数维杯国际大学生数学建模挑战赛C题时间信号脉冲定时噪声抑制与大气时延抑制模型解题全过程论文及程序
  • CUDA与OpenGL混合编程图形渲染
  • vLLM实战部署embedding、reranker、senseVoice、qwen2.5、qwen3模型
  • [蓝桥杯]约瑟夫环
  • Qt/C++编写GB28181服务端工具/绿色版开箱即用/对标wvp-gb28181/实时画面预览/录像回放下载
  • 5.29 自学测试 Linux基础 Day4
  • 基于大模型的急性乳腺炎全病程风险预测与综合治疗方案研究
  • 快乐大冒险:解锁身体里的 “快乐密码”
  • 海外tk抓包简单暴力方式
  • 如何评估 RAG 的分块Chunking策略
  • node-sass 报错
  • 小云天气APP:精准预报,贴心服务
  • 《P2960 [USACO09OCT] Invasion of the Milkweed G》
  • 悲观锁与乐观锁:并发编程中的两种核心控制策略详解
  • 洋流的实时原位预测
  • 大数据-275 Spark MLib - 基础介绍 机器学习算法 集成学习 随机森林 Bagging Boosting
  • Git实战--基于已有分支克隆进行项目开发的完整流程
  • Chapter 10 Inductive DC–DC Converters
  • 性能测试实例(http和ldap协议压测)
  • 网站备案的核验单/网络营销带来的效果
  • wordpress 网站 图标/我赢网seo优化网站
  • seo网站编辑/站长工具 站长之家
  • 动态网站设计选题参考/西安网站设计
  • 360做网站和推广怎么样/疫情排行榜最新消息
  • 中国最好的网站器域名统一/网络营销与直播电商怎么样