How do I perform a motif search on a specific fasta sequence?

I have an annotated genome file, and I have a potential ncRNA that could be an eRNA.

I have extracted this sequence to a fasta file, and want to use homer for scanning transcription factor binding sites (i.e. to work out whether it is a potential enhancer).

How can I do this? I have tried with this command but it comes up with the error:

 perl findMotifs.pl input.fa fasta motifResults/ -fasta genome.fa 

Scanning input files... !!! Something is wrong... are you sure you
chose the right length for motif finding? !!! i.e. also check your
sequence file!!!

Scanning input files... !!! Something is wrong... are you sure you
chose the right length for motif finding? !!! i.e. also check your
sequence file!!!

homer enhancer

• 4,337 views • link

updated 7.9 years ago by Ram

• written 7.9 years ago by a.rex

Can you post an extract of the fasta file(s) you're using?

• 1 views • link

written 7.9 years ago by lieven.sterck

If you already have the fasta file then I would just use MEME-ChIP; http://meme-suite.org/tools/meme-chip

• 1 views • link

written 7.9 years ago by steve

0 answers

No answers yet.

More posts like this

homer findMotifsGenome.pl
written 2.6 years ago by biology_inform •

Hi all, I have an ATAC-seq data and from differential expressed peaks I would like to find any motif that can be useful for me. …
Problem with homer in findMotifs.pl when using input and bg fasta
written 3.4 years ago by Ilario •

I am using the following command: findMotifs.pl input.fa fasta ./Output -fastaBg bg.fa -len 8,10,12 -norevopp The input and bg fasta have this structure and the …
HOMER and .fa files error
written 4.5 years ago by c.eskiw •

Good day all, I have been trying HOMER to identify txn factor motifs enriched in my ATACseq data sets. I have been using the narrow …
question about homer
written 5.5 years ago by 13719756608 •

hello,everyone i was looking for motif for my m5c-seq bed file,and i encounter some problem this is my bed file: MACS_peak_18032 Chr22 2938644 2938798 . …
HOMER: ./findMotifsGenome.pl
written 5.9 years ago by radicav06 •

Hello all! I want find motifs FOXA1 in the complete human genome. I have a file with peaks 10_FOXA1.bed. Format file: chr1 629136 629408 peak1 …
Can't get known motif enrichment result using findMotifs.pl (Homer)
written 6.6 years ago by youllae •

Hi. I'm currently using HOMER to see known motif enrichment of the list of DEGs I have. Although I followed all the instruction provided in …
Discrepancy in HOMER output
written 7.3 years ago by doan.dinh •

Hello, I need some help understanding how HOMER motif finding works. I have a ChIP-seq data for a transcription factor that I want to look …
problem with HOMER findMotifsGenome.pl
written 7.8 years ago by ste.lu

Hi All, I did and ATAC-Seq experiment in different cell lines and I was curious to see if they have different motifs in the open …
Homer error (findMotifs.pl) : are you sure you chose the right length for motif finding?
written 8.0 years ago by salamandra

I did a ChIP-seq analysis with a genome GRCh38.92 from ensembl and also used the corresponding gene anotation file (.gtf) from ensembl. Genome, annotation and …
findGenomeMotif.pl in Homer couldn't work properly
written 8.9 years ago by xuenjun1 •

./Homer/bin/findMotifsGenome.pl myc1_vs_input1.txt tair10 motif_based_on_myc1_vs_input1_cleaned_peaks/ -size 200 -mask Position file = myc1_vs_input1.txt Genome = tair10 Output Directory = motif_based_on_myc1_vs_input1_cleaned_peaks/ Fragment size set to 200 Will use …