NCBI官方真核生物全基因組注釋流程首次開(kāi)源降允!嘗鮮體驗(yàn)踩坑貼奉上~

萌哥碎碎念:

NCBI的官方注釋流程一直都是出了名的效果拔群準(zhǔn)確度極高恩闻,但可惜一直都是非開(kāi)源的,看到這么好用的工具卻用不到自己的項(xiàng)目里剧董,心里可太刺撓了幢尚。

終于!最近NCBI開(kāi)源了他們的流程翅楼,下面的內(nèi)容由小伙伴阿良供稿大家如果有好的稿子也可以投稿給我呀

項(xiàng)目地址:

Eukaryotic Genome Annotation Pipeline - External (EGAPx)

使用體驗(yàn)

總的使用感受:最好能找到國(guó)外的好友幫忙尉剩,因?yàn)橄螺d依賴文件時(shí),網(wǎng)絡(luò)速度是個(gè)問(wèn)題

運(yùn)行說(shuō)明

本文采用singularity鏡像方式運(yùn)行毅臊,以示例數(shù)據(jù)為例理茎。

原因:
采用singularity鏡像可以避免aws S3部分的設(shè)置(這里會(huì)要求你開(kāi)通亞馬遜帳號(hào))。

整個(gè)流程要先配置nextflow管嬉,可以參考:

https://www.nextflow.io/docs/latest/install.html

軟件安裝

1. 先安裝nextflow

默認(rèn)你已經(jīng)裝了java, 沒(méi)安裝的話記得安裝皂林。

curl -s https://get.nextflow.io | bash
chmod +x nextflow
sudo mv nextflow /usr/local/bin
nextflow info

如果你沒(méi)有sudo權(quán)限,試試在conda環(huán)境下安裝nextflow后運(yùn)行蚯撩。

2. 下載singularity鏡像

可直接跳到github頁(yè)面offline mode部分础倍。

選擇適合的方式獲得鏡像(考驗(yàn)?zāi)愕木W(wǎng)絡(luò))

singularity pull docker://ncbi/egapx:0.2-alpha 

3. 下載本地緩存文件

下載“一些”注釋所需要的文件(考驗(yàn)?zāi)憔W(wǎng)速的時(shí)刻到了)

git clone https://github.com/ncbi/egapx.git
cd egapx
python3 ui/egapx.py -dl -lc ../local_cache

4. 下載示例數(shù)據(jù)

下載3中的緩存文件的同時(shí),下載示例所需數(shù)據(jù)

# 方法不限胎挎,能下載下來(lái)數(shù)據(jù)就行
prefetch SRR8506572
prefetch SRR9005248
fasterq-dump --skip-technical --threads 6 --split-files --seq-defline ">\$ac.\$si.\$ri" --fasta -O sradir/  ./SRR8506572
fasterq-dump --skip-technical --threads 6 --split-files --seq-defline ">\$ac.\$si.\$ri" --fasta -O sradir/  ./SRR9005248

5. 修改示例數(shù)據(jù)的配置文件內(nèi)容

vim examples/input_D_farinae_small.yaml

6. 避坑指南 (本文最重要的地方)

biowulf_cluster.config默認(rèn)的是NIH的biowulf集群沟启,配置文件雖然已經(jīng)設(shè)好,但我無(wú)法選用犹菇,因此放棄美浦。
考慮到自己用的本地singularity鏡像,因此在后續(xù)命令中采用-e singularity參數(shù)项栏。

避坑點(diǎn):直接使用會(huì)在運(yùn)行過(guò)程中再次下載singularity鏡像,既浪費(fèi)了前面已下好的蹬竖,又會(huì)因?yàn)榫W(wǎng)絡(luò)問(wèn)題失敗沼沈。
所以記得修改egapx_config/singularity.config文件:

echo "process.container = '/path_to_/egapx_0.2-alpha.sif'" >> egapx_config/singularity_config
#或者直接vim編輯 egapx_config/singularity_config文件

7. 開(kāi)始運(yùn)行示例數(shù)據(jù)

ui/egapx.py examples/input_D_farinae_small.yaml -e singularity -w dfs_work -o dfs_out -lc ../local_cache

目前正常運(yùn)行中!有后續(xù)再更新~


運(yùn)行中.png
?著作權(quán)歸作者所有,轉(zhuǎn)載或內(nèi)容合作請(qǐng)聯(lián)系作者
  • 序言:七十年代末币厕,一起剝皮案震驚了整個(gè)濱河市列另,隨后出現(xiàn)的幾起案子,更是在濱河造成了極大的恐慌旦装,老刑警劉巖页衙,帶你破解...
    沈念sama閱讀 206,839評(píng)論 6 482
  • 序言:濱河連續(xù)發(fā)生了三起死亡事件,死亡現(xiàn)場(chǎng)離奇詭異,居然都是意外死亡店乐,警方通過(guò)查閱死者的電腦和手機(jī)艰躺,發(fā)現(xiàn)死者居然都...
    沈念sama閱讀 88,543評(píng)論 2 382
  • 文/潘曉璐 我一進(jìn)店門(mén),熙熙樓的掌柜王于貴愁眉苦臉地迎上來(lái)眨八,“玉大人腺兴,你說(shuō)我怎么就攤上這事×啵” “怎么了页响?”我有些...
    開(kāi)封第一講書(shū)人閱讀 153,116評(píng)論 0 344
  • 文/不壞的土叔 我叫張陵,是天一觀的道長(zhǎng)段誊。 經(jīng)常有香客問(wèn)我闰蚕,道長(zhǎng),這世上最難降的妖魔是什么连舍? 我笑而不...
    開(kāi)封第一講書(shū)人閱讀 55,371評(píng)論 1 279
  • 正文 為了忘掉前任没陡,我火速辦了婚禮,結(jié)果婚禮上烟瞧,老公的妹妹穿的比我還像新娘诗鸭。我一直安慰自己,他們只是感情好参滴,可當(dāng)我...
    茶點(diǎn)故事閱讀 64,384評(píng)論 5 374
  • 文/花漫 我一把揭開(kāi)白布强岸。 她就那樣靜靜地躺著,像睡著了一般砾赔。 火紅的嫁衣襯著肌膚如雪蝌箍。 梳的紋絲不亂的頭發(fā)上,一...
    開(kāi)封第一講書(shū)人閱讀 49,111評(píng)論 1 285
  • 那天暴心,我揣著相機(jī)與錄音妓盲,去河邊找鬼。 笑死专普,一個(gè)胖子當(dāng)著我的面吹牛悯衬,可吹牛的內(nèi)容都是我干的。 我是一名探鬼主播檀夹,決...
    沈念sama閱讀 38,416評(píng)論 3 400
  • 文/蒼蘭香墨 我猛地睜開(kāi)眼筋粗,長(zhǎng)吁一口氣:“原來(lái)是場(chǎng)噩夢(mèng)啊……” “哼!你這毒婦竟也來(lái)了炸渡?” 一聲冷哼從身側(cè)響起娜亿,我...
    開(kāi)封第一講書(shū)人閱讀 37,053評(píng)論 0 259
  • 序言:老撾萬(wàn)榮一對(duì)情侶失蹤,失蹤者是張志新(化名)和其女友劉穎蚌堵,沒(méi)想到半個(gè)月后买决,有當(dāng)?shù)厝嗽跇?shù)林里發(fā)現(xiàn)了一具尸體沛婴,經(jīng)...
    沈念sama閱讀 43,558評(píng)論 1 300
  • 正文 獨(dú)居荒郊野嶺守林人離奇死亡,尸身上長(zhǎng)有42處帶血的膿包…… 初始之章·張勛 以下內(nèi)容為張勛視角 年9月15日...
    茶點(diǎn)故事閱讀 36,007評(píng)論 2 325
  • 正文 我和宋清朗相戀三年督赤,在試婚紗的時(shí)候發(fā)現(xiàn)自己被綠了嘁灯。 大學(xué)時(shí)的朋友給我發(fā)了我未婚夫和他白月光在一起吃飯的照片。...
    茶點(diǎn)故事閱讀 38,117評(píng)論 1 334
  • 序言:一個(gè)原本活蹦亂跳的男人離奇死亡够挂,死狀恐怖旁仿,靈堂內(nèi)的尸體忽然破棺而出,到底是詐尸還是另有隱情孽糖,我是刑警寧澤枯冈,帶...
    沈念sama閱讀 33,756評(píng)論 4 324
  • 正文 年R本政府宣布,位于F島的核電站办悟,受9級(jí)特大地震影響尘奏,放射性物質(zhì)發(fā)生泄漏。R本人自食惡果不足惜病蛉,卻給世界環(huán)境...
    茶點(diǎn)故事閱讀 39,324評(píng)論 3 307
  • 文/蒙蒙 一炫加、第九天 我趴在偏房一處隱蔽的房頂上張望。 院中可真熱鬧铺然,春花似錦俗孝、人聲如沸。這莊子的主人今日做“春日...
    開(kāi)封第一講書(shū)人閱讀 30,315評(píng)論 0 19
  • 文/蒼蘭香墨 我抬頭看了看天上的太陽(yáng)。三九已至沽瘦,卻和暖如春革骨,著一層夾襖步出監(jiān)牢的瞬間,已是汗流浹背析恋。 一陣腳步聲響...
    開(kāi)封第一講書(shū)人閱讀 31,539評(píng)論 1 262
  • 我被黑心中介騙來(lái)泰國(guó)打工良哲, 沒(méi)想到剛下飛機(jī)就差點(diǎn)兒被人妖公主榨干…… 1. 我叫王不留,地道東北人助隧。 一個(gè)月前我還...
    沈念sama閱讀 45,578評(píng)論 2 355
  • 正文 我出身青樓筑凫,卻偏偏與公主長(zhǎng)得像,于是被迫代替她去往敵國(guó)和親并村。 傳聞我的和親對(duì)象是個(gè)殘疾皇子漏健,可洞房花燭夜當(dāng)晚...
    茶點(diǎn)故事閱讀 42,877評(píng)論 2 345

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