// File: /workspace/github_cpp_makefile_star2_O2/chhylp123[P]hifiasm/Process_Read.cpp
// Function: append_ul_t_back(all_ul_t*, unsigned long*, char*, long, char*, long, ul_ov_t*, long, float)
void append_ul_t_back(all_ul_t *x, uint64_t *rid, char* id, int64_t id_l, char* str, int64_t str_l, ul_ov_t *o, int64_t on, float p_chain_rate) {
	int64_t i, mine, maxs, ovlp, st, et, bl = 0, pc = 0;
	uint32_t o_l, o_r;
	ul_vec_t *p = NULL;
	nid_t *np = NULL;
	ul_ov_t *z = NULL;
	uc_block_t *b = NULL, tt;

	if(id) {
		kv_pushp(nid_t, x->nid, &np);
		np->n = id_l; MALLOC(np->a, np->n+1); memcpy(np->a, id, id_l); np->a[id_l] = '\0';
	}

	if(str||str_l) {
		if(rid == NULL) {
			kv_pushp(ul_vec_t, *x, &p);
			memset(p, 0, sizeof(*p));
		} else {
			if((*rid) >= x->m) kv_resize(ul_vec_t, *x, (*rid) + 1);
			if((*rid) >= x->n) {
				memset(x->a+x->n, 0, sizeof(*p)*((*rid) + 1 - x->n));
				x->n = (*rid) + 1;
			}
			p = &(x->a[(*rid)]);
		}
		// if((*rid) == 23) fprintf(stderr, "#rid->%lu, on->%ld\n", *rid, on);

		p->bb.n = p->N_site.n = p->r_base.n = 0; p->dd = 0;
		p->rlen = str_l; 
		
		if(o == NULL || on == 0) on = 0;
		for (i = on-1, st = et = str_l; i >= 0; i--) {
			z = &(o[i]);
			if(z->el) {
				mine = MIN(et, ((int64_t)z->qe)); maxs = MAX(st, ((int64_t)z->qs));
				ovlp = mine - maxs; 
				
				if(ovlp < 0) {///push original bases
					kv_pushp(uc_block_t, p->bb, &b);
					b->hid = 0/**x->mm**/; b->rev = 0; b->base = 1; b->pchain = 0; b->el = 0;
					b->qe = maxs; b->qs = b->qe + ovlp; bl += (b->qe-b->qs);
					o_l = (b->qs >= UL_FLANK?UL_FLANK:b->qs);
					o_r = ((str_l-b->qe)>=UL_FLANK?UL_FLANK:(str_l-b->qe));
					b->hid |= (o_l<<15); b->hid |= o_r;
					b->qs -= o_l; b->qe += o_r;
					b->ts = p->r_base.n; b->te = b->ts + B4L(b->qe-b->qs);
					kv_resize(uint8_t, p->r_base, b->te); p->r_base.n = b->te;
					// if(!str) fprintf(stderr, "+rid->%lu\n", *rid);
					ha_encode_base(p->r_base.a+b->ts, str+b->qs, b->qe-b->qs, &(p->N_site), b->qs);
				} 

				st = MIN(st, z->qs);
			}

			///push ovlp bases
			kv_pushp(uc_block_t, p->bb, &b);
			b->hid = (z->tn<<1)>>1; b->rev = z->rev; b->base = 0; b->el = z->el;
			b->pchain = ((z->tn&((uint32_t)(0x80000000)))?1:0);
			b->qs = z->qs; b->qe = z->qe;
			b->ts = z->ts; b->te = z->te;
			if(b->pchain) pc++;
			
			// st = MIN(st, z->qs);
		}
		
		if(st > 0) {///push original bases
			kv_pushp(uc_block_t, p->bb, &b);
			b->hid = 0/**x->mm**/; b->rev = 0; b->base = 1; b->pchain = 0; b->el = 0;
			b->qe = st; b->qs = 0; bl += (b->qe-b->qs);
			o_l = (b->qs >= UL_FLANK?UL_FLANK:b->qs);
			o_r = ((str_l-b->qe)>=UL_FLANK?UL_FLANK:(str_l-b->qe));
			b->hid |= (o_l<<15); b->hid |= o_r;
			b->qs -= o_l; b->qe += o_r;
			b->ts = p->r_base.n; b->te = b->ts + B4L(b->qe-b->qs);
			kv_resize(uint8_t, p->r_base, b->te); p->r_base.n = b->te;
			// if(!str) fprintf(stderr, "-rid->%lu, st->%ld, str_l->%ld\n", *rid, st, str_l);
			ha_encode_base(p->r_base.a+b->ts, str+b->qs, b->qe-b->qs, &(p->N_site), b->qs);	
			// push_subblock_original_bases(str, x, p, end, str_l, 321);//for debug
		}

		if(pc > 0) p->dd = 3;
		if((pc == on) && ((str_l-bl) > (str_l*p_chain_rate))) p->dd = 2;
		if((pc == on) && (bl == 0)) p->dd = 1;
		// debug_append_ul_t(o, on, p);
		// char *sst = NULL; CALLOC(sst, str_l);//for debug
		// retrieve_ul_t(NULL, sst, x, rid?*rid:x->n-1, 0, 0, -1);
		// if(memcmp(sst, str, str_l)) {
		// 	fprintf(stderr, "ap-Wrong read, id: %ld, [%d, %ld)\n", (int64_t)(rid?*rid:x->n-1), 0, str_l);
		// 	for (i = 0; i < str_l; i++) {
		// 		if(sst[i] != str[i]) fprintf(stderr, "[%ld] input:%c, decompress:%c\n", i, str[i], sst[i]);
		// 	}
		// }
		// free(sst);
		ovlp = p->bb.n>>1;
		for (i = 0; i < ovlp; i++) {
			tt = p->bb.a[i]; 
			p->bb.a[i] = p->bb.a[p->bb.n-i-1];
			p->bb.a[p->bb.n-i-1] = tt;
		}
		
	}
}

// Assembly:
pushq %rbp
pushq %r15
pushq %r14
pushq %r13
pushq %r12
pushq %rbx
subq $0x68, %rsp
movss %xmm0, 0x14(%rsp)
movq %r9, (%rsp)
movq %r8, 0x8(%rsp)
movq %rsi, %r12
movq %rdi, %r15
testq %rdx, %rdx
je 0x87dd
movq %rcx, %r13
movq %rdx, %rbx
movq (%r15), %rbp
cmpq 0x8(%r15), %rbp
jne 0x87a5
leaq (,%rbp,2), %rsi
testq %rbp, %rbp
pushq $0x2
popq %rax
cmoveq %rax, %rsi
movq %rsi, 0x8(%r15)
movq 0x10(%r15), %rdi
shlq $0x4, %rsi
callq 0x2370
movq %rax, %r14
movq %rax, 0x10(%r15)
movq (%r15), %rbp
jmp 0x87a9
movq 0x10(%r15), %r14
leaq 0x1(%rbp), %rax
movq %rax, (%r15)
shlq $0x4, %rbp
movl %r13d, 0x8(%r14,%rbp)
leal 0x1(%r13), %edi
callq 0x2280
movq %rax, (%r14,%rbp)
movq %rax, %rdi
movq %rbx, %rsi
movq %r13, %rdx
callq 0x2160
movq (%r14,%rbp), %rax
movb $0x0, (%rax,%r13)
movq 0x8(%rsp), %rax
orq (%rsp), %rax
je 0x8db5
movq 0xa8(%rsp), %rbp
movq 0xa0(%rsp), %rbx
testq %r12, %r12
je 0x88aa
movq (%r12), %rax
movq 0x58(%r15), %rcx
cmpq %rcx, %rax
setb %dl
leaq 0x1(%rax), %rsi
cmpq %rsi, %rcx
setae %cl
orb %dl, %cl
jne 0x886e
movq %rax, %rcx
shrq %rcx
orq %rax, %rcx
movq %rcx, %rax
shrq $0x2, %rax
orq %rcx, %rax
movq %rax, %rcx
shrq $0x4, %rcx
orq %rax, %rcx
movq %rcx, %rax
shrq $0x8, %rax
orq %rcx, %rax
movq %rax, %rcx
shrq $0x10, %rcx
orq %rax, %rcx
incq %rcx
movq %rcx, 0x58(%r15)
movq 0x48(%r15), %rdi
imulq $0x58, %rcx, %rsi
callq 0x2370
movq %rax, 0x48(%r15)
movq (%r12), %rax
movq 0x50(%r15), %rdx
movq %rax, %rcx
subq %rdx, %rcx
jb 0x88a0
imulq $0x58, %rdx, %rdi
addq 0x48(%r15), %rdi
imulq $0x58, %rcx, %rdx
addq $0x58, %rdx
xorl %esi, %esi
callq 0x20e0
movq (%r12), %rax
incq %rax
movq %rax, 0x50(%r15)
movq (%r12), %rax
imulq $0x58, %rax, %r13
addq 0x48(%r15), %r13
jmp 0x8912
movq 0x50(%r15), %rcx
cmpq 0x58(%r15), %rcx
jne 0x88dd
leaq (%rcx,%rcx), %rax
testq %rcx, %rcx
pushq $0x2
popq %rcx
cmoveq %rcx, %rax
movq %rax, 0x58(%r15)
movq 0x48(%r15), %rdi
imulq $0x58, %rax, %rsi
callq 0x2370
movq %rax, 0x48(%r15)
movq 0x50(%r15), %rcx
jmp 0x88e1
movq 0x48(%r15), %rax
leaq 0x1(%rcx), %rdx
movq %rdx, 0x50(%r15)
imulq $0x58, %rcx, %rcx
leaq (%rax,%rcx), %r13
xorps %xmm0, %xmm0
movups %xmm0, 0x40(%rax,%rcx)
movups %xmm0, 0x30(%rax,%rcx)
movups %xmm0, 0x20(%rax,%rcx)
movups %xmm0, 0x10(%rax,%rcx)
movups %xmm0, (%rax,%rcx)
andq $0x0, 0x50(%rax,%rcx)
andq $0x0, (%r13)
leaq 0x38(%r13), %rax
movq %rax, 0x20(%rsp)
andq $0x0, 0x40(%r13)
andq $0x0, 0x20(%r13)
movb $0x0, 0x50(%r13)
movq (%rsp), %rcx
movl %ecx, 0x18(%r13)
testq %rbx, %rbx
cmoveq %rbx, %rbp
imulq $0x1c, %rbp, %rax
leaq (%rax,%rbx), %r15
addq $-0x4, %r15
xorl %r12d, %r12d
movq %rbp, 0x18(%rsp)
movq %rcx, %r8
xorl %ebx, %ebx
testq %rbp, %rbp
jle 0x8baf
testb $0x40, 0x3(%r15)
je 0x8af6
movq %rbx, 0x38(%rsp)
movl -0x14(%r15), %eax
movl -0x10(%r15), %r14d
movq (%rsp), %rcx
cmpq %rcx, %r14
cmovgeq %rcx, %r14
cmpq %rax, %r8
movq %rax, %rbx
cmovgq %r8, %rbx
cmpq %r14, %rbx
jle 0x8aea
movq %r12, 0x30(%rsp)
movq 0x20(%r13), %rax
movq 0x30(%r13), %r12
cmpq 0x28(%r13), %rax
movq %r8, 0x28(%rsp)
jne 0x89d7
leaq (%rax,%rax), %rcx
testq %rax, %rax
pushq $0x2
popq %rax
cmoveq %rax, %rcx
movq %rcx, 0x28(%r13)
imulq $0x24, %rcx, %rsi
movq %r12, %rdi
callq 0x2370
movq %rax, %r12
movq %rax, 0x30(%r13)
movq 0x20(%r13), %rax
leaq 0x1(%rax), %rcx
movq %rcx, 0x20(%r13)
imulq $0x24, %rax, %rdi
movl $0x200, %esi            # imm = 0x200
cmpl %esi, %r14d
movl $0x200, %eax            # imm = 0x200
cmovbl %r14d, %eax
movl %ebx, %edx
movq (%rsp), %rcx
subq %rdx, %rcx
cmpq $0x200, %rcx            # imm = 0x200
cmovgel %esi, %ecx
movb $0x40, 0x20(%r12,%rdi)
movl %r14d, %edx
subl %eax, %edx
shll $0xf, %eax
orl %ecx, %eax
movl %eax, (%r12,%rdi)
movl %edx, 0x4(%r12,%rdi)
addl %ebx, %ecx
movl %ecx, 0x8(%r12,%rdi)
movl (%r13), %r8d
movl %r8d, 0xc(%r12,%rdi)
subl %edx, %ecx
movl %ecx, %esi
shrl $0x2, %esi
movl %ecx, %eax
andl $0x3, %eax
cmpl $0x1, %eax
sbbl $-0x1, %esi
addl %r8d, %esi
movl %esi, 0x10(%r12,%rdi)
movq 0x10(%r13), %rax
cmpq %rsi, 0x8(%r13)
jae 0x8ab4
addq %rdi, %r12
decq %rsi
movq %rsi, %rcx
shrq %rcx
orq %rsi, %rcx
movq %rcx, %rdx
shrq $0x2, %rdx
orq %rcx, %rdx
movq %rdx, %rcx
shrq $0x4, %rcx
orq %rdx, %rcx
movq %rcx, %rdx
shrq $0x8, %rdx
orq %rcx, %rdx
movq %rdx, %rsi
shrq $0x10, %rsi
orq %rdx, %rsi
incq %rsi
movq %rsi, 0x8(%r13)
movq %rax, %rdi
callq 0x2370
movq %rax, 0x10(%r13)
movl 0x10(%r12), %esi
movl 0xc(%r12), %r8d
movl 0x4(%r12), %edx
movl 0x8(%r12), %ecx
subl %edx, %ecx
subl %r14d, %ebx
movq 0x30(%rsp), %r12
addq %rbx, %r12
movq %rsi, (%r13)
movl %r8d, %edi
addq %rax, %rdi
movl %edx, %r8d
movq 0x8(%rsp), %rax
leaq (%rax,%r8), %rsi
movl %ecx, %edx
movq 0x20(%rsp), %rcx
callq 0x7dde
movl -0x14(%r15), %eax
movq 0x28(%rsp), %r8
cmpq %rax, %r8
cmovgeq %rax, %r8
movq 0x38(%rsp), %rbx
movq 0x20(%r13), %rcx
movq 0x30(%r13), %rax
cmpq 0x28(%r13), %rcx
jne 0x8b30
leaq (%rcx,%rcx), %rdx
testq %rcx, %rcx
pushq $0x2
popq %rcx
cmoveq %rcx, %rdx
movq %rdx, 0x28(%r13)
imulq $0x24, %rdx, %rsi
movq %rax, %rdi
movq %r8, %r14
callq 0x2370
movq %r14, %r8
movq %rax, 0x30(%r13)
movq 0x20(%r13), %rcx
leaq 0x1(%rcx), %rdx
movq %rdx, 0x20(%r13)
imulq $0x24, %rcx, %rcx
movl -0xc(%r15), %edx
movl $0x7fffffff, %esi       # imm = 0x7FFFFFFF
andl %esi, %edx
movl %edx, (%rax,%rcx)
movl (%r15), %edx
movb 0x20(%rax,%rcx), %sil
shrl $0x1a, %edx
andb $0x20, %dl
andb $-0x61, %sil
orb %dl, %sil
movb %sil, 0x20(%rax,%rcx)
movl (%r15), %edi
shrl $0x17, %edi
andb $-0x80, %dil
andb $0x3f, %sil
orb %dil, %sil
movb %sil, 0x20(%rax,%rcx)
movl -0xc(%r15), %esi
shrl $0x1f, %esi
addq %rsi, %rbx
orb %dl, %sil
orb %dil, %sil
movb %sil, 0x20(%rax,%rcx)
movsd -0x8(%r15), %xmm0
movsd -0x14(%r15), %xmm1
movlhps %xmm0, %xmm1            # xmm1 = xmm1[0],xmm0[0]
movups %xmm1, 0x4(%rax,%rcx)
addq $-0x1c, %r15
decq %rbp
jmp 0x8957
testq %r8, %r8
jle 0x8cf3
movq %r12, %r15
movq 0x20(%r13), %rax
movq 0x30(%r13), %r12
cmpq 0x28(%r13), %rax
movq 0x18(%rsp), %r14
movq %rbx, %rbp
jne 0x8c00
leaq (%rax,%rax), %rcx
testq %rax, %rax
pushq $0x2
popq %rax
cmoveq %rax, %rcx
movq %rcx, 0x28(%r13)
imulq $0x24, %rcx, %rsi
movq %r12, %rdi
movq %r8, %rbx
callq 0x2370
movq %rbx, %r8
movq %rax, %r12
movq %rax, 0x30(%r13)
movq 0x20(%r13), %rax
leaq 0x1(%rax), %rcx
movq %rcx, 0x20(%r13)
imulq $0x24, %rax, %rdi
movb $0x40, 0x20(%r12,%rdi)
andl $0x0, 0x4(%r12,%rdi)
movl %r8d, %ebx
movq (%rsp), %rax
subq %rbx, %rax
cmpq $0x200, %rax            # imm = 0x200
movl $0x200, %edx            # imm = 0x200
cmovll %eax, %edx
movl %edx, (%r12,%rdi)
addl %r8d, %edx
movl %edx, 0x8(%r12,%rdi)
movl (%r13), %esi
movl %esi, 0xc(%r12,%rdi)
movl %edx, %ecx
shrl $0x2, %ecx
movl %edx, %eax
andl $0x3, %eax
cmpl $0x1, %eax
sbbl $-0x1, %ecx
addl %esi, %ecx
movl %ecx, 0x10(%r12,%rdi)
movq 0x10(%r13), %rax
xorl %r9d, %r9d
cmpq %rcx, 0x8(%r13)
jae 0x8cc7
addq %rdi, %r12
decq %rcx
movq %rcx, %rdx
shrq %rdx
orq %rcx, %rdx
movq %rdx, %rcx
shrq $0x2, %rcx
orq %rdx, %rcx
movq %rcx, %rdx
shrq $0x4, %rdx
orq %rcx, %rdx
movq %rdx, %rcx
shrq $0x8, %rcx
orq %rdx, %rcx
movq %rcx, %rsi
shrq $0x10, %rsi
orq %rcx, %rsi
incq %rsi
movq %rsi, 0x8(%r13)
movq %rax, %rdi
callq 0x2370
movq %rax, 0x10(%r13)
movl 0x10(%r12), %ecx
movl 0xc(%r12), %esi
movl 0x4(%r12), %r9d
movl 0x8(%r12), %edx
movq %r15, %r12
addq %rbx, %r12
movq %rcx, (%r13)
movl %esi, %edi
addq %rax, %rdi
movl %r9d, %r8d
movq 0x8(%rsp), %rsi
addq %r8, %rsi
subl %r9d, %edx
movq 0x20(%rsp), %rcx
callq 0x7dde
movq %rbp, %rbx
jmp 0x8cf8
movq 0x18(%rsp), %r14
testq %rbx, %rbx
je 0x8d02
movb $0x3, 0x50(%r13)
cmpq %r14, %rbx
jne 0x8d34
testq %r12, %r12
je 0x8d28
movq (%rsp), %rax
cvtsi2ss %rax, %xmm0
subq %r12, %rax
cvtsi2ss %rax, %xmm1
mulss 0x14(%rsp), %xmm0
ucomiss %xmm0, %xmm1
jbe 0x8d34
cmpq $0x1, %r12
movb $0x2, %al
sbbb $0x0, %al
movb %al, 0x50(%r13)
movq 0x20(%r13), %rax
shrq %rax
pushq $-0x24
popq %rcx
xorl %edx, %edx
subq $0x1, %rax
jb 0x8db5
movq 0x30(%r13), %rsi
movl 0x20(%rsi,%rdx), %edi
movl %edi, 0x60(%rsp)
movups (%rsi,%rdx), %xmm0
movups 0x10(%rsi,%rdx), %xmm1
movaps %xmm1, 0x50(%rsp)
movaps %xmm0, 0x40(%rsp)
imulq $0x24, 0x20(%r13), %rdi
addq %rsi, %rdi
movups (%rcx,%rdi), %xmm0
movups 0x10(%rcx,%rdi), %xmm1
movups %xmm0, (%rsi,%rdx)
movups %xmm1, 0x10(%rsi,%rdx)
movl 0x20(%rcx,%rdi), %edi
movl %edi, 0x20(%rsi,%rdx)
imulq $0x24, 0x20(%r13), %rsi
addq 0x30(%r13), %rsi
movl 0x60(%rsp), %edi
movl %edi, 0x20(%rcx,%rsi)
movaps 0x40(%rsp), %xmm0
movaps 0x50(%rsp), %xmm1
movups %xmm1, 0x10(%rcx,%rsi)
movups %xmm0, (%rcx,%rsi)
addq $-0x24, %rcx
addq $0x24, %rdx
jmp 0x8d40
addq $0x68, %rsp
popq %rbx
popq %r12
popq %r13
popq %r14
popq %r15
popq %rbp
retq