INNER CODE UNIT · Rust

errors

google/best · src/stats.rs:531

        let errors = res.mismatches + res.non_hp_ins + res.non_hp_del + res.hp_ins + res.hp_del;
        res.read_len = res.matches + res.mismatches + res.non_hp_del + res.hp_del;
        res.ref_cov = (res.read_len as f64) / (curr_ref_seq.len() as f64);
        res.gc_content /= res.read_len as f64;
        res.concordance = (res.matches as f64) / ((res.matches + errors) as f64);
        res.concordance_gc = (res.matches as f64)
            / ((res.matches + res.mismatches + res.gc_ins + res.gc_del) as f64);
        res.concordance_qv = concordance_qv(res.concordance, errors > 0);

        for (i, &has_error) in intervals.iter().zip(&interval_has_error) {
            if !has_error {
                res.feature_stats
                    .get_mut(i.val.as_str())
                    .unwrap()
                    .identical_overlaps += 1;
            }
        }

View source record →

📰 Research Paper
Loading…
⏳ Fetching content…