|
|
@@ -503,6 +503,7 @@ impl Genome {
|
|
|
aligner: impl Fn(String) -> Result<Vec<Mapping>>,
|
|
|
) -> Result<()> {
|
|
|
let mappings = aligner(String::from_utf8(sequence.to_vec())?)?;
|
|
|
+ println!("{mappings:?}");
|
|
|
self.add_contig(name, mappings, None, String::from_utf8(sequence.to_vec())?)?;
|
|
|
Ok(())
|
|
|
}
|
|
|
@@ -995,12 +996,12 @@ mod tests {
|
|
|
use dashmap::DashSet;
|
|
|
use indicatif::MultiProgress;
|
|
|
use indicatif_log_bridge::LogWrapper;
|
|
|
- use libVariant::{in_out::dict_reader::read_dict, utils::new_pg_speed, variants::Variant};
|
|
|
+ use libVariant::{in_out::dict_reader::read_dict, variants::Variant};
|
|
|
|
|
|
use super::*;
|
|
|
use crate::{
|
|
|
genomic_graph::GenomicGraph,
|
|
|
- phase::{merge_phases, variants_phasing, write_phases_bed, HeteroVar, Phase},
|
|
|
+ phase::{variants_phasing, write_phases_bed},
|
|
|
};
|
|
|
|
|
|
fn init() {
|
|
|
@@ -1039,11 +1040,11 @@ mod tests {
|
|
|
fn test_graph() -> Result<()> {
|
|
|
init();
|
|
|
|
|
|
- let case = "MERY";
|
|
|
- let chrom = vec!["chr7"];
|
|
|
+ let case = "SPINATO";
|
|
|
+ let chrom = vec!["chr5", "chr11", "chr14"];
|
|
|
info!("This record will be captured by `cargo test`");
|
|
|
|
|
|
- let dir = format!("/data/longreads_basic_pipe/{case}/diag/neo_contigs");
|
|
|
+ let dir = format!("/data/longreads_basic_pipe/{case}/diag/asm_contigs");
|
|
|
|
|
|
// Load from fasta in dir.
|
|
|
let genome = Genome::from_contigs_sequences(&dir)?;
|