@@ -252,3 +252,15 @@ pub fn vcf_header(dict: &str) -> anyhow::Result<Vec<String>> {
Ok(header)
}
+#[cfg(test)]
+mod tests {
+ use crate::config::Config;
+
+ use super::*;
+ #[test]
+ fn vcf_dbsnp() {
+ let config = Config::default();
+ fetch_vcf(conf.db_snp, &GenomeRange::new("chr10", , end))
+ }
+}
@@ -581,6 +581,7 @@ pub fn somatic_depth_quality_ranges(
.filter(|(ctg, _)| ctg != "chrM")
.map(|(sn, _ln)| sn)
.collect();
+ println!("{contigs:#?}");
let cfg = config;