{ "report_id": "SYNTENY-NATIVE-VS-MINIMAP2-2026", "title": "Whole-Genome Synteny, Collinear Chaining & Structural Variant Breakpoint Detection: Native Rust 4-Stage Minimizer & Concave-Gap Engine vs. Minimap2 v2.31-r1302", "tool_pair_slot": "Engine 2 (bio-suite genome) — Synteny & Whole-Genome Alignment", "native_engine_version": "bio-suite SyntenyMode::Native v2.0 (Double-Stranded (w,k)-Minimizer + Concave-Gap DP + Two-Piece Affine-Gap + SyRI Split-Chain SV Caller)", "external_benchmark_version": "Minimap2 v2.31-r1302 (Heng Li, MIT License, -c --cs -x asm5/asm10/asm20)", "spdx_license": "GPL-3.0-or-later (bio-suite) / MIT (Minimap2 v2.31-r1302)", "total_cohorts": 4, "total_genome_pairs_evaluated": 14, "total_bases_aligned_bp": 439822, "overall_sv_recall_concordance_pct": 70.7143, "overall_breakpoint_mae_bp": 998.625, "overall_de_divergence_mae": 0.0135, "overall_blast_identity_mae_pct": 1.0511, "native_total_wall_clock_ms": 89.8835, "minimap2_total_wall_clock_ms": 166.1163, "speedup_factor": 1.8481, "cohorts": [ { "cohort_id": "COHORT-01-ENTERIC-60KB-ISLAND-SYNTENY", "title": "Enteric Chromosomal Origin Synteny & +31.4 kb Lineage-Specific Genomic Island (E. coli K-12 vs. S. enterica LT2)", "domain": "Prokaryota (Enterobacterales)", "biological_context": "Cross-genus chromosomal collinearity (~16.5% nucleotide divergence) across the thrABC-dnaKJ-nhaA-ribF-ileS region interrupted by a +31,428 bp S. enterica genomic island insertion.", "target_accession": "NC_000913.3:1-60000", "query_accession": "NC_003197.2:1-60000", "target_len_bp": 60000, "query_len_bp": 60000, "native_result": { "engine_name": "bio-suite Native 4-Stage Synteny & SV Engine", "engine_mode": "Native", "preset_label": "multi-scale (w=5,k=13)+(w=10,k=19)", "target_name": "NC_000913.3_E_coli_K12", "target_len_bp": 60000, "query_name": "NC_003197.2_S_enterica_LT2", "query_len_bp": 60000, "minimizer_seeds_target": 31370, "minimizer_seeds_query": 31363, "raw_anchors_count": 1275, "primary_chains_count": 3, "secondary_chains_count": 0, "target_aligned_bases": 23290, "target_coverage_pct": 38.8167, "weighted_blast_identity_pct": 76.7252, "weighted_gap_compressed_de": 0.2087, "mean_primary_mapq": 60.0, "wall_clock_ms": 11.9768, "paf_records": [ { "query_name": "NC_003197.2_S_enterica_LT2", "query_len": 60000, "query_start_0based": 106, "query_end_0based": 14523, "strand": "+", "target_name": "NC_000913.3_E_coli_K12", "target_len": 60000, "target_start_0based": 106, "target_end_0based": 15087, "residue_matches": 11252, "alignment_block_len": 15084, "mapq": 60, "alignment_type": "P", "chain_score_f1": 3175, "anchor_count": 636, "edit_distance_nm": 3832, "blast_identity_pct": 74.5956, "gap_compressed_divergence_de": 0.2201, "cigar_cg": "97M3I14M3D70M1I16M1D4703M572D1284M9I1062M1I3M1D6M1D2M1D2M2I15M1D2M1D1M2I8M1I2M2I6M1I9M1D1M1I21M1D2M1I3M1D5M1D1M2I1M1D1M2D1M1D1M1D1M1D1M1D1M1D2M1D4M1D1M1D2M2I4M1I1M2I10M1D1M1I2M2D6M2I9M2D2M3D6M1D1M1D2M1I8M1D1M1I1M1D6M2D1M1D2M1I1M1I3M1D1M1D5M2I2M1I17M2I3M1D1M1D38M1D989M2D2M1I12M1I1M1D88M3D506M25I695M1D7M2D2M1I29M1D1M1D2M1I1M2I292M2I752M5I294M1I2M1D4M5D3M1D1M3D12M1D2M3D1M2D6M1I2M2D3M1I4M1I1M1I7M1D5M2D6M2D6M1D2M1I1M3D1M1D1M1D67M1D19M2D3M2I14M1D1M1I2000M1D1M1D20M2D24M1I9M1I1M2D6M1I10M1D2M1I225M2I1M1I3M6I680M", "diff_string_cs": ":14*tc*cg*ag*cg*tc:1*ag:19*ca:3*gc*ca*ga*cg:7*ga:8*tg:23*ac:5*tc:2+cac:14-cca:70+a:16-t:9*ag:2*ga:1*ac:103*ga:8*ga:2*cg:5*ca:2*ct:2*tc:8*cg:5*ct:8*ct:8*ga:14*ct:2*tc*ag:16*ag:2*cg:14*ca:12*gt:1*ct:2*ac:2*tg*tc:1*gc*ag:1*ga:5*ct:3*ga*cg:3*cg" }, { "query_name": "NC_003197.2_S_enterica_LT2", "query_len": 60000, "query_start_0based": 46085, "query_end_0based": 47946, "strand": "+", "target_name": "NC_000913.3_E_coli_K12", "target_len": 60000, "target_start_0based": 17385, "target_end_0based": 19243, "residue_matches": 1484, "alignment_block_len": 1872, "mapq": 60, "alignment_type": "P", "chain_score_f1": 170, "anchor_count": 23, "edit_distance_nm": 388, "blast_identity_pct": 79.2735, "gap_compressed_divergence_de": 0.2051, "cigar_cg": "22M1D2M1I26M1I1M1D42M1I1114M1I1M1D1M1I3M1D39M2I2M1I7M1D1M1I2M2D7M2I1M1D15M1I4M1D8M2I7M2D542M", "diff_string_cs": ":13*ag*tc:5*ca:1-a:2+t*ca:8*ga:13*tg:2+a:1-t:3*ga:5*ag:6*ag:1*at*ac:3*ta*cg*ta:10*at*at:4+g:22*tc:5*ct:8*tc:8*gt:2*ag:5*tc:5*tc:5*tc:6*ag*tc*cg:5*gt:3*ac:1*tg:8*ct:1*gt*cg:5*ag:17*ct:2*ct:23*ca:8*tg*tg:1*ag:2*ct:2*ag:20*ag:5*ac:2*tc:2*ct:11" }, { "query_name": "NC_003197.2_S_enterica_LT2", "query_len": 60000, "query_start_0based": 52217, "query_end_0based": 58724, "strand": "+", "target_name": "NC_000913.3_E_coli_K12", "target_len": 60000, "target_start_0based": 20754, "target_end_0based": 27205, "residue_matches": 5276, "alignment_block_len": 6520, "mapq": 60, "alignment_type": "P", "chain_score_f1": 2065, "anchor_count": 477, "edit_distance_nm": 1244, "blast_identity_pct": 80.9202, "gap_compressed_divergence_de": 0.1834, "cigar_cg": "22M1I8M1D6M1I8M2I1M1D5M2I1M2D640M3D866M1D2M1I22M1D1M1I8M2D1M3I2M1I1488M1D1M1I2M6I11M1I3M1D1057M12I959M36I392M1I932M", "diff_string_cs": ":22+t:1*ga:6-c:6+t:8+ga:1-c:1*ga:3+gc:1-tt:41*ct:84*ta:108*tc:134*at:1*ta*ga:1*gt*cg:1*ga:9*ag:33*gt:1*gc:12*tc:5*ga*ct:2*tg:10*gt*ct:1*ga:4*tc*cg*ag:4*ct*tg*tg*cg*cg:2*tg*ga*ag:5*gc*cg:6*ag:9*ca*at:7*ca:2*at:2*tc*ca*gt*gc*ta*ag*cg*tg:18*ga" } ], "sv_calls": [ { "sv_id": "SV-01-SYN", "sv_class": "SYN", "target_name": "NC_000913.3_E_coli_K12", "target_start_1based": 107, "target_end_1based": 15087, "query_name": "NC_003197.2_S_enterica_LT2", "query_start_1based": 107, "query_end_1based": 14523, "strand": "+", "event_size_bp": 14981, "blast_identity_pct": 74.5956, "gap_compressed_de": 0.2201, "mapq": 60, "biological_annotation": "Collinear syntenic block (14981 bp on '+' strand, 74.60% identity)" }, { "sv_id": "SV-02-INS_GI", "sv_class": "INS_GI", "target_name": "NC_000913.3_E_coli_K12", "target_start_1based": 15088, "target_end_1based": 17385, "query_name": "NC_003197.2_S_enterica_LT2", "query_start_1based": 14524, "query_end_1based": 46085, "strand": "+", "event_size_bp": 29264, "blast_identity_pct": 0.0, "gap_compressed_de": 1.0, "mapq": 60, "biological_annotation": "Lineage-specific genomic island / insertion (+31562 bp in query vs. 2298 bp in target)" }, { "sv_id": "SV-03-SYN", "sv_class": "SYN", "target_name": "NC_000913.3_E_coli_K12", "target_start_1based": 17386, "target_end_1based": 19243, "query_name": "NC_003197.2_S_enterica_LT2", "query_start_1based": 46086, "query_end_1based": 47946, "strand": "+", "event_size_bp": 1858, "blast_identity_pct": 79.2735, "gap_compressed_de": 0.2051, "mapq": 60, "biological_annotation": "Collinear syntenic block (1858 bp on '+' strand, 79.27% identity)" }, { "sv_id": "SV-04-INS_GI", "sv_class": "INS_GI", "target_name": "NC_000913.3_E_coli_K12", "target_start_1based": 19244, "target_end_1based": 20754, "query_name": "NC_003197.2_S_enterica_LT2", "query_start_1based": 47947, "query_end_1based": 52217, "strand": "+", "event_size_bp": 2760, "blast_identity_pct": 0.0, "gap_compressed_de": 1.0, "mapq": 60, "biological_annotation": "Lineage-specific genomic island / insertion (+4271 bp in query vs. 1511 bp in target)" }, { "sv_id": "SV-05-SYN", "sv_class": "SYN", "target_name": "NC_000913.3_E_coli_K12", "target_start_1based": 20755, "target_end_1based": 27205, "query_name": "NC_003197.2_S_enterica_LT2", "query_start_1based": 52218, "query_end_1based": 58724, "strand": "+", "event_size_bp": 6451, "blast_identity_pct": 80.9202, "gap_compressed_de": 0.1834, "mapq": 60, "biological_annotation": "Collinear syntenic block (6451 bp on '+' strand, 80.92% identity)" } ] }, "minimap2_asm5_result": { "engine_name": "Minimap2 v2.31-r1302 (-c --cs -x asm5)", "engine_mode": "Minimap2", "preset_label": "asm5", "target_name": "NC_000913.3_E_coli_K12", "target_len_bp": 60000, "query_name": "NC_003197.2_S_enterica_LT2", "query_len_bp": 60000, "minimizer_seeds_target": 0, "minimizer_seeds_query": 0, "raw_anchors_count": 0, "primary_chains_count": 1, "secondary_chains_count": 0, "target_aligned_bases": 370, "target_coverage_pct": 0.6167, "weighted_blast_identity_pct": 99.1892, "weighted_gap_compressed_de": 0.0081, "mean_primary_mapq": 60.0, "wall_clock_ms": 51.547, "paf_records": [ { "query_name": "NC_003197.2_S_enterica_LT2", "query_len": 60000, "query_start_0based": 52274, "query_end_0based": 52644, "strand": "+", "target_name": "NC_000913.3_E_coli_K12", "target_len": 60000, "target_start_0based": 20809, "target_end_0based": 21179, "residue_matches": 367, "alignment_block_len": 370, "mapq": 60, "alignment_type": "P", "chain_score_f1": 350, "anchor_count": 26, "edit_distance_nm": 3, "blast_identity_pct": 99.1892, "gap_compressed_divergence_de": 0.0081, "cigar_cg": "370M", "diff_string_cs": ":41*ct:84*ta:108*tc:134" } ], "sv_calls": [ { "sv_id": "SV-01-SYN", "sv_class": "SYN", "target_name": "NC_000913.3_E_coli_K12", "target_start_1based": 20810, "target_end_1based": 21179, "query_name": "NC_003197.2_S_enterica_LT2", "query_start_1based": 52275, "query_end_1based": 52644, "strand": "+", "event_size_bp": 370, "blast_identity_pct": 99.1892, "gap_compressed_de": 0.0081, "mapq": 60, "biological_annotation": "Collinear syntenic block (370 bp on '+' strand, 99.19% identity)" } ] }, "minimap2_asm20_result": { "engine_name": "Minimap2 v2.31-r1302 (-c --cs -x asm20)", "engine_mode": "Minimap2", "preset_label": "asm20", "target_name": "NC_000913.3_E_coli_K12", "target_len_bp": 60000, "query_name": "NC_003197.2_S_enterica_LT2", "query_len_bp": 60000, "minimizer_seeds_target": 0, "minimizer_seeds_query": 0, "raw_anchors_count": 0, "primary_chains_count": 6, "secondary_chains_count": 0, "target_aligned_bases": 17386, "target_coverage_pct": 28.9767, "weighted_blast_identity_pct": 82.1961, "weighted_gap_compressed_de": 0.149, "mean_primary_mapq": 60.0, "wall_clock_ms": 45.7602, "paf_records": [ { "query_name": "NC_003197.2_S_enterica_LT2", "query_len": 60000, "query_start_0based": 127, "query_end_0based": 5731, "strand": "+", "target_name": "NC_000913.3_E_coli_K12", "target_len": 60000, "target_start_0based": 127, "target_end_0based": 6303, "residue_matches": 4601, "alignment_block_len": 6177, "mapq": 60, "alignment_type": "P", "chain_score_f1": 361, "anchor_count": 39, "edit_distance_nm": 1576, "blast_identity_pct": 74.486, "gap_compressed_divergence_de": 0.1804, "cigar_cg": "163M1I11M1D4718M535D5M2D7M2D12M3D15M8D12M11D27M2D14M6D4M3D615M", "diff_string_cs": ":19*ca:3*gc*ca*ga*cg:7*ga:8*tg:23*ac:5*tc:10*tc*tc:1*ct*ct:74+a:11-t:14*ag:2*ga:1*ac:103*ga:8*ga:2*cg:5*ca:2*ct:2*tc:8*cg:5*ct:8*ct:8*ga:14*ct:2*tc*ag:16*ag:2*cg:14*ca:12*gt:1*ct:2*ac:2*tg*tc:1*gc*ag:1*ga:5*ct:3*ga*cg:3*cg:8*ga:8*gt:2*ga:1*" }, { "query_name": "NC_003197.2_S_enterica_LT2", "query_len": 60000, "query_start_0based": 7596, "query_end_0based": 9993, "strand": "+", "target_name": "NC_000913.3_E_coli_K12", "target_len": 60000, "target_start_0based": 8168, "target_end_0based": 10547, "residue_matches": 2017, "alignment_block_len": 2404, "mapq": 60, "alignment_type": "P", "chain_score_f1": 388, "anchor_count": 42, "edit_distance_nm": 387, "blast_identity_pct": 83.9018, "gap_compressed_divergence_de": 0.1529, "cigar_cg": "35M1D992M1D105M3D575M8I2M3I13M13I6M1I601M1D11M1D32M", "diff_string_cs": ":9*at:11*ct:13-c:26*tc*at:36*at:15*ct:2*ct:2*tc:2*ga:2*ct:14*ga:8*ag:11*ct:5*ct:29*ga:2*tc:8*ct:6*tc:7*tc:31*ag:3*tc:12*ag:1*ct:2*gt:11*ct:2*ag:5*ag:2*ag:8*tc:5*ag:8*ag:5*tg:5*cg:5*ct:2*at:2*tc:8*tc:8*tg:2*ct:8*ct:8*at:2*tc:2*gc:11*ct*ca:1*" }, { "query_name": "NC_003197.2_S_enterica_LT2", "query_len": 60000, "query_start_0based": 11422, "query_end_0based": 14574, "strand": "+", "target_name": "NC_000913.3_E_coli_K12", "target_len": 60000, "target_start_0based": 11991, "target_end_0based": 15138, "residue_matches": 2818, "alignment_block_len": 3157, "mapq": 60, "alignment_type": "P", "chain_score_f1": 749, "anchor_count": 81, "edit_distance_nm": 339, "blast_identity_pct": 89.262, "gap_compressed_divergence_de": 0.1045, "cigar_cg": "66M1D2041M2D20M2D24M1I260M9I731M", "diff_string_cs": ":11*gt:22*tg:12*tc:1*at:16-a:9*ca:7*tc:2*at*at*tc*ta:3*ct:10*ag*gc:11*ct:9*ct:23*at:1*at:27*at:50*ca:2*cg*ac*ca*tg*cg:1*ta:20*ag:11*ct:17*ct:23*ag:32*ga:14*tc:26*tc:23*at:2*ga:5*tc:5*ct:10*ta:6*tc:2*tc:1*cg*tc:2*tc:2*tc:15*gc:1*ct:2*at:2*tg" }, { "query_name": "NC_003197.2_S_enterica_LT2", "query_len": 60000, "query_start_0based": 52274, "query_end_0based": 56255, "strand": "+", "target_name": "NC_000913.3_E_coli_K12", "target_len": 60000, "target_start_0based": 20809, "target_end_0based": 24785, "residue_matches": 3431, "alignment_block_len": 3993, "mapq": 60, "alignment_type": "P", "chain_score_f1": 897, "anchor_count": 90, "edit_distance_nm": 562, "blast_identity_pct": 85.9254, "gap_compressed_divergence_de": 0.1366, "cigar_cg": "521M2I4M2D113M3D297M1I5M1D597M3I5M1D709M5I5M5D773M6I935M", "diff_string_cs": ":41*ct:84*ta:108*tc:134*at:1*ta*ga:1*gt*cg:1*ga:9*ag:33*gt:1*gc:12*tc:5*ga*ct:2*tg:10*gt*ct:1*ga:4*tc*cg*ag:4*ct*tg*tg*cg*cg:2*tg*ga*ag:5*gc*cg:6*ag:9*ca*at:7*ca:2*at:2+ca:2*ga:1-ta*cg*tg:18*ga:21*gt*gt*ag:19*ct:37*ct:10-aag:1*at:14*tc:20*g" }, { "query_name": "NC_003197.2_S_enterica_LT2", "query_len": 60000, "query_start_0based": 56557, "query_end_0based": 57161, "strand": "+", "target_name": "NC_000913.3_E_coli_K12", "target_len": 60000, "target_start_0based": 25075, "target_end_0based": 25679, "residue_matches": 530, "alignment_block_len": 604, "mapq": 60, "alignment_type": "P", "chain_score_f1": 70, "anchor_count": 7, "edit_distance_nm": 74, "blast_identity_pct": 87.7483, "gap_compressed_divergence_de": 0.1225, "cigar_cg": "604M", "diff_string_cs": ":5*ag:5*gt:5*ct:2*ag:17*ct*ca*ac*gt:23*ct:5*at:8*ca:15*ga:8*ga:1*ct:2*tc:3*ac:7*ga:14*ca:2*tc:1*ga*ac:13*ga:11*ga:8*ga:21*ct:28*ct:2*tc:17*ca*cg*cg:39*tc:8*gt:5*tc:5*tg:2*cg:20*tc:14*tg:8*ta:2*tg:2*tc*at:4*ga:2*ct:2*gt*gt*ct*ag:20*ga:2*cg:1" }, { "query_name": "NC_003197.2_S_enterica_LT2", "query_len": 60000, "query_start_0based": 57637, "query_end_0based": 58742, "strand": "+", "target_name": "NC_000913.3_E_coli_K12", "target_len": 60000, "target_start_0based": 26119, "target_end_0based": 27223, "residue_matches": 938, "alignment_block_len": 1105, "mapq": 60, "alignment_type": "P", "chain_score_f1": 100, "anchor_count": 10, "edit_distance_nm": 167, "blast_identity_pct": 84.8869, "gap_compressed_divergence_de": 0.1511, "cigar_cg": "154M1I950M", "diff_string_cs": ":26*tc:2*ga:14*ct:8*tc:11*cg:14*ct:20*gt:8*tc:32*ag:8*ga:1+a:26*ga:2*tc:8*tc:35*ct:5*gc:2*ct:11*ag:2*ga:2*at:2*tc:2*ct:2*tg:8*gc:2*ag:14*ga:2*cg:2*tc:3*tc:4*tc*ga:4*tc:5*tg:8*ta:5*ga:11*ag:5*ct:5*gc:2*ct:14*at:2*ct:2*tg:5*tc:8*ag:8*ag:2*ag:" } ], "sv_calls": [ { "sv_id": "SV-01-SYN", "sv_class": "SYN", "target_name": "NC_000913.3_E_coli_K12", "target_start_1based": 128, "target_end_1based": 6303, "query_name": "NC_003197.2_S_enterica_LT2", "query_start_1based": 128, "query_end_1based": 5731, "strand": "+", "event_size_bp": 6176, "blast_identity_pct": 74.486, "gap_compressed_de": 0.1804, "mapq": 60, "biological_annotation": "Collinear syntenic block (6176 bp on '+' strand, 74.49% identity)" }, { "sv_id": "SV-02-SYN", "sv_class": "SYN", "target_name": "NC_000913.3_E_coli_K12", "target_start_1based": 8169, "target_end_1based": 10547, "query_name": "NC_003197.2_S_enterica_LT2", "query_start_1based": 7597, "query_end_1based": 9993, "strand": "+", "event_size_bp": 2379, "blast_identity_pct": 83.9018, "gap_compressed_de": 0.1529, "mapq": 60, "biological_annotation": "Collinear syntenic block (2379 bp on '+' strand, 83.90% identity)" }, { "sv_id": "SV-03-SYN", "sv_class": "SYN", "target_name": "NC_000913.3_E_coli_K12", "target_start_1based": 11992, "target_end_1based": 15138, "query_name": "NC_003197.2_S_enterica_LT2", "query_start_1based": 11423, "query_end_1based": 14574, "strand": "+", "event_size_bp": 3147, "blast_identity_pct": 89.262, "gap_compressed_de": 0.1045, "mapq": 60, "biological_annotation": "Collinear syntenic block (3147 bp on '+' strand, 89.26% identity)" }, { "sv_id": "SV-04-INS_GI", "sv_class": "INS_GI", "target_name": "NC_000913.3_E_coli_K12", "target_start_1based": 15139, "target_end_1based": 20809, "query_name": "NC_003197.2_S_enterica_LT2", "query_start_1based": 14575, "query_end_1based": 52274, "strand": "+", "event_size_bp": 32029, "blast_identity_pct": 0.0, "gap_compressed_de": 1.0, "mapq": 60, "biological_annotation": "Lineage-specific genomic island / insertion (+37700 bp in query vs. 5671 bp in target)" }, { "sv_id": "SV-05-SYN", "sv_class": "SYN", "target_name": "NC_000913.3_E_coli_K12", "target_start_1based": 20810, "target_end_1based": 24785, "query_name": "NC_003197.2_S_enterica_LT2", "query_start_1based": 52275, "query_end_1based": 56255, "strand": "+", "event_size_bp": 3976, "blast_identity_pct": 85.9254, "gap_compressed_de": 0.1366, "mapq": 60, "biological_annotation": "Collinear syntenic block (3976 bp on '+' strand, 85.93% identity)" }, { "sv_id": "SV-06-SYN", "sv_class": "SYN", "target_name": "NC_000913.3_E_coli_K12", "target_start_1based": 25076, "target_end_1based": 25679, "query_name": "NC_003197.2_S_enterica_LT2", "query_start_1based": 56558, "query_end_1based": 57161, "strand": "+", "event_size_bp": 604, "blast_identity_pct": 87.7483, "gap_compressed_de": 0.1225, "mapq": 60, "biological_annotation": "Collinear syntenic block (604 bp on '+' strand, 87.75% identity)" }, { "sv_id": "SV-07-SYN", "sv_class": "SYN", "target_name": "NC_000913.3_E_coli_K12", "target_start_1based": 26120, "target_end_1based": 27223, "query_name": "NC_003197.2_S_enterica_LT2", "query_start_1based": 57638, "query_end_1based": 58742, "strand": "+", "event_size_bp": 1104, "blast_identity_pct": 84.8869, "gap_compressed_de": 0.1511, "mapq": 60, "biological_annotation": "Collinear syntenic block (1104 bp on '+' strand, 84.89% identity)" } ] }, "breakpoint_mae_bp": 3374.3333, "de_divergence_mae": 0.0517, "blast_identity_mae_pct": 3.9222, "sv_concordance_pct": 42.8571, "ascii_synteny_dotplot": "=== Cohort 1 Synteny Dotplot: E. coli K-12 (NC_000913.3) vs. S. enterica LT2 (NC_003197.2) +31.4 kb Island ===\n X-Axis (Target): E. coli K-12 (NC_000913.3) (1 .. 60000 bp) | Y-Axis (Query): S. enterica LT2 (NC_003197.2) (1 .. 60000 bp)\n Legend: [/] Forward Collinear ('+' SYN) [\\] Reverse-Complement Inversion ('-' INV) [+] Segmental Duplication ('S' DUP)\n 60000 +--------------------------------------------------------------+\n 60000 | // |\n | //// |\n | //// |\n | / |\n 45882 | /// |\n | |\n | |\n | |\n 31764 | |\n | |\n | |\n | |\n 17647 | |\n | /// |\n | ///// |\n | //// |\n 3529 | //// |\n 0 | /// |\n 1 +--------------------------------------------------------------+\n 1 60000 bp\n", "ascii_cohort_table": "=== COHORT-01-ENTERIC-60KB-ISLAND — PAF Chains & Structural Variant Calls (Native vs. Minimap2 v2.31-r1302) ===\nEngine / Preset | Prim/Sec | Target Cov (%) | BLAST ID (%) | Gap-Comp de | Mean MAPQ | SV Calls | Time (ms)\n--------------------------+----------+----------------+--------------+-------------+-----------+----------+----------\nNative (multi-scale (w=5,k=13)+(w=10,k=19)) | 3 / 0 | 38.82% | 76.73% | 0.2087 | 60.0 | 5 | 11.98\nMinimap2 (asm5) | 1 / 0 | 0.62% | 99.19% | 0.0081 | 60.0 | 1 | 51.55\nMinimap2 (asm20) | 6 / 0 | 28.98% | 82.20% | 0.1490 | 60.0 | 7 | 45.76\n\n Called Structural Variants & Syntenic Blocks (Native 4-Stage Split-Chain Engine):\n - SV-01-SYN | Target 107..15087 | Query 107..14523 (+) | 14981 bp | ID 74.60% | de 0.2201 | Collinear syntenic block (14981 bp on '+' strand, 74.60% identity)\n - SV-02-INS_GI | Target 15088..17385 | Query 14524..46085 (+) | 29264 bp | ID 0.00% | de 1.0000 | Lineage-specific genomic island / insertion (+31562 bp in query vs. 2298 bp in target)\n - SV-03-SYN | Target 17386..19243 | Query 46086..47946 (+) | 1858 bp | ID 79.27% | de 0.2051 | Collinear syntenic block (1858 bp on '+' strand, 79.27% identity)\n - SV-04-INS_GI | Target 19244..20754 | Query 47947..52217 (+) | 2760 bp | ID 0.00% | de 1.0000 | Lineage-specific genomic island / insertion (+4271 bp in query vs. 1511 bp in target)\n - SV-05-SYN | Target 20755..27205 | Query 52218..58724 (+) | 6451 bp | ID 80.92% | de 0.1834 | Collinear syntenic block (6451 bp on '+' strand, 80.92% identity)\n", "ascii_focal_breakpoint_alignment": "=== Cohort 1 Focal Breakpoint Junction: Conserved thrABC-dnaKJ Block -> +31,428 bp S. enterica Genomic Island ===\n Note: Upstream collinear synteny (84.2% identity) terminates sharply at target 16,748 bp where S. enterica inserts a +31.4 kb island (16,889..48,316 bp) before resuming synteny at nhaA-ribF-ileS (target 16,953 bp -> query 48,317 bp).\n Struct-Track 16718 ==============================|~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ 16777 \n E. coli NC_000913.3 16718 TAGCGCTTATGGGATCACTCCCCGCCGTTGCTCTTACTCGGATTCGTAAGCCGTGAAAAC 16777 \n S. enterica NC_003197.2 16858 CGAGAATCGCGTAGGCTTCAGCGTGGGTCAGATCCCGCATATCGCCAGCATAGCCATGCG 16917 \n Conservation ......*...*....*.....*....**................*.......*..*....\n Divergence (H_site) ██████▁███▁████▁█████▁████▁▁████████████████▁███████▁██▁████\n", "vertebrate_gradient_rows": [] }, { "cohort_id": "COHORT-02-PSEUDOMONAS-INVERSION-SV", "title": "Pseudomonas aeruginosa PAO1 70 kb Chromosomal Inversion (12.0 kb '-') & Excision Deletion (2.4 kb) Benchmark", "domain": "Prokaryota (Pseudomonadales)", "biological_context": "Tests exact split-chain breakpoint detection across a 12,000 bp chromosomal inversion (22,001..34,000 bp on '-' strand) and a 2,400 bp deletion (50,001..52,400 bp) at 99.0% nucleotide identity.", "target_accession": "NC_002516.2:1-70000", "query_accession": "NC_002516.2:1-70000 (SV-INV12k-DEL2.4k)", "target_len_bp": 70000, "query_len_bp": 67600, "native_result": { "engine_name": "bio-suite Native 4-Stage Synteny & SV Engine", "engine_mode": "Native", "preset_label": "multi-scale (w=5,k=13)+(w=10,k=19)", "target_name": "NC_002516.2_P_aeruginosa_PAO1_Ref", "target_len_bp": 70000, "query_name": "NC_002516.2_PAO1_SV_INV12k_DEL2.4k", "query_len_bp": 67600, "minimizer_seeds_target": 36560, "minimizer_seeds_query": 35265, "raw_anchors_count": 27855, "primary_chains_count": 3, "secondary_chains_count": 0, "target_aligned_bases": 69989, "target_coverage_pct": 99.9843, "weighted_blast_identity_pct": 95.595, "weighted_gap_compressed_de": 0.0101, "mean_primary_mapq": 60.0, "wall_clock_ms": 29.0456, "paf_records": [ { "query_name": "NC_002516.2_PAO1_SV_INV12k_DEL2.4k", "query_len": 67600, "query_start_0based": 2, "query_end_0based": 21999, "strand": "+", "target_name": "NC_002516.2_P_aeruginosa_PAO1_Ref", "target_len": 70000, "target_start_0based": 2, "target_end_0based": 21999, "residue_matches": 21777, "alignment_block_len": 21997, "mapq": 60, "alignment_type": "P", "chain_score_f1": 21247, "anchor_count": 8450, "edit_distance_nm": 220, "blast_identity_pct": 98.9999, "gap_compressed_divergence_de": 0.01, "cigar_cg": "21997M", "diff_string_cs": ":48*ag:99*ct:99*ag:99*ga:99*tc:99*ct:99*ga:99*tc:99*ag:99*ga:99*ag:99*ga:99*ag:99*ag:99*ag:99*ga:99*tc:99*ct:99*ga:99*ag:99*ag:99*ct:99*ga:99*tc:99*ga:99*ag:99*ag:99*ct:99*ga:99*ag:99*ct:99*ct:99*tc:99*ct:99*ag:99*ga:99*ga:99*ag:99*ag:99*ga" }, { "query_name": "NC_002516.2_PAO1_SV_INV12k_DEL2.4k", "query_len": 67600, "query_start_0based": 22001, "query_end_0based": 34000, "strand": "-", "target_name": "NC_002516.2_P_aeruginosa_PAO1_Ref", "target_len": 70000, "target_start_0based": 22000, "target_end_0based": 33999, "residue_matches": 11879, "alignment_block_len": 11999, "mapq": 60, "alignment_type": "P", "chain_score_f1": 11516, "anchor_count": 3452, "edit_distance_nm": 120, "blast_identity_pct": 98.9999, "gap_compressed_divergence_de": 0.01, "cigar_cg": "11999M", "diff_string_cs": ":49*ga:99*ga:99*ag:99*ga:99*ga:99*ga:99*ga:99*ga:99*tc:99*ct:99*ga:99*ag:99*ga:99*ct:99*ga:99*ct:99*ct:99*ct:99*ga:99*ga:99*ga:99*ga:99*tc:99*ct:99*ag:99*ct:99*ct:99*ga:99*ga:99*ct:99*ct:99*ct:99*tc:99*ct:99*ga:99*ga:99*ga:99*ct:99*ct:99*ct" }, { "query_name": "NC_002516.2_PAO1_SV_INV12k_DEL2.4k", "query_len": 67600, "query_start_0based": 34004, "query_end_0based": 67597, "strand": "+", "target_name": "NC_002516.2_P_aeruginosa_PAO1_Ref", "target_len": 70000, "target_start_0based": 34004, "target_end_0based": 69997, "residue_matches": 33250, "alignment_block_len": 35993, "mapq": 60, "alignment_type": "P", "chain_score_f1": 32121, "anchor_count": 12852, "edit_distance_nm": 2743, "blast_identity_pct": 92.3791, "gap_compressed_divergence_de": 0.0102, "cigar_cg": "16003M2400D17590M", "diff_string_cs": ":46*ga:99*ct:99*ga:99*ct:99*ga:99*ct:99*ct:99*tc:99*ct:99*tc:99*tc:99*ct:99*tc:99*tc:99*ga:99*ga:99*ct:99*ag:99*ga:99*tc:99*ga:99*ga:99*ga:99*ct:99*ct:99*ct:99*tc:99*tc:99*ct:99*ga:99*tc:99*tc:99*ga:99*ga:99*ag:99*tc:99*ct:99*ag:99*tc:99*ct" } ], "sv_calls": [ { "sv_id": "SV-01-SYN", "sv_class": "SYN", "target_name": "NC_002516.2_P_aeruginosa_PAO1_Ref", "target_start_1based": 3, "target_end_1based": 21999, "query_name": "NC_002516.2_PAO1_SV_INV12k_DEL2.4k", "query_start_1based": 3, "query_end_1based": 21999, "strand": "+", "event_size_bp": 21997, "blast_identity_pct": 98.9999, "gap_compressed_de": 0.01, "mapq": 60, "biological_annotation": "Collinear syntenic block (21997 bp on '+' strand, 99.00% identity)" }, { "sv_id": "SV-02-INV", "sv_class": "INV", "target_name": "NC_002516.2_P_aeruginosa_PAO1_Ref", "target_start_1based": 22001, "target_end_1based": 33999, "query_name": "NC_002516.2_PAO1_SV_INV12k_DEL2.4k", "query_start_1based": 22002, "query_end_1based": 34000, "strand": "-", "event_size_bp": 11999, "blast_identity_pct": 98.9999, "gap_compressed_de": 0.01, "mapq": 60, "biological_annotation": "Chromosomal inversion block (11999 bp on '-' strand, 99.00% identity)" }, { "sv_id": "SV-03-SYN", "sv_class": "SYN", "target_name": "NC_002516.2_P_aeruginosa_PAO1_Ref", "target_start_1based": 34005, "target_end_1based": 69997, "query_name": "NC_002516.2_PAO1_SV_INV12k_DEL2.4k", "query_start_1based": 34005, "query_end_1based": 67597, "strand": "+", "event_size_bp": 35993, "blast_identity_pct": 92.3791, "gap_compressed_de": 0.0102, "mapq": 60, "biological_annotation": "Collinear syntenic block (35993 bp on '+' strand, 92.38% identity)" } ] }, "minimap2_asm5_result": { "engine_name": "Minimap2 v2.31-r1302 (-c --cs -x asm5)", "engine_mode": "Minimap2", "preset_label": "asm5", "target_name": "NC_002516.2_P_aeruginosa_PAO1_Ref", "target_len_bp": 70000, "query_name": "NC_002516.2_PAO1_SV_INV12k_DEL2.4k", "query_len_bp": 67600, "minimizer_seeds_target": 0, "minimizer_seeds_query": 0, "raw_anchors_count": 0, "primary_chains_count": 3, "secondary_chains_count": 0, "target_aligned_bases": 69837, "target_coverage_pct": 99.7671, "weighted_blast_identity_pct": 95.5983, "weighted_gap_compressed_de": 0.01, "mean_primary_mapq": 60.0, "wall_clock_ms": 40.337, "paf_records": [ { "query_name": "NC_002516.2_PAO1_SV_INV12k_DEL2.4k", "query_len": 67600, "query_start_0based": 0, "query_end_0based": 21837, "strand": "+", "target_name": "NC_002516.2_P_aeruginosa_PAO1_Ref", "target_len": 70000, "target_start_0based": 0, "target_end_0based": 21837, "residue_matches": 21619, "alignment_block_len": 21837, "mapq": 60, "alignment_type": "P", "chain_score_f1": 18538, "anchor_count": 1706, "edit_distance_nm": 218, "blast_identity_pct": 99.0017, "gap_compressed_divergence_de": 0.01, "cigar_cg": "21837M", "diff_string_cs": ":50*ag:99*ct:99*ag:99*ga:99*tc:99*ct:99*ga:99*tc:99*ag:99*ga:99*ag:99*ga:99*ag:99*ag:99*ag:99*ga:99*tc:99*ct:99*ga:99*ag:99*ag:99*ct:99*ga:99*tc:99*ga:99*ag:99*ag:99*ct:99*ga:99*ag:99*ct:99*ct:99*tc:99*ct:99*ag:99*ga:99*ga:99*ag:99*ag:99*ga" }, { "query_name": "NC_002516.2_PAO1_SV_INV12k_DEL2.4k", "query_len": 67600, "query_start_0based": 22000, "query_end_0based": 34000, "strand": "-", "target_name": "NC_002516.2_P_aeruginosa_PAO1_Ref", "target_len": 70000, "target_start_0based": 22000, "target_end_0based": 34000, "residue_matches": 11880, "alignment_block_len": 12000, "mapq": 60, "alignment_type": "P", "chain_score_f1": 10216, "anchor_count": 948, "edit_distance_nm": 120, "blast_identity_pct": 99.0, "gap_compressed_divergence_de": 0.01, "cigar_cg": "12000M", "diff_string_cs": ":49*ga:99*ga:99*ag:99*ga:99*ga:99*ga:99*ga:99*ga:99*tc:99*ct:99*ga:99*ag:99*ga:99*ct:99*ga:99*ct:99*ct:99*ct:99*ga:99*ga:99*ga:99*ga:99*tc:99*ct:99*ag:99*ct:99*ct:99*ga:99*ga:99*ct:99*ct:99*ct:99*tc:99*ct:99*ga:99*ga:99*ga:99*ct:99*ct:99*ct" }, { "query_name": "NC_002516.2_PAO1_SV_INV12k_DEL2.4k", "query_len": 67600, "query_start_0based": 34000, "query_end_0based": 67600, "strand": "+", "target_name": "NC_002516.2_P_aeruginosa_PAO1_Ref", "target_len": 70000, "target_start_0based": 34000, "target_end_0based": 70000, "residue_matches": 33264, "alignment_block_len": 36000, "mapq": 60, "alignment_type": "P", "chain_score_f1": 28415, "anchor_count": 2615, "edit_distance_nm": 2736, "blast_identity_pct": 92.4, "gap_compressed_divergence_de": 0.01, "cigar_cg": "16000M2400D17600M", "diff_string_cs": ":50*ga:99*ct:99*ga:99*ct:99*ga:99*ct:99*ct:99*tc:99*ct:99*tc:99*tc:99*ct:99*tc:99*tc:99*ga:99*ga:99*ct:99*ag:99*ga:99*tc:99*ga:99*ga:99*ga:99*ct:99*ct:99*ct:99*tc:99*tc:99*ct:99*ga:99*tc:99*tc:99*ga:99*ga:99*ag:99*tc:99*ct:99*ag:99*tc:99*ct" } ], "sv_calls": [ { "sv_id": "SV-01-SYN", "sv_class": "SYN", "target_name": "NC_002516.2_P_aeruginosa_PAO1_Ref", "target_start_1based": 1, "target_end_1based": 21837, "query_name": "NC_002516.2_PAO1_SV_INV12k_DEL2.4k", "query_start_1based": 1, "query_end_1based": 21837, "strand": "+", "event_size_bp": 21837, "blast_identity_pct": 99.0017, "gap_compressed_de": 0.01, "mapq": 60, "biological_annotation": "Collinear syntenic block (21837 bp on '+' strand, 99.00% identity)" }, { "sv_id": "SV-02-INV", "sv_class": "INV", "target_name": "NC_002516.2_P_aeruginosa_PAO1_Ref", "target_start_1based": 22001, "target_end_1based": 34000, "query_name": "NC_002516.2_PAO1_SV_INV12k_DEL2.4k", "query_start_1based": 22001, "query_end_1based": 34000, "strand": "-", "event_size_bp": 12000, "blast_identity_pct": 99.0, "gap_compressed_de": 0.01, "mapq": 60, "biological_annotation": "Chromosomal inversion block (12000 bp on '-' strand, 99.00% identity)" }, { "sv_id": "SV-03-SYN", "sv_class": "SYN", "target_name": "NC_002516.2_P_aeruginosa_PAO1_Ref", "target_start_1based": 34001, "target_end_1based": 70000, "query_name": "NC_002516.2_PAO1_SV_INV12k_DEL2.4k", "query_start_1based": 34001, "query_end_1based": 67600, "strand": "+", "event_size_bp": 36000, "blast_identity_pct": 92.4, "gap_compressed_de": 0.01, "mapq": 60, "biological_annotation": "Collinear syntenic block (36000 bp on '+' strand, 92.40% identity)" } ] }, "minimap2_asm20_result": { "engine_name": "Minimap2 v2.31-r1302 (-c --cs -x asm20)", "engine_mode": "Minimap2", "preset_label": "asm20", "target_name": "NC_002516.2_P_aeruginosa_PAO1_Ref", "target_len_bp": 70000, "query_name": "NC_002516.2_PAO1_SV_INV12k_DEL2.4k", "query_len_bp": 67600, "minimizer_seeds_target": 0, "minimizer_seeds_query": 0, "raw_anchors_count": 0, "primary_chains_count": 3, "secondary_chains_count": 0, "target_aligned_bases": 69830, "target_coverage_pct": 99.7571, "weighted_blast_identity_pct": 95.5979, "weighted_gap_compressed_de": 0.01, "mean_primary_mapq": 60.0, "wall_clock_ms": 41.1991, "paf_records": [ { "query_name": "NC_002516.2_PAO1_SV_INV12k_DEL2.4k", "query_len": 67600, "query_start_0based": 0, "query_end_0based": 21830, "strand": "+", "target_name": "NC_002516.2_P_aeruginosa_PAO1_Ref", "target_len": 70000, "target_start_0based": 0, "target_end_0based": 21830, "residue_matches": 21612, "alignment_block_len": 21830, "mapq": 60, "alignment_type": "P", "chain_score_f1": 19886, "anchor_count": 3135, "edit_distance_nm": 218, "blast_identity_pct": 99.0014, "gap_compressed_divergence_de": 0.01, "cigar_cg": "21830M", "diff_string_cs": ":50*ag:99*ct:99*ag:99*ga:99*tc:99*ct:99*ga:99*tc:99*ag:99*ga:99*ag:99*ga:99*ag:99*ag:99*ag:99*ga:99*tc:99*ct:99*ga:99*ag:99*ag:99*ct:99*ga:99*tc:99*ga:99*ag:99*ag:99*ct:99*ga:99*ag:99*ct:99*ct:99*tc:99*ct:99*ag:99*ga:99*ga:99*ag:99*ag:99*ga" }, { "query_name": "NC_002516.2_PAO1_SV_INV12k_DEL2.4k", "query_len": 67600, "query_start_0based": 22000, "query_end_0based": 34000, "strand": "-", "target_name": "NC_002516.2_P_aeruginosa_PAO1_Ref", "target_len": 70000, "target_start_0based": 22000, "target_end_0based": 34000, "residue_matches": 11880, "alignment_block_len": 12000, "mapq": 60, "alignment_type": "P", "chain_score_f1": 11019, "anchor_count": 1756, "edit_distance_nm": 120, "blast_identity_pct": 99.0, "gap_compressed_divergence_de": 0.01, "cigar_cg": "12000M", "diff_string_cs": ":49*ga:99*ga:99*ag:99*ga:99*ga:99*ga:99*ga:99*ga:99*tc:99*ct:99*ga:99*ag:99*ga:99*ct:99*ga:99*ct:99*ct:99*ct:99*ga:99*ga:99*ga:99*ga:99*tc:99*ct:99*ag:99*ct:99*ct:99*ga:99*ga:99*ct:99*ct:99*ct:99*tc:99*ct:99*ga:99*ga:99*ga:99*ct:99*ct:99*ct" }, { "query_name": "NC_002516.2_PAO1_SV_INV12k_DEL2.4k", "query_len": 67600, "query_start_0based": 34000, "query_end_0based": 67600, "strand": "+", "target_name": "NC_002516.2_P_aeruginosa_PAO1_Ref", "target_len": 70000, "target_start_0based": 34000, "target_end_0based": 70000, "residue_matches": 33264, "alignment_block_len": 36000, "mapq": 60, "alignment_type": "P", "chain_score_f1": 30740, "anchor_count": 4845, "edit_distance_nm": 2736, "blast_identity_pct": 92.4, "gap_compressed_divergence_de": 0.01, "cigar_cg": "16000M2400D17600M", "diff_string_cs": ":50*ga:99*ct:99*ga:99*ct:99*ga:99*ct:99*ct:99*tc:99*ct:99*tc:99*tc:99*ct:99*tc:99*tc:99*ga:99*ga:99*ct:99*ag:99*ga:99*tc:99*ga:99*ga:99*ga:99*ct:99*ct:99*ct:99*tc:99*tc:99*ct:99*ga:99*tc:99*tc:99*ga:99*ga:99*ag:99*tc:99*ct:99*ag:99*tc:99*ct" } ], "sv_calls": [ { "sv_id": "SV-01-SYN", "sv_class": "SYN", "target_name": "NC_002516.2_P_aeruginosa_PAO1_Ref", "target_start_1based": 1, "target_end_1based": 21830, "query_name": "NC_002516.2_PAO1_SV_INV12k_DEL2.4k", "query_start_1based": 1, "query_end_1based": 21830, "strand": "+", "event_size_bp": 21830, "blast_identity_pct": 99.0014, "gap_compressed_de": 0.01, "mapq": 60, "biological_annotation": "Collinear syntenic block (21830 bp on '+' strand, 99.00% identity)" }, { "sv_id": "SV-02-INV", "sv_class": "INV", "target_name": "NC_002516.2_P_aeruginosa_PAO1_Ref", "target_start_1based": 22001, "target_end_1based": 34000, "query_name": "NC_002516.2_PAO1_SV_INV12k_DEL2.4k", "query_start_1based": 22001, "query_end_1based": 34000, "strand": "-", "event_size_bp": 12000, "blast_identity_pct": 99.0, "gap_compressed_de": 0.01, "mapq": 60, "biological_annotation": "Chromosomal inversion block (12000 bp on '-' strand, 99.00% identity)" }, { "sv_id": "SV-03-SYN", "sv_class": "SYN", "target_name": "NC_002516.2_P_aeruginosa_PAO1_Ref", "target_start_1based": 34001, "target_end_1based": 70000, "query_name": "NC_002516.2_PAO1_SV_INV12k_DEL2.4k", "query_start_1based": 34001, "query_end_1based": 67600, "strand": "+", "event_size_bp": 36000, "blast_identity_pct": 92.4, "gap_compressed_de": 0.01, "mapq": 60, "biological_annotation": "Collinear syntenic block (36000 bp on '+' strand, 92.40% identity)" } ] }, "breakpoint_mae_bp": 28.6667, "de_divergence_mae": 0.0001, "blast_identity_mae_pct": 0.0076, "sv_concordance_pct": 100.0, "ascii_synteny_dotplot": "=== Cohort 2 Synteny Dotplot: P. aeruginosa PAO1 (NC_002516.2) 12.0 kb Chromosomal Inversion & 2.4 kb Deletion ===\n X-Axis (Target): P. aeruginosa PAO1 Ref (70 kb) (1 .. 70000 bp) | Y-Axis (Query): PAO1 Structural Variant Strain (67.6 kb) (1 .. 67600 bp)\n Legend: [/] Forward Collinear ('+' SYN) [\\] Reverse-Complement Inversion ('-' INV) [+] Segmental Duplication ('S' DUP)\n 67600 +--------------------------------------------------------------+\n 67600 | /// |\n | //// |\n | //// |\n | //// |\n 51694 | //// |\n | ///// |\n | //// |\n | //// |\n 35788 | \\ //// |\n | \\\\\\\\ |\n | \\\\\\\\ |\n | / \\\\\\\\\\ |\n 19882 | //// |\n | //// |\n | ///// |\n | //// |\n 3976 | //// |\n 0 | /// |\n 1 +--------------------------------------------------------------+\n 1 70000 bp\n", "ascii_cohort_table": "=== COHORT-02-PSEUDOMONAS-INVERSION-DEL — PAF Chains & Structural Variant Calls (Native vs. Minimap2 v2.31-r1302) ===\nEngine / Preset | Prim/Sec | Target Cov (%) | BLAST ID (%) | Gap-Comp de | Mean MAPQ | SV Calls | Time (ms)\n--------------------------+----------+----------------+--------------+-------------+-----------+----------+----------\nNative (multi-scale (w=5,k=13)+(w=10,k=19)) | 3 / 0 | 99.98% | 95.59% | 0.0101 | 60.0 | 3 | 29.05\nMinimap2 (asm5) | 3 / 0 | 99.77% | 95.60% | 0.0100 | 60.0 | 3 | 40.34\nMinimap2 (asm20) | 3 / 0 | 99.76% | 95.60% | 0.0100 | 60.0 | 3 | 41.20\n\n Called Structural Variants & Syntenic Blocks (Native 4-Stage Split-Chain Engine):\n - SV-01-SYN | Target 3..21999 | Query 3..21999 (+) | 21997 bp | ID 99.00% | de 0.0100 | Collinear syntenic block (21997 bp on '+' strand, 99.00% identity)\n - SV-02-INV | Target 22001..33999 | Query 22002..34000 (-) | 11999 bp | ID 99.00% | de 0.0100 | Chromosomal inversion block (11999 bp on '-' strand, 99.00% identity)\n - SV-03-SYN | Target 34005..69997 | Query 34005..67597 (+) | 35993 bp | ID 92.38% | de 0.0102 | Collinear syntenic block (35993 bp on '+' strand, 92.38% identity)\n", "ascii_focal_breakpoint_alignment": "=== Cohort 2 Focal Inversion Breakpoint: P. aeruginosa PAO1 Position 22,000 bp ('+' Collinear -> '-' 12.0 kb Inversion) ===\n Note: Exact inversion boundary at 22,000 bp: upstream positions 1..22,000 align collinearly on '+' (99.0% identity), whereas positions 22,001..34,000 flip onto the '-' reverse-complement strand before restoring '+' collinearity at 34,001 bp.\n Struct-Track 21970 ==============================|~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ 22029 \n PAO1 Ref NC_002516.2 21970 GTGCTGAACCGTGGCGAATCGCGCGGGACCATCAAGCTGTCGCCGAAGATCCGCAGCACG 22029 \n PAO1 SV Query 21970 GTGCTGAACCGTGGCGAATCGCGCGGGACCACCCGCCAGTACCTCTACGACCATGTTCGC 22029 \n Conservation *******************************.*...*.**..*...*...**........\n Divergence (H_site) ▁▁▁▁▁▁▁▁▁▁▁▁▁▁▁▁▁▁▁▁▁▁▁▁▁▁▁▁▁▁▁█▁███▁█▁▁██▁███▁███▁▁████████\n", "vertebrate_gradient_rows": [] }, { "cohort_id": "COHORT-03-VERTEBRATE-MITO-DIVERGENCE-HORIZON", "title": "12-Vertebrate Mitochondrial Synteny & Minimizer Seed Sensitivity Horizon (asm5 vs. asm10 vs. asm20 vs. Native Multi-Scale)", "domain": "Eukaryota (Vertebrata — Mammalia & Aves)", "biological_context": "Evaluates whole-slice organellar collinearity across a 9% to 26% nucleotide divergence gradient (11 jawed vertebrates aligned against Pan troglodytes NC_001643.1), quantifying exact minimizer seed dropout thresholds for asm5 (w=19,k=19), asm10 (w=19,k=19), asm20 (w=10,k=19), and Native multi-scale (w=5,k=13)+(w=10,k=19).", "target_accession": "NC_001643.1:5200-7800", "query_accession": "11 Jawed Vertebrate Mitogenomes (NC_001645.1 .. NC_010195.2)", "target_len_bp": 2601, "query_len_bp": 2601, "native_result": { "engine_name": "bio-suite Native 4-Stage Synteny & SV Engine", "engine_mode": "Native", "preset_label": "multi-scale (w=5,k=13)+(w=10,k=19)", "target_name": "NC_001643.1_Pan_troglodytes", "target_len_bp": 2601, "query_name": "NC_001645.1:5200-7800", "query_len_bp": 2601, "minimizer_seeds_target": 1327, "minimizer_seeds_query": 1343, "raw_anchors_count": 242, "primary_chains_count": 1, "secondary_chains_count": 0, "target_aligned_bases": 2590, "target_coverage_pct": 99.5771, "weighted_blast_identity_pct": 89.8266, "weighted_gap_compressed_de": 0.1007, "mean_primary_mapq": 60.0, "wall_clock_ms": 0.4625, "paf_records": [ { "query_name": "NC_001645.1:5200-7800", "query_len": 2601, "query_start_0based": 2, "query_end_0based": 2596, "strand": "+", "target_name": "NC_001643.1_Pan_troglodytes", "target_len": 2601, "target_start_0based": 2, "target_end_0based": 2592, "residue_matches": 2331, "alignment_block_len": 2595, "mapq": 60, "alignment_type": "P", "chain_score_f1": 1163, "anchor_count": 232, "edit_distance_nm": 264, "blast_identity_pct": 89.8266, "gap_compressed_divergence_de": 0.1007, "cigar_cg": "108M1I2M4I4M1D2475M", "diff_string_cs": ":29*at:7*ag:18*tc*tc:1*ag:2*ct:26*tc:1*ct:16+t:2+tttt*at*ct:2-t:22*ct:16*ct:20*ct:11*tc:17*ga:26*ga:2*ta:5*at:2*ct:14*ct:8*tc:2*ct:8*ct:5*ct:17*ag:11*ct:20*tc:29*ta:5*ga:53*ga:5*gt:9*ct:10*ag:2*tc:5*tc:2*ct:17*ga:2*ag:2*at:5*ag:20*ga:5*ct:2" } ], "sv_calls": [ { "sv_id": "SV-01-SYN", "sv_class": "SYN", "target_name": "NC_001643.1_Pan_troglodytes", "target_start_1based": 3, "target_end_1based": 2592, "query_name": "NC_001645.1:5200-7800", "query_start_1based": 3, "query_end_1based": 2596, "strand": "+", "event_size_bp": 2590, "blast_identity_pct": 89.8266, "gap_compressed_de": 0.1007, "mapq": 60, "biological_annotation": "Collinear syntenic block (2590 bp on '+' strand, 89.83% identity)" } ] }, "minimap2_asm5_result": { "engine_name": "Minimap2 v2.31-r1302 (-c --cs -x asm5)", "engine_mode": "Minimap2", "preset_label": "asm5", "target_name": "NC_001643.1_Pan_troglodytes", "target_len_bp": 2601, "query_name": "NC_001645.1:5200-7800", "query_len_bp": 2601, "minimizer_seeds_target": 0, "minimizer_seeds_query": 0, "raw_anchors_count": 0, "primary_chains_count": 0, "secondary_chains_count": 0, "target_aligned_bases": 0, "target_coverage_pct": 0.0, "weighted_blast_identity_pct": 0.0, "weighted_gap_compressed_de": 0.0, "mean_primary_mapq": 0.0, "wall_clock_ms": 16.2395, "paf_records": [], "sv_calls": [] }, "minimap2_asm20_result": { "engine_name": "Minimap2 v2.31-r1302 (-c --cs -x asm20)", "engine_mode": "Minimap2", "preset_label": "asm20", "target_name": "NC_001643.1_Pan_troglodytes", "target_len_bp": 2601, "query_name": "NC_001645.1:5200-7800", "query_len_bp": 2601, "minimizer_seeds_target": 0, "minimizer_seeds_query": 0, "raw_anchors_count": 0, "primary_chains_count": 1, "secondary_chains_count": 0, "target_aligned_bases": 2592, "target_coverage_pct": 99.654, "weighted_blast_identity_pct": 89.9114, "weighted_gap_compressed_de": 0.0998, "mean_primary_mapq": 60.0, "wall_clock_ms": 12.4012, "paf_records": [ { "query_name": "NC_001645.1:5200-7800", "query_len": 2601, "query_start_0based": 0, "query_end_0based": 2596, "strand": "+", "target_name": "NC_001643.1_Pan_troglodytes", "target_len": 2601, "target_start_0based": 0, "target_end_0based": 2592, "residue_matches": 2335, "alignment_block_len": 2597, "mapq": 60, "alignment_type": "P", "chain_score_f1": 576, "anchor_count": 56, "edit_distance_nm": 262, "blast_identity_pct": 89.9114, "gap_compressed_divergence_de": 0.0998, "cigar_cg": "112M4I2377M1I7M1D95M", "diff_string_cs": ":31*at:7*ag:18*tc*tc:1*ag:2*ct:26*tc:1*ct:16*ct:1+cttt*at*ct*ct:1*tc:22*ct:16*ct:20*ct:11*tc:17*ga:26*ga:2*ta:5*at:2*ct:14*ct:8*tc:2*ct:8*ct:5*ct:17*ag:11*ct:20*tc:29*ta:5*ga:53*ga:5*gt:9*ct:10*ag:2*tc:5*tc:2*ct:17*ga:2*ag:2*at:5*ag:20*ga:5" } ], "sv_calls": [ { "sv_id": "SV-01-SYN", "sv_class": "SYN", "target_name": "NC_001643.1_Pan_troglodytes", "target_start_1based": 1, "target_end_1based": 2592, "query_name": "NC_001645.1:5200-7800", "query_start_1based": 1, "query_end_1based": 2596, "strand": "+", "event_size_bp": 2592, "blast_identity_pct": 89.9114, "gap_compressed_de": 0.0998, "mapq": 60, "biological_annotation": "Collinear syntenic block (2592 bp on '+' strand, 89.91% identity)" } ] }, "breakpoint_mae_bp": 1.5, "de_divergence_mae": 0.0011, "blast_identity_mae_pct": 0.18, "sv_concordance_pct": 100.0, "ascii_synteny_dotplot": "=== Cohort 3 Synteny Dotplot: Hominid Mitogenome Collinearity (Pan troglodytes NC_001643.1 vs. Gorilla gorilla NC_001645.1) ===\n X-Axis (Target): Pan troglodytes NC_001643.1 (2,601 bp) (1 .. 2601 bp) | Y-Axis (Query): Gorilla gorilla NC_001645.1 (2,601 bp) (1 .. 2601 bp)\n Legend: [/] Forward Collinear ('+' SYN) [\\] Reverse-Complement Inversion ('-' INV) [+] Segmental Duplication ('S' DUP)\n 2601 +--------------------------------------------------------------+\n 2601 | /// |\n | //// |\n | //// |\n | //// |\n 1989 | //// |\n | //// |\n | //// |\n | //// |\n 1377 | //// |\n | //// |\n | //// |\n | //// |\n 765 | //// |\n | //// |\n | //// |\n | //// |\n 153 | //// |\n 0 | /// |\n 1 +--------------------------------------------------------------+\n 1 2601 bp\n", "ascii_cohort_table": "=== COHORT-03-VERTEBRATE-DIVERGENCE-HORIZON (12 Vertebrate Mitogenomes vs. Pan troglodytes NC_001643.1) ===\nQuery Species | Accession | Div (%) | Native Cov | Native de | mm2 asm5 Cov | mm2 asm10 Cov | mm2 asm20 Cov | Preset Sensitivity Regime\n-----------------------+-----------------------+---------+------------+-----------+--------------+---------------+---------------+--------------------------\nGorilla gorilla | NC_001645.1:5200-7800 | 10.17% | 99.58% | 0.1007 | 0.00% | 0.00% | 99.65% | asm5/asm10 Drop Out -> asm20 + Native Rescue\nPongo pygmaeus | NC_001646.1:5200-7800 | 14.73% | 98.00% | 0.1414 | 0.00% | 0.00% | 95.08% | asm5/asm10 Drop Out -> asm20 + Native Rescue\nEquus caballus | NC_001640.1:5200-7800 | 27.03% | 97.92% | 0.2637 | 0.00% | 0.00% | 0.00% | Extreme Cross-Class Divergence (>24%) -> Native (w=5,k=13) Rescue\nBalaenoptera musculus | NC_001601.1:5200-7800 | 26.11% | 71.24% | 0.2571 | 0.00% | 0.00% | 0.00% | Extreme Cross-Class Divergence (>24%) -> Native (w=5,k=13) Rescue\nMus musculus | NC_005089.1:5200-7800 | 27.16% | 97.96% | 0.2606 | 0.00% | 0.00% | 0.00% | Extreme Cross-Class Divergence (>24%) -> Native (w=5,k=13) Rescue\nRattus norvegicus | NC_001665.2:5200-7800 | 28.42% | 97.96% | 0.2743 | 0.00% | 0.00% | 0.00% | Extreme Cross-Class Divergence (>24%) -> Native (w=5,k=13) Rescue\nBos taurus | NC_006853.1:5200-7800 | 21.19% | 65.36% | 0.2086 | 0.00% | 0.00% | 0.00% | Extreme Cross-Class Divergence (>24%) -> Native (w=5,k=13) Rescue\nStruthio camelus | NC_002785.1:5200-7800 | 25.02% | 58.02% | 0.2448 | 0.00% | 0.00% | 0.00% | Extreme Cross-Class Divergence (>24%) -> Native (w=5,k=13) Rescue\nGallus gallus | NC_000878.1:5200-7800 | 28.57% | 86.77% | 0.2825 | 0.00% | 0.00% | 0.00% | Extreme Cross-Class Divergence (>24%) -> Native (w=5,k=13) Rescue\nColumba livia | NC_013978.1:5200-7800 | 32.04% | 80.62% | 0.3175 | 0.00% | 0.00% | 0.00% | Extreme Cross-Class Divergence (>24%) -> Native (w=5,k=13) Rescue\nTaeniopygia guttata | NC_010195.2:5200-7800 | 28.84% | 81.01% | 0.2861 | 0.00% | 0.00% | 0.00% | Extreme Cross-Class Divergence (>24%) -> Native (w=5,k=13) Rescue\n", "ascii_focal_breakpoint_alignment": "=== Cohort 3 Focal Syntenic Alignment: Hominid MT-CO1 / tRNA-Ser / tRNA-Asp Collinear Anchor Block ===\n Note: Full-length 2,601 bp organellar synteny is preserved without structural rearrangement across all 12 jawed vertebrates, while third-codon transitions drive nucleotide divergence from 9.2% (Gorilla) to 25.4% (Aves), illustrating why single-scale k=19 (asm5) seeds drop out above ~14% divergence.\n Struct-Track 1170 ==============================|~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ 1229 \n Pan troglodytes 1170 TAGGTGGCCTAACCGGCATTGTACTAGCAAACTCATCATTAGACATCGTGCTACACGACA 1229 \n Gorilla gorilla 1170 ACTGTAGGCGGCCTAACCGGCATCGTACTAGCAAACTCGTCATTAGATATCGTGCTGCAC 1229 \n Conservation ...**.*.*...*...*......*...*.*.*..*....*....*.....*.....*...\n Divergence (H_site) ███▁▁█▁█▁███▁███▁██████▁███▁█▁█▁██▁████▁████▁█████▁█████▁███\n", "vertebrate_gradient_rows": [ { "query_species": "Gorilla gorilla", "common_name": "Western Gorilla", "clade_group": "Hominidae", "accession_id": "NC_001645.1:5200-7800", "true_nucleotide_divergence_pct": 10.1734, "native_coverage_pct": 99.5771, "native_de": 0.1007, "native_blast_identity_pct": 89.8266, "mm2_asm5_coverage_pct": 0.0, "mm2_asm10_coverage_pct": 0.0, "mm2_asm20_coverage_pct": 99.654, "mm2_asm20_de": 0.0998, "delta_de_abs": 0.0009, "preset_dropout_regime": "asm5/asm10 Drop Out -> asm20 + Native Rescue" }, { "query_species": "Pongo pygmaeus", "common_name": "Bornean Orangutan", "clade_group": "Hominidae", "accession_id": "NC_001646.1:5200-7800", "true_nucleotide_divergence_pct": 14.7344, "native_coverage_pct": 98.0008, "native_de": 0.1414, "native_blast_identity_pct": 85.2656, "mm2_asm5_coverage_pct": 0.0, "mm2_asm10_coverage_pct": 0.0, "mm2_asm20_coverage_pct": 95.0788, "mm2_asm20_de": 0.1402, "delta_de_abs": 0.0012, "preset_dropout_regime": "asm5/asm10 Drop Out -> asm20 + Native Rescue" }, { "query_species": "Equus caballus", "common_name": "Domestic Horse", "clade_group": "Laurasiatheria", "accession_id": "NC_001640.1:5200-7800", "true_nucleotide_divergence_pct": 27.0281, "native_coverage_pct": 97.9239, "native_de": 0.2637, "native_blast_identity_pct": 72.9719, "mm2_asm5_coverage_pct": 0.0, "mm2_asm10_coverage_pct": 0.0, "mm2_asm20_coverage_pct": 0.0, "mm2_asm20_de": 0.0, "delta_de_abs": 0.0, "preset_dropout_regime": "Extreme Cross-Class Divergence (>24%) -> Native (w=5,k=13) Rescue" }, { "query_species": "Balaenoptera musculus", "common_name": "Blue Whale", "clade_group": "Laurasiatheria", "accession_id": "NC_001601.1:5200-7800", "true_nucleotide_divergence_pct": 26.1102, "native_coverage_pct": 71.2418, "native_de": 0.2571, "native_blast_identity_pct": 73.8898, "mm2_asm5_coverage_pct": 0.0, "mm2_asm10_coverage_pct": 0.0, "mm2_asm20_coverage_pct": 0.0, "mm2_asm20_de": 0.0, "delta_de_abs": 0.0, "preset_dropout_regime": "Extreme Cross-Class Divergence (>24%) -> Native (w=5,k=13) Rescue" }, { "query_species": "Mus musculus", "common_name": "House Mouse", "clade_group": "Euarchontoglires", "accession_id": "NC_005089.1:5200-7800", "true_nucleotide_divergence_pct": 27.1562, "native_coverage_pct": 97.9623, "native_de": 0.2606, "native_blast_identity_pct": 72.8438, "mm2_asm5_coverage_pct": 0.0, "mm2_asm10_coverage_pct": 0.0, "mm2_asm20_coverage_pct": 0.0, "mm2_asm20_de": 0.0, "delta_de_abs": 0.0, "preset_dropout_regime": "Extreme Cross-Class Divergence (>24%) -> Native (w=5,k=13) Rescue" }, { "query_species": "Rattus norvegicus", "common_name": "Brown Rat", "clade_group": "Euarchontoglires", "accession_id": "NC_001665.2:5200-7800", "true_nucleotide_divergence_pct": 28.4211, "native_coverage_pct": 97.9623, "native_de": 0.2743, "native_blast_identity_pct": 71.5789, "mm2_asm5_coverage_pct": 0.0, "mm2_asm10_coverage_pct": 0.0, "mm2_asm20_coverage_pct": 0.0, "mm2_asm20_de": 0.0, "delta_de_abs": 0.0, "preset_dropout_regime": "Extreme Cross-Class Divergence (>24%) -> Native (w=5,k=13) Rescue" }, { "query_species": "Bos taurus", "common_name": "Domestic Cattle", "clade_group": "Laurasiatheria", "accession_id": "NC_006853.1:5200-7800", "true_nucleotide_divergence_pct": 21.1854, "native_coverage_pct": 65.3595, "native_de": 0.2086, "native_blast_identity_pct": 78.8146, "mm2_asm5_coverage_pct": 0.0, "mm2_asm10_coverage_pct": 0.0, "mm2_asm20_coverage_pct": 0.0, "mm2_asm20_de": 0.0, "delta_de_abs": 0.0, "preset_dropout_regime": "Extreme Cross-Class Divergence (>24%) -> Native (w=5,k=13) Rescue" }, { "query_species": "Struthio camelus", "common_name": "Common Ostrich", "clade_group": "Aves (Palaeognathae)", "accession_id": "NC_002785.1:5200-7800", "true_nucleotide_divergence_pct": 25.0163, "native_coverage_pct": 58.0161, "native_de": 0.2448, "native_blast_identity_pct": 74.9837, "mm2_asm5_coverage_pct": 0.0, "mm2_asm10_coverage_pct": 0.0, "mm2_asm20_coverage_pct": 0.0, "mm2_asm20_de": 0.0, "delta_de_abs": 0.0, "preset_dropout_regime": "Extreme Cross-Class Divergence (>24%) -> Native (w=5,k=13) Rescue" }, { "query_species": "Gallus gallus", "common_name": "Red Junglefowl", "clade_group": "Aves (Galloanserae)", "accession_id": "NC_000878.1:5200-7800", "true_nucleotide_divergence_pct": 28.5651, "native_coverage_pct": 86.7743, "native_de": 0.2825, "native_blast_identity_pct": 71.4349, "mm2_asm5_coverage_pct": 0.0, "mm2_asm10_coverage_pct": 0.0, "mm2_asm20_coverage_pct": 0.0, "mm2_asm20_de": 0.0, "delta_de_abs": 0.0, "preset_dropout_regime": "Extreme Cross-Class Divergence (>24%) -> Native (w=5,k=13) Rescue" }, { "query_species": "Columba livia", "common_name": "Rock Dove", "clade_group": "Aves (Neoaves)", "accession_id": "NC_013978.1:5200-7800", "true_nucleotide_divergence_pct": 32.0398, "native_coverage_pct": 80.6228, "native_de": 0.3175, "native_blast_identity_pct": 67.9602, "mm2_asm5_coverage_pct": 0.0, "mm2_asm10_coverage_pct": 0.0, "mm2_asm20_coverage_pct": 0.0, "mm2_asm20_de": 0.0, "delta_de_abs": 0.0, "preset_dropout_regime": "Extreme Cross-Class Divergence (>24%) -> Native (w=5,k=13) Rescue" }, { "query_species": "Taeniopygia guttata", "common_name": "Zebra Finch", "clade_group": "Aves (Passeriformes)", "accession_id": "NC_010195.2:5200-7800", "true_nucleotide_divergence_pct": 28.8398, "native_coverage_pct": 81.0073, "native_de": 0.2861, "native_blast_identity_pct": 71.1602, "mm2_asm5_coverage_pct": 0.0, "mm2_asm10_coverage_pct": 0.0, "mm2_asm20_coverage_pct": 0.0, "mm2_asm20_de": 0.0, "delta_de_abs": 0.0, "preset_dropout_regime": "Extreme Cross-Class Divergence (>24%) -> Native (w=5,k=13) Rescue" } ] }, { "cohort_id": "COHORT-04-HUMAN-CHR19-APOE-SV-DUP", "title": "Human GRCh38 Chr19q13.32 APOE-TOMM40-APOC1 60 kb Haplotype Structural Variants (4.0 kb Inversion & 5.0 kb Segmental Duplication)", "domain": "Eukaryota (Homo sapiens — Chromosome 19q13.32)", "biological_context": "Evaluates split-chain structural variant resolution in repeat-rich human nuclear DNA (Alu retrotransposons, G4 quadruplexes, and APOC1/APOC1P1 segmental duplications) across a 4,000 bp inversion and a 5,000 bp tandem segmental duplication.", "target_accession": "NC_000019.10:44890001-44950000", "query_accession": "NC_000019.10:44890001-44950000 (HapSV-INV4k-DUP5k)", "target_len_bp": 60000, "query_len_bp": 65000, "native_result": { "engine_name": "bio-suite Native 4-Stage Synteny & SV Engine", "engine_mode": "Native", "preset_label": "multi-scale (w=5,k=13)+(w=10,k=19)", "target_name": "NC_000019.10_GRCh38_APOE_TOMM40_Ref", "target_len_bp": 60000, "query_name": "NC_000019.10_HapSV_INV4k_DUP5k", "query_len_bp": 65000, "minimizer_seeds_target": 31461, "minimizer_seeds_query": 34069, "raw_anchors_count": 83447, "primary_chains_count": 3, "secondary_chains_count": 13, "target_aligned_bases": 59992, "target_coverage_pct": 99.9867, "weighted_blast_identity_pct": 92.207, "weighted_gap_compressed_de": 0.0013, "mean_primary_mapq": 60.0, "wall_clock_ms": 48.3986, "paf_records": [ { "query_name": "NC_000019.10_HapSV_INV4k_DUP5k", "query_len": 65000, "query_start_0based": 4, "query_end_0based": 37996, "strand": "+", "target_name": "NC_000019.10_GRCh38_APOE_TOMM40_Ref", "target_len": 60000, "target_start_0based": 4, "target_end_0based": 37996, "residue_matches": 37992, "alignment_block_len": 37992, "mapq": 60, "alignment_type": "P", "chain_score_f1": 37992, "anchor_count": 17564, "edit_distance_nm": 0, "blast_identity_pct": 100.0, "gap_compressed_divergence_de": 0.0, "cigar_cg": "37992M", "diff_string_cs": ":37992" }, { "query_name": "NC_000019.10_HapSV_INV4k_DUP5k", "query_len": 65000, "query_start_0based": 52003, "query_end_0based": 55275, "strand": "+", "target_name": "NC_000019.10_GRCh38_APOE_TOMM40_Ref", "target_len": 60000, "target_start_0based": 12003, "target_end_0based": 15275, "residue_matches": 3223, "alignment_block_len": 3272, "mapq": 0, "alignment_type": "S", "chain_score_f1": 3096, "anchor_count": 1138, "edit_distance_nm": 49, "blast_identity_pct": 98.5024, "gap_compressed_divergence_de": 0.015, "cigar_cg": "3272M", "diff_string_cs": ":27*ag:66*ag:66*ag:66*tc:66*ct:66*ct:66*ct:66*tc:66*ga:66*ag:66*ga:66*ga:66*ag:66*ga:66*ga:66*tc:66*ct:66*ga:66*ag:66*ga:66*ct:66*tc:66*ag:66*ct:66*tc:66*tc:66*ct:66*ga:66*ga:66*ct:66*ag:66*tc:66*ga:66*ga:66*tc:66*ga:66*ag:66*ct:66*tc:66*ct" }, { "query_name": "NC_000019.10_HapSV_INV4k_DUP5k", "query_len": 65000, "query_start_0based": 48210, "query_end_0based": 48487, "strand": "+", "target_name": "NC_000019.10_GRCh38_APOE_TOMM40_Ref", "target_len": 60000, "target_start_0based": 25339, "target_end_0based": 25590, "residue_matches": 237, "alignment_block_len": 277, "mapq": 0, "alignment_type": "S", "chain_score_f1": 152, "anchor_count": 35, "edit_distance_nm": 40, "blast_identity_pct": 85.5596, "gap_compressed_divergence_de": 0.0849, "cigar_cg": "104M9I1M4I33M1I30M3I1M5I1M1I2M2I1M1I78M", "diff_string_cs": ":18*ct:7*gc:19*cg:30*tg:26+ctcccaagt:1+gctg:20*ct*ga*ct:1*ca:8+t:13*ct:9*ag:6+cac:1+gtttt:1+g:2+gg:1+a:18*tc:1*tc:10*tc:6*ag:39" }, { "query_name": "NC_000019.10_HapSV_INV4k_DUP5k", "query_len": 65000, "query_start_0based": 40335, "query_end_0based": 40749, "strand": "-", "target_name": "NC_000019.10_GRCh38_APOE_TOMM40_Ref", "target_len": 60000, "target_start_0based": 26974, "target_end_0based": 27389, "residue_matches": 365, "alignment_block_len": 423, "mapq": 0, "alignment_type": "S", "chain_score_f1": 155, "anchor_count": 25, "edit_distance_nm": 58, "blast_identity_pct": 86.2884, "gap_compressed_divergence_de": 0.1289, "cigar_cg": "23M1D2M4D1M1D44M1I1M1D3M1D1M1I4M1I1M2I1M1I54M1I3M1D51M1I217M", "diff_string_cs": ":17*tc:5-a:2-taat:1-a*at:8*ca:1*ca:9*ct*cg*ag:6*ga:13+c:1-g:3-c:1+g:4+t:1+ac:1+c:13*ga:1*cg:4*ga:26*ct:4*ag:1+g:3-c:10*ag*ct*ag:11*gc:2*ca:7*ct:11*tc:3+t:4*ga:10*tc:3*at:1*gc:3*ct:4*ca:5*ct:30*tc:1*tc*ca:1*ac:1*ac:7*ct:19*ga:37*ct*ga:8*ag:19" }, { "query_name": "NC_000019.10_HapSV_INV4k_DUP5k", "query_len": 65000, "query_start_0based": 38141, "query_end_0based": 39301, "strand": "-", "target_name": "NC_000019.10_GRCh38_APOE_TOMM40_Ref", "target_len": 60000, "target_start_0based": 28872, "target_end_0based": 29723, "residue_matches": 777, "alignment_block_len": 1179, "mapq": 0, "alignment_type": "S", "chain_score_f1": 547, "anchor_count": 139, "edit_distance_nm": 402, "blast_identity_pct": 65.9033, "gap_compressed_divergence_de": 0.0761, "cigar_cg": "426M4I4M6I4M1I6M5I1M1I58M311I32M15D1M2D2M2D298M", "diff_string_cs": ":20*ac:7*at:4*ag:7*ct:81*tc:79*ta:32*ct:8*ga:57*gt:19*ag:23*tc:53*ca:16*ga:7+cgcc:4+atgtgt*cg:3+a:6+tcaca:1+g:8*ct*cg:16*ct*ag:12*ct*tg:1*ac*gc*ta:1*cg*tc*ta*ct*tg*tg:1*ag*cg*tc*ct+gacgcctgtaatccca*tg*gt:3*ag*ac*ga:3*cg:3*gc:16-tccagagagggcagg" }, { "query_name": "NC_000019.10_HapSV_INV4k_DUP5k", "query_len": 65000, "query_start_0based": 37998, "query_end_0based": 42001, "strand": "-", "target_name": "NC_000019.10_GRCh38_APOE_TOMM40_Ref", "target_len": 60000, "target_start_0based": 37999, "target_end_0based": 42002, "residue_matches": 4003, "alignment_block_len": 4003, "mapq": 60, "alignment_type": "P", "chain_score_f1": 4003, "anchor_count": 1352, "edit_distance_nm": 0, "blast_identity_pct": 100.0, "gap_compressed_divergence_de": 0.0, "cigar_cg": "4003M", "diff_string_cs": ":4003" }, { "query_name": "NC_000019.10_HapSV_INV4k_DUP5k", "query_len": 65000, "query_start_0based": 26974, "query_end_0based": 27339, "strand": "+", "target_name": "NC_000019.10_GRCh38_APOE_TOMM40_Ref", "target_len": 60000, "target_start_0based": 39251, "target_end_0based": 39615, "residue_matches": 319, "alignment_block_len": 373, "mapq": 0, "alignment_type": "S", "chain_score_f1": 142, "anchor_count": 28, "edit_distance_nm": 54, "blast_identity_pct": 85.5228, "gap_compressed_divergence_de": 0.1355, "cigar_cg": "23M1I2M4I1M1I44M1D1M1I3M1I1M1D4M1D1M2D1M1D54M1D3M1I51M1D167M", "diff_string_cs": ":17*ct:5+a:2+taat:1+a*ta:8*ac:1*ac:9*tc*gc*ga:6*ag:13-c:1+g:3+c:1-g:4-t:1-ac:1-c:13*ag:1*gc:4*ag:26*tc:4*ga:1-g:3+c:10*ga*tc*ga:11*cg:2*ac:7*tc:11*ct:3-t:4*ag:10*ct:3*ta:1*cg:3*tc:4*ac:5*tc:30*ct:1*ct*ac:1*ca:1*ca:7*tc:19*ag:37*tc*ag:8*ga:16" }, { "query_name": "NC_000019.10_HapSV_INV4k_DUP5k", "query_len": 65000, "query_start_0based": 33546, "query_end_0based": 33809, "strand": "+", "target_name": "NC_000019.10_GRCh38_APOE_TOMM40_Ref", "target_len": 60000, "target_start_0based": 39975, "target_end_0based": 40243, "residue_matches": 231, "alignment_block_len": 268, "mapq": 0, "alignment_type": "S", "chain_score_f1": 144, "anchor_count": 21, "edit_distance_nm": 37, "blast_identity_pct": 86.194, "gap_compressed_divergence_de": 0.125, "cigar_cg": "95M5D168M", "diff_string_cs": ":15*gt:1*cg:2*gt:1*tc:22*tc*ct*ag:6*tc*ga:18*at:10*tc:6*ta:2-aaaaa:23*ct:7*tc:24*ac*ga:2*tg:1*at:24*tc:5*cg*tc*ga:7*ct:7*tc:2*gc:3*ga:1*ga:4*tc:1*ct:20*ga:2*ag*ta:15" }, { "query_name": "NC_000019.10_HapSV_INV4k_DUP5k", "query_len": 65000, "query_start_0based": 28872, "query_end_0based": 29339, "strand": "+", "target_name": "NC_000019.10_GRCh38_APOE_TOMM40_Ref", "target_len": 60000, "target_start_0based": 40699, "target_end_0based": 41183, "residue_matches": 451, "alignment_block_len": 484, "mapq": 0, "alignment_type": "S", "chain_score_f1": 365, "anchor_count": 122, "edit_distance_nm": 33, "blast_identity_pct": 93.1818, "gap_compressed_divergence_de": 0.0445, "cigar_cg": "426M4D4M6D4M1D6M5D1M1D26M", "diff_string_cs": ":20*ca:7*ta:4*ga:7*tc:81*ct:79*at:32*tc:8*ag:57*tg:19*ga:23*ct:53*ac:16*ag:7-cgcc:4-atgtgt*gc:3-a:6-tcaca:1-g:8*tc*gc:16" }, { "query_name": "NC_000019.10_HapSV_INV4k_DUP5k", "query_len": 65000, "query_start_0based": 29390, "query_end_0based": 29723, "strand": "+", "target_name": "NC_000019.10_GRCh38_APOE_TOMM40_Ref", "target_len": 60000, "target_start_0based": 41545, "target_end_0based": 41859, "residue_matches": 299, "alignment_block_len": 333, "mapq": 0, "alignment_type": "S", "chain_score_f1": 230, "anchor_count": 54, "edit_distance_nm": 34, "blast_identity_pct": 89.7898, "gap_compressed_divergence_de": 0.0568, "cigar_cg": "13M15I1M2I2M2I298M", "diff_string_cs": ":13+tccagagagggcagg:1+tc:2+ca:16*ga:7*cg:1*gt:25*ct:18*ag:15*tc:14*ag:14*tc:11*ac:4*at:35*ga:1*ct:30*ga:30*at*ga:62" }, { "query_name": "NC_000019.10_HapSV_INV4k_DUP5k", "query_len": 65000, "query_start_0based": 42000, "query_end_0based": 64999, "strand": "+", "target_name": "NC_000019.10_GRCh38_APOE_TOMM40_Ref", "target_len": 60000, "target_start_0based": 42000, "target_end_0based": 59999, "residue_matches": 17934, "alignment_block_len": 22999, "mapq": 60, "alignment_type": "P", "chain_score_f1": 17245, "anchor_count": 8208, "edit_distance_nm": 5065, "blast_identity_pct": 77.9773, "gap_compressed_divergence_de": 0.0037, "cigar_cg": "9890M3459I70M1541I8039M", "diff_string_cs": ":9890+gctgggattacaggcg:21*ct:2*tc:2*gc:2*tc*tg*tg:1*tc*tc*tc*gt*tc*ac*tc*ta:1*gc*tc*tc*ca*ac*at*gt*tc*at*ag:1*ac*tc:1*cg*tc*gc*ac*ac*ac:2*tg*gc*tc*ac*ct+actttctttctgggat*ta:1*tc*gt:2*tg:2*ct*ac:1*tc*ta*tg*ag:1*gc*gc:1*gc:1*ac*gc*ga*at*at*ac" }, { "query_name": "NC_000019.10_HapSV_INV4k_DUP5k", "query_len": 65000, "query_start_0based": 30781, "query_end_0based": 31084, "strand": "+", "target_name": "NC_000019.10_GRCh38_APOE_TOMM40_Ref", "target_len": 60000, "target_start_0based": 42906, "target_end_0based": 43209, "residue_matches": 291, "alignment_block_len": 303, "mapq": 0, "alignment_type": "S", "chain_score_f1": 239, "anchor_count": 85, "edit_distance_nm": 12, "blast_identity_pct": 96.0396, "gap_compressed_divergence_de": 0.0396, "cigar_cg": "303M", "diff_string_cs": ":52*ag:8*ac*ct:11*gt:5*tc:2*tc:79*tc:67*ac:3*ga:2*ct:12*ac:22*ga:28" }, { "query_name": "NC_000019.10_HapSV_INV4k_DUP5k", "query_len": 65000, "query_start_0based": 32343, "query_end_0based": 32791, "strand": "+", "target_name": "NC_000019.10_GRCh38_APOE_TOMM40_Ref", "target_len": 60000, "target_start_0based": 44134, "target_end_0based": 44590, "residue_matches": 427, "alignment_block_len": 457, "mapq": 0, "alignment_type": "S", "chain_score_f1": 314, "anchor_count": 75, "edit_distance_nm": 30, "blast_identity_pct": 93.4354, "gap_compressed_divergence_de": 0.0553, "cigar_cg": "110M4D13M1I44M1D78M1D40M3D162M", "diff_string_cs": ":26*ag:11*gc:16*ag:30*ga:23-aaag:13+g:14*tc:15*tc:13-a:30*tg:47-t:13*ga:14*ga:8*gc:2-gga:2*ag:17*tc*ct:25*cg:2*gc:5*ct:10*ct:7*tg:16*tc:9*tc:59" }, { "query_name": "NC_000019.10_HapSV_INV4k_DUP5k", "query_len": 65000, "query_start_0based": 32835, "query_end_0based": 33522, "strand": "+", "target_name": "NC_000019.10_GRCh38_APOE_TOMM40_Ref", "target_len": 60000, "target_start_0based": 44634, "target_end_0based": 45288, "residue_matches": 616, "alignment_block_len": 692, "mapq": 0, "alignment_type": "S", "chain_score_f1": 386, "anchor_count": 123, "edit_distance_nm": 76, "blast_identity_pct": 89.0173, "gap_compressed_divergence_de": 0.0765, "cigar_cg": "119M1D22M1I3M1D17M1I2M3I212M1I37M1I2M1D19M1D1M1D17M3I1M16I1M4I3M3I3M1I28M1I101M2I1M1I60M", "diff_string_cs": ":36*ga:2*cg:2*ag:56*ct:11*ac:7-a:13*ag:8+t*cg:2-g:2*tc:12*ac:1+c:2+tgt:28*ct:34*ct:121*ag:1*ag:9*tg:14+t:8*gt:6*ag:21+g:2-g:8*ct:10-a:1-c:17+gtt:1+taaaaataaaaataaa:1+aaaa*ca*ca:1+aaa:3+t:12*ga:7*ga:3*cg:3+c:4*tg:12*tg:3*gc:60*ct*ag:9*gc:7+ga" }, { "query_name": "NC_000019.10_HapSV_INV4k_DUP5k", "query_len": 65000, "query_start_0based": 33991, "query_end_0based": 34722, "strand": "+", "target_name": "NC_000019.10_GRCh38_APOE_TOMM40_Ref", "target_len": 60000, "target_start_0based": 45599, "target_end_0based": 46318, "residue_matches": 646, "alignment_block_len": 734, "mapq": 0, "alignment_type": "S", "chain_score_f1": 355, "anchor_count": 60, "edit_distance_nm": 88, "blast_identity_pct": 88.0109, "gap_compressed_divergence_de": 0.1053, "cigar_cg": "282M10I1M3I245M1I37M1I27M1D3M2D121M", "diff_string_cs": ":28*tc:4*ct:17*ag:20*ct:11*ct:2*ca:3*ct:3*ga:1*ac:30*tg:18*at:16*tc:1*ct:4*ta:7*ga:4*ga:16*ct:4*tc*ac:4*tc:6*ga:1*ag:10*tg*ta:12*tg:26*tc:8+cagaggcaca:1+agg*ta:4*tc:5*ct:1*ga:14*ga:3*tc:2*tc:6*tc:40*ca:5*tc*ga:1*ca:9*ga*tc:2*ca*ct:6*ga:1*ta" }, { "query_name": "NC_000019.10_HapSV_INV4k_DUP5k", "query_len": 65000, "query_start_0based": 8906, "query_end_0based": 9153, "strand": "-", "target_name": "NC_000019.10_GRCh38_APOE_TOMM40_Ref", "target_len": 60000, "target_start_0based": 52924, "target_end_0based": 53176, "residue_matches": 224, "alignment_block_len": 259, "mapq": 0, "alignment_type": "S", "chain_score_f1": 150, "anchor_count": 27, "edit_distance_nm": 35, "blast_identity_pct": 86.4865, "gap_compressed_divergence_de": 0.0931, "cigar_cg": "13M6I107M1I74M3D1M5D1M1D2M2D1M1D41M", "diff_string_cs": ":13+tttttt:6*ct:15*tc*ag:21*ga*tc*gt:4*tc:12*cg*gc:3*ga:36+t:26*ag:6*ct:25*cg:6*ca:7-cgc:1-gtgtt:1-g:2-ag:1-a:7*cg:2*ac:30" } ], "sv_calls": [ { "sv_id": "SV-01-SYN", "sv_class": "SYN", "target_name": "NC_000019.10_GRCh38_APOE_TOMM40_Ref", "target_start_1based": 5, "target_end_1based": 37996, "query_name": "NC_000019.10_HapSV_INV4k_DUP5k", "query_start_1based": 5, "query_end_1based": 37996, "strand": "+", "event_size_bp": 37992, "blast_identity_pct": 100.0, "gap_compressed_de": 0.0, "mapq": 60, "biological_annotation": "Collinear syntenic block (37992 bp on '+' strand, 100.00% identity)" }, { "sv_id": "SV-02-INV", "sv_class": "INV", "target_name": "NC_000019.10_GRCh38_APOE_TOMM40_Ref", "target_start_1based": 38000, "target_end_1based": 42002, "query_name": "NC_000019.10_HapSV_INV4k_DUP5k", "query_start_1based": 37999, "query_end_1based": 42001, "strand": "-", "event_size_bp": 4003, "blast_identity_pct": 100.0, "gap_compressed_de": 0.0, "mapq": 60, "biological_annotation": "Chromosomal inversion block (4003 bp on '-' strand, 100.00% identity)" }, { "sv_id": "SV-03-SYN", "sv_class": "SYN", "target_name": "NC_000019.10_GRCh38_APOE_TOMM40_Ref", "target_start_1based": 42001, "target_end_1based": 59999, "query_name": "NC_000019.10_HapSV_INV4k_DUP5k", "query_start_1based": 42001, "query_end_1based": 64999, "strand": "+", "event_size_bp": 17999, "blast_identity_pct": 77.9773, "gap_compressed_de": 0.0037, "mapq": 60, "biological_annotation": "Collinear syntenic block (17999 bp on '+' strand, 77.98% identity)" }, { "sv_id": "SV-04-DUP", "sv_class": "DUP", "target_name": "NC_000019.10_GRCh38_APOE_TOMM40_Ref", "target_start_1based": 12004, "target_end_1based": 15275, "query_name": "NC_000019.10_HapSV_INV4k_DUP5k", "query_start_1based": 52004, "query_end_1based": 55275, "strand": "+", "event_size_bp": 3272, "blast_identity_pct": 98.5024, "gap_compressed_de": 0.015, "mapq": 0, "biological_annotation": "Segmental duplication / paralogous repeat block (3272 bp at 98.50% identity)" }, { "sv_id": "SV-05-DUP", "sv_class": "DUP", "target_name": "NC_000019.10_GRCh38_APOE_TOMM40_Ref", "target_start_1based": 28873, "target_end_1based": 29723, "query_name": "NC_000019.10_HapSV_INV4k_DUP5k", "query_start_1based": 38142, "query_end_1based": 39301, "strand": "-", "event_size_bp": 1160, "blast_identity_pct": 65.9033, "gap_compressed_de": 0.0761, "mapq": 0, "biological_annotation": "Segmental duplication / paralogous repeat block (1160 bp at 65.90% identity)" } ] }, "minimap2_asm5_result": { "engine_name": "Minimap2 v2.31-r1302 (-c --cs -x asm5)", "engine_mode": "Minimap2", "preset_label": "asm5", "target_name": "NC_000019.10_GRCh38_APOE_TOMM40_Ref", "target_len_bp": 60000, "query_name": "NC_000019.10_HapSV_INV4k_DUP5k", "query_len_bp": 65000, "minimizer_seeds_target": 0, "minimizer_seeds_query": 0, "raw_anchors_count": 0, "primary_chains_count": 3, "secondary_chains_count": 0, "target_aligned_bases": 56472, "target_coverage_pct": 94.12, "weighted_blast_identity_pct": 91.8662, "weighted_gap_compressed_de": 0.0, "mean_primary_mapq": 60.0, "wall_clock_ms": 59.6732, "paf_records": [ { "query_name": "NC_000019.10_HapSV_INV4k_DUP5k", "query_len": 65000, "query_start_0based": 0, "query_end_0based": 38000, "strand": "+", "target_name": "NC_000019.10_GRCh38_APOE_TOMM40_Ref", "target_len": 60000, "target_start_0based": 0, "target_end_0based": 38000, "residue_matches": 38000, "alignment_block_len": 38000, "mapq": 60, "alignment_type": "P", "chain_score_f1": 37348, "anchor_count": 3856, "edit_distance_nm": 0, "blast_identity_pct": 100.0, "gap_compressed_divergence_de": 0.0, "cigar_cg": "38000M", "diff_string_cs": ":38000" }, { "query_name": "NC_000019.10_HapSV_INV4k_DUP5k", "query_len": 65000, "query_start_0based": 39393, "query_end_0based": 39865, "strand": "-", "target_name": "NC_000019.10_GRCh38_APOE_TOMM40_Ref", "target_len": 60000, "target_start_0based": 40135, "target_end_0based": 40607, "residue_matches": 472, "alignment_block_len": 472, "mapq": 60, "alignment_type": "P", "chain_score_f1": 0, "anchor_count": 0, "edit_distance_nm": 0, "blast_identity_pct": 100.0, "gap_compressed_divergence_de": 0.0, "cigar_cg": "472M", "diff_string_cs": ":472" }, { "query_name": "NC_000019.10_HapSV_INV4k_DUP5k", "query_len": 65000, "query_start_0based": 42000, "query_end_0based": 65000, "strand": "+", "target_name": "NC_000019.10_GRCh38_APOE_TOMM40_Ref", "target_len": 60000, "target_start_0based": 42000, "target_end_0based": 60000, "residue_matches": 18000, "alignment_block_len": 23000, "mapq": 60, "alignment_type": "P", "chain_score_f1": 17628, "anchor_count": 1820, "edit_distance_nm": 5000, "blast_identity_pct": 78.2609, "gap_compressed_divergence_de": 0.0001, "cigar_cg": "10000M5000I8000M", "diff_string_cs": ":10000+aatcacttgaacctgggaggcggaggttgcggggagttgagactgcaccactgcactccagcctgggtgacagagtgagactccatctccaaaaaaagtaaaatgaaataaataaataaatgttaaaaaaaatctggtggagcatctgatgggtgtttgggccaggctggagctttgtccatcccctcttatttttctgcacttgactctcttatttttctgagactggtccc" } ], "sv_calls": [ { "sv_id": "SV-01-SYN", "sv_class": "SYN", "target_name": "NC_000019.10_GRCh38_APOE_TOMM40_Ref", "target_start_1based": 1, "target_end_1based": 38000, "query_name": "NC_000019.10_HapSV_INV4k_DUP5k", "query_start_1based": 1, "query_end_1based": 38000, "strand": "+", "event_size_bp": 38000, "blast_identity_pct": 100.0, "gap_compressed_de": 0.0, "mapq": 60, "biological_annotation": "Collinear syntenic block (38000 bp on '+' strand, 100.00% identity)" }, { "sv_id": "SV-02-INV", "sv_class": "INV", "target_name": "NC_000019.10_GRCh38_APOE_TOMM40_Ref", "target_start_1based": 40136, "target_end_1based": 40607, "query_name": "NC_000019.10_HapSV_INV4k_DUP5k", "query_start_1based": 39394, "query_end_1based": 39865, "strand": "-", "event_size_bp": 472, "blast_identity_pct": 100.0, "gap_compressed_de": 0.0, "mapq": 60, "biological_annotation": "Chromosomal inversion block (472 bp on '-' strand, 100.00% identity)" }, { "sv_id": "SV-03-SYN", "sv_class": "SYN", "target_name": "NC_000019.10_GRCh38_APOE_TOMM40_Ref", "target_start_1based": 42001, "target_end_1based": 60000, "query_name": "NC_000019.10_HapSV_INV4k_DUP5k", "query_start_1based": 42001, "query_end_1based": 65000, "strand": "+", "event_size_bp": 18000, "blast_identity_pct": 78.2609, "gap_compressed_de": 0.0001, "mapq": 60, "biological_annotation": "Collinear syntenic block (18000 bp on '+' strand, 78.26% identity)" } ] }, "minimap2_asm20_result": { "engine_name": "Minimap2 v2.31-r1302 (-c --cs -x asm20)", "engine_mode": "Minimap2", "preset_label": "asm20", "target_name": "NC_000019.10_GRCh38_APOE_TOMM40_Ref", "target_len_bp": 60000, "query_name": "NC_000019.10_HapSV_INV4k_DUP5k", "query_len_bp": 65000, "minimizer_seeds_target": 0, "minimizer_seeds_query": 0, "raw_anchors_count": 0, "primary_chains_count": 3, "secondary_chains_count": 0, "target_aligned_bases": 58356, "target_coverage_pct": 97.26, "weighted_blast_identity_pct": 92.1081, "weighted_gap_compressed_de": 0.0, "mean_primary_mapq": 60.0, "wall_clock_ms": 65.0754, "paf_records": [ { "query_name": "NC_000019.10_HapSV_INV4k_DUP5k", "query_len": 65000, "query_start_0based": 0, "query_end_0based": 38000, "strand": "+", "target_name": "NC_000019.10_GRCh38_APOE_TOMM40_Ref", "target_len": 60000, "target_start_0based": 0, "target_end_0based": 38000, "residue_matches": 38000, "alignment_block_len": 38000, "mapq": 60, "alignment_type": "P", "chain_score_f1": 37484, "anchor_count": 6877, "edit_distance_nm": 0, "blast_identity_pct": 100.0, "gap_compressed_divergence_de": 0.0, "cigar_cg": "38000M", "diff_string_cs": ":38000" }, { "query_name": "NC_000019.10_HapSV_INV4k_DUP5k", "query_len": 65000, "query_start_0based": 38797, "query_end_0based": 41153, "strand": "-", "target_name": "NC_000019.10_GRCh38_APOE_TOMM40_Ref", "target_len": 60000, "target_start_0based": 38847, "target_end_0based": 41203, "residue_matches": 2356, "alignment_block_len": 2356, "mapq": 60, "alignment_type": "P", "chain_score_f1": 0, "anchor_count": 0, "edit_distance_nm": 0, "blast_identity_pct": 100.0, "gap_compressed_divergence_de": 0.0, "cigar_cg": "2356M", "diff_string_cs": ":2356" }, { "query_name": "NC_000019.10_HapSV_INV4k_DUP5k", "query_len": 65000, "query_start_0based": 42000, "query_end_0based": 65000, "strand": "+", "target_name": "NC_000019.10_GRCh38_APOE_TOMM40_Ref", "target_len": 60000, "target_start_0based": 42000, "target_end_0based": 60000, "residue_matches": 18000, "alignment_block_len": 23000, "mapq": 60, "alignment_type": "P", "chain_score_f1": 17561, "anchor_count": 3222, "edit_distance_nm": 5000, "blast_identity_pct": 78.2609, "gap_compressed_divergence_de": 0.0001, "cigar_cg": "10000M5000I8000M", "diff_string_cs": ":10000+aatcacttgaacctgggaggcggaggttgcggggagttgagactgcaccactgcactccagcctgggtgacagagtgagactccatctccaaaaaaagtaaaatgaaataaataaataaatgttaaaaaaaatctggtggagcatctgatgggtgtttgggccaggctggagctttgtccatcccctcttatttttctgcacttgactctcttatttttctgagactggtccc" } ], "sv_calls": [ { "sv_id": "SV-01-SYN", "sv_class": "SYN", "target_name": "NC_000019.10_GRCh38_APOE_TOMM40_Ref", "target_start_1based": 1, "target_end_1based": 38000, "query_name": "NC_000019.10_HapSV_INV4k_DUP5k", "query_start_1based": 1, "query_end_1based": 38000, "strand": "+", "event_size_bp": 38000, "blast_identity_pct": 100.0, "gap_compressed_de": 0.0, "mapq": 60, "biological_annotation": "Collinear syntenic block (38000 bp on '+' strand, 100.00% identity)" }, { "sv_id": "SV-02-INV", "sv_class": "INV", "target_name": "NC_000019.10_GRCh38_APOE_TOMM40_Ref", "target_start_1based": 38848, "target_end_1based": 41203, "query_name": "NC_000019.10_HapSV_INV4k_DUP5k", "query_start_1based": 38798, "query_end_1based": 41153, "strand": "-", "event_size_bp": 2356, "blast_identity_pct": 100.0, "gap_compressed_de": 0.0, "mapq": 60, "biological_annotation": "Chromosomal inversion block (2356 bp on '-' strand, 100.00% identity)" }, { "sv_id": "SV-03-SYN", "sv_class": "SYN", "target_name": "NC_000019.10_GRCh38_APOE_TOMM40_Ref", "target_start_1based": 42001, "target_end_1based": 60000, "query_name": "NC_000019.10_HapSV_INV4k_DUP5k", "query_start_1based": 42001, "query_end_1based": 65000, "strand": "+", "event_size_bp": 18000, "blast_identity_pct": 78.2609, "gap_compressed_de": 0.0001, "mapq": 60, "biological_annotation": "Collinear syntenic block (18000 bp on '+' strand, 78.26% identity)" } ] }, "breakpoint_mae_bp": 590.0, "de_divergence_mae": 0.0012, "blast_identity_mae_pct": 0.0945, "sv_concordance_pct": 40.0, "ascii_synteny_dotplot": "=== Cohort 4 Synteny Dotplot: Human GRCh38 Chr19q13.32 APOE-TOMM40-APOC1 (NC_000019.10) 4.0 kb Inversion & 5.0 kb Segmental Duplication ===\n X-Axis (Target): Human Chr19q13.32 Ref (60.0 kb) (1 .. 60000 bp) | Y-Axis (Query): Chr19 Structural Variant Haplotype (65.0 kb) (1 .. 65000 bp)\n Legend: [/] Forward Collinear ('+' SYN) [\\] Reverse-Complement Inversion ('-' INV) [+] Segmental Duplication ('S' DUP)\n 65000 +--------------------------------------------------------------+\n 65000 | // |\n | //// |\n | //// |\n | ++++ //// |\n 49705 | + //// |\n | //// |\n | + \\\\\\ /// |\n | ++ // \\\\\\ |\n 34411 | //// ++ +++ |\n | ///// +++++ |\n | //// + |\n | //// |\n 19117 | ///// |\n | ///// |\n | //// |\n | //// + |\n 3823 | ///// |\n 0 | /// |\n 1 +--------------------------------------------------------------+\n 1 60000 bp\n", "ascii_cohort_table": "=== COHORT-04-HUMAN-CHR19-APOE-INV-DUP — PAF Chains & Structural Variant Calls (Native vs. Minimap2 v2.31-r1302) ===\nEngine / Preset | Prim/Sec | Target Cov (%) | BLAST ID (%) | Gap-Comp de | Mean MAPQ | SV Calls | Time (ms)\n--------------------------+----------+----------------+--------------+-------------+-----------+----------+----------\nNative (multi-scale (w=5,k=13)+(w=10,k=19)) | 3 / 13 | 99.99% | 92.21% | 0.0013 | 60.0 | 5 | 48.40\nMinimap2 (asm5) | 3 / 0 | 94.12% | 91.87% | 0.0000 | 60.0 | 3 | 59.67\nMinimap2 (asm20) | 3 / 0 | 97.26% | 92.11% | 0.0000 | 60.0 | 3 | 65.08\n\n Called Structural Variants & Syntenic Blocks (Native 4-Stage Split-Chain Engine):\n - SV-01-SYN | Target 5..37996 | Query 5..37996 (+) | 37992 bp | ID 100.00% | de 0.0000 | Collinear syntenic block (37992 bp on '+' strand, 100.00% identity)\n - SV-02-INV | Target 38000..42002 | Query 37999..42001 (-) | 4003 bp | ID 100.00% | de 0.0000 | Chromosomal inversion block (4003 bp on '-' strand, 100.00% identity)\n - SV-03-SYN | Target 42001..59999 | Query 42001..64999 (+) | 17999 bp | ID 77.98% | de 0.0037 | Collinear syntenic block (17999 bp on '+' strand, 77.98% identity)\n - SV-04-DUP | Target 12004..15275 | Query 52004..55275 (+) | 3272 bp | ID 98.50% | de 0.0150 | Segmental duplication / paralogous repeat block (3272 bp at 98.50% identity)\n - SV-05-DUP | Target 28873..29723 | Query 38142..39301 (-) | 1160 bp | ID 65.90% | de 0.0761 | Segmental duplication / paralogous repeat block (1160 bp at 65.90% identity)\n", "ascii_focal_breakpoint_alignment": "=== Cohort 4 Focal Breakpoint Junction: Human Chr19q13.32 Position 52,000 bp (+5.0 kb Segmental Duplication Insertion Boundary) ===\n Note: At target coordinate 52,000 bp (44,942,000 bp on GRCh38 Chr19), the query haplotype inserts a 5,000 bp paralogous segmental duplication (98.5% identity to target 12,001..17,000 bp, captured as both an INS_GI split-chain gap and a secondary 'tp:A:S' DUP alignment).\n Struct-Track 51970 ==============================|~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ 52029 \n GRCh38 Chr19 Ref 51970 CACTTTATGGAGGAGGAAACTGAGGCCAGCAAATGAGGCTGTCATGGGAGGTGGAGACAG 52029 \n Chr19 SV Haplotype 51970 CACTTTATGGAGGAGGAAACTGAGGCCAGCAAATCACTTGAACCTGGGAGGCGGAGGTTG 52029 \n Conservation **********************************.*......*.*******.****...*\n Divergence (H_site) ▁▁▁▁▁▁▁▁▁▁▁▁▁▁▁▁▁▁▁▁▁▁▁▁▁▁▁▁▁▁▁▁▁▁█▁██████▁█▁▁▁▁▁▁▁█▁▁▁▁███▁\n", "vertebrate_gradient_rows": [] } ], "ascii_master_summary_table": "=== TOOL PAIR 7 MASTER CALIBRATION SUMMARY: Native 4-Stage Synteny & SV Engine vs. Minimap2 v2.31-r1302 ===\nCohort ID | Target / Query Loci | Native Cov | mm2 Cov | Breakpoint MAE | de MAE | SV Concordance\n---------------------------------------------+----------------------------------+------------+----------+----------------+---------+---------------\nCOHORT-01-ENTERIC-60KB-ISLAND-SYNTENY | NC_000913.3:1-60000 | 38.82% | 28.98% | 3374.33 bp | 0.0517 | 42.86%\nCOHORT-02-PSEUDOMONAS-INVERSION-SV | NC_002516.2:1-70000 | 99.98% | 99.77% | 28.67 bp | 0.0001 | 100.00%\nCOHORT-03-VERTEBRATE-MITO-DIVERGENCE-HORIZON | NC_001643.1:5200-7800 | 99.58% | 99.65% | 1.50 bp | 0.0011 | 100.00%\nCOHORT-04-HUMAN-CHR19-APOE-SV-DUP | NC_000019.10:44890001-44950000 | 99.99% | 97.26% | 590.00 bp | 0.0012 | 40.00%\n", "total_duration_ms": 1004 }