seal 0.1.6

Implementation of Needleman-Wunsch & Smith-Waterman sequence alignment.
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

use crate::pair::step_mask::StepMask;

#[derive(Copy, Clone, Debug, PartialEq)]
pub struct Cursor {
    pub x: usize,
    pub y: usize,
}

impl Cursor {
    pub fn apply_forwards_step(&mut self, step_mask: StepMask) {
        self.apply_step(step_mask, true)
    }

    pub fn apply_backwards_step(&mut self, step_mask: StepMask) {
        self.apply_step(step_mask, false)
    }

    fn apply_step(&mut self, step_mask: StepMask, forward: bool) {
        let delta = match step_mask {
            StepMask::ALIGN => (1, 1),
            StepMask::INSERT => (0, 1),
            StepMask::DELETE => (1, 0),
            StepMask::STOP => (0, 0),
            _ => {
                panic!("Invalid step_mask.");
            }
        };
        if forward {
            self.x += delta.0;
            self.y += delta.1;
        } else {
            self.x -= delta.0;
            self.y -= delta.1;
        }
    }
}

#[cfg(test)]
mod tests {
    use super::*;

    #[test]
    fn apply_works() {
        {
            let mut cursor = Cursor { x: 10, y: 10 };
            cursor.apply_backwards_step(StepMask::STOP);
            assert_eq!(cursor, Cursor { x: 10, y: 10 });
        }
        {
            let mut cursor = Cursor { x: 10, y: 10 };
            cursor.apply_backwards_step(StepMask::ALIGN);
            assert_eq!(cursor, Cursor { x: 9, y: 9 });
        }
        {
            let mut cursor = Cursor { x: 10, y: 10 };
            cursor.apply_backwards_step(StepMask::INSERT);
            assert_eq!(cursor, Cursor { x: 10, y: 9 });
        }
        {
            let mut cursor = Cursor { x: 10, y: 10 };
            cursor.apply_backwards_step(StepMask::DELETE);
            assert_eq!(cursor, Cursor { x: 9, y: 10 });
        }
    }
}