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%

Hisat2 Bowtie2比对结果解读相关推荐

【bioinfo】hisat2/bowtie2比对结果summary文件解读
hisat2/bowtie2比对后,比对结果基本信息统计的summary文件解读. 参考: https://www.cnblogs.com/leezx/p/8540862.html seqanswer ...
linux质控命令,RNA-seq摸索：2.sra下载数据→fastqc质控→hisat2/bowtie2/STAR/salmon比对→Samtools格式转换→IGV可视化结果...
我觉得这个作者这样的分类特别清晰,值得学 new_workflow/ │ └── annotation/ │ │ └── genome/ │ │ └── input/ │ │ └── output/ ...
2019/2/19_生信软件trim_galore和hisat2软件的结果解读
每日工作 2019/2/19 trim_galore软件的结果解读 trim_galore报告: SUMMARISING RUN PARAMETERS ======================== ...
conda安装bowtie2的报错:undefined symbol
使用conda安装bowtie2遇到undefined symbol报错 /public/home/xuzhougeng/miniconda3/envs/bsseq/bin/bowtie2-build ...
Linux做RNA-seq上游分析基本流程
本机处理器两核心两线程(多核心多线程在比对时占优),操作系统是Linux,发行版是deepin20.1社区版参考 http://www.bio-info-trainee.com/2218.html ...
c 语言头文件seqlsit,2020-10-24 RNAseq 从fq开始分析全流程
分析流程 image.png 1.上传四个样本原始就文件到服务器 Step 1. Analysing Sequence Quality with FastQC Description "Fa ...
RNA-seq | 转录组标准分析流程和常用工具软件介绍
笔记内容摘要:RNA-seq转录组基础知识与标准分析流程,简单记录学习过程. 转录组分析是对样本转录产物RNA的深入挖掘研究.通常情况下,植物的表型差异可能由许多因素控制,其中包括基因的转录环节,不同 ...
bwa、bowtie2、tophat、hisat2 比对软件学习中的笔记整理
对常用的比对软件学习进行用法整理记录.记录的内容相对简单,详细说明及用法还得参考软件使用说明书 bwa.bowtie2.tophat.hisatbwabwa(Burrows-Wheeler Align ...
RNA-seq流程学习笔记（7）-使用Hisat2进行序列比对
参考文章: RNAseq(4)–Hisat2进行序列比对及Samtools格式转化 RNA-seq(5):序列比对:Hisat2 hisat2比对软件将reads比对到参考基因组 hisat2比对 R ...

Hisat2 Bowtie2比对结果解读