bismark和bsmap比对原理

news/2024/12/29 16:34:33/

1,bismark比对方法比较简单粗暴,它制作两种类型的基因组:1),将所有的C转化为T的基因组;2),将所有的G转化为A的基因组。它将bisulfate(重chong硫酸盐)处理后的测序reads也进行上述转化,这样也得到两种类型的reads。将两种类型的reads分别比对到两种基因组上,就得到4中比对结果,选择得分最高的比对结果作为最终结果。

示意图如下:

在这里插入图片描述

为什么这样做可以?

先看下下面的图:
在这里插入图片描述
重硫酸盐处理将未甲基化的C转化为U,PCR后U变成A,进而变成T。任意取一条比对后的reads,经过PCR后得到BSW,BSWR,BSC,BSCR四种reads,因为BSW和BSWR完全互补,BSC和BSCR完全互补,我们只需要取BSW和BSCR去查看就好了(如果BSW能比对上,那么BSWR一定也可以比对上)。

参考基因组往往是单链,我们选择Watson链作为参考基因组,它可以制作两种类型的参考基因组:
1,ATGTTTGTTTGAG
2,ACATTCACTTAAA

BSW也可以制作两种:
1,ATGTTTGTTTGAG
2,ACATTTATTTAAA

BSCR也可以制作两种:
1,ATGTTTATTTAAA
2,ACATTCACTTAAA

可以看到BSW的1和BSCR的2分别与参考基因组的1和2完美匹配。

可能导致错误的地方:1,SNP位点;2,PCR引入的错误;3,重硫酸盐导致的突变;4,测序错误。

2,bsmap容许有T比对到C或者T上,而C只能比对到C上,用容许错配的方法代替了上述制作多个基因组的问题。

具体比对过程可以通过下面两种方式实现(bsmap文章中未写明)
1,watson链和crick链都用来做参考基因组,这样只要容许有T比对到C或者T上,就可以实现所有的比对;
2,只用watson链做参考基因组,既需要容许有T比对到C或者T上,又需要容许A比对到A或G上。

在这里插入图片描述
在这里插入图片描述


http://www.ppmy.cn/news/54349.html

相关文章

Springboot Mybatis使用pageHelper实现分页查询

以下介绍实战中数据库框架使用的是mybatis,对整合mybatis此处不做介绍。 使用pageHelper实现分页查询其实非常简单,共两步: 一、导入依赖; 二、添加配置; 那么开始, 第一步: pom.xml添加依…

antDesignPro6: 如何设置环境变量,取值自动根据不同环境,动态修改(3步)。

官网文档:环境变量 - Ant Design Pro Pro 脚手架默认使用 Umi 作为底层框架,在 Umi 内可通过指定 UMI_ENV 环境变量来区分不同环境的配置文件,UMI_ENV 需要在 package.json 内配置。当 UMI_ENV 为 test 时,则必须在 config 目录下…

CSS3 grid网格布局

文章目录 CSS3 grid网格布局概述grid属性说明使用grid-template-rows & grid-template-columns 定义行高和列宽grid-auto-flow 定义项目的排列顺序grid-auto-rows & grid-auto-columns 定义多余网格的行高和列宽row-gap & column-gap 设置行间距和列间距gap 简写形…

ACM MM23 Workshop|多媒体+无人机

摘要: 无人驾驶飞行器 (UAV),也称为无人机,由于能够从空中捕获高质量的多媒体数据,近年来变得越来越流行。 随着航空摄影、电影摄影和测绘等多媒体应用的兴起,无人机已成为收集丰富多样的多媒体内容的强大工具。 本次…

机器视觉软件开发的多样性

当前四大主流的机器视觉 一.Halcon机器视觉软件 HALCON是德国MVtec公司开发的一套完善的标准的机器视觉算法包(商业用途),在欧洲以及日本的工业界已经是公认具有最佳效能的Machine Vision软件,HALCON支持Windows,Linux和Mac OS X操作环境;编程接口支持 c,c++,cpp.net,delphi…

vue webpack 打包优化方法

两种方法: 一、gzip压缩 前言 很多Vue项目,都有一个通病,随着代码越来越庞大,首页加载越来越慢,有时候出去演示,网络不好的时候,转半天,登录页面还没有出来。这主要是一些大的js文…

MySQL更新表的记录详解

文章目录 前言一、更新数据记录 1.特定数据记录2.所有数据记录总结 前言 更新数据记录是数据操作中常见的操作,可以更新表中已经存在数据记录中的值。在MySQL中可以通过UPDATE语句来实现更新数据记录,该SQL语句可以通过如下几种方式使用:更新…

HTML 状态消息

HTML 错误消息 当浏览器向 Web 服务器请求服务时,可能会发生错误,并且服务器可能会返回错误代码,例如 "404 Not Found"。 通常这些错误被称为 HTML 错误消息。 但是这些消息应称为 HTTP 状态消息。实际上,服务器总会为每个请求返回一条消息。最常见的消息是 2…