aardvark-bio 0.10.2

Aardvark - A tool for sniffing out the differences in vari-Ants
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


use std::path::{Path, PathBuf};

/// Wrapper that indexes a VCF file using noodles indexer, generating a .tbi file.
/// # Arguments
/// * `vcf_fn` - the filename to index
pub fn index_vcf(vcf_fn: &Path) -> anyhow::Result<()> {
    // first, build the index
    let index = noodles::vcf::fs::index(vcf_fn)?;

    // add the .tbi extension to the filename
    let mut tbi_fn = vcf_fn.to_owned().into_os_string();
    tbi_fn.push(".tbi");
    let tbi_fn = PathBuf::from(tbi_fn);
    
    // write the index out to file
    noodles::tabix::fs::write(&tbi_fn, &index)?;

    Ok(())
}

/// Wrapper that indexes a BED file using noodles indexer, generating a .tbi file.
/// # Arguments
/// * `bed_fn` - the filename to index
pub fn index_bed(bed_fn: &Path) -> anyhow::Result<()> {
    // first, build the index
    let index = noodles::bed::fs::index(bed_fn)?;

    // add the .tbi extension to the filename
    let mut tbi_fn = bed_fn.to_owned().into_os_string();
    tbi_fn.push(".tbi");
    let tbi_fn = PathBuf::from(tbi_fn);

    // write the index out to file
    noodles::tabix::fs::write(&tbi_fn, &index)?;

    Ok(())
}