gen 0.1.23

A sequence graph and version control system.
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
78
79
80
81
82
83
84
85

use std::{
    fs::File,
    io::{self, BufWriter, Write},
};

use convert_case::{Case, Casing};
use gen_core::{HashId, Strand};

#[derive(Clone, Debug, Eq, Hash, PartialEq, Ord, PartialOrd)]
pub struct Segment {
    pub sequence: String,
    pub node_id: HashId,
    pub sequence_start: i64,
    pub sequence_end: i64,
    pub strand: Strand,
}

#[derive(Clone, Debug, Eq, Hash, PartialEq, Ord, PartialOrd)]
pub struct Link {
    pub source_segment_id: String,
    pub source_strand: Strand,
    pub target_segment_id: String,
    pub target_strand: Strand,
}

#[derive(Clone, Debug, Eq, Hash, PartialEq)]
pub struct Path {
    pub name: String,
    pub segment_ids: Vec<String>,
    pub node_strands: Vec<Strand>,
}

impl Segment {
    pub fn segment_id(&self) -> String {
        format!(
            "{}.{}.{}",
            self.node_id, self.sequence_start, self.sequence_end
        )
    }
}

fn segment_line(segment: &Segment) -> String {
    // NOTE: We encode the node ID and start coordinate in the segment ID
    format!("S\t{}\t{}\n", segment.segment_id(), segment.sequence)
}

fn link_line(link: &Link) -> String {
    format!(
        "L\t{}\t{}\t{}\t{}\t0M\n",
        link.source_segment_id, link.source_strand, link.target_segment_id, link.target_strand
    )
}

pub fn path_line(path: &Path) -> String {
    let segments = path
        .segment_ids
        .iter()
        .zip(path.node_strands.iter())
        .map(|(segment_id, node_strand)| format!("{segment_id}{node_strand}"))
        .collect::<Vec<String>>()
        .join(",");
    format!("P\t{}\t{}\t*\n", path.name.to_case(Case::Train), segments)
}

pub fn write_segments(writer: &mut BufWriter<File>, segments: &[&Segment]) -> io::Result<()> {
    for segment in segments {
        writer.write_all(&segment_line(segment).into_bytes())?;
    }
    Ok(())
}

pub fn write_links(writer: &mut BufWriter<File>, links: &[&Link]) -> io::Result<()> {
    for link in links {
        writer.write_all(&link_line(link).into_bytes())?;
    }
    Ok(())
}

pub fn bool_to_strand(direction: bool) -> Strand {
    if direction {
        Strand::Forward
    } else {
        Strand::Reverse
    }
}