用HOMER里的scanMotifGenomeWide.pl脚本在全基因组范围找已知的motif

直接使用conda安装就可以
安装好以后scanMotifGenomeWide.pl脚本会在虚拟环境的bin目录下
使用命令是

scanMotifGenomeWide.pl at_Ebox.motif at.col0.chr.fna > at_Ebox.at.output.bed
第一个位置参数是motif文件
motif文件的格式是

序列对应的是出现概率最大的碱基
下面数字每列对应的碱基顺序应该是 ACGT
每行对应的是序列每个位置4个不同的碱基出现的概率
第二个位置参数对应的是基因组序列
最后接输出文件
输出是bed文件

这个运行速度很快
只要是自己感兴趣的序列,整理成这种motif的格式,可以直接利用这个脚本定位在基因组上的位置

声明:文中观点不代表本站立场。本文传送门:https://eyangzhen.com/424450.html

联系我们
联系我们
分享本页
返回顶部