New answers tagged sequence-analysis
1
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 ...
0
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.
Top 50 recent answers are included