Hisat2 Bowtie2比对结果解读

article/2025/11/5 10:17:59

Bowtie的中文意思是:领结,蝴蝶结

在这里插入图片描述

Bowtie2用户手册:

http://bowtie-bio.sourceforge.net/bowtie2/manual.shtml

在看比对结果前需要了解三个概念:

1.Aligned concordantly 合理比对

在这里插入图片描述

主要和比对参数:–fr/–rf/–ff 有关,默认是:–fr
在参数–fr下:前两种比对就是合理比对,也即:concordantly,后两种比对会成为不合理比对,也即:aligned discordantly,如果两个reads距离太远也会成为不合理比对。

2.Multiple alignments 多重比对

定义:一对reads/一个pair比对上genome的多个位置,也即:aligned >1 times 或是 aligned concordantly >1 times,这对reads/这个pair就会成为多重比对。相对应的就是 unique比对,也即:aligned exactly 1 time 或是:aligned concordantly exactly 1 time。

这里 bowtie2的默认参数是:
Default mode: search for multiple alignments, report the best one, 也即是:多重比对,只取最好的那个alignment。

3.Bowtie2在双端比对下:

如果一对reads/一个pair没有比对上任何地方,也即 aligned concordantly 0 times,那么就将这个pair拆分成两个single-end reads来进行比对。

双端比对:

For datasets consisting of pairs, the summary might look like this:10000 reads; of these:10000 (100.00%) were paired; of these:650 (6.50%) aligned concordantly 0 times8823 (88.23%) aligned concordantly exactly 1 time527 (5.27%) aligned concordantly >1 times----650 pairs aligned concordantly 0 times; of these:34 (5.23%) aligned discordantly 1 time----616 pairs aligned 0 times concordantly or discordantly; of these:1232 mates make up the pairs; of these:660 (53.57%) aligned 0 times571 (46.35%) aligned exactly 1 time1 (0.08%) aligned >1 times96.70% overall alignment rate

翻译下即:

检测到10000个pairs,也即20000条reads;其中,10000个pairs都是paired也即双端数据;其中,650个pairs没有合理比对上;8823个pairs只合理比对上一次;527个pairs合理比对多于一次;在650个pairs没有合理比对上的pairs中:34个pairs可以不合理比对上一次;剩下的616个pairs,也即1232条reads中:660条reads没有比对上;571条reads可以比对上一次;1条read比对多于1次

总共的比对率是:(88232+5272+34*2+571+1)/20000=96.7%


http://chatgpt.dhexx.cn/article/B2R9NohX.shtml

相关文章

hisat2构建转录组索引的问题

hisat2构建转录组索引的问题 目前最新版为为hisat2. 安装过程如下 wget ftp://ftp.ccb.jhu.edu/pub/infphilo/hisat2/downloads/hisat2-2.1.0-Linux_x86_64.zip unzip hisat2-2.1.0-Linux_x86_64.zip 下载解压缩即可。 在进行比对前,首先需要对参考基因组建立索引…

hisat2比对结果图文解析

hisat2比对结果图文解析 看到hisat2输出的这个比对结果统计信息,是不是有点懵? 下面用图的方式进行简单剖析: 以上参考自https://blog.csdn.net/weixin_34191845/article/details/86453897 参考文献 Hisat2 bowtie2比对结果解读&#xff08…

RNA-seq——三、使用Hisat2进行序列比对

步骤 1. 下载对应的index2. 序列比对3. samtools:将sam文件转为bam文件4. 将bam文件载入IGV 为什么要比对:https://www.jianshu.com/p/681e02e7f9af Jimmy老师主要演示了四种比对工具,分别为hisat2、subjunc、bowtie2、bwa。除了subjunc能够直…

hisat2的index差别

下载的时候发下hisat2 主页中有多个index文件,一时间不解,搜索后发现如下评价。 目录 1.下载三个index:2.重命名为:3.hisat2比对命令:4.比对率:结论 1.下载三个index: 2.重命名为: …

Hisat2下载

功能: 将测序结果比对到参考基因组上 网站: http://ccb.jhu.edu/software/hisat2/index.shtml 安装:mkdir ~/biosoft && cd ~/biosoft wget ftp://ftp.ccb.jhu.edu/pub/infphilo/hisat2/downloads/hisat2-2.1.0-Linux_x86_64.zip …

安装Hisat2

一、(MobaXterm_Personal)安装aspera 首先进行预编译解压安装: mkdir Biosofts unzip hisat2-2.2.1-Linux_x86_64.zip -d ~/Biosofts/ cd ~/Biosofts ll ###zip文件,unzip解压,-d制指定目录 安装完成: 设置…

RNA-seq流程学习笔记(7)-使用Hisat2进行序列比对

参考文章: RNAseq(4)–Hisat2进行序列比对及Samtools格式转化 RNA-seq(5):序列比对:Hisat2 hisat2比对软件将reads比对到参考基因组 hisat2比对 RNA-seq数据分析使用方法(陈建国 译) 转录组分析2——比对基因组 RNA-seq练习 第二部…

RNA-seq分析htseq-count的使用

HTSeq作为一款可以处理高通量数据的python包,由Simon Anders, Paul Theodor Pyl, Wolfgang Huber等人携手推出HTSeq — A Python framework to work with high-throughput sequencing data。自发布以来就备受广大分析人员青睐,其提供了许多功能给那些熟悉python的大佬们去自信…

转录组学习之序列比对(Hisat2)[学习笔记通俗易懂版]

转录组学习之序列比对(hisat2)[学习笔记通俗易懂版] data :2023.7.25 recorder :CYH-BI 特别注意:本文为我自己学习的学习记录,没有任何权威,只能仅供初学者提供思路与参考。 本文知乎地址:https://zhua…

Hisat2安装及比对

Hisat2和STAR是目前转录组分析过程中用来做比对的两款主要工具,记得有一篇好像是2017年的文章专门比较了几款转录组比对工具对结果的影响,结论中认为两款软件在实际使用过程中对结果影响及耗时区别不大,我认为选一款就可以,之前总…

Hisat2 比对到参考基因组

比对的流程:建立索引→比对到参考基因组→SAM转BAM文件→BAM建立索引 1.准备参考基因组、建立索引 ## 参考基因组准备:注意参考基因组版本信息 # 下载,Ensembl:http://asia.ensembl.org/index.html # http://ftp.ensembl.org/pub/release-…

数据项组成数据元素,数据元素组成数据

数据元素:是组成数据的、有一定意义的基本单位。 数据项:一个数据可以由若干个数据项组成。数据项是数据不可分割的最小单位。 数据元素:字段、域、属性 数据项:元素、结点、顶点、记录 数据项组成数据元素,数据元…

数据结构 基本概念(数据项--数据元素--数据对象-数据类型-抽象数据类型)

//数据结构基本概念 #include<iostream> using namespace std;/* 数据 – 程序的操作对象&#xff0c;用于描述客观事物 数据的特点&#xff1a; 可以输入到计算机 可以被计算机程序处理 数据是一个抽象的概念&#xff0c;将其进行分类后得到程序设计语言中的类型。如&am…

根据结构体数组中某一数据项对结构体数组排序

/* *copyright(c) 2018,HH *All rights reserved. *作 者&#xff1a;HH *完成日期&#xff1a;2018年8月17日 *版本号&#xff1a;v1.0 * *问题描述:输入结构体数组&#xff0c;并根据结构体中的某一数据项对整个结构体数组进行排序 *输入描述&#xff1a;&#xff1b; *程序输…

点击echarts柱状图动态改变数据项颜色样式

首先附上参考文章连接&#xff1a;https://blog.csdn.net/weixin_42870683/article/details/103528254添加链接描述 今天来实现点击echarts柱状图&#xff0c;动态改变柱状图数据项颜色样式的案例。只要认真做&#xff0c;很容易学会~ 首先引入ECharts.js文件 <!-- 引入 …

数据结构考研:数据、数据元素、数据项、数据对象、数据结构的区别/详细解释(计算机/软件工程/王道论坛)

一、问题背景 博主最近在准备2020年春招复习数据结构这门功课时&#xff0c;采用了王道论坛的《2020年数据结构考研复习指导》这本书&#xff0c;该书的第一章节便是数据结构的基本概念&#xff1a;数据、数据元素、数据项、数据对象、数据结构。但该书对这几个容易混淆的概念…

数据,数据元素 数据项,数据对象的详细理解

1.数据(Data)&#xff1a;数据就是用户输入到计算机被计算机程序处理的一些符号&#xff0c;比如图片还有声音等.... 2.数据元素(Data Element)&#xff1a;是数据的基本单位&#xff0c;数据元素用于完整的描述一个对象&#xff0c;比如一个学生表&#xff0c;学生表也是由 数…

element el-transfer穿梭框的使用,自定义穿梭框的数据项

本文概述基本穿梭框效果自定义数据项实现右侧列表元素变化时触发方法 本文概述 查看element官方文档可以了解到 el-transfer穿梭框的基本使用方法&#xff0c;本文主要介绍数据项的自定义渲染和右侧列表元素变化时触发方法的实现。 自定义数据项&#xff0c;除了基本的el-tra…

数据结构之数据、数据元素、数据项、数据对象之间的关系

首图 # 前言 博主介绍&#xff1a; -- 本人是小奇&#xff0c;一个喜欢读书、交友、技术分享的Java程序员&#xff0c;更喜欢整理分享一些Java面试题&#xff0c;如果大家有需要可以关注我的公众号『小奇Java面试』领取资料。 文章目录 首图 前言 文章目录 关键字定义 1、数…

数据、数据元素、数据项、数据对象的理解

理解数据结构与算法时&#xff0c;觉得大多数的书上讲的定义特别抽象&#xff0c; 数据、数据元素、数据项、数据对象这四个概念特别难懂 本人理解如下&#xff1a; 假设有两张表&#xff0c;A表为人员表&#xff0c;B表为课程表&#xff0c;这两张表就是数据 、 而单独的一…