Take the 2-minute tour ×
Biology Stack Exchange is a question and answer site for biology researchers, academics, and students. It's 100% free, no registration required.

I have a vector sequence and would like to find the following nucleotide sequence in it.

AASYWSRA

This query sequence uses several degenerate symbols, defined as:

  • S = C or G (Strong)
  • Y = C or T (pYrimidine)
  • W = A or T (Weak)
  • R = A or G (puRine)

(See this page for a complete description of the IUPAC nucleotide codes.)

Is there a tool to align this sequence with the vector and find all the possible sites?

EDIT: I found this website but couldn't find any of the resources helpful. Maybe one of these will work but I found most of them too complicated. If you know which one will work and explain how to use it I will be grateful.

share|improve this question
    
I could've sworn NCBI BLAST can do that? Have you tried it? –  Armatus Apr 30 '13 at 22:12
add comment

3 Answers

If you're looking for an exact match, you don't really need a complex aligner. Perl regular expressions are pretty powerful at string transformations or conditional matching of substrings. For example, to find all matches of AASYWSRA in a nucleotide sequence $seq, you can do:

@matches = $seq =~ m/AA[CG][CT][AT][CG][AG]A/g;

The [] brackets are known as Perl character classes, which allow different characters in a position.

share|improve this answer
add comment

I've never utilized the feature, but I know Vmatch allows you to define your own alphabets and symbol mappings, which should allow you to handle the IUPAC ambiguous nucleotide symbols. The default alphabet/mapping may even include these symbols.

share|improve this answer
add comment

I have used MAST to search for particular motifs in a sequence.

From the MAST intro page, it looks like you'll have to reformat your motif a little, but it should suffice for your purposes:

Motifs are represented as position-dependent scoring matrices that describe the score of each possible letter at each position in the pattern.

Here is more info on formatting, which also indicates that the alphabet you use will be interpreted correctly.

share|improve this answer
add comment

Your Answer

 
discard

By posting your answer, you agree to the privacy policy and terms of service.

Not the answer you're looking for? Browse other questions tagged or ask your own question.