Pairwise DNA Alignment in Seqqio: Compare Constructs, Gaps, and Sequence Coverage

This tutorial compares one 16-base synthetic reference with an unchanged copy, a one-base substitution, a one-base insertion, and an eight-base fragment. The sequences are short so every result can be inspected directly. The same decisions matter when you compare longer local constructs: what span was assessed, how gaps were scored, and which original sequence each result belongs to.

Open Pairwise DNA, select Single pair, and put the reference in input A. Put the unchanged copy in input B for the first run. Select Global · end-to-end, set Gap opening to 10 and Gap extension to 0.5, confirm the displayed matrix is NUC.4.4, then select Align pair.

The two FASTA entries above describe the inputs; for Single pair, paste or load one record into each side. Both sequences retain their supplied orientation. Seqqio does not automatically try a reverse complement. If a query was supplied from the opposite strand, resolve its orientation explicitly before judging differences.

Gap costs in the form are positive penalties. A gap of length k costs opening + (k - 1) × extension, so a one-base gap costs 10 and a four-base gap costs 11.5. Global mode also penalizes terminal gaps. The official Biopython pairwise-alignment tutorial explains how global and local objectives and affine gap parameters affect an alignment. Keep the settings fixed when comparing this tutorial's scores.

Keep A unchanged and replace B with each variant below. The substitution changes reference base 12 from A to C. The insertion adds an A to the reference string. These are original software controls, not experimentally characterized constructs.

Both input coverages are 100% for these three comparisons. Seqqio's identity is identical columns divided by all alignment columns, including gap columns. The insertion therefore gives 16/17, rounded to 94.12%, rather than 100%. Its alignment was:

The gap preserves all bases from both original sequences. Nearby repeated A bases can allow equivalent gap placements, so do not turn one displayed optimum into proof of a uniquely resolved insertion coordinate. The score evaluates a path under the chosen model; it is not a probability of relatedness, an E-value, or experimental evidence.

Replace B with TGCACTGA, the eight-base segment at reference bases 5–12. Run Global first, then switch only Alignment mode to Local · best matching region and align again.

Global includes every original base, padding the shorter B sequence with four terminal gaps on each side. Eight matches contribute 40, and the two four-base gap runs cost 11.5 each, giving 17. Local selects the eight matching bases without those flanks, giving 40. The Local result answers where the fragment fits; it says nothing about the unaligned half of A.

Coverage measures the original source span included in the selected alignment divided by that source's full length. It is not the fraction of alignment columns containing a match. Read the two source coverages separately when input lengths differ. Comparing full candidate constructs usually calls for an end-to-end question; locating a known partial fragment calls for a regional question.

Pairwise DNA accepts supported IUPAC DNA symbols, but a literal symbol match is not necessarily a confirmed base match. In our control, Global alignment of NNNN with NNNN had 100% literal identity and score −4 under NUC.4.4. Those four N symbols do not tell you which concrete bases are present.

For the same N-only inputs, Local returned score 0, an empty alignment, zero coverage, and no defined identity percentage. An empty best region is a valid outcome under this profile. Preserve that distinction instead of reporting “0% identity,” substituting concrete bases, or treating the result as a calculation failure.

Once the single-pair controls make sense, choose Batch. Reference A vs every query in B compares one selected A record independently with each B record. A1 vs B1, A2 vs B2, … pairs two equally sized FASTA inputs by original record order. Matching titles do not define paired-batch partners. Confirm order and identifiers before running.

Save the selected pair's aligned FASTA and the dataset report. Keep the original inputs, matrix, mode, gap costs, method profile Seqqio-PairwiseDna-Affine-v1, source coordinates, and both coverage values with your review. Batch pairing organizes separate comparisons; it does not create a multiple sequence alignment.

Use the existing Global Alignment guide and Local Protein Alignment guide for the conceptual background. If your question is instead “where does this full DNA motif occur with a bounded number of edits?”, the Seqqio Fuzzy DNA Search tutorial addresses that different result set.

Pairwise DNA fits local construct review, variant demonstrations, and comparisons between a reference and supplied fragments when an explicit scoring model is appropriate. Local mode returns one optimal region, not all possible matching sites. Use a dedicated reference-search or read-mapping workflow when you need indexed genome search, mapping quality, database significance, or large collections of reads. Per-tool work limits still apply.

Pairwise DNA is included in Seqqio's 39-application Windows 64-bit toolkit. The current offer is a US$99 one-time purchase with no activation key. Review the platform requirements and final checkout total on the product page.

Explore Seqqio and the US$99 toolkit