how to map a motif in a series of sequences?

0

Entering edit mode

6.5 years ago

haris.zafr • 0

Hi everyone,

I have about 700 sequences of COI gene and I would like to see a certain motif that I also already have, if and if yes, where it hits in each of them. What approaches and tools could I use?

Thank you very much!!

gene • 1.2k views

ADD COMMENT • link 6.5 years ago by haris.zafr • 0

0

Entering edit mode

What is your motif? Amino acids? DNA? Any redunant/ambiguous characters?

ADD REPLY • link 6.5 years ago by Joe 21k

0

Entering edit mode

It is DNA and we have ambiguous characters. Thanks again!! :)

ADD REPLY • link 6.5 years ago by haris.zafr • 0

0

Entering edit mode

The easiest way to do it I would think would be to create a multiple sequence alignment of all of the genes and the motif of interest. If the sequences are too long for this (MSA scales pretty horribly), you may want to consider HMM approaches, or simply making a BLAST database of your 700 sequences and querying your motif against them all - I think BLAST tolerates certain ambiguous characters but I think you're going to have to give us a bit more information to get any further (some example data for instance).

ADD REPLY • link 6.5 years ago by Joe 21k

Login before adding your answer.