Fuzzy DNA Search in Seqqio: Substitutions, Indels, and Strand Hits
Find approximate DNA matches with substitutions alone or substitutions plus indels, then interpret reverse-strand coordinates, alignments, and CIGAR operations.
The two modes are separate scientific profiles: Seqqio-FuzzyDna-Hamming-v1 and Seqqio-FuzzyDna-Edit-v1. Both can search the forward target, its reverse complement, or both. Hits may overlap. The tool reports sequence similarity under the chosen edit rule; it does not provide an E-value, genomic uniqueness score, primer-specificity verdict, or local-alignment significance.
Choose substitutions or full edit distance
| Mode | Allowed differences | Target span |
|---|---|---|
| Substitutions | Base replacements only | Always the same length as the query |
| Edit distance | Replacements, target insertions, and query deletions | Can be shorter, equal to, or longer than the query |
Open Fuzzy Search DNA, paste a query from 1 to 4,096 bases, select the mode, and set zero to 64 edits, always fewer than the query length. Choose Forward, Reverse complement, or Both, then paste one target or an ordered FASTA batch. Only A, C, G, and T are accepted; ambiguity is not treated as a free mismatch.
Reproduce a two-strand Hamming search
Use query ACG, target TTACGTT, Substitutions mode, zero edits, and Both strands. The executed Seqqio run returned two exact hits. Hit 1 is forward at bases 3-5. Hit 2 is reverse-complement at bases 4-6. Each has CIGAR 3=, zero substitutions, and an aligned query and target of ACG.
This small example demonstrates why strand selection belongs in the result. Coordinates are always one-based and refer to the original target record, even when the aligned target is shown in reverse-complement orientation. A palindromic occurrence can legitimately appear as a separate hit on each requested strand.
See how indels change the target span
Now use query ACGT, target TTACGGTAA, Edit distance, one edit, and Both strands. The engine returned four hits: three forward and one reverse. The forward hits cover bases 3-5 with 3=1D, bases 3-6 with 3=1X, and bases 3-7 with 2=1I2=. The reverse hit covers bases 3-6 with 1X3=.
| Operation | Meaning | Effect on target span |
|---|---|---|
= | Equal query and target base | Consumes both |
X | Substituted base | Consumes both |
I | Base inserted in the target relative to the query | Consumes target only |
D | Query base deleted relative to the target | Consumes query only |
For the five-base forward hit, 2=1I2= aligns query AC-GT to target ACGGT. The edit count is the number of nucleotide operations, not the number of gap runs. Select a hit in Seqqio to inspect the complete aligned strings, difference highlights, original coordinates, and record identity.
Understand the edit-distance result set
Edit-distance mode uses a full-query Levenshtein objective with a free target prefix and reports one stable best nonempty suffix for each endpoint. It is not an enumeration of every tied start interval, a Smith-Waterman local alignment, or a database-search significance model. A zero-hit run is a successful completed search, while a work-budget error is reported separately and never presented as zero evidence.
Seqqio's independent validation compared edit scores and paths with Edlib, whose official project documentation describes global, prefix, and infix edit-distance modes and path recovery. Reverse-complement behavior was also checked against the Biopython sequence API. The Seqqio profile keeps its own documented endpoint and tie-breaking rules.
When fuzzy DNA search is a good fit
Use this tool to inspect a motif with a bounded number of substitutions or indels, demonstrate strand-aware approximate matching, audit coordinate-specific hits, or process local FASTA records with a reproducible rule. Use a genome-scale aligner or dedicated off-target service when you need indexed reference search, affine gaps, mapping quality, statistical significance, or domain-specific specificity decisions.
Fuzzy Search DNA is one of 39 applications in the Seqqio Windows 64-bit workspace. The Seqqio overview explains its local record, history, batch, and export model. The complete toolkit is offered as a US$99 one-time purchase with no activation key.
References
- Sosic M, Sikic M and contributors. Edlib Official Edlib project documentation Independent exact edit-distance and alignment-path implementation used during validation.
- Biopython contributors. Bio.Seq module Official Biopython API documentation Independent reverse-complement sequence behavior used during validation.