vcf-reformatter 0.7.6

Fast VCF file parser and reformatter with VEP and SnpEff annotation support which can output to MAF
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

use regex::Regex;

/// Extract a `TAG=...` annotation value from the INFO field (parsed_lines[7]),
/// stripping it (and any resulting stray `;`) out of the remaining INFO string.
pub(crate) fn extract_tag_regex(parsed_lines: &mut [String], tag: &str) -> Option<String> {
    if parsed_lines.len() < 8 {
        return None;
    }

    let info_field = &parsed_lines[7];
    // Anchor to an INFO key boundary, or `MYCSQ=` would match as `CSQ=`.
    let re = Regex::new(&format!(r"(?:^|;){tag}=([^;]+)")).ok()?;

    let captures = re.captures(info_field)?;
    let value = captures.get(1)?.as_str().to_string();

    let new_info = re.replace(info_field, "").to_string();
    let cleaned_info = new_info
        .replace(";;", ";")
        .trim_start_matches(';')
        .trim_end_matches(';')
        .to_string();

    parsed_lines[7] = cleaned_info;

    Some(value)
}

pub fn extract_csq_regex(parsed_lines: &mut [String]) -> Option<String> {
    extract_tag_regex(parsed_lines, "CSQ")
}