qbix 0.1.1

Random access to BAM records by read name using a compact .qbi index
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80

#include <stdint.h>
#include <stddef.h>
#include <stdlib.h>
#include <htslib/bgzf.h>
#include <htslib/hts.h>
#include <htslib/kstring.h>
#include <htslib/sam.h>
#include <htslib/tbx.h>

int64_t qbix_hts_shim_bgzf_tell(htsFile *fp)
{
    BGZF *bgzf = hts_get_bgzfp(fp);
    return bgzf == 0 ? -1 : bgzf_tell(bgzf);
}

int64_t qbix_hts_shim_bgzf_seek(htsFile *fp, int64_t offset)
{
    BGZF *bgzf = hts_get_bgzfp(fp);
    return bgzf == 0 ? -1 : bgzf_seek(bgzf, offset, SEEK_SET);
}

int qbix_hts_shim_bgzf_set_cache_size(htsFile *fp, int size)
{
    BGZF *bgzf = hts_get_bgzfp(fp);
    if (bgzf == 0) return -1;
    bgzf_set_cache_size(bgzf, size);
    return 0;
}

const char *qbix_hts_shim_bam_qname(bam1_t *b)
{
    return bam_get_qname(b);
}

const char *qbix_hts_shim_sam_hdr_text(const sam_hdr_t *h)
{
    return sam_hdr_str((sam_hdr_t *)h);
}

size_t qbix_hts_shim_sam_hdr_text_len(const sam_hdr_t *h)
{
    return sam_hdr_length((sam_hdr_t *)h);
}

int qbix_hts_shim_set_unknown_sort_order(sam_hdr_t *h)
{
    if (sam_hdr_count_lines(h, "HD") == 0) {
        if (sam_hdr_add_line(h, "HD", "VN", "1.6", "SO", "unknown", NULL) < 0)
            return -1;
    } else if (sam_hdr_update_hd(h, "SO", "unknown") < 0) {
        return -1;
    }
    if (sam_hdr_remove_tag_hd(h, "SS") < 0) return -1;
    return sam_hdr_remove_tag_hd(h, "GO") < 0 ? -1 : 0;
}

int qbix_hts_shim_add_pg(sam_hdr_t *h, const char *version, const char *command_line)
{
    return sam_hdr_add_pg(h, "qbix", "VN", version, "CL", command_line, NULL);
}

int qbix_hts_shim_sam_format1(const sam_hdr_t *h, const bam1_t *b, char **out, size_t *len)
{
    kstring_t str = {0, 0, 0};
    if (out == 0 || len == 0) return -1;
    *out = 0;
    *len = 0;
    if (sam_format1(h, b, &str) < 0) {
        free(str.s);
        return -1;
    }
    *out = str.s;
    *len = str.l;
    return 0;
}

void qbix_hts_shim_free(void *ptr)
{
    free(ptr);
}