openvax / varcode / 25761723707
88%

DEFAULT BRANCH: main
Ran
Jobs 0
Files 0
Run time –
Badge
Embed ▾
README BADGES
x

If you need to use a raster PNG badge, change the '.svg' to '.png' in the link

Markdown

Textile

RDoc

HTML

Rst
25761723707

push

github

web-flow 
varcode.transforms.left_align_indels (#369, planning + scaffold) (#374)

* Scaffold varcode.transforms.left_align_indels signature (#369)

Planning PR — implementation tracked in the description. This commit
adds only:

* left_align_indels function signature with the full contract
  docstring (cardinality, provenance, metadata behavior, idempotence,
  reference-tier coverage including the partial-shift case for
  no-FASTA exon-boundary crossings).
* Module-docstring promotion from "planned" to "shipped" alongside
  pair_breakends.
* __all__ entry.

Body raises NotImplementedError pointing at #369. Unlike #364's
pair_breakends scaffold, no new core attributes — left_align_indels
reuses Variant.source_variants and the reference helpers from #372.

* Implement left_align_indels: algorithm, tests, docs (#369)

Implementation, test matrix, and docs for the left-alignment indel
transform. Builds on the genome-FASTA infrastructure from #372 — no
new dependencies, no new public API surface beyond the function
itself.

## Algorithm

Standard variant-normalization left-shift (Tan/Abecasis/Kang 2015),
adapted to varcode's normalized indel representation. The shift
check position differs between deletions (reference[start-1]) and
insertions (reference[start]) because varcode's start convention
differs for the two cases — deletion start points at the first
deleted base, insertion start at the anchor base before the inserted
sequence.

Shift terminates on natural disagreement (the reference base doesn't
match the rotating payload) OR coverage exhaustion (reference_range
returns empty — chromosome start, off-end-of-FASTA, exon boundary
without a FASTA). The latter case sets info["left_align_partial"]
on the shifted variant so downstream tools know a fuller shift might
have been possible with a chromosome FASTA.

## Tests (18 cases, no real FASTA dependency)

Uses the _DictBackedFasta mock established in test_genome_sequence.py:

* Pass-throughs: SNV, MNV, complex va... (continued)
Source Files on build 25761723707
Detailed source file information is not available for this build.