生信小白——初學(xué) bowtie2

本文就以生信小白的角度碧浊,簡述bowtie2的下載、安裝及使用必孤。關(guān)于bowtie2的使用教程非常豐富,具體的概念本文不再贅述兴想,只對操作過程進(jìn)行簡(xi)簡(zhi)單(ru)單(wei)的介紹赡勘,希望對同樣的初學(xué)者提供一些幫助闸与。

我是使用Mac終端連接服務(wù)器環(huán)境下使用的。


下載及安裝

起初我是在服務(wù)器環(huán)境下使用 wget+鏈接 下載軟件厂画,運(yùn)行時發(fā)現(xiàn)安裝失敗拷邢。

就是這樣瞭稼,安裝失敗

于是在服務(wù)器里下載miniconda,并添加了bioconda鏡像欲虚。(conda安裝教程http://www.reibang.com/p/edaa744ea47d)悔雹。conda環(huán)境下輸入conda iinstall gatk bowtie2荠商。檢索結(jié)果輸入y。

檢索后輸入“y”

在miniconda/bin目錄下輸入ls初肉,可以看到bowtie2牙咏,輸入 ./bowtie2 -help嘹裂,出現(xiàn)help說明軟件就可以用了。

安裝成功

使用

序列數(shù)據(jù)

我已經(jīng)把需要比對的序列和參考基因組上傳服務(wù)器了。read 目錄下是我的測序文件a_0.fq伊磺、 a_20.fq 和 a_60.fq,ref 目錄下是我的參考基因組 Tas.V1.genome.fa

read目錄下是我測序獲得的三條序列豪筝,ref目錄下是我的參考基因組

建索引

進(jìn)入ref目錄续崖。輸入bowtie2-build Tsa.v1.genome.fa(參考基因組) Tas(輸出文件名前綴)

建索引

需要一段時間严望,我170M的基因組大概1分鐘多時間莺丑。完成后會獲得6個 .bt2后綴的文件墩蔓,如下,說明索引建立成功昏名。

獲得索引文件

序列比對

這里的參數(shù)較多轻局,選取主要的使用即可样刷。通常測序文件是雙端的置鼻,這里我的是序列是單端,所以使用的參數(shù)也不同储藐。首先是雙端序列(我沒跑雙端嘶是,結(jié)果就不呈現(xiàn)了),單端序列結(jié)果如下辖源。

進(jìn)入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(指定輸出文件名)

單端序列

初步結(jié)果顯示,a_0.fp 文件中有172333條序列洲尊,其中只有8.79%的序列比對到參考基因組上坞嘀。

比對結(jié)果

以此類推,我獲得了其余兩個序列的sam文件棺滞。

獲得三個序列的sam文件


參考及其感謝

hsy_hzauer??http://www.reibang.com/u/49bcf3035565

賣萌哥????conda的安裝與使用(最新版):http://www.reibang.com/p/edaa744ea47d

郭威 ? ?序列比對軟件bowtie2--快速入門:https://v.qq.com/x/page/m0846mwp8cl.html

最后編輯于
?著作權(quán)歸作者所有,轉(zhuǎn)載或內(nèi)容合作請聯(lián)系作者
  • 序言:七十年代末,一起剝皮案震驚了整個濱河市矮男,隨后出現(xiàn)的幾起案子毡鉴,更是在濱河造成了極大的恐慌,老刑警劉巖憎瘸,帶你破解...
    沈念sama閱讀 222,681評論 6 517
  • 序言:濱河連續(xù)發(fā)生了三起死亡事件幌甘,死亡現(xiàn)場離奇詭異甘晤,居然都是意外死亡,警方通過查閱死者的電腦和手機(jī)遏弱,發(fā)現(xiàn)死者居然都...
    沈念sama閱讀 95,205評論 3 399
  • 文/潘曉璐 我一進(jìn)店門,熙熙樓的掌柜王于貴愁眉苦臉地迎上來塞弊,“玉大人漱逸,你說我怎么就攤上這事泪姨。” “怎么了饰抒?”我有些...
    開封第一講書人閱讀 169,421評論 0 362
  • 文/不壞的土叔 我叫張陵肮砾,是天一觀的道長。 經(jīng)常有香客問我袋坑,道長仗处,這世上最難降的妖魔是什么? 我笑而不...
    開封第一講書人閱讀 60,114評論 1 300
  • 正文 為了忘掉前任枣宫,我火速辦了婚禮,結(jié)果婚禮上也颤,老公的妹妹穿的比我還像新娘洋幻。我一直安慰自己,他們只是感情好翅娶,可當(dāng)我...
    茶點(diǎn)故事閱讀 69,116評論 6 398
  • 文/花漫 我一把揭開白布文留。 她就那樣靜靜地躺著,像睡著了一般竭沫。 火紅的嫁衣襯著肌膚如雪燥翅。 梳的紋絲不亂的頭發(fā)上,一...
    開封第一講書人閱讀 52,713評論 1 312
  • 那天输吏,我揣著相機(jī)與錄音权旷,去河邊找鬼替蛉。 笑死贯溅,一個胖子當(dāng)著我的面吹牛,可吹牛的內(nèi)容都是我干的躲查。 我是一名探鬼主播它浅,決...
    沈念sama閱讀 41,170評論 3 422
  • 文/蒼蘭香墨 我猛地睜開眼,長吁一口氣:“原來是場噩夢啊……” “哼镣煮!你這毒婦竟也來了姐霍?” 一聲冷哼從身側(cè)響起,我...
    開封第一講書人閱讀 40,116評論 0 277
  • 序言:老撾萬榮一對情侶失蹤典唇,失蹤者是張志新(化名)和其女友劉穎镊折,沒想到半個月后,有當(dāng)?shù)厝嗽跇淞掷锇l(fā)現(xiàn)了一具尸體介衔,經(jīng)...
    沈念sama閱讀 46,651評論 1 320
  • 正文 獨(dú)居荒郊野嶺守林人離奇死亡恨胚,尸身上長有42處帶血的膿包…… 初始之章·張勛 以下內(nèi)容為張勛視角 年9月15日...
    茶點(diǎn)故事閱讀 38,714評論 3 342
  • 正文 我和宋清朗相戀三年,在試婚紗的時候發(fā)現(xiàn)自己被綠了炎咖。 大學(xué)時的朋友給我發(fā)了我未婚夫和他白月光在一起吃飯的照片赃泡。...
    茶點(diǎn)故事閱讀 40,865評論 1 353
  • 序言:一個原本活蹦亂跳的男人離奇死亡寒波,死狀恐怖,靈堂內(nèi)的尸體忽然破棺而出升熊,到底是詐尸還是另有隱情俄烁,我是刑警寧澤,帶...
    沈念sama閱讀 36,527評論 5 351
  • 正文 年R本政府宣布级野,位于F島的核電站页屠,受9級特大地震影響,放射性物質(zhì)發(fā)生泄漏蓖柔。R本人自食惡果不足惜卷中,卻給世界環(huán)境...
    茶點(diǎn)故事閱讀 42,211評論 3 336
  • 文/蒙蒙 一、第九天 我趴在偏房一處隱蔽的房頂上張望渊抽。 院中可真熱鬧蟆豫,春花似錦、人聲如沸懒闷。這莊子的主人今日做“春日...
    開封第一講書人閱讀 32,699評論 0 25
  • 文/蒼蘭香墨 我抬頭看了看天上的太陽愤估。三九已至帮辟,卻和暖如春,著一層夾襖步出監(jiān)牢的瞬間玩焰,已是汗流浹背由驹。 一陣腳步聲響...
    開封第一講書人閱讀 33,814評論 1 274
  • 我被黑心中介騙來泰國打工, 沒想到剛下飛機(jī)就差點(diǎn)兒被人妖公主榨干…… 1. 我叫王不留昔园,地道東北人蔓榄。 一個月前我還...
    沈念sama閱讀 49,299評論 3 379
  • 正文 我出身青樓,卻偏偏與公主長得像默刚,于是被迫代替她去往敵國和親甥郑。 傳聞我的和親對象是個殘疾皇子,可洞房花燭夜當(dāng)晚...
    茶點(diǎn)故事閱讀 45,870評論 2 361

推薦閱讀更多精彩內(nèi)容