Pairwise Alignment
pairwise_alignmentGlobal (Needleman-Wunsch), local (Smith-Waterman) or semi-global/fitting pairwise alignment of two sequences, with match/mismatch scoring and affine gap costs (Gotoh).
Input Schema
| Name | Required | Description | Default |
|---|---|---|---|
| gap | No | Affine gap EXTEND penalty, charged per gap position (including the first). | |
| mode | No | "global" penalises end gaps in both sequences; "local" returns the best-scoring subalignment; "semiglobal" is a fitting alignment — seqB is consumed end to end while seqA's terminal overhangs are free and are not emitted, so a partial read placed on a longer reference is not smeared across it. | global |
| seqA | Yes | First sequence (raw or FASTA; nucleotide or protein). | |
| seqB | Yes | Second sequence (raw or FASTA; nucleotide or protein). | |
| match | No | Match score. | |
| gapOpen | No | Extra one-off penalty charged on top of gap for a gap's first position. Defaults to 1.5 * gap, so a k-base gap costs gap * (k + 1.5) and one contiguous k-base indel is cheaper than k scattered 1-base gaps. Pass 0 for a purely linear penalty. | |
| mismatch | No | Mismatch penalty. |