Augustus

在基因組注釋中,預(yù)測基因結(jié)構(gòu)是最核心的一環(huán),在眾多的從頭預(yù)測軟件中第煮,augustus應(yīng)用較廣。但是由于安裝過程依賴軟件較多抑党,如果沒有管理員權(quán)限一般不要輕易嘗試包警,有興趣可以參考 基因預(yù)測軟件中的奧古斯都

1. 腳本運行

$ augustus --strand=both --genemode=partial --singlestrand=false --hintsfile=hints.gff --extrinsicCfgFile=extrinsic.cfg --protein=on --introns=on --start=on --stop=on --cds=on --codingseq=on --alternatives-from-evidence=true --gff3=on --UTR=on ----outfile=out.gff --species=human genome.fa

$ augustus --noprediction=true --species=SPECIES sequences.gb

參數(shù)解釋:

--strand=both,--strand=forward or --strand=backward 表示注釋基因在兩條鏈還是其中一條底靠;

--genemodel=partial, --genemodel=intronless, --genemodel=complete,--genemodel=atleastone or --genemodel=exactlyone

partial :?允許在序列邊界預(yù)測不完整的基因(默認值) ? ?intronless :?只預(yù)測單外顯子基因害晦,如在原核生物和一些真核生物中? ? complete :?只預(yù)測完整基因 ? atleastone :?預(yù)測至少一個完整的基因 ? ?exactlyone :?準確預(yù)測一個完整的基因;

--singlestrand=true?獨立預(yù)測每條鏈上的基因暑中,允許在相反的鏈上有重疊的基因壹瘟。這個選項在默認情況下是關(guān)閉的;

--hintsfile=hintsfilename?當使用此選項時鳄逾,將打開考慮提示(外部信息)的預(yù)測稻轨。hintsfilename包含gff格式的提示;

--extrinsicCfgFile=cfgfilename?可選的雕凹。此文件包含用于提示及其boni和mali的使用源列表殴俱。如果沒有指定配置目錄中的extrinsic.cfg文件政冻,則使用$AUGUSTUS配置路徑;(不太懂)

--maxDNAPieceSize=n?此值指定是為了運行核心算法(Viterbi)而將序列分割成的片段的最大長度线欲。默認的是——maxDNAPieceSize = 200000明场;

--protein=on/off

--introns=on/off

--start=on/off

--stop=on/off

--cds=on/off

--codingseq=on/off

輸出選項。輸出預(yù)測蛋白序列李丰,內(nèi)含子苦锨,起始密碼子,終止密碼子趴泌≈凼妫或者在“初始”、“內(nèi)部”嗜憔、“終端”和“單外顯子”之外使用“cds”魏蔗。cds不包括停止密碼子(除非stopCodonExcludedFromCDS=false),而終端和單個外顯子包含停止密碼子;

--AUGUSTUS_CONFIG_PATH=path?配置config的路徑(如果沒有就指定為環(huán)境變量)痹筛;

--gff3=on/off 輸出gff3格式莺治;

--proteinprofile=filename 讀取蛋白序列;

--predictionStart=A, --predictionEnd=B 針對某一特定區(qū)域進行注釋帚稠,可以快速注釋某一感興趣的區(qū)域谣旁;

--UTR=on/off?預(yù)測除編碼序列外的未翻譯區(qū)域。目前滋早,這只適用于人類榄审,galdieria, toxopl asma和caenorhabditis;

--outfile=filename?打印輸出到文件名而不是標準輸出杆麸;

--noInFrameStop=true/false?不要返回帶有終止密碼子的轉(zhuǎn)錄本搁进。否則,可能會發(fā)生內(nèi)含子跨越終止密碼子昔头。默認值:假饼问;(不太懂)

--noprediction=true/false?Augustus也可以以genebank格式文件為輸入文件,進行基因預(yù)測揭斧,并將預(yù)測結(jié)果和genebank的結(jié)果進行比較后得出一個精確性的統(tǒng)計結(jié)果莱革。

?著作權(quán)歸作者所有,轉(zhuǎn)載或內(nèi)容合作請聯(lián)系作者
  • 序言:七十年代末,一起剝皮案震驚了整個濱河市讹开,隨后出現(xiàn)的幾起案子盅视,更是在濱河造成了極大的恐慌,老刑警劉巖旦万,帶你破解...
    沈念sama閱讀 206,839評論 6 482
  • 序言:濱河連續(xù)發(fā)生了三起死亡事件闹击,死亡現(xiàn)場離奇詭異,居然都是意外死亡成艘,警方通過查閱死者的電腦和手機赏半,發(fā)現(xiàn)死者居然都...
    沈念sama閱讀 88,543評論 2 382
  • 文/潘曉璐 我一進店門梅忌,熙熙樓的掌柜王于貴愁眉苦臉地迎上來,“玉大人除破,你說我怎么就攤上這事∏砬唬” “怎么了瑰枫?”我有些...
    開封第一講書人閱讀 153,116評論 0 344
  • 文/不壞的土叔 我叫張陵,是天一觀的道長丹莲。 經(jīng)常有香客問我光坝,道長,這世上最難降的妖魔是什么甥材? 我笑而不...
    開封第一講書人閱讀 55,371評論 1 279
  • 正文 為了忘掉前任盯另,我火速辦了婚禮,結(jié)果婚禮上洲赵,老公的妹妹穿的比我還像新娘鸳惯。我一直安慰自己,他們只是感情好叠萍,可當我...
    茶點故事閱讀 64,384評論 5 374
  • 文/花漫 我一把揭開白布芝发。 她就那樣靜靜地躺著,像睡著了一般苛谷。 火紅的嫁衣襯著肌膚如雪辅鲸。 梳的紋絲不亂的頭發(fā)上,一...
    開封第一講書人閱讀 49,111評論 1 285
  • 那天腹殿,我揣著相機與錄音独悴,去河邊找鬼。 笑死锣尉,一個胖子當著我的面吹牛刻炒,可吹牛的內(nèi)容都是我干的。 我是一名探鬼主播自沧,決...
    沈念sama閱讀 38,416評論 3 400
  • 文/蒼蘭香墨 我猛地睜開眼落蝙,長吁一口氣:“原來是場噩夢啊……” “哼!你這毒婦竟也來了暂幼?” 一聲冷哼從身側(cè)響起筏勒,我...
    開封第一講書人閱讀 37,053評論 0 259
  • 序言:老撾萬榮一對情侶失蹤,失蹤者是張志新(化名)和其女友劉穎旺嬉,沒想到半個月后管行,有當?shù)厝嗽跇淞掷锇l(fā)現(xiàn)了一具尸體,經(jīng)...
    沈念sama閱讀 43,558評論 1 300
  • 正文 獨居荒郊野嶺守林人離奇死亡邪媳,尸身上長有42處帶血的膿包…… 初始之章·張勛 以下內(nèi)容為張勛視角 年9月15日...
    茶點故事閱讀 36,007評論 2 325
  • 正文 我和宋清朗相戀三年捐顷,在試婚紗的時候發(fā)現(xiàn)自己被綠了荡陷。 大學(xué)時的朋友給我發(fā)了我未婚夫和他白月光在一起吃飯的照片。...
    茶點故事閱讀 38,117評論 1 334
  • 序言:一個原本活蹦亂跳的男人離奇死亡迅涮,死狀恐怖废赞,靈堂內(nèi)的尸體忽然破棺而出,到底是詐尸還是另有隱情叮姑,我是刑警寧澤唉地,帶...
    沈念sama閱讀 33,756評論 4 324
  • 正文 年R本政府宣布,位于F島的核電站传透,受9級特大地震影響耘沼,放射性物質(zhì)發(fā)生泄漏。R本人自食惡果不足惜朱盐,卻給世界環(huán)境...
    茶點故事閱讀 39,324評論 3 307
  • 文/蒙蒙 一群嗤、第九天 我趴在偏房一處隱蔽的房頂上張望。 院中可真熱鬧兵琳,春花似錦狂秘、人聲如沸。這莊子的主人今日做“春日...
    開封第一講書人閱讀 30,315評論 0 19
  • 文/蒼蘭香墨 我抬頭看了看天上的太陽。三九已至羡榴,卻和暖如春碧查,著一層夾襖步出監(jiān)牢的瞬間,已是汗流浹背校仑。 一陣腳步聲響...
    開封第一講書人閱讀 31,539評論 1 262
  • 我被黑心中介騙來泰國打工忠售, 沒想到剛下飛機就差點兒被人妖公主榨干…… 1. 我叫王不留,地道東北人迄沫。 一個月前我還...
    沈念sama閱讀 45,578評論 2 355
  • 正文 我出身青樓稻扬,卻偏偏與公主長得像,于是被迫代替她去往敵國和親羊瘩。 傳聞我的和親對象是個殘疾皇子泰佳,可洞房花燭夜當晚...
    茶點故事閱讀 42,877評論 2 345

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