ex1/
ex1.rs

1use bam_builder::{bam_order::BamSortOrder, BamBuilder};
2
3fn main() {
4    // Create a builder with all defaults except the read_len is 100
5    let mut builder = BamBuilder::new(
6        100,
7        30,
8        "Pair".to_owned(),
9        None,
10        BamSortOrder::Unsorted,
11        None,
12        None,
13    );
14
15    // Create a builder for read pair spec
16    let records = builder
17        .pair_builder()
18        .contig(0)
19        .start1(0)
20        .start2(200)
21        .unmapped1(false)
22        .unmapped2(false)
23        .bases1("A".repeat(100))
24        .bases2("C".repeat(100))
25        .build()
26        .unwrap();
27    println!("{:?}", records);
28
29    // Add the pair to bam builder
30    builder.add_pair(records);
31    println!("{:?}", builder);
32    // Write records to a file
33    builder
34        .to_path(std::path::Path::new(&String::from("./test.bam")))
35        .unwrap();
36}