1 /** htslib-1.9 tbx.h as D module 2 * 3 * Changes include: 4 * Removed if(n)defs 5 * Change numeric #defines to enum int 6 * Changed ^typedef struct {...} <name>$ to ^struct <name> {...}$ 7 * extern const to __gshared 8 * made #define function macros into inline functions (tbx_itr* -> hts_itr*) 9 * In D, const on either LHS or RHS of function declaration applies to the function, not return value, unless parents included: 10 * changed ^const <type> <fnname> to ^const(<type>) <fnname> 11 */ 12 module dhtslib.htslib.tbx; 13 14 import std.stdint : int32_t; 15 16 import dhtslib.htslib.hts; 17 import dhtslib.htslib.bgzf; 18 19 extern (C): 20 /// @file htslib/tbx.h 21 /// Tabix API functions. 22 /* 23 Copyright (C) 2009, 2012-2015 Genome Research Ltd. 24 Copyright (C) 2010, 2012 Broad Institute. 25 Author: Heng Li <lh3@sanger.ac.uk> 26 Permission is hereby granted, free of charge, to any person obtaining a copy 27 of this software and associated documentation files (the "Software"), to deal 28 in the Software without restriction, including without limitation the rights 29 to use, copy, modify, merge, publish, distribute, sublicense, and/or sell 30 copies of the Software, and to permit persons to whom the Software is 31 furnished to do so, subject to the following conditions: 32 The above copyright notice and this permission notice shall be included in 33 all copies or substantial portions of the Software. 34 THE SOFTWARE IS PROVIDED "AS IS", WITHOUT WARRANTY OF ANY KIND, EXPRESS OR 35 IMPLIED, INCLUDING BUT NOT LIMITED TO THE WARRANTIES OF MERCHANTABILITY, 36 FITNESS FOR A PARTICULAR PURPOSE AND NONINFRINGEMENT. IN NO EVENT SHALL 37 THE AUTHORS OR COPYRIGHT HOLDERS BE LIABLE FOR ANY CLAIM, DAMAGES OR OTHER 38 LIABILITY, WHETHER IN AN ACTION OF CONTRACT, TORT OR OTHERWISE, ARISING 39 FROM, OUT OF OR IN CONNECTION WITH THE SOFTWARE OR THE USE OR OTHER 40 DEALINGS IN THE SOFTWARE. */ 41 42 enum int TBX_MAX_SHIFT = 31;/// ??? 43 44 enum int TBX_GENERIC = 0; /// generic flat file 45 enum int TBX_SAM = 1; /// SAM 46 enum int TBX_VCF = 2; /// VCF 47 enum int TBX_UCSC = 0x10000; /// ?UCSC flat file? 48 49 /// tabix config 50 struct tbx_conf_t { 51 int32_t preset; /// ? 52 /// seq col., beg col. and end col. 53 int32_t sc, bc, ec; // seq col., beg col. and end col. 54 /// ? 55 int32_t meta_char, line_skip; 56 } 57 58 /// tabix data 59 struct tbx_t { 60 tbx_conf_t conf; /// tabix config 61 hts_idx_t *idx; /// index data 62 void *dict; /// ?dictionary 63 } 64 65 //extern const tbx_conf_t tbx_conf_gff, tbx_conf_bed, tbx_conf_psltbl, tbx_conf_sam, tbx_conf_vcf; 66 /// prebaked TABIX config data for GFF3, BED, PSL table, SAM, VCF 67 extern (C) extern __gshared const tbx_conf_t tbx_conf_gff, tbx_conf_bed, tbx_conf_psltbl, tbx_conf_sam, tbx_conf_vcf; 68 69 alias tbx_itr_destroy = hts_itr_destroy; 70 71 /* hts_itr_t *hts_itr_query(const hts_idx_t *idx, int tid, int beg, int end, hts_readrec_func *readrec); */ 72 //#define tbx_itr_queryi(tbx, tid, beg, end) hts_itr_query((tbx)->idx, (tid), (beg), (end), tbx_readrec) 73 /// tabix query by integer based tid(contig)/start/end 74 pragma(inline, true) 75 auto tbx_itr_queryi(const tbx_t *tbx, int tid, int beg, int end) 76 { return hts_itr_query(tbx.idx, tid, beg, end, &tbx_readrec); } 77 78 /* hts_itr_t *hts_itr_querys(const hts_idx_t *idx, const char *reg, hts_name2id_f getid, void *hdr, hts_itr_query_func *itr_query, hts_readrec_func *readrec); */ 79 //#define tbx_itr_querys(tbx, s) hts_itr_querys((tbx)->idx, (s), (hts_name2id_f)(tbx_name2id), (tbx), hts_itr_query, tbx_readrec) 80 /// tabix query by string "chr:start-end" 81 pragma(inline, true) 82 auto tbx_itr_querys(const tbx_t *tbx, const char *s) 83 { 84 return hts_itr_querys(tbx.idx, s, 85 cast(hts_name2id_f)(&tbx_name2id), 86 cast(void*)tbx, 87 &hts_itr_query, 88 &tbx_readrec); 89 } 90 91 /* int hts_itr_next(BGZF *fp, hts_itr_t *iter, void *r, void *data) HTS_RESULT_USED; */ 92 //#define tbx_itr_next(htsfp, tbx, itr, r) hts_itr_next(hts_get_bgzfp(htsfp), (itr), (r), (tbx)) 93 /// advance tabix iterator 94 pragma(inline, true) 95 auto tbx_itr_next(htsFile *htsfp, tbx_t *tbx, hts_itr_t *itr, void *r) 96 { return hts_itr_next(hts_get_bgzfp(htsfp), itr, r, tbx); } 97 98 /* int hts_itr_next(BGZF *fp, hts_itr_t *iter, void *r, void *data) HTS_RESULT_USED; */ 99 //#define tbx_bgzf_itr_next(bgzfp, tbx, itr, r) hts_itr_next((bgzfp), (itr), (r), (tbx)) 100 /// advance tabix iterator 101 pragma(inline, true) 102 auto tbx_bgzf_itr_next(BGZF *bgzfp, tbx_t *tbx, hts_itr_t *itr, void *r) 103 { return hts_itr_next(bgzfp, itr, r, tbx); } 104 105 /// contig name to integer id 106 int tbx_name2id(tbx_t *tbx, const char *ss); 107 108 /** Internal helper function used by tbx_itr_next() defined in hts.c -- do not use directly */ 109 BGZF *hts_get_bgzfp(htsFile *fp); 110 /** Called by tabix iterator to read the next record */ 111 int tbx_readrec(BGZF *fp, void *tbxv, void *sv, int *tid, int *beg, int *end); 112 113 /// build tabix index 114 tbx_t *tbx_index(BGZF *fp, int min_shift, const tbx_conf_t *conf); 115 /// build tabix index 116 int tbx_index_build(const char *fn, int min_shift, const tbx_conf_t *conf); 117 /// ditto 118 int tbx_index_build2(const char *fn, const char *fnidx, int min_shift, const tbx_conf_t *conf); 119 /// ditto 120 int tbx_index_build3(const char *fn, const char *fnidx, int min_shift, int n_threads, const tbx_conf_t *conf); 121 /// load tabix index 122 tbx_t *tbx_index_load(const char *fn); 123 /// ditto 124 tbx_t *tbx_index_load2(const char *fn, const char *fnidx); 125 /// return C-style array of sequence names (NB: free the array but not the values) 126 const(char **) tbx_seqnames(tbx_t *tbx, int *n); // free the array but not the values 127 /// destroy/dealloc tabix data 128 void tbx_destroy(tbx_t *tbx);