Fimo扫描基因组序列
2020-02-18 本文已影响0人
生信编程日常
在做motif分析时,经常用Fimo扫描基因组序列得到motif对应的序列位置,进而进行下一步的分析。说明文档可参考:http://meme-suite.org/doc/fimo.html
![](https://img.haomeiwen.com/i20297934/db30b8e5af72a240.png)
fimo的基本用法为:
fimo [options] <motif file> <sequence file>
motif file就是motif文件,提供MEME and DREME输出的文件就可以,也可以从CISBP,JASPAR或者HOCOMOCO等数据库下载得到。比如如下格式:
![](https://img.haomeiwen.com/i20297934/b9f81ff72645c584.png)
sequence file是序列文件,用全基因组还是提出来的基因组片段都可以。
其他参数还有:
![](https://img.haomeiwen.com/i20297934/a842febadef77c64.png)
用到的比较多的还有--thresh,即p值的阈值,默认是1e-4;--qv--thresh可以修改q值的阈值。
默认会输出fimo_out/的文件夹,里边有fimo.html, fimo.tsv, fimo.gff, cisml.xml和fimo.xml五种不同格式的结果。如果选择--text则只会出现tsv格式的结果。
欢迎关注~
![](https://img.haomeiwen.com/i20297934/51f7609ac7b56bc3.jpg)