Find Protein Motifs in Seqqio: Patterns, Regular Expressions, and FASTA Batches
Find an exact amino-acid string or a bounded sequence pattern while retaining every FASTA record. Worked examples distinguish hit count, covered residues and biological interpretation.
A useful motif search keeps three things separate: the definition of the query, the positions that match it, and the interpretation of those matches. This tutorial uses short artificial sequences so that you can inspect every result directly before applying the same settings to a larger batch.
Run an exact-string control first
>two_pattern_regions
MAGDAAGKGDVVG
>no_match_control
MACDEFGHIK
Open Protein Pattern Find, enter the FASTA batch, set Pattern syntax to Literal motif, and use GDAAG as the query pattern. The first record contains one match at residues 3–7. The second record completes successfully with no matches.
A successful no-hit record is a result, not a failed calculation. Keep it in the batch report. Dropping no-hit rows changes the denominator when you later describe how many sequences contain the pattern. Conversely, a record that cannot be evaluated should retain its error status rather than being counted as a successful negative.
Allow a defined amount of variation with bounded regex
Change Pattern syntax to Bounded regex and use GD[A-Z]{2}G. The query specifies G, then D, then exactly two amino-acid letters, then G. Because the accepted input alphabet is the canonical 20 amino acids, this example searches five-residue strings within that alphabet.
| Record | Matches, one-based inclusive | Hit count | Covered residues |
|---|---|---|---|
| two_pattern_regions | 3–7: GDAAG; 9–13: GDVVG | 2 | 10 |
| no_match_control | None | 0 | 0 |
The fixed letters and bounded length are part of your hypothesis. Broadening a query can increase matches without increasing biological specificity. Document the exact pattern in the exported analysis rather than summarizing it later as “the GD motif.” A reader should be able to reconstruct which strings were eligible.
Decide whether overlapping starts belong in the analysis
Use AAAAA as a separate sequence and AAA as a literal motif. With Include overlapping starts enabled, the executed result contains hits at 1–3, 2–4 and 3–5. With it disabled, this example contains only the first hit at 1–3.
| Setting | Hits | Covered residues |
|---|---|---|
| Overlapping starts enabled | 3 | 5 |
| Overlapping starts disabled | 1 | 3 |
Three matches of length three cover five distinct residues, not nine. Use the appropriate quantity for the question: hit count describes reported occurrences, while coverage describes positions included in at least one occurrence. For variable-length expressions, including overlapping starts does not promise that every possible alternative parse at a single start is enumerated.
Keep query errors separate from biological negatives
The bounded-regex mode deliberately rejects unbounded patterns such as GD.*G. The executed example returned unbounded_protein_pattern. The query A{0,2} can match an empty string and returned empty_protein_pattern_match. Rewrite these queries with a finite, positive length appropriate to the question before interpreting any result.
The protein input accepts the canonical 20 amino-acid letters. An X or another unsupported residue symbol is not an automatic wildcard. Resolve ambiguous inputs from their source or retain their exclusions explicitly. Do not silently delete letters, because doing so also changes downstream coordinates.
Know when a sequence pattern is the wrong model
ScanProsite is a separate resource for searching PROSITE patterns and profiles. PROSITE notation is not automatically interchangeable with Seqqio’s bounded regular expressions. If you need an established family signature or profile-based analysis, use the appropriate documented resource and inspect its evidence and scoring model.
Seqqio’s local pattern search is useful when you already have a clearly defined literal string or bounded sequence rule and want to inspect its positions across your own records. It does not perform a profile-HMM search, establish structural context, or turn a short shared string into evidence of common function.
Save a complete and reproducible search
Retain the input FASTA, exact query, syntax mode, overlap setting, Seqqio version and exported TSV. Reconcile successful hit records, successful no-hit records and failed records with the original record count. Use stable identifiers when joining the report to annotations or an experimental sample table.
Protein Pattern Find is part of Seqqio’s Windows 64-bit sequence workspace. It fits a researcher who wants a local, inspectable batch search with retained settings and exports. The product page describes the available toolkit and purchase options.
References
- SIB Swiss Institute of Bioinformatics. ScanProsite PROSITE A separate tool for PROSITE patterns/profiles and user-defined protein searches.