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
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127

//! Variant record samples key.

/// Read depth for each allele (`AD`).
pub const READ_DEPTHS: &str = "AD";

/// Read depth for each allele on the forward strand (`ADF`).
pub const FORWARD_STRAND_READ_DEPTHS: &str = "ADF";

/// Read depth for each allele on the reverse strand (`ADR`).
pub const REVERSE_STRAND_READ_DEPTHS: &str = "ADR";

/// Read depth (`DP`).
pub const READ_DEPTH: &str = "DP";

/// Expected alternate allele counts (`EC`).
pub const EXPECTED_ALTERNATE_ALLELE_COUNTS: &str = "EC";

/// Filter indicating if this genotype was "called" (`FT`).
pub const FILTER: &str = "FT";

/// Genotype likelihoods (`GL`).
pub const GENOTYPE_LIKELIHOODS: &str = "GL";

/// Genotype posterior probabilities (`GP`).
pub const GENOTYPE_POSTERIOR_PROBABILITIES: &str = "GP";

/// Conditional genotype quality (`GQ`).
pub const CONDITIONAL_GENOTYPE_QUALITY: &str = "GQ";

/// Genotype (`GT`).
pub const GENOTYPE: &str = "GT";

/// Haplotype quality (`HQ`).
pub const HAPLOTYPE_QUALITY: &str = "HQ";

/// RMS mapping quality (`MQ`).
pub const MAPPING_QUALITY: &str = "MQ";

/// Phred-scaled genotype likelihoods rounded to the closest integer (`PL`).
pub const ROUNDED_GENOTYPE_LIKELIHOODS: &str = "PL";

/// Phred-scaled genotype posterior probabilities rounded to the closest integer (`PP`).
pub const ROUNDED_GENOTYPE_POSTERIOR_PROBABILITIES: &str = "PP";

/// Phasing quality (`PQ`).
pub const PHASING_QUALITY: &str = "PQ";

/// Phase set (`PS`).
pub const PHASE_SET: &str = "PS";

/// Phase set list (`PSL`).
///
/// Added in VCF 4.4.
pub const PHASE_SET_LIST: &str = "PSL";

/// Phase set list ordinal (`PSO`).
///
/// Added in VCF 4.4.
pub const PHASE_SET_LIST_ORDINALS: &str = "PSO";

/// Phase set list quality (`PSQ`).
///
/// Added in VCF 4.4.
pub const PHASE_SET_LIST_QUALITIES: &str = "PSQ";

/// Copy number genotype for imprecise events (`CN`).
pub const GENOTYPE_COPY_NUMBER: &str = "CN";

/// Confidence interval around copy number (`CICN`).
///
/// Added in VCF 4.4.
pub const COPY_NUMBER_CONFIDENCE_INTERVAL: &str = "CICN";

/// Copy number genotype quality for imprecise events (`CNQ`).
pub const GENOTYPE_COPY_NUMBER_QUALITY: &str = "CNQ";

/// Copy number genotype likelihood for imprecise events (`CNL`).
pub const GENOTYPE_COPY_NUMBER_LIKELIHOODS: &str = "CNL";

/// Copy number posterior probabilities (`CNP`).
pub const GENOTYPE_COPY_NUMBER_POSTERIOR_PROBABILITIES: &str = "CNP";

/// Phred style probability score that the variant is novel (`NQ`).
pub const NOVEL_VARIANT_QUALITY_SCORE: &str = "NQ";

/// Unique haplotype identifier (`HAP`).
pub const HAPLOTYPE_ID: &str = "HAP";

/// Unique identifier of ancestral haplotype (`AHAP`).
pub const ANCESTRAL_HAPLOTYPE_ID: &str = "AHAP";

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

    #[test]
    fn test_fmt() {
        assert_eq!(READ_DEPTHS, "AD");
        assert_eq!(FORWARD_STRAND_READ_DEPTHS, "ADF");
        assert_eq!(REVERSE_STRAND_READ_DEPTHS, "ADR");
        assert_eq!(READ_DEPTH, "DP");
        assert_eq!(EXPECTED_ALTERNATE_ALLELE_COUNTS, "EC");
        assert_eq!(FILTER, "FT");
        assert_eq!(GENOTYPE_LIKELIHOODS, "GL");
        assert_eq!(GENOTYPE_POSTERIOR_PROBABILITIES, "GP");
        assert_eq!(CONDITIONAL_GENOTYPE_QUALITY, "GQ");
        assert_eq!(GENOTYPE, "GT");
        assert_eq!(HAPLOTYPE_QUALITY, "HQ");
        assert_eq!(MAPPING_QUALITY, "MQ");
        assert_eq!(ROUNDED_GENOTYPE_LIKELIHOODS, "PL");
        assert_eq!(ROUNDED_GENOTYPE_POSTERIOR_PROBABILITIES, "PP");
        assert_eq!(PHASING_QUALITY, "PQ");
        assert_eq!(PHASE_SET, "PS");
        assert_eq!(PHASE_SET_LIST, "PSL");
        assert_eq!(PHASE_SET_LIST_ORDINALS, "PSO");
        assert_eq!(PHASE_SET_LIST_QUALITIES, "PSQ");

        assert_eq!(GENOTYPE_COPY_NUMBER, "CN");
        assert_eq!(COPY_NUMBER_CONFIDENCE_INTERVAL, "CICN");
        assert_eq!(GENOTYPE_COPY_NUMBER_QUALITY, "CNQ");
        assert_eq!(GENOTYPE_COPY_NUMBER_LIKELIHOODS, "CNL");
        assert_eq!(GENOTYPE_COPY_NUMBER_POSTERIOR_PROBABILITIES, "CNP");
        assert_eq!(NOVEL_VARIANT_QUALITY_SCORE, "NQ");
        assert_eq!(HAPLOTYPE_ID, "HAP");
        assert_eq!(ANCESTRAL_HAPLOTYPE_ID, "AHAP");
    }
}