kmergenie 安裝及使用

由于需要做微生物的De novo無(wú)參組裝志鹃,查了一些資料發(fā)現(xiàn)大多軟件需要先確定K值肠槽,而kmergenie能夠根據(jù)數(shù)據(jù)的得到最佳的K值悍缠,下邊記錄一下安裝過(guò)程

依賴(lài):

python >= 2.7?

R (r-base) (如果沒(méi)有直接使用 conda install r-base 進(jìn)行安裝)

1雁乡、使用conda安裝(失敗)

conda 是生物信息學(xué)學(xué)習(xí)的必須要掌握的軟件=滞>氩住!

首先在anaconda 中查到kmergenie它匕,搜索結(jié)果如下:

然后使用命令:

conda install?kmergenie? ?進(jìn)行安裝

安裝后使用時(shí)發(fā)現(xiàn)出現(xiàn)錯(cuò)誤:

ModuleNotFoundError: No module named 'readfq'

然后我繼續(xù)用conda安裝readfq:

conda install readfq

仍然沒(méi)有解決這個(gè)問(wèn)題展融,查找了資料說(shuō)是要把readfq的依賴(lài)放到kmergenie中:

kmergenie missing readfq · Issue #37856 · bioconda/bioconda-recipes · GitHub

果斷放棄,選擇自己安裝:

2豫柬、編譯安裝

首先從官網(wǎng)(KmerGenie)下載安裝包告希,安裝步驟如下:

tar -xzvf?kmergenie-1.7051.tar.gz

cd?kmergenie-1.7051/

make

python setup.py install --user? (--user 沒(méi)有管理員權(quán)限的用戶要加上,不然會(huì)報(bào)錯(cuò))

完成后執(zhí)行命令試下: kergenie -h

這就安裝成功了烧给。下邊將其添加進(jìn)環(huán)境變量就行了燕偶,一般我是更改.bashrc 文件在最后一行添加下邊內(nèi)容:

export PATH=$PATH:/public2/home/Software/kmergenie-1.7051

然后 source .bashrc 就ok啦!

3础嫡、使用

kergenie? sequence.txt? -o? kergenie_output/sample? -k 140 -l 15 -s 10 -t 10

sequence.txt 為分析文件(fastq)的絕對(duì)路徑指么,多個(gè)文件要換行

-k 為選擇分析最大的K值

-l 為選擇分析最小的K值

-s 最小值到最大值的步長(zhǎng)

-t 線程數(shù)

注意:第一次分析是s可以設(shè)大點(diǎn),可以根據(jù)分析結(jié)果給出的最佳K值再縮小再進(jìn)一步的分析(第一次:I:15,k:140,s:10? 第二次I:15,k:80,s:6)

參考:

基因組組裝----SOAPdenovo2 - 簡(jiǎn)書(shū)

http://kmergenie.bx.psu.edu/README

?著作權(quán)歸作者所有,轉(zhuǎn)載或內(nèi)容合作請(qǐng)聯(lián)系作者
  • 序言:七十年代末榴鼎,一起剝皮案震驚了整個(gè)濱河市伯诬,隨后出現(xiàn)的幾起案子,更是在濱河造成了極大的恐慌巫财,老刑警劉巖盗似,帶你破解...
    沈念sama閱讀 219,589評(píng)論 6 508
  • 序言:濱河連續(xù)發(fā)生了三起死亡事件,死亡現(xiàn)場(chǎng)離奇詭異平项,居然都是意外死亡赫舒,警方通過(guò)查閱死者的電腦和手機(jī),發(fā)現(xiàn)死者居然都...
    沈念sama閱讀 93,615評(píng)論 3 396
  • 文/潘曉璐 我一進(jìn)店門(mén)闽瓢,熙熙樓的掌柜王于貴愁眉苦臉地迎上來(lái)接癌,“玉大人,你說(shuō)我怎么就攤上這事扣讼∪泵停” “怎么了?”我有些...
    開(kāi)封第一講書(shū)人閱讀 165,933評(píng)論 0 356
  • 文/不壞的土叔 我叫張陵,是天一觀的道長(zhǎng)枯夜。 經(jīng)常有香客問(wèn)我弯汰,道長(zhǎng)艰山,這世上最難降的妖魔是什么湖雹? 我笑而不...
    開(kāi)封第一講書(shū)人閱讀 58,976評(píng)論 1 295
  • 正文 為了忘掉前任,我火速辦了婚禮曙搬,結(jié)果婚禮上摔吏,老公的妹妹穿的比我還像新娘。我一直安慰自己纵装,他們只是感情好征讲,可當(dāng)我...
    茶點(diǎn)故事閱讀 67,999評(píng)論 6 393
  • 文/花漫 我一把揭開(kāi)白布。 她就那樣靜靜地躺著橡娄,像睡著了一般诗箍。 火紅的嫁衣襯著肌膚如雪。 梳的紋絲不亂的頭發(fā)上挽唉,一...
    開(kāi)封第一講書(shū)人閱讀 51,775評(píng)論 1 307
  • 那天滤祖,我揣著相機(jī)與錄音,去河邊找鬼瓶籽。 笑死匠童,一個(gè)胖子當(dāng)著我的面吹牛,可吹牛的內(nèi)容都是我干的塑顺。 我是一名探鬼主播汤求,決...
    沈念sama閱讀 40,474評(píng)論 3 420
  • 文/蒼蘭香墨 我猛地睜開(kāi)眼,長(zhǎng)吁一口氣:“原來(lái)是場(chǎng)噩夢(mèng)啊……” “哼严拒!你這毒婦竟也來(lái)了扬绪?” 一聲冷哼從身側(cè)響起,我...
    開(kāi)封第一講書(shū)人閱讀 39,359評(píng)論 0 276
  • 序言:老撾萬(wàn)榮一對(duì)情侶失蹤裤唠,失蹤者是張志新(化名)和其女友劉穎勒奇,沒(méi)想到半個(gè)月后,有當(dāng)?shù)厝嗽跇?shù)林里發(fā)現(xiàn)了一具尸體巧骚,經(jīng)...
    沈念sama閱讀 45,854評(píng)論 1 317
  • 正文 獨(dú)居荒郊野嶺守林人離奇死亡赊颠,尸身上長(zhǎng)有42處帶血的膿包…… 初始之章·張勛 以下內(nèi)容為張勛視角 年9月15日...
    茶點(diǎn)故事閱讀 38,007評(píng)論 3 338
  • 正文 我和宋清朗相戀三年,在試婚紗的時(shí)候發(fā)現(xiàn)自己被綠了劈彪。 大學(xué)時(shí)的朋友給我發(fā)了我未婚夫和他白月光在一起吃飯的照片竣蹦。...
    茶點(diǎn)故事閱讀 40,146評(píng)論 1 351
  • 序言:一個(gè)原本活蹦亂跳的男人離奇死亡,死狀恐怖沧奴,靈堂內(nèi)的尸體忽然破棺而出痘括,到底是詐尸還是另有隱情,我是刑警寧澤,帶...
    沈念sama閱讀 35,826評(píng)論 5 346
  • 正文 年R本政府宣布纲菌,位于F島的核電站挠日,受9級(jí)特大地震影響,放射性物質(zhì)發(fā)生泄漏翰舌。R本人自食惡果不足惜嚣潜,卻給世界環(huán)境...
    茶點(diǎn)故事閱讀 41,484評(píng)論 3 331
  • 文/蒙蒙 一、第九天 我趴在偏房一處隱蔽的房頂上張望椅贱。 院中可真熱鬧懂算,春花似錦、人聲如沸庇麦。這莊子的主人今日做“春日...
    開(kāi)封第一講書(shū)人閱讀 32,029評(píng)論 0 22
  • 文/蒼蘭香墨 我抬頭看了看天上的太陽(yáng)山橄。三九已至垮媒,卻和暖如春,著一層夾襖步出監(jiān)牢的瞬間航棱,已是汗流浹背睡雇。 一陣腳步聲響...
    開(kāi)封第一講書(shū)人閱讀 33,153評(píng)論 1 272
  • 我被黑心中介騙來(lái)泰國(guó)打工, 沒(méi)想到剛下飛機(jī)就差點(diǎn)兒被人妖公主榨干…… 1. 我叫王不留丧诺,地道東北人入桂。 一個(gè)月前我還...
    沈念sama閱讀 48,420評(píng)論 3 373
  • 正文 我出身青樓,卻偏偏與公主長(zhǎng)得像驳阎,于是被迫代替她去往敵國(guó)和親抗愁。 傳聞我的和親對(duì)象是個(gè)殘疾皇子,可洞房花燭夜當(dāng)晚...
    茶點(diǎn)故事閱讀 45,107評(píng)論 2 356

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