比对程序 Bowtie 的简单使用

主要简述 bowtie 的入门用法，完成从原始序列到最后mapping 之后的 sam格式文件。

---

参考文档

最重要的参考文档还是官方文档：

Bowtie-mannual

Getting started with Bowtie

github-bowtie

文献

---

下载，安装

下载：
bowtie-1.1.2
下载已经预编译好的软件即可。

安装：
先解压到一个文件夹

export PATH=$PATH:/home/user/bowtie/
如果需要永久设置，需要把上面一行加到 ./bashrc 文件中即可

测试：
在任何目录下：

$ bowtie
# 会出现各种信息，用法等等

---

以下是比较简单的例子：

建立索引

$ bowtie-build tigr7.fq /media/文档/BowtieIndex/tigr7
# bowtie-build : 是建立 索引
# tigr7.fa : 是基因组序列
# /media/文档/BowtieIndex/tigr7 : 指的是输出的索引文件 放在 /media/文档/BowtieIndex/ 目录下，并且所有index文件的前缀名为 tigr7,

mapping

把原始 fq 序列文件 mapping 到 参考基因组上，生成 sam 格式文件

/home/bowtie/bowtie -p 4 -x /media/文档/BowtieIndex/tigr7 -1 /media/文档/1.fq -2 /media/文档/2.fq -S /media/文档/eg2.sam
# -p 4 : 指的是用 4 个线程去运行程序
# -x : 指的是 The basename of the index for the reference genome，即后跟 参考基因组的位置
# -1 : /media/文档/1.fq 指的是 paired-end seq 中 的一端测序结果
# -2 : /media/文档/2.fq 指的是 paired-end seq 中 的另一端测序结果
# -S : /media/文档/eg2.sam 指的是 File to write SAM alignments to

---

这样， 就生成了所谓的 sam 格式文件， 以便于后续用 samtools 去处理。