生信小白——初学 bowtie2

本文就以生信小白的角度,简述bowtie2的下载、安装及使用。关于bowtie2的使用教程非常丰富,具体的概念本文不再赘述,只对操作过程进行简(xi)简(zhi)单(ru)单(wei)的介绍,希望对同样的初学者提供一些帮助。

我是使用Mac终端连接服务器环境下使用的。


下载及安装

起初我是在服务器环境下使用 wget+链接 下载软件,运行时发现安装失败。

就是这样,安装失败

于是在服务器里下载miniconda,并添加了bioconda镜像。(conda安装教程https://www.jianshu.com/p/edaa744ea47d)。conda环境下输入conda iinstall gatk bowtie2。检索结果输入y。

检索后输入“y”

在miniconda/bin目录下输入ls,可以看到bowtie2,输入 ./bowtie2 -help,出现help说明软件就可以用了。

安装成功

使用

序列数据

我已经把需要比对的序列和参考基因组上传服务器了。read 目录下是我的测序文件a_0.fq、 a_20.fq 和 a_60.fq,ref 目录下是我的参考基因组 Tas.V1.genome.fa

read目录下是我测序获得的三条序列,ref目录下是我的参考基因组

建索引

进入ref目录。输入bowtie2-build Tsa.v1.genome.fa(参考基因组) Tas(输出文件名前缀)

建索引

需要一段时间,我170M的基因组大概1分钟多时间。完成后会获得6个 .bt2后缀的文件,如下,说明索引建立成功。

获得索引文件

序列比对

这里的参数较多,选取主要的使用即可。通常测序文件是双端的,这里我的是序列是单端,所以使用的参数也不同。首先是双端序列(我没跑双端,结果就不呈现了),单端序列结果如下。

进入read目录下输入:

双端:$ bowtie2 -x ../ref/Tsa(索引所在目录,因为在上一级,所以使用../) -1 read1_1.fq(第一端) -2 read1_2.fq(第二端) -S read1.sam(指定输出文件名)

双端序列

单端:$ bowtie2 -x ../ref/Tsa(索引所在目录) -U a_0.fq(第一端) -S a_0_Tsa.sam(指定输出文件名)

单端序列

初步结果显示,a_0.fp 文件中有172333条序列,其中只有8.79%的序列比对到参考基因组上。

比对结果

以此类推,我获得了其余两个序列的sam文件。

获得三个序列的sam文件


参考及其感谢

hsy_hzauer  https://www.jianshu.com/u/49bcf3035565

卖萌哥    conda的安装与使用(最新版):https://www.jianshu.com/p/edaa744ea47d

郭威    序列比对软件bowtie2--快速入门:https://v.qq.com/x/page/m0846mwp8cl.html

最后编辑于
©著作权归作者所有,转载或内容合作请联系作者
  • 序言:七十年代末,一起剥皮案震惊了整个滨河市,随后出现的几起案子,更是在滨河造成了极大的恐慌,老刑警刘岩,带你破解...
    沈念sama阅读 217,657评论 6 505
  • 序言:滨河连续发生了三起死亡事件,死亡现场离奇诡异,居然都是意外死亡,警方通过查阅死者的电脑和手机,发现死者居然都...
    沈念sama阅读 92,889评论 3 394
  • 文/潘晓璐 我一进店门,熙熙楼的掌柜王于贵愁眉苦脸地迎上来,“玉大人,你说我怎么就摊上这事。” “怎么了?”我有些...
    开封第一讲书人阅读 164,057评论 0 354
  • 文/不坏的土叔 我叫张陵,是天一观的道长。 经常有香客问我,道长,这世上最难降的妖魔是什么? 我笑而不...
    开封第一讲书人阅读 58,509评论 1 293
  • 正文 为了忘掉前任,我火速办了婚礼,结果婚礼上,老公的妹妹穿的比我还像新娘。我一直安慰自己,他们只是感情好,可当我...
    茶点故事阅读 67,562评论 6 392
  • 文/花漫 我一把揭开白布。 她就那样静静地躺着,像睡着了一般。 火红的嫁衣衬着肌肤如雪。 梳的纹丝不乱的头发上,一...
    开封第一讲书人阅读 51,443评论 1 302
  • 那天,我揣着相机与录音,去河边找鬼。 笑死,一个胖子当着我的面吹牛,可吹牛的内容都是我干的。 我是一名探鬼主播,决...
    沈念sama阅读 40,251评论 3 418
  • 文/苍兰香墨 我猛地睁开眼,长吁一口气:“原来是场噩梦啊……” “哼!你这毒妇竟也来了?” 一声冷哼从身侧响起,我...
    开封第一讲书人阅读 39,129评论 0 276
  • 序言:老挝万荣一对情侣失踪,失踪者是张志新(化名)和其女友刘颖,没想到半个月后,有当地人在树林里发现了一具尸体,经...
    沈念sama阅读 45,561评论 1 314
  • 正文 独居荒郊野岭守林人离奇死亡,尸身上长有42处带血的脓包…… 初始之章·张勋 以下内容为张勋视角 年9月15日...
    茶点故事阅读 37,779评论 3 335
  • 正文 我和宋清朗相恋三年,在试婚纱的时候发现自己被绿了。 大学时的朋友给我发了我未婚夫和他白月光在一起吃饭的照片。...
    茶点故事阅读 39,902评论 1 348
  • 序言:一个原本活蹦乱跳的男人离奇死亡,死状恐怖,灵堂内的尸体忽然破棺而出,到底是诈尸还是另有隐情,我是刑警宁泽,带...
    沈念sama阅读 35,621评论 5 345
  • 正文 年R本政府宣布,位于F岛的核电站,受9级特大地震影响,放射性物质发生泄漏。R本人自食恶果不足惜,却给世界环境...
    茶点故事阅读 41,220评论 3 328
  • 文/蒙蒙 一、第九天 我趴在偏房一处隐蔽的房顶上张望。 院中可真热闹,春花似锦、人声如沸。这庄子的主人今日做“春日...
    开封第一讲书人阅读 31,838评论 0 22
  • 文/苍兰香墨 我抬头看了看天上的太阳。三九已至,却和暖如春,着一层夹袄步出监牢的瞬间,已是汗流浃背。 一阵脚步声响...
    开封第一讲书人阅读 32,971评论 1 269
  • 我被黑心中介骗来泰国打工, 没想到刚下飞机就差点儿被人妖公主榨干…… 1. 我叫王不留,地道东北人。 一个月前我还...
    沈念sama阅读 48,025评论 2 370
  • 正文 我出身青楼,却偏偏与公主长得像,于是被迫代替她去往敌国和亲。 传闻我的和亲对象是个残疾皇子,可洞房花烛夜当晚...
    茶点故事阅读 44,843评论 2 354

推荐阅读更多精彩内容