1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
//! Lightweight HGVS variant parser.
//!
//! `tinyhgvs` parses a HGVS variant into explicit Rust structs and enums
//! that describe:
//!
//! - the reference sequence context such as `NM_004006.2` or `NP_003997.1`
//! - the coordinate type such as coding DNA (`c.`), genomic DNA (`g.`), RNA
//! (`r.`), or protein (`p.`)
//! - the biological description itself, represented as either a nucleotide
//! variant, a nucleotide allele, or a protein consequence
//!
//! The crate is intentionally small. It aims to represent common, high-value
//! HGVS syntax clearly, while returning structured errors for syntax families
//! tracked in the unsupported inventory.
//!
//! The main entry points are:
//!
//! - [`parse_hgvs`] to parse a string into [`HgvsVariant`]
//! - [`ParseHgvsError`] to inspect invalid or unsupported input
//!
//! # Reading the Parsed Model
//!
//! The [`HgvsVariant`] separates a HGVS syntax into three top-level parts:
//!
//! - `reference`: the reference source for a variant.
//! - `coordinate_system`: the one-letter HGVS coordinate type.
//! - `description`: the nucleotide or protein variant description, including
//! location and base edits or effects.
//!
//! # Examples
//!
//! A substitution crossing exon/intron border (intronic):
//!
//! ```rust
//! use tinyhgvs::{NucleotideAnchor, NucleotideEdit, VariantDescription, parse_hgvs};
//!
//! let variant = parse_hgvs("NM_004006.2:c.357+1G>A").unwrap();
//! let description = variant.description;
//!
//! match description {
//! VariantDescription::Nucleotide(nucleotide) => {
//! assert_eq!(nucleotide.location.start().unwrap().anchor().unwrap(), NucleotideAnchor::Absolute);
//! assert_eq!(nucleotide.location.start().unwrap().coordinate().unwrap(), 357);
//! assert_eq!(nucleotide.location.start().unwrap().offset().unwrap(), 1);
//! assert!(matches!(
//! nucleotide.edit,
//! NucleotideEdit::Substitution { ref reference, ref alternate }
//! if reference == "G" && alternate == "A"
//! ));
//! }
//! _ => unreachable!("expected nucleotide variant"),
//! }
//! ```
//!
//! A nucleotide allele keeps one first allele, an optional second
//! established allele, and any later unphased additions:
//!
//! ```rust
//! use tinyhgvs::{AllelePhase, VariantDescription, parse_hgvs};
//!
//! let variant = parse_hgvs("NM_004006.2:c.[2376G>C];[2376=]").unwrap();
//!
//! match variant.description {
//! VariantDescription::NucleotideAllele(allele) => {
//! assert_eq!(allele.allele_one.variants.len(), 1);
//! assert!(allele.allele_two.is_some());
//! assert_eq!(allele.phase, Some(AllelePhase::Trans));
//! }
//! _ => unreachable!("expected nucleotide allele"),
//! }
//! ```
//!
//! A nonsense mutation leading to an early termination at protein-level:
//!
//! ```rust
//! use tinyhgvs::{CoordinateSystem, ProteinEffect, VariantDescription, parse_hgvs};
//!
//! let variant = parse_hgvs("NP_003997.1:p.Trp24Ter").unwrap();
//! assert_eq!(variant.coordinate_system, CoordinateSystem::Protein);
//!
//! match variant.description {
//! VariantDescription::Protein(protein) => {
//! assert!(!protein.is_predicted);
//! assert!(matches!(protein.effect, ProteinEffect::Edit { .. }));
//! }
//! _ => unreachable!("expected protein variant"),
//! }
//! ```
//!
//! Unsupported syntax is reported with a stable diagnostic code:
//!
//! ```rust
//! use tinyhgvs::parse_hgvs;
//!
//! let error = parse_hgvs("NM_004006.2:c.[2376G>C];[?]").unwrap_err();
//! assert_eq!(error.code(), "unsupported.allele_unknown_variant");
//! ```
pub use ;
pub use ;
pub use parse_hgvs;