從零開始學生信-meme的使用


上一章講了meme的安裝,本章進行使用的講解,需要快捷食用本章節的請關注公眾號:生物信息分析學習

 

 

meme的使用,網上有很多例子。

1. 網頁版:快捷方便,一台電腦同時可以申請運行四個任務。

先進入網站:http://meme-suite.org/tools/streme

 

 想用什么套件選什么,一目了然。

不懂可以參考手冊:

 

 這是上一章為什么一定想裝streme的原因(因為meme只願意支持50個sequences,雖然多少都能跑,但是好像不太對)。

streme提交的參數如下:

  • streme --verbosity 1 --oc . --dna --p cLN.bed_3.fa --minw 5 --maxw 30 --pvt 0.05 --dfile description --totallength 1e7 --time 18000

只有minw和maxw是自己設置的,minw=5 因為我的文件最短的sequence是5(其實一般設置6或者8),maxw=30 因為網頁版最大支持30.

命令行版本如下:

meme $1 -dna -oc . -nostatus -time 18000 -mod zoops -nmotifs 1000 -minw 5 -maxw 50 -objfun classic -revcomp -markov_order 0 -evt 0.05 -o $1.meme

其中$1為你的序列文件,注意不要出現ATCGN以外的字符(網頁版的會報錯,命令行的沒有嘗試)。

E-value設置為0.05以內,否則結果也不可信,還浪費時間。

=======================================================================================

現在看一下meme和streme的結果:

meme命令行的結果:

 

 streme的網頁版結果:

 

 明顯,meme的結果是異常的,同時我手動搜索過motif序列,meme中只有兩個sites的motif,實際上存在於很多序列中。

 


免責聲明!

本站轉載的文章為個人學習借鑒使用,本站對版權不負任何法律責任。如果侵犯了您的隱私權益,請聯系本站郵箱yoyou2525@163.com刪除。



 
粵ICP備18138465號   © 2018-2025 CODEPRJ.COM