// File: /workspace/github_c_makefile_star2_O2/icebert[P]pblat/htslib/sam.c
// Function: bam_hdr_read
bam_hdr_t *bam_hdr_read(BGZF *fp)
{
    bam_hdr_t *h;
    char buf[4];
    int magic_len;
    int32_t i, name_len, num_names = 0;
    size_t bufsize;
    ssize_t bytes;
#ifndef UCSC_CRAM
    // check EOF
    int has_EOF = bgzf_check_EOF(fp);
    if (has_EOF < 0) {
        perror("[W::bam_hdr_read] bgzf_check_EOF");
    } else if (has_EOF == 0 && hts_verbose >= 2)
        fprintf(stderr, "[W::%s] EOF marker is absent. The input is probably truncated.\n", __func__);
#endif
    // read "BAM1"
    magic_len = bgzf_read(fp, buf, 4);
    if (magic_len != 4 || strncmp(buf, "BAM\1", 4)) {
        if (hts_verbose >= 1) fprintf(stderr, "[E::%s] invalid BAM binary header\n", __func__);
        return 0;
    }
    h = bam_hdr_init();
    if (!h) goto nomem;

    // read plain text and the number of reference sequences
    bytes = bgzf_read(fp, &h->l_text, 4);
    if (bytes != 4) goto read_err;
    if (fp->is_be) ed_swap_4p(&h->l_text);

    bufsize = ((size_t) h->l_text) + 1;
    if (bufsize < h->l_text) goto nomem; // so large that adding 1 overflowed
    h->text = (char*)malloc(bufsize);
    if (!h->text) goto nomem;
    h->text[h->l_text] = 0; // make sure it is NULL terminated
    bytes = bgzf_read(fp, h->text, h->l_text);
    if (bytes != h->l_text) goto read_err;

    bytes = bgzf_read(fp, &h->n_targets, 4);
    if (bytes != 4) goto read_err;
    if (fp->is_be) ed_swap_4p(&h->n_targets);

    if (h->n_targets < 0) goto invalid;

    // read reference sequence names and lengths
    if (h->n_targets > 0) {
        h->target_name = (char**)calloc(h->n_targets, sizeof(char*));
        if (!h->target_name) goto nomem;
        h->target_len = (uint32_t*)calloc(h->n_targets, sizeof(uint32_t));
        if (!h->target_len) goto nomem;
    }
    else {
        h->target_name = NULL;
        h->target_len = NULL;
    }

    for (i = 0; i != h->n_targets; ++i) {
        bytes = bgzf_read(fp, &name_len, 4);
        if (bytes != 4) goto read_err;
        if (fp->is_be) ed_swap_4p(&name_len);
        if (name_len <= 0) goto invalid;

        h->target_name[i] = (char*)malloc(name_len);
        if (!h->target_name[i]) goto nomem;
        num_names++;

        bytes = bgzf_read(fp, h->target_name[i], name_len);
        if (bytes != name_len) goto read_err;

        if (h->target_name[i][name_len - 1] != '\0') {
            /* Fix missing NUL-termination.  Is this being too nice?
               We could alternatively bail out with an error. */
            char *new_name;
            if (name_len == INT32_MAX) goto invalid;
            new_name = realloc(h->target_name[i], name_len + 1);
            if (new_name == NULL) goto nomem;
            h->target_name[i] = new_name;
            h->target_name[i][name_len] = '\0';
        }

        bytes = bgzf_read(fp, &h->target_len[i], 4);
        if (bytes != 4) goto read_err;
        if (fp->is_be) ed_swap_4p(&h->target_len[i]);
    }
    return h;

 nomem:
    if (hts_verbose >= 1) fprintf(stderr, "[E::%s] out of memory\n", __func__);
    goto clean;

 read_err:
    if (hts_verbose >= 1) {
        if (bytes < 0) {
            fprintf(stderr, "[E::%s] error reading BGZF stream\n", __func__);
        } else {
            fprintf(stderr, "[E::%s] truncated bam header\n", __func__);
        }
    }
    goto clean;

 invalid:
    if (hts_verbose >= 1) {
        fprintf(stderr, "[E::%s] invalid BAM binary header\n", __func__);
    }

 clean:
    if (h != NULL) {
        h->n_targets = num_names; // ensure we free only allocated target_names
        bam_hdr_destroy(h);
    }
    return NULL;
}

// Assembly:
pushq %rbp
pushq %r15
pushq %r14
pushq %r13
pushq %r12
pushq %rbx
pushq %rax
movq %rdi, %r14
callq 0x5d323
testl %eax, %eax
js 0x6146d
jne 0x613f0
leaq 0x53a01(%rip), %rax     # 0xb4de8
cmpl $0x2, (%rax)
jge 0x61485
leaq 0x4(%rsp), %rsi
pushq $0x4
popq %rdx
movq %r14, %rdi
callq 0x5ca56
cmpl $0x4, %eax
jne 0x6140f
cmpl $0x14d4142, 0x4(%rsp)   # imm = 0x14D4142
je 0x61433
leaq 0x539d2(%rip), %rax     # 0xb4de8
cmpl $0x0, (%rax)
jg 0x614a9
xorl %ebx, %ebx
movq %rbx, %rax
addq $0x8, %rsp
popq %rbx
popq %r12
popq %r13
popq %r14
popq %r15
popq %rbp
retq
callq 0x61280
movq %rax, %rbx
testq %rax, %rax
je 0x61533
leaq 0x8(%rbx), %r15
pushq $0x4
popq %rdx
movq %r14, %rdi
movq %r15, %rsi
callq 0x5ca56
cmpq $0x4, %rax
jne 0x6147e
testb $0xc, 0x2(%r14)
je 0x614cf
movl (%r15), %eax
bswapl %eax
movl %eax, (%r15)
jmp 0x614d2
leaq 0x27dfc(%rip), %rdi     # 0x89270
callq 0x6590
jmp 0x613f0
xorl %ebp, %ebp
jmp 0x616d6
movq 0x51b4c(%rip), %rax     # 0xb2fd8
movq (%rax), %rdi
leaq 0x2614e(%rip), %rsi     # 0x875e4
leaq 0x27df4(%rip), %rdx     # 0x89291
xorl %eax, %eax
callq 0x6a90
jmp 0x613f0
movq 0x51b28(%rip), %rax     # 0xb2fd8
movq (%rax), %rdi
leaq 0x27de9(%rip), %rsi     # 0x892a3
leaq 0x27dd0(%rip), %rdx     # 0x89291
xorl %ebx, %ebx
xorl %eax, %eax
callq 0x6a90
jmp 0x61421
movl (%r15), %eax
movl %eax, %r12d
leaq 0x1(%r12), %rdi
callq 0x69d0
movq %rax, 0x28(%rbx)
testq %rax, %rax
je 0x61533
movb $0x0, (%rax,%r12)
movq %r14, %rdi
movq %rax, %rsi
movq %r12, %rdx
callq 0x5ca56
movl (%r15), %ecx
xorl %ebp, %ebp
cmpq %rcx, %rax
jne 0x616d6
pushq $0x4
popq %rdx
movq %r14, %rdi
movq %rbx, %rsi
callq 0x5ca56
cmpq $0x4, %rax
jne 0x616d6
testb $0xc, 0x2(%r14)
je 0x61570
movl (%rbx), %r13d
bswapl %r13d
movl %r13d, (%rbx)
jmp 0x61573
xorl %ebp, %ebp
leaq 0x538ac(%rip), %rax     # 0xb4de8
cmpl $0x0, (%rax)
jg 0x6154f
testq %rbx, %rbx
jne 0x616e2
jmp 0x6141f
movq 0x51a82(%rip), %rax     # 0xb2fd8
movq (%rax), %rdi
leaq 0x27d66(%rip), %rsi     # 0x892c6
leaq 0x27d2a(%rip), %rdx     # 0x89291
xorl %eax, %eax
callq 0x6a90
jmp 0x61541
movl (%rbx), %r13d
testl %r13d, %r13d
js 0x615a9
je 0x615d1
movl %r13d, %r15d
pushq $0x8
popq %rsi
movq %r15, %rdi
callq 0x6470
movq %rax, 0x20(%rbx)
xorl %ebp, %ebp
testq %rax, %rax
je 0x61535
pushq $0x4
popq %rsi
movq %r15, %rdi
callq 0x6470
movq %rax, 0x10(%rbx)
testq %rax, %rax
jne 0x615db
jmp 0x61535
xorl %ebp, %ebp
leaq 0x53836(%rip), %rax     # 0xb4de8
cmpl $0x0, (%rax)
jle 0x616e2
movq 0x51a16(%rip), %rax     # 0xb2fd8
movq (%rax), %rdi
leaq 0x27cd7(%rip), %rsi     # 0x892a3
jmp 0x61710
andq $0x0, 0x20(%rbx)
andq $0x0, 0x10(%rbx)
xorl %ebp, %ebp
pushq $0x4
popq %r12
movl %r13d, %eax
cmpq %rax, %rbp
je 0x61421
movq %r14, %rdi
movq %rsp, %rsi
movq %r12, %rdx
callq 0x5ca56
cmpq $0x4, %rax
jne 0x616d6
testb $0xc, 0x2(%r14)
movl (%rsp), %eax
je 0x61614
bswapl %eax
movl %eax, (%rsp)
testl %eax, %eax
jle 0x615ab
movl %eax, %r13d
movq %r13, %rdi
callq 0x69d0
movq 0x20(%rbx), %rcx
movq %rax, (%rcx,%rbp,8)
movq 0x20(%rbx), %rax
movq (%rax,%rbp,8), %rsi
testq %rsi, %rsi
je 0x61535
leaq 0x1(%rbp), %r15
movq %r14, %rdi
movq %r13, %rdx
callq 0x5ca56
movslq (%rsp), %rsi
cmpq %rsi, %rax
jne 0x616d3
movq 0x20(%rbx), %rcx
movq (%rcx,%rbp,8), %rdi
cmpb $0x0, -0x1(%rdi,%rax)
je 0x61698
cmpl $0x7fffffff, %esi       # imm = 0x7FFFFFFF
je 0x61720
incq %rsi
callq 0x6c80
testq %rax, %rax
je 0x61728
movq 0x20(%rbx), %rcx
movq %rax, (%rcx,%rbp,8)
movq 0x20(%rbx), %rax
movq (%rax,%rbp,8), %rax
movslq (%rsp), %rcx
movb $0x0, (%rax,%rcx)
leaq (,%rbp,4), %rsi
addq 0x10(%rbx), %rsi
movq %r14, %rdi
movq %r12, %rdx
callq 0x5ca56
cmpq $0x4, %rax
jne 0x616d3
testb $0xc, 0x2(%r14)
je 0x616c8
movq 0x10(%rbx), %rax
movl (%rax,%rbp,4), %ecx
bswapl %ecx
movl %ecx, (%rax,%rbp,4)
movl (%rbx), %r13d
movq %r15, %rbp
jmp 0x615e1
movq %r15, %rbp
leaq 0x5370b(%rip), %rcx     # 0xb4de8
cmpl $0x0, (%rcx)
jg 0x616f1
movl %ebp, (%rbx)
movq %rbx, %rdi
callq 0x6128b
jmp 0x6141f
movq 0x518e0(%rip), %rcx     # 0xb2fd8
movq (%rcx), %rdi
testq %rax, %rax
jns 0x61709
leaq 0x27bd6(%rip), %rsi     # 0x892dd
jmp 0x61710
leaq 0x27bf0(%rip), %rsi     # 0x89300
leaq 0x27b7a(%rip), %rdx     # 0x89291
xorl %eax, %eax
callq 0x6a90
jmp 0x616e2
movq %r15, %rbp
jmp 0x615ab
movq %r15, %rbp
jmp 0x61535