audit_cds_annotations
Check every annotated CDS of one or more GenBank maps against its own bases and genetic code: a frame that does not close (frameshift), a stop inside the reading frame, a first codon the code cannot start with, a last codon that is not a stop, and a /translation the bases do not give. A partial end (< or >) is not judged, and a /pseudo CDS carries no flag. Where an annotation ends without a stop it reads on in frame (round the origin of a circular map) to the first stop and says how many codons that adds, whether it runs into another annotated CDS (a fusion annotated in pieces) and the read-through a random backbone of the same composition would give. Extraction and translation match Biopython on all 1,871 CDS of 695 NCBI vector records.
Input Schema
| Name | Required | Description | Default |
|---|---|---|---|
| _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. | |
| genbank | Yes | the GenBank files' contents, one string per file (a file may hold several records). | |
| _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. |
Output Schema
| Name | Required | Description | Default |
|---|---|---|---|
| flags | Yes | per flag, how many CDS carry it and in how many records | |
| result | Yes | the full result | |
| flagged | Yes | every CDS with a flag: record, feature, label, location and the flags | |
| readThrough | Yes | every annotation that ends without a stop, with where reading on stops |