chhylp123 / hifiasm

Hifiasm: a haplotype-resolved assembler for accurate Hifi reads
MIT License
534 stars 87 forks source link

uc_block_qse_cutoff assertion fail #441

Closed olekto closed 1 year ago

olekto commented 1 year ago

Hi Haoyu, I get this assertion failure with both 0.19.4 and 0.19.0: hifiasm: gfa_ut.cpp:4917: void uc_block_qse_cutoff(uint32_t, uint32_t, uc_block_t, uint32_t, uint32_t): Assertion `(nqs) >= x->qs && (*nqe) <= x->qe' failed.

I reused some bin files from 0.18.9 ec.bin, ovlp.reverse.bin, ovlp.source.bin , and ul.ovlp.bin. Do you know what goes wrong? I've attached the full log below. The species is a salmon.

Thank you.

Ole

Full log:

> [M::ha_analyze_count] lowest: count[5] = 10180782
> [M::ha_analyze_count] highest: count[29] = 57795668
> [M::ha_hist_line]     2: *************************************** 22315710
> [M::ha_hist_line]     3: ****************** 10450141
> [M::ha_hist_line]     4: **************** 9516574
> [M::ha_hist_line]     5: ****************** 10180782
> [M::ha_hist_line]     6: ******************** 11827351
> [M::ha_hist_line]     7: ************************ 13654493
> [M::ha_hist_line]     8: *************************** 15780579
> [M::ha_hist_line]     9: ******************************* 18026773
> [M::ha_hist_line]    10: *********************************** 20396887
> [M::ha_hist_line]    11: *************************************** 22601527
> [M::ha_hist_line]    12: ******************************************* 24823568
> [M::ha_hist_line]    13: ********************************************** 26696094
> [M::ha_hist_line]    14: ************************************************** 28627814
> [M::ha_hist_line]    15: **************************************************** 30006952
> [M::ha_hist_line]    16: ****************************************************** 31384688
> [M::ha_hist_line]    17: ******************************************************** 32544899
> [M::ha_hist_line]    18: *********************************************************** 34192213
> [M::ha_hist_line]    19: ************************************************************** 35966809
> [M::ha_hist_line]    20: ****************************************************************** 37927182
> [M::ha_hist_line]    21: ********************************************************************** 40290475
> [M::ha_hist_line]    22: ************************************************************************** 42924748
> [M::ha_hist_line]    23: ******************************************************************************* 45928415
> [M::ha_hist_line]    24: ************************************************************************************* 48893939
> [M::ha_hist_line]    25: ***************************************************************************************** 51632336
> [M::ha_hist_line]    26: ********************************************************************************************** 54218156
> [M::ha_hist_line]    27: ************************************************************************************************** 56358880
> [M::ha_hist_line]    28: **************************************************************************************************** 57590182
> [M::ha_hist_line]    29: **************************************************************************************************** 57795668
> [M::ha_hist_line]    30: *************************************************************************************************** 57459775
> [M::ha_hist_line]    31: ************************************************************************************************ 55590501
> [M::ha_hist_line]    32: ******************************************************************************************* 52864290
> [M::ha_hist_line]    33: ************************************************************************************* 49354709
> [M::ha_hist_line]    34: ****************************************************************************** 45283335
> [M::ha_hist_line]    35: ********************************************************************** 40741856
> [M::ha_hist_line]    36: ************************************************************** 35946567
> [M::ha_hist_line]    37: ****************************************************** 31047059
> [M::ha_hist_line]    38: ********************************************** 26464075
> [M::ha_hist_line]    39: ************************************** 22115917
> [M::ha_hist_line]    40: ******************************** 18304186
> [M::ha_hist_line]    41: ************************** 14879701
> [M::ha_hist_line]    42: ********************* 11970202
> [M::ha_hist_line]    43: **************** 9467196
> [M::ha_hist_line]    44: ************* 7645373
> [M::ha_hist_line]    45: *********** 6118316
> [M::ha_hist_line]    46: ********* 4957786
> [M::ha_hist_line]    47: ******* 4058191
> [M::ha_hist_line]    48: ****** 3430742
> [M::ha_hist_line]    49: ***** 2926025
> [M::ha_hist_line]    50: ***** 2610654
> [M::ha_hist_line]    51: **** 2355391
> [M::ha_hist_line]    52: **** 2202660
> [M::ha_hist_line]    53: **** 2057439
> [M::ha_hist_line]    54: *** 1986575
> [M::ha_hist_line]    55: *** 1880970
> [M::ha_hist_line]    56: *** 1850757
> [M::ha_hist_line]    57: *** 1787708
> [M::ha_hist_line]    58: *** 1743432
> [M::ha_hist_line]    59: *** 1663175
> [M::ha_hist_line]    60: *** 1630869
> [M::ha_hist_line]    61: *** 1552109
> [M::ha_hist_line]    62: *** 1519524
> [M::ha_hist_line]    63: *** 1455889
> [M::ha_hist_line]    64: ** 1404575
> [M::ha_hist_line]    65: ** 1322413
> [M::ha_hist_line]    66: ** 1270187
> [M::ha_hist_line]    67: ** 1188612
> [M::ha_hist_line]    68: ** 1129429
> [M::ha_hist_line]    69: ** 1056399
> [M::ha_hist_line]    70: ** 1001915
> [M::ha_hist_line]    71: ** 938552
> [M::ha_hist_line]    72: ** 892939
> [M::ha_hist_line]    73: * 842465
> [M::ha_hist_line]    74: * 799335
> [M::ha_hist_line]    75: * 753467
> [M::ha_hist_line]    76: * 713958
> [M::ha_hist_line]    77: * 674452
> [M::ha_hist_line]    78: * 651155
> [M::ha_hist_line]    79: * 621715
> [M::ha_hist_line]    80: * 594499
> [M::ha_hist_line]    81: * 570188
> [M::ha_hist_line]    82: * 551561
> [M::ha_hist_line]    83: * 533443
> [M::ha_hist_line]    84: * 524029
> [M::ha_hist_line]    85: * 496820
> [M::ha_hist_line]    86: * 483086
> [M::ha_hist_line]    87: * 476582
> [M::ha_hist_line]    88: * 461984
> [M::ha_hist_line]    89: * 450542
> [M::ha_hist_line]    90: * 443350
> [M::ha_hist_line]    91: * 432203
> [M::ha_hist_line]    92: * 425832
> [M::ha_hist_line]    93: * 416998
> [M::ha_hist_line]    94: * 406989
> [M::ha_hist_line]    95: * 394719
> [M::ha_hist_line]    96: * 385162
> [M::ha_hist_line]    97: * 373520
> [M::ha_hist_line]    98: * 365359
> [M::ha_hist_line]    99: * 358934
> [M::ha_hist_line]   100: * 352159
> [M::ha_hist_line]   101: * 343358
> [M::ha_hist_line]   102: * 337238
> [M::ha_hist_line]   103: * 325839
> [M::ha_hist_line]   104: * 320692
> [M::ha_hist_line]   105: * 315515
> [M::ha_hist_line]   106: * 303257
> [M::ha_hist_line]   107: * 297011
> [M::ha_hist_line]   108: * 294470
> [M::ha_hist_line]  rest: ************************************************ 27967432
> [M::ha_analyze_count] left: none
> [M::ha_analyze_count] right: none
> [M::ha_ft_gen] peak_hom: 29; peak_het: -1
> [M::ha_ct_shrink::1663.885*6.29] ==> counted 20326276 distinct minimizer k-mers
> [M::ha_ft_gen::1670.724*6.26@32.967GB] ==> filtered out 20326276 k-mers occurring 145 or more times
> [M::ha_opt_update_cov] updated max_n_chain to 145
> [M::yak_count] collected 2150418983 minimizers
> [M::ha_pt_gen::2995.337*5.88] ==> counted 120121427 distinct minimizer k-mers
> [M::ha_pt_gen] count[4095] = 0 (for sanity check)
> [M::ha_analyze_count] lowest: count[5] = 699590
> [M::ha_analyze_count] highest: count[28] = 2669743
> [M::ha_hist_line]     1: ****************************************************************************************************> 44092346
> [M::ha_hist_line]     2: ***************************************************************************************** 2372838
> [M::ha_hist_line]     3: *********************************** 947720
> [M::ha_hist_line]     4: *************************** 724263
> [M::ha_hist_line]     5: ************************** 699590
> [M::ha_hist_line]     6: ***************************** 763913
> [M::ha_hist_line]     7: ******************************** 850664
> [M::ha_hist_line]     8: ************************************ 962740
> [M::ha_hist_line]     9: **************************************** 1079764
> [M::ha_hist_line]    10: ********************************************* 1201636
> [M::ha_hist_line]    11: ************************************************* 1314451
> [M::ha_hist_line]    12: ***************************************************** 1420466
> [M::ha_hist_line]    13: ******************************************************** 1505708
> [M::ha_hist_line]    14: *********************************************************** 1580926
> [M::ha_hist_line]    15: ************************************************************* 1628444
> [M::ha_hist_line]    16: *************************************************************** 1669090
> [M::ha_hist_line]    17: **************************************************************** 1702246
> [M::ha_hist_line]    18: ****************************************************************** 1750713
> [M::ha_hist_line]    19: ******************************************************************** 1808358
> [M::ha_hist_line]    20: ********************************************************************** 1881949
> [M::ha_hist_line]    21: ************************************************************************** 1968072
> [M::ha_hist_line]    22: ****************************************************************************** 2075345
> [M::ha_hist_line]    23: ********************************************************************************** 2198868
> [M::ha_hist_line]    24: *************************************************************************************** 2323091
> [M::ha_hist_line]    25: ******************************************************************************************* 2434162
> [M::ha_hist_line]    26: *********************************************************************************************** 2546474
> [M::ha_hist_line]    27: ************************************************************************************************** 2622430
> [M::ha_hist_line]    28: **************************************************************************************************** 2669743
> [M::ha_hist_line]    29: **************************************************************************************************** 2667515
> [M::ha_hist_line]    30: *************************************************************************************************** 2636346
> [M::ha_hist_line]    31: *********************************************************************************************** 2543994
> [M::ha_hist_line]    32: ****************************************************************************************** 2414606
> [M::ha_hist_line]    33: ************************************************************************************ 2250991
> [M::ha_hist_line]    34: ***************************************************************************** 2058813
> [M::ha_hist_line]    35: ********************************************************************* 1850756
> [M::ha_hist_line]    36: ************************************************************* 1636319
> [M::ha_hist_line]    37: ***************************************************** 1416480
> [M::ha_hist_line]    38: ********************************************* 1214385
> [M::ha_hist_line]    39: ************************************** 1021510
> [M::ha_hist_line]    40: ******************************** 856937
> [M::ha_hist_line]    41: *************************** 708149
> [M::ha_hist_line]    42: ********************** 582262
> [M::ha_hist_line]    43: ****************** 474447
> [M::ha_hist_line]    44: *************** 397199
> [M::ha_hist_line]    45: ************ 332480
> [M::ha_hist_line]    46: *********** 281659
> [M::ha_hist_line]    47: ********* 244014
> [M::ha_hist_line]    48: ******** 217091
> [M::ha_hist_line]    49: ******* 194942
> [M::ha_hist_line]    50: ******* 181860
> [M::ha_hist_line]    51: ****** 169901
> [M::ha_hist_line]    52: ****** 161931
> [M::ha_hist_line]    53: ****** 154685
> [M::ha_hist_line]    54: ****** 151556
> [M::ha_hist_line]    55: ***** 146448
> [M::ha_hist_line]    56: ***** 143574
> [M::ha_hist_line]    57: ***** 139804
> [M::ha_hist_line]    58: ***** 135932
> [M::ha_hist_line]    59: ***** 131313
> [M::ha_hist_line]    60: ***** 127854
> [M::ha_hist_line]    61: ***** 123441
> [M::ha_hist_line]    62: **** 119948
> [M::ha_hist_line]    63: **** 115618
> [M::ha_hist_line]    64: **** 112469
> [M::ha_hist_line]    65: **** 106949
> [M::ha_hist_line]    66: **** 102565
> [M::ha_hist_line]    67: **** 97244
> [M::ha_hist_line]    68: **** 94257
> [M::ha_hist_line]    69: *** 89733
> [M::ha_hist_line]    70: *** 86094
> [M::ha_hist_line]    71: *** 81774
> [M::ha_hist_line]    72: *** 79531
> [M::ha_hist_line]    73: *** 75193
> [M::ha_hist_line]    74: *** 72208
> [M::ha_hist_line]    75: *** 68996
> [M::ha_hist_line]    76: ** 66133
> [M::ha_hist_line]    77: ** 64113
> [M::ha_hist_line]    78: ** 62224
> [M::ha_hist_line]    79: ** 58881
> [M::ha_hist_line]    80: ** 58149
> [M::ha_hist_line]    81: ** 55945
> [M::ha_hist_line]    82: ** 54409
> [M::ha_hist_line]    83: ** 53588
> [M::ha_hist_line]    84: ** 51900
> [M::ha_hist_line]    85: ** 50104
> [M::ha_hist_line]    86: ** 48970
> [M::ha_hist_line]    87: ** 47214
> [M::ha_hist_line]    88: ** 46050
> [M::ha_hist_line]    89: ** 45296
> [M::ha_hist_line]    90: ** 44596
> [M::ha_hist_line]    91: ** 43491
> [M::ha_hist_line]    92: ** 42213
> [M::ha_hist_line]    93: ** 41311
> [M::ha_hist_line]    94: ** 40277
> [M::ha_hist_line]    95: * 39125
> [M::ha_hist_line]    96: * 38532
> [M::ha_hist_line]    97: * 37495
> [M::ha_hist_line]    98: * 36206
> [M::ha_hist_line]    99: * 36011
> [M::ha_hist_line]   100: * 34712
> [M::ha_hist_line]   101: * 33939
> [M::ha_hist_line]   102: * 33383
> [M::ha_hist_line]   103: * 32552
> [M::ha_hist_line]   104: * 31495
> [M::ha_hist_line]   105: * 30816
> [M::ha_hist_line]   106: * 30166
> [M::ha_hist_line]   107: * 29423
> [M::ha_hist_line]   108: * 29273
> [M::ha_hist_line]   109: * 28500
> [M::ha_hist_line]   110: * 27951
> [M::ha_hist_line]   111: * 26712
> [M::ha_hist_line]   112: * 26227
> [M::ha_hist_line]   113: * 25598
> [M::ha_hist_line]   114: * 25443
> [M::ha_hist_line]   115: * 24614
> [M::ha_hist_line]   116: * 23690
> [M::ha_hist_line]   117: * 23774
> [M::ha_hist_line]   118: * 23079
> [M::ha_hist_line]   119: * 22745
> [M::ha_hist_line]   120: * 22174
> [M::ha_hist_line]   121: * 21640
> [M::ha_hist_line]   122: * 21386
> [M::ha_hist_line]   123: * 20758
> [M::ha_hist_line]   124: * 20635
> [M::ha_hist_line]   125: * 20016
> [M::ha_hist_line]   126: * 19658
> [M::ha_hist_line]   127: * 19031
> [M::ha_hist_line]   128: * 18651
> [M::ha_hist_line]   129: * 18247
> [M::ha_hist_line]   130: * 17888
> [M::ha_hist_line]   131: * 17761
> [M::ha_hist_line]   132: * 17304
> [M::ha_hist_line]   133: * 16782
> [M::ha_hist_line]   134: * 16440
> [M::ha_hist_line]   135: * 16011
> [M::ha_hist_line]   136: * 15551
> [M::ha_hist_line]   137: * 15295
> [M::ha_hist_line]   138: * 15192
> [M::ha_hist_line]   139: * 14909
> [M::ha_hist_line]   140: * 14648
> [M::ha_hist_line]   141: * 13844
> [M::ha_hist_line]  rest: ***** 133498
> [M::ha_analyze_count] left: none
> [M::ha_analyze_count] right: none
> [M::ha_pt_gen] peak_hom: 28; peak_het: -1
> [M::ha_ct_shrink::2995.668*5.88] ==> counted 76029081 distinct minimizer k-mers
> [M::ha_pt_gen::] counting in normal mode
> [M::yak_count] collected 2150418983 minimizers
> [M::ha_pt_gen::3507.182*6.77] ==> indexed 2106326637 positions, counted 76029081 distinct minimizer k-mers
> [M::ha_assemble::16440.053*13.20@41.637GB] ==> corrected reads for round 1
> [M::ha_assemble] # bases: 78526930206; # corrected bases: 164481046; # recorrected bases: 209675
> [M::ha_assemble] size of buffer: 1.046GB
> [M::yak_count] collected 2136684417 minimizers
> [M::ha_pt_gen::16862.096*13.22] ==> counted 77519250 distinct minimizer k-mers
> [M::ha_pt_gen] count[4095] = 0 (for sanity check)
> [M::ha_analyze_count] lowest: count[5] = 571346
> [M::ha_analyze_count] highest: count[30] = 2646196
> [M::ha_hist_line]     1: ****************************************************************************************************> 4230298
> [M::ha_hist_line]     2: ************************ 626066
> [M::ha_hist_line]     3: ****************** 474525
> [M::ha_hist_line]     4: ******************* 509503
> [M::ha_hist_line]     5: ********************** 571346
> [M::ha_hist_line]     6: ************************* 670384
> [M::ha_hist_line]     7: ***************************** 774030
> [M::ha_hist_line]     8: ********************************** 892466
> [M::ha_hist_line]     9: ************************************** 1014431
> [M::ha_hist_line]    10: ******************************************* 1144261
> [M::ha_hist_line]    11: *********************************************** 1256844
> [M::ha_hist_line]    12: **************************************************** 1367388
> [M::ha_hist_line]    13: ******************************************************* 1451159
> [M::ha_hist_line]    14: ********************************************************** 1541607
> [M::ha_hist_line]    15: ************************************************************ 1583148
> [M::ha_hist_line]    16: ************************************************************** 1629208
> [M::ha_hist_line]    17: *************************************************************** 1654801
> [M::ha_hist_line]    18: **************************************************************** 1703465
> [M::ha_hist_line]    19: ****************************************************************** 1749465
> [M::ha_hist_line]    20: ********************************************************************* 1820264
> [M::ha_hist_line]    21: ************************************************************************ 1897111
> [M::ha_hist_line]    22: **************************************************************************** 2000591
> [M::ha_hist_line]    23: ******************************************************************************** 2112589
> [M::ha_hist_line]    24: ************************************************************************************* 2244267
> [M::ha_hist_line]    25: ***************************************************************************************** 2353974
> [M::ha_hist_line]    26: ********************************************************************************************* 2472891
> [M::ha_hist_line]    27: ************************************************************************************************* 2566430
> [M::ha_hist_line]    28: *************************************************************************************************** 2632270
> [M::ha_hist_line]    29: **************************************************************************************************** 2640793
> [M::ha_hist_line]    30: **************************************************************************************************** 2646196
> [M::ha_hist_line]    31: ************************************************************************************************* 2573269
> [M::ha_hist_line]    32: ********************************************************************************************* 2463738
> [M::ha_hist_line]    33: **************************************************************************************** 2317130
> [M::ha_hist_line]    34: ********************************************************************************* 2144355
> [M::ha_hist_line]    35: ************************************************************************** 1947321
> [M::ha_hist_line]    36: ****************************************************************** 1740567
> [M::ha_hist_line]    37: ********************************************************* 1516001
> [M::ha_hist_line]    38: ************************************************** 1313483
> [M::ha_hist_line]    39: ****************************************** 1116506
> [M::ha_hist_line]    40: *********************************** 938578
> [M::ha_hist_line]    41: ***************************** 780041
> [M::ha_hist_line]    42: ************************* 648919
> [M::ha_hist_line]    43: ******************** 522695
> [M::ha_hist_line]    44: ***************** 438976
> [M::ha_hist_line]    45: ************** 362583
> [M::ha_hist_line]    46: ************ 307816
> [M::ha_hist_line]    47: ********** 259091
> [M::ha_hist_line]    48: ********* 231969
> [M::ha_hist_line]    49: ******** 202960
> [M::ha_hist_line]    50: ******* 187926
> [M::ha_hist_line]    51: ******* 172521
> [M::ha_hist_line]    52: ****** 167129
> [M::ha_hist_line]    53: ****** 154981
> [M::ha_hist_line]    54: ****** 154227
> [M::ha_hist_line]    55: ***** 144950
> [M::ha_hist_line]    56: ***** 145281
> [M::ha_hist_line]    57: ***** 139987
> [M::ha_hist_line]    58: ***** 137710
> [M::ha_hist_line]    59: ***** 131797
> [M::ha_hist_line]    60: ***** 131304
> [M::ha_hist_line]    61: ***** 123125
> [M::ha_hist_line]    62: ***** 121546
> [M::ha_hist_line]    63: **** 118159
> [M::ha_hist_line]    64: **** 115301
> [M::ha_hist_line]    65: **** 108427
> [M::ha_hist_line]    66: **** 107907
> [M::ha_hist_line]    67: **** 99271
> [M::ha_hist_line]    68: **** 96986
> [M::ha_hist_line]    69: *** 91546
> [M::ha_hist_line]    70: *** 88914
> [M::ha_hist_line]    71: *** 83177
> [M::ha_hist_line]    72: *** 83352
> [M::ha_hist_line]    73: *** 76886
> [M::ha_hist_line]    74: *** 74822
> [M::ha_hist_line]    75: *** 70636
> [M::ha_hist_line]    76: *** 68393
> [M::ha_hist_line]    77: ** 64750
> [M::ha_hist_line]    78: ** 64672
> [M::ha_hist_line]    79: ** 60045
> [M::ha_hist_line]    80: ** 59396
> [M::ha_hist_line]    81: ** 57197
> [M::ha_hist_line]    82: ** 55180
> [M::ha_hist_line]    83: ** 53168
> [M::ha_hist_line]    84: ** 53563
> [M::ha_hist_line]    85: ** 51361
> [M::ha_hist_line]    86: ** 49601
> [M::ha_hist_line]    87: ** 48374
> [M::ha_hist_line]    88: ** 46915
> [M::ha_hist_line]    89: ** 45245
> [M::ha_hist_line]    90: ** 45769
> [M::ha_hist_line]    91: ** 43392
> [M::ha_hist_line]    92: ** 42890
> [M::ha_hist_line]    93: ** 41660
> [M::ha_hist_line]    94: ** 40874
> [M::ha_hist_line]    95: ** 39881
> [M::ha_hist_line]    96: ** 40275
> [M::ha_hist_line]    97: * 37623
> [M::ha_hist_line]    98: * 36846
> [M::ha_hist_line]    99: * 36220
> [M::ha_hist_line]   100: * 35907
> [M::ha_hist_line]   101: * 34616
> [M::ha_hist_line]   102: * 34537
> [M::ha_hist_line]   103: * 32648
> [M::ha_hist_line]   104: * 32477
> [M::ha_hist_line]   105: * 31747
> [M::ha_hist_line]   106: * 31048
> [M::ha_hist_line]   107: * 29362
> [M::ha_hist_line]   108: * 30265
> [M::ha_hist_line]   109: * 28735
> [M::ha_hist_line]   110: * 28250
> [M::ha_hist_line]   111: * 27356
> [M::ha_hist_line]   112: * 26959
> [M::ha_hist_line]   113: * 26267
> [M::ha_hist_line]   114: * 25858
> [M::ha_hist_line]   115: * 25234
> [M::ha_hist_line]   116: * 24820
> [M::ha_hist_line]   117: * 24038
> [M::ha_hist_line]   118: * 23155
> [M::ha_hist_line]   119: * 23052
> [M::ha_hist_line]   120: * 23323
> [M::ha_hist_line]   121: * 22111
> [M::ha_hist_line]   122: * 20899
> [M::ha_hist_line]   123: * 20793
> [M::ha_hist_line]   124: * 20981
> [M::ha_hist_line]   125: * 20514
> [M::ha_hist_line]   126: * 20169
> [M::ha_hist_line]   127: * 19107
> [M::ha_hist_line]   128: * 19270
> [M::ha_hist_line]   129: * 18598
> [M::ha_hist_line]   130: * 18270
> [M::ha_hist_line]   131: * 17941
> [M::ha_hist_line]   132: * 18132
> [M::ha_hist_line]   133: * 17465
> [M::ha_hist_line]   134: * 16651
> [M::ha_hist_line]   135: * 16172
> [M::ha_hist_line]   136: * 16342
> [M::ha_hist_line]   137: * 15252
> [M::ha_hist_line]   138: * 15321
> [M::ha_hist_line]   139: * 14647
> [M::ha_hist_line]   140: * 14930
> [M::ha_hist_line]   141: * 14081
> [M::ha_hist_line]   142: * 14002
> [M::ha_hist_line]   143: * 13365
> [M::ha_hist_line]  rest: ****** 145656
> [M::ha_analyze_count] left: none
> [M::ha_analyze_count] right: none
> [M::ha_pt_gen] peak_hom: 30; peak_het: -1
> [M::ha_ct_shrink::16862.287*13.22] ==> counted 73288952 distinct minimizer k-mers
> [M::ha_pt_gen::] counting in normal mode
> [M::yak_count] collected 2136684417 minimizers
> [M::ha_pt_gen::17336.859*13.21] ==> indexed 2132454119 positions, counted 73288952 distinct minimizer k-mers
> [M::ha_assemble::28986.522*13.92@51.053GB] ==> corrected reads for round 2
> [M::ha_assemble] # bases: 78518923696; # corrected bases: 9193685; # recorrected bases: 345225
> [M::ha_assemble] size of buffer: 0.937GB
> [M::yak_count] collected 2136213806 minimizers
> [M::ha_pt_gen::29416.997*13.91] ==> counted 75997726 distinct minimizer k-mers
> [M::ha_pt_gen] count[4095] = 0 (for sanity check)
> [M::ha_analyze_count] lowest: count[5] = 561030
> [M::ha_analyze_count] highest: count[30] = 2646503
> [M::ha_hist_line]     1: ****************************************************************************************************> 2857987
> [M::ha_hist_line]     2: ******************** 541877
> [M::ha_hist_line]     3: ***************** 443957
> [M::ha_hist_line]     4: ******************* 495480
> [M::ha_hist_line]     5: ********************* 561030
> [M::ha_hist_line]     6: ************************* 663441
> [M::ha_hist_line]     7: ***************************** 769194
> [M::ha_hist_line]     8: ********************************** 888047
> [M::ha_hist_line]     9: ************************************** 1011549
> [M::ha_hist_line]    10: ******************************************* 1141019
> [M::ha_hist_line]    11: *********************************************** 1254855
> [M::ha_hist_line]    12: **************************************************** 1366210
> [M::ha_hist_line]    13: ******************************************************* 1450001
> [M::ha_hist_line]    14: ********************************************************** 1539983
> [M::ha_hist_line]    15: ************************************************************ 1581573
> [M::ha_hist_line]    16: ************************************************************** 1628711
> [M::ha_hist_line]    17: *************************************************************** 1654551
> [M::ha_hist_line]    18: **************************************************************** 1702464
> [M::ha_hist_line]    19: ****************************************************************** 1747999
> [M::ha_hist_line]    20: ********************************************************************* 1819369
> [M::ha_hist_line]    21: ************************************************************************ 1896413
> [M::ha_hist_line]    22: **************************************************************************** 2000222
> [M::ha_hist_line]    23: ******************************************************************************** 2112130
> [M::ha_hist_line]    24: ************************************************************************************* 2244015
> [M::ha_hist_line]    25: ***************************************************************************************** 2352145
> [M::ha_hist_line]    26: ********************************************************************************************* 2472580
> [M::ha_hist_line]    27: ************************************************************************************************* 2566246
> [M::ha_hist_line]    28: *************************************************************************************************** 2632687
> [M::ha_hist_line]    29: **************************************************************************************************** 2640475
> [M::ha_hist_line]    30: **************************************************************************************************** 2646503
> [M::ha_hist_line]    31: ************************************************************************************************* 2574777
> [M::ha_hist_line]    32: ********************************************************************************************* 2464287
> [M::ha_hist_line]    33: **************************************************************************************** 2317808
> [M::ha_hist_line]    34: ********************************************************************************* 2146876
> [M::ha_hist_line]    35: ************************************************************************** 1949072
> [M::ha_hist_line]    36: ****************************************************************** 1742391
> [M::ha_hist_line]    37: ********************************************************* 1518121
> [M::ha_hist_line]    38: ************************************************** 1314949
> [M::ha_hist_line]    39: ****************************************** 1118248
> [M::ha_hist_line]    40: ************************************ 939839
> [M::ha_hist_line]    41: ****************************** 781371
> [M::ha_hist_line]    42: ************************* 649577
> [M::ha_hist_line]    43: ******************** 522994
> [M::ha_hist_line]    44: ***************** 439662
> [M::ha_hist_line]    45: ************** 363000
> [M::ha_hist_line]    46: ************ 308502
> [M::ha_hist_line]    47: ********** 259523
> [M::ha_hist_line]    48: ********* 232401
> [M::ha_hist_line]    49: ******** 202562
> [M::ha_hist_line]    50: ******* 188584
> [M::ha_hist_line]    51: ******* 172744
> [M::ha_hist_line]    52: ****** 167045
> [M::ha_hist_line]    53: ****** 154730
> [M::ha_hist_line]    54: ****** 154422
> [M::ha_hist_line]    55: ***** 144973
> [M::ha_hist_line]    56: ****** 145678
> [M::ha_hist_line]    57: ***** 139653
> [M::ha_hist_line]    58: ***** 137657
> [M::ha_hist_line]    59: ***** 131624
> [M::ha_hist_line]    60: ***** 131824
> [M::ha_hist_line]    61: ***** 123245
> [M::ha_hist_line]    62: ***** 121655
> [M::ha_hist_line]    63: **** 118005
> [M::ha_hist_line]    64: **** 115633
> [M::ha_hist_line]    65: **** 108503
> [M::ha_hist_line]    66: **** 107841
> [M::ha_hist_line]    67: **** 99230
> [M::ha_hist_line]    68: **** 97238
> [M::ha_hist_line]    69: *** 91582
> [M::ha_hist_line]    70: *** 89312
> [M::ha_hist_line]    71: *** 83485
> [M::ha_hist_line]    72: *** 83448
> [M::ha_hist_line]    73: *** 76663
> [M::ha_hist_line]    74: *** 75159
> [M::ha_hist_line]    75: *** 70681
> [M::ha_hist_line]    76: *** 68205
> [M::ha_hist_line]    77: ** 64700
> [M::ha_hist_line]    78: ** 64477
> [M::ha_hist_line]    79: ** 60303
> [M::ha_hist_line]    80: ** 59637
> [M::ha_hist_line]    81: ** 57264
> [M::ha_hist_line]    82: ** 55086
> [M::ha_hist_line]    83: ** 53314
> [M::ha_hist_line]    84: ** 53439
> [M::ha_hist_line]    85: ** 51136
> [M::ha_hist_line]    86: ** 49977
> [M::ha_hist_line]    87: ** 48594
> [M::ha_hist_line]    88: ** 47044
> [M::ha_hist_line]    89: ** 45072
> [M::ha_hist_line]    90: ** 45730
> [M::ha_hist_line]    91: ** 43259
> [M::ha_hist_line]    92: ** 43064
> [M::ha_hist_line]    93: ** 41833
> [M::ha_hist_line]    94: ** 40866
> [M::ha_hist_line]    95: ** 39939
> [M::ha_hist_line]    96: ** 40355
> [M::ha_hist_line]    97: * 37536
> [M::ha_hist_line]    98: * 36676
> [M::ha_hist_line]    99: * 36145
> [M::ha_hist_line]   100: * 35890
> [M::ha_hist_line]   101: * 34728
> [M::ha_hist_line]   102: * 34457
> [M::ha_hist_line]   103: * 32761
> [M::ha_hist_line]   104: * 32568
> [M::ha_hist_line]   105: * 31827
> [M::ha_hist_line]   106: * 30941
> [M::ha_hist_line]   107: * 29539
> [M::ha_hist_line]   108: * 30249
> [M::ha_hist_line]   109: * 28631
> [M::ha_hist_line]   110: * 28099
> [M::ha_hist_line]   111: * 27098
> [M::ha_hist_line]   112: * 27121
> [M::ha_hist_line]   113: * 26235
> [M::ha_hist_line]   114: * 26023
> [M::ha_hist_line]   115: * 25231
> [M::ha_hist_line]   116: * 24741
> [M::ha_hist_line]   117: * 24163
> [M::ha_hist_line]   118: * 23139
> [M::ha_hist_line]   119: * 23091
> [M::ha_hist_line]   120: * 23546
> [M::ha_hist_line]   121: * 22070
> [M::ha_hist_line]   122: * 21113
> [M::ha_hist_line]   123: * 20745
> [M::ha_hist_line]   124: * 20979
> [M::ha_hist_line]   125: * 20382
> [M::ha_hist_line]   126: * 20286
> [M::ha_hist_line]   127: * 19128
> [M::ha_hist_line]   128: * 19470
> [M::ha_hist_line]   129: * 18473
> [M::ha_hist_line]   130: * 18208
> [M::ha_hist_line]   131: * 17983
> [M::ha_hist_line]   132: * 18059
> [M::ha_hist_line]   133: * 17438
> [M::ha_hist_line]   134: * 16719
> [M::ha_hist_line]   135: * 16372
> [M::ha_hist_line]   136: * 16256
> [M::ha_hist_line]   137: * 15425
> [M::ha_hist_line]   138: * 15333
> [M::ha_hist_line]   139: * 14777
> [M::ha_hist_line]   140: * 14834
> [M::ha_hist_line]   141: * 14142
> [M::ha_hist_line]   142: * 13893
> [M::ha_hist_line]   143: * 13449
> [M::ha_hist_line]  rest: ****** 149176
> [M::ha_analyze_count] left: none
> [M::ha_analyze_count] right: none
> [M::ha_pt_gen] peak_hom: 30; peak_het: -1
> [M::ha_ct_shrink::29417.178*13.91] ==> counted 73139739 distinct minimizer k-mers
> [M::ha_pt_gen::] counting in normal mode
> [M::yak_count] collected 2136213806 minimizers
> [M::ha_pt_gen::29864.608*13.91] ==> indexed 2133355819 positions, counted 73139739 distinct minimizer k-mers
> [M::ha_assemble::41441.128*14.20@70.208GB] ==> corrected reads for round 3
> [M::ha_assemble] # bases: 78518213515; # corrected bases: 3081352; # recorrected bases: 261603
> [M::ha_assemble] size of buffer: 0.913GB
> [M::yak_count] collected 2136097998 minimizers
> [M::ha_pt_gen::41854.126*14.20] ==> counted 75694455 distinct minimizer k-mers
> [M::ha_pt_gen] count[4095] = 0 (for sanity check)
> [M::ha_analyze_count] lowest: count[5] = 557291
> [M::ha_analyze_count] highest: count[30] = 2646697
> [M::ha_hist_line]     1: ************************************************************************************************** 2601316
> [M::ha_hist_line]     2: ******************** 518709
> [M::ha_hist_line]     3: **************** 433798
> [M::ha_hist_line]     4: ******************* 491216
> [M::ha_hist_line]     5: ********************* 557291
> [M::ha_hist_line]     6: ************************* 661019
> [M::ha_hist_line]     7: ***************************** 767948
> [M::ha_hist_line]     8: ********************************** 887247
> [M::ha_hist_line]     9: ************************************** 1010664
> [M::ha_hist_line]    10: ******************************************* 1140414
> [M::ha_hist_line]    11: *********************************************** 1253583
> [M::ha_hist_line]    12: **************************************************** 1365947
> [M::ha_hist_line]    13: ******************************************************* 1449267
> [M::ha_hist_line]    14: ********************************************************** 1540107
> [M::ha_hist_line]    15: ************************************************************ 1581119
> [M::ha_hist_line]    16: ************************************************************** 1628696
> [M::ha_hist_line]    17: *************************************************************** 1654341
> [M::ha_hist_line]    18: **************************************************************** 1702347
> [M::ha_hist_line]    19: ****************************************************************** 1747856
> [M::ha_hist_line]    20: ********************************************************************* 1819527
> [M::ha_hist_line]    21: ************************************************************************ 1896109
> [M::ha_hist_line]    22: **************************************************************************** 2000023
> [M::ha_hist_line]    23: ******************************************************************************** 2111607
> [M::ha_hist_line]    24: ************************************************************************************* 2243833
> [M::ha_hist_line]    25: ***************************************************************************************** 2352191
> [M::ha_hist_line]    26: ********************************************************************************************* 2472759
> [M::ha_hist_line]    27: ************************************************************************************************* 2566252
> [M::ha_hist_line]    28: *************************************************************************************************** 2632924
> [M::ha_hist_line]    29: **************************************************************************************************** 2640463
> [M::ha_hist_line]    30: **************************************************************************************************** 2646697
> [M::ha_hist_line]    31: ************************************************************************************************* 2575026
> [M::ha_hist_line]    32: ********************************************************************************************* 2464759
> [M::ha_hist_line]    33: **************************************************************************************** 2318129
> [M::ha_hist_line]    34: ********************************************************************************* 2147273
> [M::ha_hist_line]    35: ************************************************************************** 1949181
> [M::ha_hist_line]    36: ****************************************************************** 1742704
> [M::ha_hist_line]    37: ********************************************************* 1518094
> [M::ha_hist_line]    38: ************************************************** 1314973
> [M::ha_hist_line]    39: ****************************************** 1118274
> [M::ha_hist_line]    40: ************************************ 940152
> [M::ha_hist_line]    41: ****************************** 781336
> [M::ha_hist_line]    42: ************************* 649681
> [M::ha_hist_line]    43: ******************** 523312
> [M::ha_hist_line]    44: ***************** 439662
> [M::ha_hist_line]    45: ************** 362884
> [M::ha_hist_line]    46: ************ 308794
> [M::ha_hist_line]    47: ********** 259357
> [M::ha_hist_line]    48: ********* 232426
> [M::ha_hist_line]    49: ******** 202533
> [M::ha_hist_line]    50: ******* 188718
> [M::ha_hist_line]    51: ******* 172434
> [M::ha_hist_line]    52: ****** 167028
> [M::ha_hist_line]    53: ****** 154614
> [M::ha_hist_line]    54: ****** 154496
> [M::ha_hist_line]    55: ***** 145104
> [M::ha_hist_line]    56: ****** 145802
> [M::ha_hist_line]    57: ***** 139634
> [M::ha_hist_line]    58: ***** 137666
> [M::ha_hist_line]    59: ***** 131711
> [M::ha_hist_line]    60: ***** 131613
> [M::ha_hist_line]    61: ***** 123123
> [M::ha_hist_line]    62: ***** 121884
> [M::ha_hist_line]    63: **** 117920
> [M::ha_hist_line]    64: **** 115589
> [M::ha_hist_line]    65: **** 108544
> [M::ha_hist_line]    66: **** 107820
> [M::ha_hist_line]    67: **** 99325
> [M::ha_hist_line]    68: **** 97266
> [M::ha_hist_line]    69: *** 91748
> [M::ha_hist_line]    70: *** 89159
> [M::ha_hist_line]    71: *** 83429
> [M::ha_hist_line]    72: *** 83501
> [M::ha_hist_line]    73: *** 76711
> [M::ha_hist_line]    74: *** 75092
> [M::ha_hist_line]    75: *** 70843
> [M::ha_hist_line]    76: *** 68310
> [M::ha_hist_line]    77: ** 64568
> [M::ha_hist_line]    78: ** 64563
> [M::ha_hist_line]    79: ** 60123
> [M::ha_hist_line]    80: ** 59723
> [M::ha_hist_line]    81: ** 57229
> [M::ha_hist_line]    82: ** 55041
> [M::ha_hist_line]    83: ** 53344
> [M::ha_hist_line]    84: ** 53534
> [M::ha_hist_line]    85: ** 51123
> [M::ha_hist_line]    86: ** 50007
> [M::ha_hist_line]    87: ** 48603
> [M::ha_hist_line]    88: ** 47072
> [M::ha_hist_line]    89: ** 45060
> [M::ha_hist_line]    90: ** 45701
> [M::ha_hist_line]    91: ** 43232
> [M::ha_hist_line]    92: ** 42935
> [M::ha_hist_line]    93: ** 41906
> [M::ha_hist_line]    94: ** 40971
> [M::ha_hist_line]    95: ** 40052
> [M::ha_hist_line]    96: ** 40364
> [M::ha_hist_line]    97: * 37496
> [M::ha_hist_line]    98: * 36714
> [M::ha_hist_line]    99: * 36121
> [M::ha_hist_line]   100: * 35861
> [M::ha_hist_line]   101: * 34724
> [M::ha_hist_line]   102: * 34647
> [M::ha_hist_line]   103: * 32580
> [M::ha_hist_line]   104: * 32548
> [M::ha_hist_line]   105: * 31839
> [M::ha_hist_line]   106: * 30816
> [M::ha_hist_line]   107: * 29595
> [M::ha_hist_line]   108: * 30157
> [M::ha_hist_line]   109: * 28630
> [M::ha_hist_line]   110: * 28210
> [M::ha_hist_line]   111: * 27087
> [M::ha_hist_line]   112: * 27163
> [M::ha_hist_line]   113: * 26278
> [M::ha_hist_line]   114: * 26011
> [M::ha_hist_line]   115: * 25202
> [M::ha_hist_line]   116: * 24687
> [M::ha_hist_line]   117: * 24254
> [M::ha_hist_line]   118: * 23105
> [M::ha_hist_line]   119: * 23166
> [M::ha_hist_line]   120: * 23419
> [M::ha_hist_line]   121: * 22081
> [M::ha_hist_line]   122: * 21130
> [M::ha_hist_line]   123: * 20870
> [M::ha_hist_line]   124: * 20950
> [M::ha_hist_line]   125: * 20305
> [M::ha_hist_line]   126: * 20238
> [M::ha_hist_line]   127: * 19169
> [M::ha_hist_line]   128: * 19517
> [M::ha_hist_line]   129: * 18574
> [M::ha_hist_line]   130: * 18121
> [M::ha_hist_line]   131: * 18010
> [M::ha_hist_line]   132: * 18162
> [M::ha_hist_line]   133: * 17415
> [M::ha_hist_line]   134: * 16694
> [M::ha_hist_line]   135: * 16344
> [M::ha_hist_line]   136: * 16229
> [M::ha_hist_line]   137: * 15416
> [M::ha_hist_line]   138: * 15346
> [M::ha_hist_line]   139: * 14805
> [M::ha_hist_line]   140: * 14897
> [M::ha_hist_line]   141: * 14097
> [M::ha_hist_line]   142: * 13900
> [M::ha_hist_line]   143: * 13509
> [M::ha_hist_line]  rest: ****** 150311
> [M::ha_analyze_count] left: none
> [M::ha_analyze_count] right: none
> [M::ha_pt_gen] peak_hom: 30; peak_het: -1
> [M::ha_ct_shrink::41854.303*14.20] ==> counted 73093139 distinct minimizer k-mers
> [M::ha_pt_gen::] counting in normal mode
> [M::yak_count] collected 2136097998 minimizers
> [M::ha_pt_gen::42292.426*14.20] ==> indexed 2133496682 positions, counted 73093139 distinct minimizer k-mers
> [M::ha_assemble::46968.151*14.27@75.620GB] ==> found overlaps for the final round
> [M::ha_print_ovlp_stat] # overlaps: 203526838
> [M::ha_print_ovlp_stat] # strong overlaps: 107525203
> [M::ha_print_ovlp_stat] # weak overlaps: 96001635
> [M::ha_print_ovlp_stat] # exact overlaps: 184993604
> [M::ha_print_ovlp_stat] # inexact overlaps: 18533234
> [M::ha_print_ovlp_stat] # overlaps without large indels: 201290843
> [M::ha_print_ovlp_stat] # reverse overlaps: 102334556
> [M::ha_opt_update_cov_min] updated max_n_chain to 150
> Writing reads to disk... 
> Reads has been written.
> Writing ma_hit_ts to disk... 
> ma_hit_ts has been written.
> Writing ma_hit_ts to disk... 
> ma_hit_ts has been written.
> bin files have been written.
> [M::ul_load::] ==> UL
> [M::ha_opt_update_cov] updated max_n_chain to 150
> [M::append_inexact_edges] # inserted inexact edges: 3395670
> [M::gen_cov_track::] # bases: 22469526191
> [M::dedup_HiFis::] # unitigs: 1300562, # edges: 7393744, # cc_num: 3471239
> [M::ha_ct_shrink::47731.025*14.09] ==> counted 10555847 distinct minimizer k-mers
> [M::ha_ft_ul_gen::47733.674*14.09@85.027GB] ==> filtered out 10555847 k-mers occurring 150 or more times
> [M::yak_count] collected 1793730484 minimizers
> [M::ha_pt_ul_gen::47902.739*14.09] ==> counted 90139273 distinct minimizer k-mers
> [M::ha_ct_shrink::47902.890*14.09] ==> counted 90139273 distinct minimizer k-mers
> [M::yak_count] collected 1793730484 minimizers
> [M::ha_pt_ul_gen::48122.755*14.08] ==> indexed 1793730484 positions
> [M::uidx_l_build] Index has been built.
> [M::uidx_write] Index has been written.
> [M::scall_ul_pipeline::38579.504] ==> Qualification
> [M::scall_ul_pipeline::] ==> # reads: 1511834, # bases: 0
> [M::scall_ul_pipeline::] ==> # fully covered reads: 1100049
> [M::write_all_ul_t] Index has been written.
> [M::print_ovlp_src_bl_stat::] ==> # HiFi reads:5754000, # covered HiFi reads:2672433, # chained HiFi reads:2208051
> [M::print_ovlp_src_bl_stat::] ==> # passed UL reads:1511834, # fully corrected UL reads:1054995, # almost fully corrected UL reads:39671, # UL reads have primary chains:187737
> [M::ha_opt_update_cov] updated max_n_chain to 150
> [M::work_ul_gchains_consensus::] # try:332054, # done:200822
> [M::clean_contain_g0::] # discard cread::183218
> [M::update_sg_uo::] ==> # gfa reads:1875234, # covered gfa reads:1740266
> [M::ul_refine_alignment::] ==> UL refinement...
> [M::ha_opt_update_cov] updated max_n_chain to 150
> [M::gen_cov_track::] # bases: 8199315173
> [M::work_ul_gchains::] # try:299276, # done:221285
> [M::update_sg_uo::] ==> # gfa reads:1784779, # covered gfa reads:1730400
> [M::ha_opt_update_cov] updated max_n_chain to 150
> [M::work_ul_gchains_consensus::] # try:21108, # done:1892
> [M::clean_contain_g0::] # discard cread::581
> [M::update_sg_uo::] ==> # gfa reads:1784198, # covered gfa reads:1729812
> [M::ul_refine_alignment::] ==> UL refinement...
> [M::ha_opt_update_cov] updated max_n_chain to 150
> [M::gen_cov_track::] # bases: 7941449228
> [M::work_ul_gchains::] # try:48391, # done:13507
> [M::update_sg_uo::] ==> # gfa reads:1782318, # covered gfa reads:1729049
> [M::ha_opt_update_cov] updated max_n_chain to 150
> [M::work_ul_gchains_consensus::] # try:18385, # done:135
> [M::clean_contain_g0::] # discard cread::14
> [M::update_sg_uo::] ==> # gfa reads:1782304, # covered gfa reads:1729014
> [M::ul_refine_alignment::] ==> UL refinement...
> [M::ha_opt_update_cov] updated max_n_chain to 150
> [M::gen_cov_track::] # bases: 7816522731
> [M::work_ul_gchains::] # try:42453, # done:8940
> [M::update_sg_uo::] ==> # gfa reads:1781269, # covered gfa reads:1729174
> [M::ha_opt_update_cov] updated max_n_chain to 150
> [M::work_ul_gchains_consensus::] # try:18205, # done:91
> [M::clean_contain_g0::] # discard cread::2
> [M::update_sg_uo::] ==> # gfa reads:1781267, # covered gfa reads:1729161
> [M::ul_refine_alignment::] ==> UL refinement...
> [M::ha_opt_update_cov] updated max_n_chain to 150
> [M::gen_cov_track::] # bases: 7735445481
> [M::work_ul_gchains::] # try:41704, # done:8649
> [M::update_sg_uo::] ==> # gfa reads:1780212, # covered gfa reads:1729241
> [M::ha_opt_update_cov] updated max_n_chain to 150
> [M::work_ul_gchains_consensus::] # try:18015, # done:83
> [M::clean_contain_g0::] # discard cread::2
> [M::update_sg_uo::] ==> # gfa reads:1780210, # covered gfa reads:1729228
> [M::ha_opt_update_cov] updated max_n_chain to 150
> [M::work_ul_gchains_consensus::] # try:17856, # done:0
> [M::clean_contain_g0::] # discard cread::2
> [M::update_sg_uo::] ==> # gfa reads:1780187, # covered gfa reads:1729205
> [M::purge_dups] homozygous read coverage threshold: 29
> [M::purge_dups] purge duplication coverage threshold: 37
> [M::ug_ext_gfa::] # tips::2078
> Writing raw unitig GFA to disk... 
> [M::update_sg_uo::] ==> # gfa reads:1770259, # covered gfa reads:1719925
> [M::ul_realignment::] ==> starting UL
> [M::ha_opt_update_cov] updated max_n_chain to 150
> [M::gen_cov_track::] # bases: 6889354816
> [M::ha_ct_shrink::87429.171*14.32] ==> counted 18918250 distinct minimizer k-mers
> [M::ha_ft_ul_gen::87432.059*14.32@98.160GB] ==> filtered out 18918250 k-mers occurring 32 or more times
> [M::yak_count] collected 528917359 minimizers
> [M::ha_pt_ul_gen::87499.187*14.32] ==> counted 91758444 distinct minimizer k-mers
> [M::ha_ct_shrink::87499.344*14.32] ==> counted 91758444 distinct minimizer k-mers
> [M::yak_count] collected 528917359 minimizers
> [M::ha_pt_ul_gen::87573.019*14.32] ==> indexed 528917359 positions
> [M::uidx_l_build] Index has been built.
> [M::uidx_write] Index has been written.
> [M::ha_opt_update_cov] updated max_n_chain to 150
> [M::cal_graph_ovlp_binning::2890.776] ==> Qualification
> [M::write_emask_t] Index has been written.
> [M::rescall_ul_pipeline::16152.305] ==> Qualification
> [M::rescall_ul_pipeline::] ==> # reads: 1511834, # bases: 0, # fully corrected reads: 1378027
> [M::write_all_ul_t] Index has been written.
> hifiasm: gfa_ut.cpp:4917: void uc_block_qse_cutoff(uint32_t, uint32_t, uc_block_t*, uint32_t*, uint32_t*): Assertion `(*nqs) >= x->qs && (*nqe) <= x->qe' failed.
chhylp123 commented 1 year ago

Ah, it looks like a bug... Is it possible that you can share the bin files with us? Thanks in advance.

olekto commented 1 year ago

Thank you Haoyu. I have sent you an e-mail with links to the bin files.

Ole

olekto commented 1 year ago

The assembly successfully completed with 0.19.5. Thank you!