This is a test version of Biostars. For the public version, visit https://www.biostars.org.
Is there a code to find consensus motif

Hello,

I've been trying to write a code to find a consensus motif in a given sequence, and for this purpose, I was only able to reach till finding a substring in a string. I want to be able to allot multiple nucleotides/amino acids at each position, and also enter N/X representing any of the nucleotides/amino acids. I would very much appreciate any help.

Thanks.

P.S. The post tags represent the languages I'm comfortable understanding.

Edit: Example of the consensus motif - A/T A A G C A A/T/G N N A

Sequence - CGATCGTG TAAGCAGCTA GTCATG

Bolded sequence is the consensus

c awk shell python
• 3,650 views • link
updated by Carlo Yague Carlo Yague • written by vinayjrao vinayjrao

1 answer

In shell using grep and regular expressions:

echo 'CGATCGTG TAAGCAGCTA GTCATG' | grep  -o "[AT]AAGCA[ATG]..A"
TAAGCAGCTA

'N' is expressed as '.', meaning that it can take any value. Multiple nucleotides at one position are put into square brackets.

• 1 views • link
written by Carlo Yague Carlo Yague

Thanks a lot. It's perfect.

• 1 views • link
written by vinayjrao vinayjrao

In the same lines of Carlo Yague

echo 'CGATCGTG TAAGCAGCTA GTCATG' | grep -Po \([AT]\)A{2}GCA[\1G].{2}A
TAAGCAGCTA
• 1 views • link
written by cpad0112 cpad0112

Thanks. This works too. I could use the .{2} when I have larger repeats of any nucleotide/amino acid. Although, I would like to know why it [\1G] and not [ATG]?

• 1 views • link
written by vinayjrao vinayjrao

The first AT is made a group and every time and anywhere you can call it by its serial number (1 here)

• 1 views • link
written by cpad0112 cpad0112

That's an extremely handy option. Thanks again :)

• 1 views • link
written by vinayjrao vinayjrao

Log in to answer this question.

More posts like this