load_trace
Read a Sanger trace (.ab1 / .scf / .ztr) from this project and hand its BASE CALLS onward as a record, with the per-base quality and the trim bounds beside them. This is what makes a sequencing read a pipeline input: verify it against a reference, align it, gate on its quality. The trace SIGNAL stays here — a record crosses a wire as bases, and a peak cannot be written as one.
Input Schema
| Name | Required | Description | Default |
|---|---|---|---|
| trim | No | emit the quality-trimmed read rather than the whole call list. The bounds are the sliding-window mean-phred bounds the chromatogram panel trims by, so panel and pipeline cut at the same base. | |
| _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. | |
| document | No | which trace file to read: its name (`read_F.ab1`) or its project file id. A file id is project-specific, so a template ships this unset. | |
| min_quality | No | the mean phred a window must reach to be kept when trimming. Only applies when trim is true. | |
| reverse_complement | No | read it on the other strand. A reverse primer's trace is the reverse complement of the strand you are verifying, and a name ending _R is the convention for it. |
Output Schema
| Name | Required | Description | Default |
|---|---|---|---|
| trace | Yes | what the FILE says about itself: format, instrument comments, call count, trim bounds. Not the signal, which cannot ride a wire | |
| length | Yes | how many bases that is | |
| record | Yes | the read as a linear record, ready for any downstream step | |
| sequence | Yes | its bases, trimmed when trimming is on | |
| ambiguous | Yes | how many calls are not A, C, G or T: the mixed positions | |
| qualities | Yes | per-base phred, parallel to the bases | |
| meanQuality | Yes | the ARITHMETIC mean of the per-base phred, which is what FastQC plots. Always the higher of the two, because a phred score is a log of an error rate: keep it for comparison with other tools, gate on readQuality | |
| qualityPlot | Yes | per-base phred as a plot along the read, with Q20 (1 error in 100) marked | |
| readQuality | Yes | the read's own quality: the phred of its MEAN ERROR PROBABILITY, which is what seqkit's --avg-qual reports. Wire THIS into a gate to stop a bad read early | |
| expectedAccuracy | Yes | the same number as a percentage of bases expected to be correct |