jean_core 0.1.0

Core package of jean
Documentation 
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

use num_enum::{FromPrimitive, IntoPrimitive};
use strum_macros::EnumCount;

use crate::{
  prelude::{Complement, Gap, Transcribe},
  sequence::rna,
};

#[repr(u8)]
#[derive(
  Debug, Clone, Copy, PartialEq, Eq, PartialOrd, Ord, IntoPrimitive, FromPrimitive, EnumCount,
)]
pub enum Base {
  A,
  C,
  G,
  T,
  #[default]
  GAP,
}

impl Gap for Base {
  fn gap() -> Self {
    Self::GAP
  }
}

impl TryFrom<char> for Base {
  type Error = String;

  fn try_from(c: char) -> Result<Self, Self::Error> {
    match c.to_ascii_uppercase() {
      'A' => Ok(Self::A),
      'C' => Ok(Self::C),
      'G' => Ok(Self::G),
      'T' => Ok(Self::T),
      '-' => Ok(Self::GAP),
      _ => Err(format!("Invalid DNA base '{c}'.")),
    }
  }
}

impl From<Base> for char {
  fn from(base: Base) -> Self {
    match base {
      Base::A => 'A',
      Base::C => 'C',
      Base::G => 'G',
      Base::T => 'T',
      Base::GAP => '-',
    }
  }
}

impl Complement for Base {
  fn complement(&self) -> Self {
    match self {
      Self::A => Self::T,
      Self::C => Self::G,
      Self::G => Self::C,
      Self::T => Self::A,
      Self::GAP => Self::GAP,
    }
  }
}

impl Transcribe<rna::Base> for Base {
  fn transcribe(&self) -> rna::Base {
    match self {
      Self::A => rna::Base::A,
      Self::C => rna::Base::C,
      Self::G => rna::Base::G,
      Self::T => rna::Base::U,
      Self::GAP => rna::Base::GAP,
    }
  }
}