RepeatModeler2.0安裝及使用

http://www.reibang.com/p/bccdc3c68ace

RepeatModeler已更新烟馅,擁有了一些新的功能尾抑,可以識(shí)別LTR的結(jié)構(gòu)

在基因組注釋中第一步就是重復(fù)序列的屏蔽均函,目前常用的從頭注釋pipeline就是RepeatModeler + RepeatMasker。

RepeatModeler的依賴工具可以conda安裝或者手動(dòng)安裝我磁,推薦conda

安裝需要的環(huán)境

perl和Text::Soundex模塊-- V5.8.8及以上版本
Python 3 和 h5py 模塊 -- RepeatMasker軟件需要
RECON -- De Novo Repeat Finder
RepeatScout -- De Novo Repeat Finder
TRF-- Tandem Repeat Finder
RepeatMasker & Libraries
RMBlast NCBI Blast的修改版本,可與RepeatMasker和RepeatModeler一起使用,可選
ABBlast 比NCBI Blast更快趴泌,靈敏度更高的搜索引擎敞葛,可選

Text::Soundex和h5py 模塊
cpan install Text::Soundex #perl模塊安裝
pip3 install numpy scipy matplotlib -i https://pypi.tuna.tsinghua.edu.cn/simple #h5py必須先裝此模塊誉察,-i表示清華鏡像
pip3 install h5py -i https://pypi.tuna.tsinghua.edu.cn/simple 

1、RepeatMasker惹谐、TRF持偏、RMblast
在正式安裝RepeatModeler主程序之前,需要提前配置好其依賴的工具氨肌。
見(jiàn)前文(RepeatMasker的安裝與使用)(http://www.reibang.com/p/ffdbedae80fa)鸿秆。

2、RECON(從頭預(yù)測(cè)功能實(shí)現(xiàn)的核心組件)

conda install RECON    #自動(dòng)配置在conda環(huán)境中
wget http://www.repeatmasker.org/RepeatModeler/RECON-1.08.tar.gz
tar -zxvf RECON-1.08.tar.gz 
cd RECON-1.08/src
make ; make install #已經(jīng)將使用軟件安裝在 RECON-1.08/bin下

3怎囚、RepeatScout(從頭預(yù)測(cè)功能實(shí)現(xiàn)的核心組件)

conda install RepeatScout


wget http://www.repeatmasker.org/RepeatScout-1.0.6.tar.gz
tar -xzvf RepeatScout-1.0.6.tar.gz
cd RepeatScout-1.0.6 
make # 會(huì)生產(chǎn)build_lmer_table和RepeatScout兩個(gè)程序

4卿叽、可選軟件,運(yùn)行LTR結(jié)構(gòu)搜索必須軟件恳守,也可以不安裝
LtrHarvest(是GenomeTools套件的一部分考婴,安裝genometools即可。)
MAFFT

conda install genometools-genometools
conda install -c bioconda mafft



wget http://genometools.org/pub/genometools-1.6.2.tar.gz
tar -pzxvf genometools-1.6.2.tar.gz
cd genometools-1.6.2
make threads=yes  #設(shè)置多線程模式
make prefix=/software/annotation/genometools/1.6.2 install



Ltr_retriever催烘,解壓即可使用

https://github.com/oushujun/LTR_retriever/archive/master.zip
unzip LTR_retriever-master

CD-HIT和Ninja沥阱,解壓后編譯一下就好了

wget https://github.com/weizhongli/cdhit/archive/master.zip
unzip cdhit-master.zip && cd cdhit-master/
sudo make

wget https://github.com/TravisWheelerLab/NINJA/archive/master.zip
unzip master.zip && cd NINJA/
sudo make

5、RepeatModeler
前面的工具配置好后伊群,源碼編譯RepeatModeler考杉,和RepeatMasker安裝過(guò)程很像,需要一步步指定所依賴的環(huán)境舰始。

#RepeatModeler
wget http://www.repeatmasker.org/RepeatModeler/RepeatModeler-2.0.1.tar.gz
tar xzvf RepeatModeler-2.0.1.tar.gz

cd RepeatModeler-2.0.1.tar.gz
chmod -R 755 *
perl ./configure

“./configure”執(zhí)行后崇棠,根據(jù)提示信息一步步來(lái)。

首先是perl環(huán)境丸卷,推薦使用/usr/bin環(huán)境下的perl枕稀,即系統(tǒng)perl,使用conda中的perl后面總出問(wèn)題∥溃回車?yán)^續(xù)范抓。
RepeatModeler安裝路徑,默認(rèn)自動(dòng)指定安裝路徑食铐,回車?yán)^續(xù)匕垫。
指定RepeatMasker的安裝路徑,回車?yán)^續(xù)虐呻。
指定RECON的安裝路徑象泵,回車?yán)^續(xù)。
指定RepeatScout的安裝路徑斟叼,回車?yán)^續(xù)偶惠。
指定NSEG的安裝路徑,回車?yán)^續(xù)朗涩。
指定TRF的安裝路徑(這個(gè)同時(shí)也是RepeatMasker安裝必需的)猪落,回車?yán)^續(xù)实撒。

最后是序列搜索引擎(這個(gè)同時(shí)也是RepeatMasker安裝必需的)蛙吏,例如我們這里選擇2洲赵,指定RMBlast主程序所在路徑,回車返回主界面后识腿,再選擇3出革,就完成了。
你也可以指定多種序列搜索引擎后渡讼,再選擇3骂束,不過(guò)實(shí)際運(yùn)行時(shí),一次只能選擇一種序列比對(duì)方式成箫。

這時(shí)會(huì)有提示:
這個(gè)版本的RepeatModeler可以檢測(cè)LTR的結(jié)構(gòu)展箱,是否配置,可以根據(jù)需要自行選擇蹬昌。若需要就按y混驰,然后根據(jù)提示選擇路徑完成配置。

image

這就完成RepeatModeler的安裝凳厢。

終于把一些列的環(huán)境添加完了账胧,最后配置環(huán)境變量竞慢。

#例如先紫,我的RepeatModeler安裝路徑是在“/home/my/software/RepeatModeler-open-1.0.11”
export PATH=/home/my/software/RepeatModeler-2.0.1:$PATH

#這時(shí)候沒(méi)啥問(wèn)題的話應(yīng)該可以看到幫助界面了
RepeatModeler -h

使用示例:

BuildDatabase -name human human_genomic.fa
RepeatModeler -pa 4 -database human -LTRStruct 
RepeatMasker -e rmblast -lib human-families.fa -pa 4 human_genomic.fa

后續(xù)結(jié)果跟RepeatModeler1.0的結(jié)果差不多
祝大家科研順利!3镏蟆遮精!

?著作權(quán)歸作者所有,轉(zhuǎn)載或內(nèi)容合作請(qǐng)聯(lián)系作者
  • 序言:七十年代末,一起剝皮案震驚了整個(gè)濱河市,隨后出現(xiàn)的幾起案子本冲,更是在濱河造成了極大的恐慌准脂,老刑警劉巖,帶你破解...
    沈念sama閱讀 218,122評(píng)論 6 505
  • 序言:濱河連續(xù)發(fā)生了三起死亡事件檬洞,死亡現(xiàn)場(chǎng)離奇詭異狸膏,居然都是意外死亡,警方通過(guò)查閱死者的電腦和手機(jī)添怔,發(fā)現(xiàn)死者居然都...
    沈念sama閱讀 93,070評(píng)論 3 395
  • 文/潘曉璐 我一進(jìn)店門湾戳,熙熙樓的掌柜王于貴愁眉苦臉地迎上來(lái),“玉大人广料,你說(shuō)我怎么就攤上這事砾脑。” “怎么了艾杏?”我有些...
    開(kāi)封第一講書(shū)人閱讀 164,491評(píng)論 0 354
  • 文/不壞的土叔 我叫張陵韧衣,是天一觀的道長(zhǎng)。 經(jīng)常有香客問(wèn)我购桑,道長(zhǎng)畅铭,這世上最難降的妖魔是什么? 我笑而不...
    開(kāi)封第一講書(shū)人閱讀 58,636評(píng)論 1 293
  • 正文 為了忘掉前任勃蜘,我火速辦了婚禮顶瞒,結(jié)果婚禮上,老公的妹妹穿的比我還像新娘元旬。我一直安慰自己榴徐,他們只是感情好,可當(dāng)我...
    茶點(diǎn)故事閱讀 67,676評(píng)論 6 392
  • 文/花漫 我一把揭開(kāi)白布匀归。 她就那樣靜靜地躺著坑资,像睡著了一般。 火紅的嫁衣襯著肌膚如雪穆端。 梳的紋絲不亂的頭發(fā)上袱贮,一...
    開(kāi)封第一講書(shū)人閱讀 51,541評(píng)論 1 305
  • 那天,我揣著相機(jī)與錄音体啰,去河邊找鬼攒巍。 笑死,一個(gè)胖子當(dāng)著我的面吹牛荒勇,可吹牛的內(nèi)容都是我干的柒莉。 我是一名探鬼主播,決...
    沈念sama閱讀 40,292評(píng)論 3 418
  • 文/蒼蘭香墨 我猛地睜開(kāi)眼沽翔,長(zhǎng)吁一口氣:“原來(lái)是場(chǎng)噩夢(mèng)啊……” “哼兢孝!你這毒婦竟也來(lái)了窿凤?” 一聲冷哼從身側(cè)響起,我...
    開(kāi)封第一講書(shū)人閱讀 39,211評(píng)論 0 276
  • 序言:老撾萬(wàn)榮一對(duì)情侶失蹤跨蟹,失蹤者是張志新(化名)和其女友劉穎雳殊,沒(méi)想到半個(gè)月后,有當(dāng)?shù)厝嗽跇?shù)林里發(fā)現(xiàn)了一具尸體窗轩,經(jīng)...
    沈念sama閱讀 45,655評(píng)論 1 314
  • 正文 獨(dú)居荒郊野嶺守林人離奇死亡夯秃,尸身上長(zhǎng)有42處帶血的膿包…… 初始之章·張勛 以下內(nèi)容為張勛視角 年9月15日...
    茶點(diǎn)故事閱讀 37,846評(píng)論 3 336
  • 正文 我和宋清朗相戀三年,在試婚紗的時(shí)候發(fā)現(xiàn)自己被綠了痢艺。 大學(xué)時(shí)的朋友給我發(fā)了我未婚夫和他白月光在一起吃飯的照片寝并。...
    茶點(diǎn)故事閱讀 39,965評(píng)論 1 348
  • 序言:一個(gè)原本活蹦亂跳的男人離奇死亡,死狀恐怖腹备,靈堂內(nèi)的尸體忽然破棺而出衬潦,到底是詐尸還是另有隱情,我是刑警寧澤植酥,帶...
    沈念sama閱讀 35,684評(píng)論 5 347
  • 正文 年R本政府宣布镀岛,位于F島的核電站,受9級(jí)特大地震影響友驮,放射性物質(zhì)發(fā)生泄漏漂羊。R本人自食惡果不足惜,卻給世界環(huán)境...
    茶點(diǎn)故事閱讀 41,295評(píng)論 3 329
  • 文/蒙蒙 一卸留、第九天 我趴在偏房一處隱蔽的房頂上張望走越。 院中可真熱鬧,春花似錦耻瑟、人聲如沸旨指。這莊子的主人今日做“春日...
    開(kāi)封第一講書(shū)人閱讀 31,894評(píng)論 0 22
  • 文/蒼蘭香墨 我抬頭看了看天上的太陽(yáng)谆构。三九已至,卻和暖如春框都,著一層夾襖步出監(jiān)牢的瞬間搬素,已是汗流浹背。 一陣腳步聲響...
    開(kāi)封第一講書(shū)人閱讀 33,012評(píng)論 1 269
  • 我被黑心中介騙來(lái)泰國(guó)打工魏保, 沒(méi)想到剛下飛機(jī)就差點(diǎn)兒被人妖公主榨干…… 1. 我叫王不留熬尺,地道東北人。 一個(gè)月前我還...
    沈念sama閱讀 48,126評(píng)論 3 370
  • 正文 我出身青樓谓罗,卻偏偏與公主長(zhǎng)得像粱哼,于是被迫代替她去往敵國(guó)和親。 傳聞我的和親對(duì)象是個(gè)殘疾皇子妥衣,可洞房花燭夜當(dāng)晚...
    茶點(diǎn)故事閱讀 44,914評(píng)論 2 355

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