[ DEBUG pbstarphase::hla::realigner] Mapping { query_name: None, query_len: Some(21500), query_start: 8376, query_end: 21332, strand: Forward, target_name: Some("HLA:HLA25943"), target_len: 12971, target_start: 0, target_end: 12971, match_len: 12950, block_len: 12973, mapq: 1, is_primary: true, is_supplementary: false, alignment: Some(Alignment { nm: 23, cigar: Some([(101, 0), (1, 2), (1826, 0), (1, 2), (2176, 0), (1, 1), (216, 0), (1, 2), (1681, 0), (1, 2), (520, 0), (1, 2), (2353, 0), (9, 2), (617, 0), (1, 2), (228, 0), (1, 2), (1227, 0), (1, 1), (1835, 0), (1, 2), (174, 0)]), cigar_str: Some("8376S101M1D1826M1D2176M1I216M1D1681M1D520M1D2353M9D617M1D228M1D1227M1I1835M1D174M168S"), md: Some("101^C1826^A1632T759^T1681^A520^T1219C1133^TTCCCGAAG225C391^C228^G2936G125^G174"), cs: Some(":101-c:1826-a:1632*tg:543+t:216-t:1681-a:520-t:1219*ct:1133-ttcccgaag:225*ct:391-c:228-g:1227+a:1709*gc:125-g:174"), alignment_score: Some(12830) }) } [ DEBUG pbstarphase::hla::realigner] MappingStats { seq_len: 12971, nm: 23, unmapped: 0, clipped_start: None, clipped_end: None } [ DEBUG pbstarphase::hla::realigner] [ DEBUG pbstarphase::hla::realigner] Mapping { query_name: None, query_len: Some(21500), query_start: 8247, query_end: 21138, strand: Forward, target_name: Some("HLA:HLA23724"), target_len: 12902, target_start: 0, target_end: 12902, match_len: 12885, block_len: 12906, mapq: 0, is_primary: false, is_supplementary: false, alignment: Some(Alignment { nm: 21, cigar: Some([(230, 0), (1, 2), (4002, 0), (1, 1), (216, 0), (1, 2), (286, 0), (2, 1), (1393, 0), (1, 2), (520, 0), (1, 2), (2353, 0), (9, 2), (617, 0), (1, 2), (228, 0), (1, 2), (1227, 0), (1, 1), (1815, 0)]), cigar_str: Some("8247S230M1D4002M1I216M1D286M2I1393M1D520M1D2353M9D617M1D228M1D1227M1I1815M362S"), md: Some("230^C4218^T1679^A520^T2118C234^TTCCCGAAG225C391^C228^G3042"), cs: Some(":230-c:4002+t:216-t:286+ca:1393-a:520-t:2118*ct:234-ttcccgaag:225*ct:391-c:228-g:1227+a:1815"), alignment_score: Some(12779) }) } [ DEBUG pbstarphase::hla::realigner] MappingStats { seq_len: 12902, nm: 21, unmapped: 0, clipped_start: None, clipped_end: None } [ DEBUG pbstarphase::hla::realigner] [ DEBUG pbstarphase::hla::realigner] Mapping { query_name: None, query_len: Some(21500), query_start: 7761, query_end: 8684, strand: Forward, target_name: Some("HLA:HLA00915"), target_len: 13445, target_start: 0, target_end: 919, match_len: 819, block_len: 936, mapq: 1, is_primary: false, is_supplementary: true, alignment: Some(Alignment { nm: 117, cigar: Some([(24, 0), (2, 2), (11, 0), (1, 1), (68, 0), (3, 1), (36, 0), (2, 2), (128, 0), (5, 1), (41, 0), (1, 1), (56, 0), (2, 2), (36, 0), (1, 1), (107, 0), (2, 2), (115, 0), (2, 2), (32, 0), (1, 1), (26, 0), (1, 2), (16, 0), (1, 2), (5, 0), (1, 2), (96, 0), (5, 1), (109, 0)]), cigar_str: Some("7761S24M2D11M1I68M3I36M2D128M5I41M1I56M2D36M1I107M2D115M2D32M1I26M1D16M1D5M1D96M5I109M12816S"), md: Some("8A15^CA0T3G23A10A10T18G0T31T10T1^TA0C11C2A31A1C7T2G1G16T4A1T3T0A1A4C6G2A33T20G0G12T8T0G0G0G0C2C5A9C3G2C0A7^CA0C6G13T0G7G10T12T17A37G17A0T13^AG12G8C3G10T7T12G5A2A15A19G3C0C4G2^TG1G20C2G7G3G6A13^G6T2C6^A1C3^T5T6C16G29A8T2A6C42A16G23G0C28G12"), cs: Some(":8*ag:15-ca*tg:3*gt:6+t:17*ag:10*ag:10*tc:18*gc*ta:8+tta:23*tc:10*tg:1-ta*ca:11*ct:2*at:31*at:1*cg:7*tc:2*ga:1*ga:16*tc:4*ac:1*tc:3*tc*ag:1*ag:4*ct:6*gt:2*ag:19+ttttc:14*tc:20*ga*ga:4+c:8*ta:8*tc*ga*gt*gc*ct:2*cg:5*ac:9*ct:3*ga:2*ct*ag:7-ca*ct:6*gt:13*tc*ga:7*ga:5+t:5*tc:12*ta:17*ag:37*ga:17*ag*ta:13-ag:12*ga:8*ct:3*gt:10*tc:7*tc:12*gc:5*ag:2*ag:15*ac:19*ga:3*ct*ct:4*ga:2-tg:1*gt:20*cg:2*ga:6+a:1*ga:3*gc:6*ag:13-g:6*tc:2*ct:6-a:1*ct:3-t:5*tc:6*cg:16*ga:29*ag:8*ta:2*ag:6*cg:17+ctttt:25*ag:16*ga:23*ga*ct:28*ga:12"), alignment_score: Some(321) }) } [ DEBUG pbstarphase::hla::realigner] MappingStats { seq_len: 13445, nm: 117, unmapped: 12526, clipped_start: None, clipped_end: None } [ DEBUG pbstarphase::hla::realigner] [ DEBUG pbstarphase::hla::realigner] m64043_200405_180950/35784429/ccs => "HLA-DRB3" "01:01:02:05" => Forward, 0.00163=(21+0)/12902, 0.00163=21/(12902-0) [ DEBUG pbstarphase::hla::realigner] 8247..21138 aligned to 01:01:02:05:0..12902