BGI/SOAPdenovo2/sparsePregraph/inc/glf.h

   1 #ifndef GLF_H_
   2 #define GLF_H_
   3
   4 typedef struct
   5 {
   6         unsigned char ref_base: 4, dummy: 4; /** "XACMGRSVTWYHKDBN"[ref_base] gives the reference base */
   7         unsigned char max_mapQ; /** maximum mapping quality */
   8         unsigned char lk[10];   /** log likelihood ratio, capped at 255 */
   9         unsigned min_lk: 8, depth: 24; /** minimum lk capped at 255, and the number of mapped reads */
  10 } glf1_t;
  11
  12 #include <stdint.h>
  13 #include "bgzf.h"
  14 typedef BGZF * glfFile;
  15
  16 #define GLF3_RTYPE_END   0
  17 #define GLF3_RTYPE_SUB   1
  18 #define GLF3_RTYPE_INDEL 2
  19
  20 typedef struct
  21 {
  22         uint8_t ref_base: 4, rtype: 4; /** "XACMGRSVTWYHKDBN"[ref_base] gives the reference base */
  23         uint8_t rms_mapQ; /** RMS mapping quality */
  24         uint8_t lk[10];   /** log likelihood ratio, capped at 255 */
  25         uint32_t min_lk: 8, depth: 24; /** minimum lk capped at 255, and the number of mapped reads */
  26         int32_t offset; /** the first base in a chromosome has offset zero. */
  27         // for indel (lkHom1, lkHom2 and lkHet are the first three elements in lk[10])
  28         int16_t indel_len[2];
  29         int32_t max_len; // maximum indel len; will be modified by glf3_read1()
  30         char * indel_seq[2];
  31 } glf3_t;
  32
  33 typedef struct
  34 {
  35         int32_t l_text;
  36         uint8_t * text;
  37 } glf3_header_t;
  38
  39 #ifdef __cplusplus
  40 extern "C" {
  41 #endif
  42
  43 #define glf3_init1() ((glf3_t*)calloc(1, sizeof(glf3_t)))
  44 #define glf3_destroy1(g3) do { free((g3)->indel_seq[0]); free((g3)->indel_seq[1]); free(g3); } while (0)
  45
  46         glf3_header_t * glf3_header_init();
  47         glf3_header_t * glf3_header_read ( glfFile fp );
  48         void glf3_header_write ( glfFile fp, const glf3_header_t * h );
  49         void glf3_header_destroy ( glf3_header_t * h );
  50         char * glf3_ref_read ( glfFile fp, int * len );
  51         void glf3_ref_write ( glfFile fp, const char * name, int len );
  52         int glf3_write1 ( glfFile fp, const glf3_t * g3 );
  53         int glf3_read1 ( glfFile fp, glf3_t * g3 );
  54
  55 #ifdef __cplusplus
  56 }
  57 #endif
  58
  59 #endif