trimal的使用

之前使用gblocks裁剪比對序列,感覺裁剪的不是很智能,大量可能有用的序列被裁剪掉,所以嘗試使用trimal。

1.下載安裝

在自己目錄下輸入

git clone https://github.com/scapella/trimal.git # 安裝trimal

導入trimal安裝包(如果服務器連不上網也可以在這個網站下載安裝包導入服務器)
然后unzip解壓安裝包之后編譯軟件

cd trimal/source
make

使用

trimal常用命令如下:

 trimal -in example1 -out output1 -automated1

-in <inputfile>多種格式的輸入文件(clustal,fasta,NBRF / PIR,nexus,phylip3.2,phylip)。
-compareset <inputfile>包含要比較的路線的文件的路徑輸入列表。
-matrix <inpufile>用戶定義的相似性矩陣的輸入文件(默認為Blosum62)。

-out <outputfile>以相同輸入格式輸出對齊(默認標準輸出)。(默認輸入格式)
-htmlout <outputfile>獲取trimal在HTML文件中工作的摘要。

-CUSTAL格式的-clustal輸出文件
-fasta FASTA格式的輸出文件
-nbrf NBRF / PIR格式的輸出文件
-nexus NEXUS格式的輸出文件
-mega MEGA格式的輸出文件
-phylip3.2 PHYLIP3.2格式的輸出文件
-phylip PHYLIP / PHYLIP4格式的輸出文件

互補獲得互補對齊。
-colnumbering獲取舊對齊和新對齊中列之間的關系。

-select {n,l,mk}選擇要從對齊中刪除的列。(請參閱用戶指南)。
-gt -gapthreshold <n> 1 - (允許間隙的序列的分數)。
-st -simthreshold <n>允許的最小平均相似度。
-ct -conthreshold <n>允許的最小一致性值。
-cons <n>保留原始路線中位置的最小百分比。

-nogaps刪除對齊中有間隙的所有位置。
-noallgaps刪除僅由間隙組成的列。

-gappyout在“gappyout”模式下使用自動選擇。該方法僅使用基于間隙分布的信息。(請參閱用戶指南)。
-strict在“嚴格”模式下使用自動選擇。(請參閱用戶指南)。
-strictplus在“strictplus”模式下使用自動選擇。(請參閱用戶指南)。

-automated1使用基于相似性統計的自動方法的啟發式選擇。(請參閱用戶指南)。

-resoverlap位置與列中其他位置的最小重疊被視為“好位置”。(請參閱用戶指南)。
-seqoverlap序列必須具有的“良好位置”的最小百分比才能被保存。(請參閱用戶指南)。

-w <n>(一半)窗口大小,位置i的得分是窗口(i-n)到(i + n)的平均值。
-gw <n>(half)窗口大小僅適用于基于Gaps的統計/方法。
-sw <n>(half)窗口大小僅適用于基于相似性的統計/方法。
-cw <n>(half)窗口大小僅適用于基于Consistency的統計信息/方法。

-sgc打印輸入對齊中列的間隙百分比計數。
-sgt打印累積的間隙百分比計數。
-scc輸入對齊中列的保留值。
-sct打印累計保護值計數。
-sfc從比較文件方法打印所選對齊中列的比較值。
-sft從比較文件方法打印所選對齊的累計比較值計數

?著作權歸作者所有,轉載或內容合作請聯系作者
平臺聲明:文章內容(如有圖片或視頻亦包括在內)由作者上傳并發布,文章內容僅代表作者本人觀點,簡書系信息發布平臺,僅提供信息存儲服務。
  • 序言:七十年代末,一起剝皮案震驚了整個濱河市,隨后出現的幾起案子,更是在濱河造成了極大的恐慌,老刑警劉巖,帶你破解...
    沈念sama閱讀 227,837評論 6 531
  • 序言:濱河連續發生了三起死亡事件,死亡現場離奇詭異,居然都是意外死亡,警方通過查閱死者的電腦和手機,發現死者居然都...
    沈念sama閱讀 98,196評論 3 414
  • 文/潘曉璐 我一進店門,熙熙樓的掌柜王于貴愁眉苦臉地迎上來,“玉大人,你說我怎么就攤上這事。” “怎么了?”我有些...
    開封第一講書人閱讀 175,688評論 0 373
  • 文/不壞的土叔 我叫張陵,是天一觀的道長。 經常有香客問我,道長,這世上最難降的妖魔是什么? 我笑而不...
    開封第一講書人閱讀 62,654評論 1 309
  • 正文 為了忘掉前任,我火速辦了婚禮,結果婚禮上,老公的妹妹穿的比我還像新娘。我一直安慰自己,他們只是感情好,可當我...
    茶點故事閱讀 71,456評論 6 406
  • 文/花漫 我一把揭開白布。 她就那樣靜靜地躺著,像睡著了一般。 火紅的嫁衣襯著肌膚如雪。 梳的紋絲不亂的頭發上,一...
    開封第一講書人閱讀 54,955評論 1 321
  • 那天,我揣著相機與錄音,去河邊找鬼。 笑死,一個胖子當著我的面吹牛,可吹牛的內容都是我干的。 我是一名探鬼主播,決...
    沈念sama閱讀 43,044評論 3 440
  • 文/蒼蘭香墨 我猛地睜開眼,長吁一口氣:“原來是場噩夢啊……” “哼!你這毒婦竟也來了?” 一聲冷哼從身側響起,我...
    開封第一講書人閱讀 42,195評論 0 287
  • 序言:老撾萬榮一對情侶失蹤,失蹤者是張志新(化名)和其女友劉穎,沒想到半個月后,有當地人在樹林里發現了一具尸體,經...
    沈念sama閱讀 48,725評論 1 333
  • 正文 獨居荒郊野嶺守林人離奇死亡,尸身上長有42處帶血的膿包…… 初始之章·張勛 以下內容為張勛視角 年9月15日...
    茶點故事閱讀 40,608評論 3 354
  • 正文 我和宋清朗相戀三年,在試婚紗的時候發現自己被綠了。 大學時的朋友給我發了我未婚夫和他白月光在一起吃飯的照片。...
    茶點故事閱讀 42,802評論 1 369
  • 序言:一個原本活蹦亂跳的男人離奇死亡,死狀恐怖,靈堂內的尸體忽然破棺而出,到底是詐尸還是另有隱情,我是刑警寧澤,帶...
    沈念sama閱讀 38,318評論 5 358
  • 正文 年R本政府宣布,位于F島的核電站,受9級特大地震影響,放射性物質發生泄漏。R本人自食惡果不足惜,卻給世界環境...
    茶點故事閱讀 44,048評論 3 347
  • 文/蒙蒙 一、第九天 我趴在偏房一處隱蔽的房頂上張望。 院中可真熱鬧,春花似錦、人聲如沸。這莊子的主人今日做“春日...
    開封第一講書人閱讀 34,422評論 0 26
  • 文/蒼蘭香墨 我抬頭看了看天上的太陽。三九已至,卻和暖如春,著一層夾襖步出監牢的瞬間,已是汗流浹背。 一陣腳步聲響...
    開封第一講書人閱讀 35,673評論 1 281
  • 我被黑心中介騙來泰國打工, 沒想到剛下飛機就差點兒被人妖公主榨干…… 1. 我叫王不留,地道東北人。 一個月前我還...
    沈念sama閱讀 51,424評論 3 390
  • 正文 我出身青樓,卻偏偏與公主長得像,于是被迫代替她去往敵國和親。 傳聞我的和親對象是個殘疾皇子,可洞房花燭夜當晚...
    茶點故事閱讀 47,762評論 2 372

推薦閱讀更多精彩內容