noodles-util 0.71.0

noodles support utilities
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
86
87
88
89

use std::io::{self, BufReader, Read};

use noodles_bcf as bcf;
use noodles_bgzf as bgzf;
use noodles_vcf::{self as vcf, variant::io::Read as _};

use crate::variant::Record;

pub(super) enum Inner<R> {
    Bcf(bcf::io::Reader<bgzf::io::Reader<BufReader<R>>>),
    BcfRaw(bcf::io::Reader<BufReader<R>>),
    Vcf(vcf::io::Reader<BufReader<R>>),
    VcfGz(vcf::io::Reader<bgzf::io::Reader<BufReader<R>>>),
}

impl<R> Inner<R>
where
    R: Read,
{
    pub(super) fn read_header(&mut self) -> io::Result<vcf::Header> {
        match self {
            Self::Bcf(reader) => reader.read_header(),
            Self::BcfRaw(reader) => reader.read_header(),
            Self::Vcf(reader) => reader.read_header(),
            Self::VcfGz(reader) => reader.read_header(),
        }
    }

    pub(super) fn read_record(&mut self, record: &mut Record) -> io::Result<usize> {
        match self {
            Inner::Bcf(reader) => {
                if !matches!(record, Record::Bcf(_)) {
                    *record = Record::Bcf(bcf::Record::default());
                }

                if let Record::Bcf(r) = record {
                    reader.read_record(r)
                } else {
                    unreachable!();
                }
            }
            Inner::BcfRaw(reader) => {
                if !matches!(record, Record::Bcf(_)) {
                    *record = Record::Bcf(bcf::Record::default());
                }

                if let Record::Bcf(r) = record {
                    reader.read_record(r)
                } else {
                    unreachable!();
                }
            }
            Inner::Vcf(reader) => {
                if !matches!(record, Record::Vcf(_)) {
                    *record = Record::Vcf(vcf::Record::default());
                }

                if let Record::Vcf(r) = record {
                    reader.read_record(r)
                } else {
                    unreachable!();
                }
            }
            Inner::VcfGz(reader) => {
                if !matches!(record, Record::Vcf(_)) {
                    *record = Record::Vcf(vcf::Record::default());
                }

                if let Record::Vcf(r) = record {
                    reader.read_record(r)
                } else {
                    unreachable!();
                }
            }
        }
    }

    pub(super) fn records<'a>(
        &'a mut self,
        header: &'a vcf::Header,
    ) -> impl Iterator<Item = io::Result<Box<dyn vcf::variant::Record>>> + 'a {
        match self {
            Inner::Bcf(reader) => reader.variant_records(header),
            Inner::BcfRaw(reader) => reader.variant_records(header),
            Inner::Vcf(reader) => reader.variant_records(header),
            Inner::VcfGz(reader) => reader.variant_records(header),
        }
    }
}