ITSx軟件的安裝(提取ITS序列)

該軟件可以自動(dòng)從序列中提取ITS1偿凭、ITS2、5.8S的序列并標(biāo)注位置,下面是這款軟件的安裝與使用方法寒屯。
第一步 拼接核糖體基因組

cd /home/monkeyflower/bioworkplace/test
conda activate getorganelle
get_organelle_from_reads.py -1test.1.fastq.gz -2 test.2.fastq.gz -F embplant_nr -o test -R 10 -t 2 -k 35,85,115
#轉(zhuǎn)到工作目錄,再激活getorganelle,在進(jìn)行核糖體基因組拼接。

第二步 下載ITSx安裝包

cd /home/monkeyflower/biosoft
wget -c https://microbiology.se/sw/ITSx_1.1.2.tar.gz
tar -zxvf ITSx_1.1.2.tar.gz
#先轉(zhuǎn)到/biosoft文件夾寡夹,再下載并解壓縮ITSx处面。解壓縮的文件夾里包含了ITSx, ITSx_db, user’s guide, README.txt , license.txt這幾個(gè)文件與文件夾。

第三步 安裝軟件

cd /home/monkeyflower/xiaodeng
mkdir /home/monkeyflower/xiaodeng/bin
mv /home/monkeyflower/biosoft/ITSx_1.1.2/ITSx_db /home/monkeyflower/xiaodeng/bin
mv /home/monkeyflower/biosoft/ITSx_1.1.2/ITSx /home/monkeyflower/xiaodeng/bin
#轉(zhuǎn)到自己的工作目錄菩掏;再創(chuàng)建一個(gè)bin文件魂角,用來存放ITSx的二進(jìn)制命令;再將ITSx_db和ITSx移動(dòng)到bin文件夾中智绸。

第四步 添加到環(huán)境變量

vim ~/.bashrc
a
#按a進(jìn)入編輯模式
export PATH=/home/monkeyflower/xiaodeng/bin:$PATH
Esc
#按Esc鍵退出編輯
:wq
#輸入:wq保存
source ~/.bashrc
#激活環(huán)境

第五步 測試

mv /home/monkeyflower/biosoft/ITSx_1.1.2/test.fasta /home/monkeyflower/bioworkplace/AAA
#將測試文件移動(dòng)到工作目錄
ITSx -i test.fasta -o test
#報(bào)錯(cuò):-bash: /home/monkeyflower/xiaodeng/bin/ITSx: 權(quán)限不夠
ls -l
#發(fā)現(xiàn)ITSx的權(quán)限為-rw-rw-r--野揪,權(quán)限較低。
chmod 775 ITSx
#更改權(quán)限為-rwxrwxr-x
ITSx -i test.fasta -o test
#再次測試成功传于。

第六步 結(jié)果文件解讀
結(jié)果會(huì)產(chǎn)生6個(gè)文件囱挑;包括匯總文件1個(gè),ITS亞區(qū)詳細(xì)位置表格1個(gè)沼溜,graph文件1個(gè)平挑,完整ITS序列文件1個(gè),ITS1和ITS2區(qū)域的FASTA文件各1個(gè)系草。此外通熄,如果發(fā)現(xiàn)不包含任何ITS區(qū)域的序列,也會(huì)寫出代表這些條目的序列ID列表找都。其中.ITS1.fasta唇辨、.ITS2.fasta就是我們所需要的序列。
也有更加簡單的安裝方式能耻,直接conda install -c bioconda itsx赏枚。但有時(shí)候就喜歡折騰折騰自己
總結(jié)起來,就是先將核糖體基因組拼接好晓猛,再用ITSx進(jìn)行提取饿幅,核心的兩條命令是

get_organelle_from_reads.py -1test.1.fastq.gz -2 test.2.fastq.gz -F embplant_nr -o test -R 10 -t 2 -k 35,85,115
#拼接核糖體基因
ITSx -i test.fasta -o test
#提取ITS序列
最后編輯于
?著作權(quán)歸作者所有,轉(zhuǎn)載或內(nèi)容合作請(qǐng)聯(lián)系作者
  • 序言:七十年代末,一起剝皮案震驚了整個(gè)濱河市戒职,隨后出現(xiàn)的幾起案子栗恩,更是在濱河造成了極大的恐慌,老刑警劉巖洪燥,帶你破解...
    沈念sama閱讀 210,978評(píng)論 6 490
  • 序言:濱河連續(xù)發(fā)生了三起死亡事件磕秤,死亡現(xiàn)場離奇詭異,居然都是意外死亡捧韵,警方通過查閱死者的電腦和手機(jī)市咆,發(fā)現(xiàn)死者居然都...
    沈念sama閱讀 89,954評(píng)論 2 384
  • 文/潘曉璐 我一進(jìn)店門,熙熙樓的掌柜王于貴愁眉苦臉地迎上來再来,“玉大人床绪,你說我怎么就攤上這事。” “怎么了癞己?”我有些...
    開封第一講書人閱讀 156,623評(píng)論 0 345
  • 文/不壞的土叔 我叫張陵,是天一觀的道長梭伐。 經(jīng)常有香客問我痹雅,道長,這世上最難降的妖魔是什么糊识? 我笑而不...
    開封第一講書人閱讀 56,324評(píng)論 1 282
  • 正文 為了忘掉前任绩社,我火速辦了婚禮,結(jié)果婚禮上赂苗,老公的妹妹穿的比我還像新娘愉耙。我一直安慰自己,他們只是感情好拌滋,可當(dāng)我...
    茶點(diǎn)故事閱讀 65,390評(píng)論 5 384
  • 文/花漫 我一把揭開白布朴沿。 她就那樣靜靜地躺著,像睡著了一般败砂。 火紅的嫁衣襯著肌膚如雪赌渣。 梳的紋絲不亂的頭發(fā)上,一...
    開封第一講書人閱讀 49,741評(píng)論 1 289
  • 那天昌犹,我揣著相機(jī)與錄音坚芜,去河邊找鬼。 笑死斜姥,一個(gè)胖子當(dāng)著我的面吹牛鸿竖,可吹牛的內(nèi)容都是我干的。 我是一名探鬼主播铸敏,決...
    沈念sama閱讀 38,892評(píng)論 3 405
  • 文/蒼蘭香墨 我猛地睜開眼缚忧,長吁一口氣:“原來是場噩夢啊……” “哼!你這毒婦竟也來了搞坝?” 一聲冷哼從身側(cè)響起搔谴,我...
    開封第一講書人閱讀 37,655評(píng)論 0 266
  • 序言:老撾萬榮一對(duì)情侶失蹤,失蹤者是張志新(化名)和其女友劉穎桩撮,沒想到半個(gè)月后敦第,有當(dāng)?shù)厝嗽跇淞掷锇l(fā)現(xiàn)了一具尸體,經(jīng)...
    沈念sama閱讀 44,104評(píng)論 1 303
  • 正文 獨(dú)居荒郊野嶺守林人離奇死亡店量,尸身上長有42處帶血的膿包…… 初始之章·張勛 以下內(nèi)容為張勛視角 年9月15日...
    茶點(diǎn)故事閱讀 36,451評(píng)論 2 325
  • 正文 我和宋清朗相戀三年芜果,在試婚紗的時(shí)候發(fā)現(xiàn)自己被綠了。 大學(xué)時(shí)的朋友給我發(fā)了我未婚夫和他白月光在一起吃飯的照片融师。...
    茶點(diǎn)故事閱讀 38,569評(píng)論 1 340
  • 序言:一個(gè)原本活蹦亂跳的男人離奇死亡右钾,死狀恐怖,靈堂內(nèi)的尸體忽然破棺而出,到底是詐尸還是另有隱情舀射,我是刑警寧澤窘茁,帶...
    沈念sama閱讀 34,254評(píng)論 4 328
  • 正文 年R本政府宣布,位于F島的核電站脆烟,受9級(jí)特大地震影響山林,放射性物質(zhì)發(fā)生泄漏。R本人自食惡果不足惜邢羔,卻給世界環(huán)境...
    茶點(diǎn)故事閱讀 39,834評(píng)論 3 312
  • 文/蒙蒙 一驼抹、第九天 我趴在偏房一處隱蔽的房頂上張望。 院中可真熱鬧拜鹤,春花似錦框冀、人聲如沸。這莊子的主人今日做“春日...
    開封第一講書人閱讀 30,725評(píng)論 0 21
  • 文/蒼蘭香墨 我抬頭看了看天上的太陽。三九已至极谊,卻和暖如春诡右,著一層夾襖步出監(jiān)牢的瞬間,已是汗流浹背轻猖。 一陣腳步聲響...
    開封第一講書人閱讀 31,950評(píng)論 1 264
  • 我被黑心中介騙來泰國打工帆吻, 沒想到剛下飛機(jī)就差點(diǎn)兒被人妖公主榨干…… 1. 我叫王不留,地道東北人咙边。 一個(gè)月前我還...
    沈念sama閱讀 46,260評(píng)論 2 360
  • 正文 我出身青樓猜煮,卻偏偏與公主長得像,于是被迫代替她去往敵國和親败许。 傳聞我的和親對(duì)象是個(gè)殘疾皇子王带,可洞房花燭夜當(dāng)晚...
    茶點(diǎn)故事閱讀 43,446評(píng)論 2 348

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