How toconvert a set of motif consensus sequences in IUPAC to meme format?
1
1
Entering edit mode
20 months ago

I am trying to run FIMO from MEME Suite to to search for occurrences of a list of motifs in my DNA sequences. The list of motifs are present in IUPAC format, given by their consensus sequences. I have a text file containing the motif name followed by its consensus sequence or Sequence logo as given below:

>AHR_HUMAN.H11MO.0.B
dKhGCGTGh
>ALX1_HUMAN.H11MO.0.B
vTRATYGnATTA
>ANDR_HUMAN.H11MO.0.A
WKThYYddbhTRTTTRYh
>ANDR_HUMAN.H11MO.1.A
dGnACWbTbWGTdCYb

and so on.

I am trying to run FIMO, which expects the motifs in meme format as input.

Can anyone tell me how to convert all the motifs to meme format in one step? I have about 100 motifs I need to scan, and I don't want to scan one motif at a time.

In the online server of Meme Suite, I can just paste the list of motifs (as given above) and it works. But that is awfully slow and I want to run it on my own server.

MEME FIMO motif transfac meme suite • 890 views
ADD COMMENT
2
Entering edit mode
20 months ago

They have a utility for that and converting from just about any other format as well..

ADD COMMENT
0
Entering edit mode

Yes there is the iupac2meme utility.

iupac2meme Usage: iupac2meme [options] <consensus_sequence>+ It works fine if I type in a single motif sequence at a time. How do I use this for multiple motifs in one go?

ADD REPLY
1
Entering edit mode

Using a for or while loop. See this stackoverflow answer for details.

ADD REPLY

Login before adding your answer.

Traffic: 1310 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6