htslib + zig

const std = @import("std");

const hts = @cImport({
    // htslib.org
    @cInclude("htslib/hts.h");
    @cInclude("htslib/sam.h");
});

pub fn main() !void {
    const stdout = std.io.getStdOut().writer();
    var v = hts.hts_version();
    try stdout.print("version, {s}!\n", .{v});

    var fh = hts.hts_open("/data/human/hg002.22.bam", "r");
    if (fh == null) {
        std.debug.warn("null file: {}", .{v});
    }
    try stdout.print("eof?:{d}\n", .{hts.hts_check_EOF(fh)});

    var hdr = hts.sam_hdr_read(fh);
    try stdout.print("null header?:{d}\n", .{hdr == null});

    var idx = hts.sam_index_load(fh, "/data/human/hg002.22.bam");
    try stdout.print("null idx?:{d}\n", .{idx == null});

    var qiter = hts.sam_itr_querys(idx, hdr, "22:450000");
    try stdout.print("null itr?:{d}\n", .{qiter == null});

    var b = hts.bam_init1();
    // ? https://github.com/ziglang/zig/issues/1580
    var ret = hts.sam_itr_next(fh, qiter, b);

    try stdout.print("next ret?:{d}\n", .{ret});
}