当前位置: 首页 > 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
http://www.lqws.cn/news/84709.html

相关文章:

  • 调试的本质:从混沌走向秩序
  • 双指针题解——反转字符串中的单词【LeetCode】
  • FastAPI安全认证:从密码到令牌的魔法之旅
  • 嵌入式Linux 期末复习指南(下)
  • Mysql水平分表(基于Mycat)及常用分片规则
  • Spring 5 响应式编程:构建高性能全栈应用的关键
  • hooks组件-useState
  • 吴恩达机器学习笔记(1)—引言
  • 设计模式——访问者设计模式(行为型)
  • wow Warlock shushia [Dreadsteed]
  • 地图 APP 和购物 APP 是最急切上 AI的地方
  • Artificial Analysis2025年Q1人工智能发展六大趋势总结
  • ThreadLocal ,底层原理,强引用,弱引用,内存泄漏
  • Vue3(watch,watchEffect,标签中ref的使用,TS,props,生命周期)
  • FastAPI+Pyomo实现线性回归解决饮食问题
  • 函数调用的机器级实现(二):栈帧的访问与切换机制
  • 极客时间:用 FAISS、LangChain 和 Google Colab 模拟 LLM 的短期与长期记忆
  • 【springcloud】快速搭建一套分布式服务springcloudalibaba(四)
  • python爬虫:Ruia的详细使用(一个基于asyncio和aiohttp的异步爬虫框架)
  • Langchian - 自定义提示词模板 提取结构化的数据
  • 【redis实战篇】第七天
  • 在 Linux 服务器上无需 sudo 权限解压/打包 .7z 的方法(实用命令)
  • 小团队如何落地 Scrum 模型:从 0 到 1 的实战指南
  • rabbitmq Direct交换机简介
  • C++——AVL平衡树
  • Java递归编程中的StackOverflowError问题分析与解决方案
  • 题目 3230: 蓝桥杯2024年第十五届省赛真题-星际旅行
  • 数字孪生智慧水利解决方案:数字化场景、智慧化模拟、精准化决策,构建数字孪生流域为核心的智慧水利体系
  • 【笔记】Windows 部署 Suna 开源项目完整流程记录
  • 前端面试宝典---前端水印