使用homer进行差异基因启动子motif富集分析中，报错，运行结果knowResults文件为空。
运行代码
findMotifs.pl gene_list.txt athaliana -len 12 -start -800 -end 200
geneLlist.txt文件为单列基因ID文件
报错信息：
Progress: Step7 - Known motif enrichment
Reading input files...
0 total sequences read
506 motifs loaded
Cache length = 11180
Using hypergeometric scoring
Checking enrichment of 506 motif(s)
|0% 50% 100%|
=================================================================================
Illegal division by zero at /home/smx-z3/bin/findKnownMotifs.pl line 152.
Progress: Step8 - De novo motif finding (HOMER)
-blen automatically set to 2
Scanning input files...
!!! Something is wrong... are you sure you chose the right length for motif finding?
!!! i.e. also check your sequence file!!!
Use of uninitialized value in numeric gt (>) at /home/bin/compareMotifs.pl line 1394.
!!! Filtered out all motifs!!!
Job finished
是哪里出问题了
阅读全文
收起全文
