align_read
Pairwise-align ONE sequencing read against a reference and hand on the two gapped rows, so the read can be graded. Needleman-Wunsch with free end gaps and the IUPAC substitution matrix — SeqTrace's aligner, the same call the trace editor makes, so a pipeline and the panel align a read identically. Leave reference unset to align against the record this step runs on.
Input Schema
| Name | Required | Description | Default |
|---|---|---|---|
| read | Yes | the read's bases — wire a trace reader's `sequence` output here, or paste them. | |
| _files | No | The project: record references shaped like _subject. Their records become library documents a tool names in sequences/parts/document; a trace or alignment is found by file name. | |
| _subject | No | The record to act on: exactly one of sequence, accession (+source), sample (an id from list_samples) or file ({name, text} or {name, base64}); record_index picks one record of several. | |
| reference | No | what to align against. Unset, it is the record this step runs on (the open plasmid, or an upstream product wired into the subject). |
Output Schema
| Name | Required | Description | Default |
|---|---|---|---|
| columns | Yes | the alignment's width, gaps included | |
| identity | Yes | % identity over the columns where both rows have a base (gaps excluded) | |
| mismatches | Yes | columns where the two rows carry different bases | |
| alignedRead | Yes | the read with gaps inserted — wire into `verify_read.aligned_read` | |
| referenceLength | Yes | the reference's ungapped length in bp — wire into `verify_read.target_len` | |
| alignedReference | Yes | the reference with gaps inserted — wire into `verify_read.aligned_ref` |