+typedef struct __bmc_aux_t {
+ int max;
+ uint32_t *info;
+ uint16_t *info16;
+ errmod_t *em;
+} bmc_aux_t;
+
+uint32_t pysam_glf_depth( glf1_t * g )
+{
+ return g->depth;
+}
+
+
+void pysam_dump_glf( glf1_t * g, bam_maqcns_t * c )
+{
+ int x = 0;
+ fprintf(stderr,
+ "glf: ref_base=%i, max_mapQ=%i, min_lk=%i, depth=%i",
+ g->ref_base,
+ g->max_mapQ,
+ g->min_lk,
+ g->depth );
+
+ for (x = 0; x < 10; ++x)
+ fprintf(stderr, ", lk%x=%i, ", x, g->lk[x]);
+
+ fprintf(stderr,
+ "maqcns: het_rate=%f, theta=%f, n_hap=%i, cap_mapQ=%i, errmod=%i, min_baseQ=%i, eta=%f, q_r=%f, aux_max=%i",
+ c->het_rate,
+ c->theta,
+ c->n_hap,
+ c->cap_mapQ,
+ c->errmod,
+ c->min_baseQ,
+ c->eta,
+ c->q_r,
+ c->aux->max);
+
+ for (x = 0; x < c->aux->max; ++x)
+ {
+ fprintf(stderr, ", info-%i=%i ", x, c->aux->info[x]);
+ if (c->aux->info[x] == 0) break;
+ }
+
+ for (x = 0; x < c->aux->max; ++x)
+ {
+ fprintf(stderr, ", info16-%i=%i ", x, c->aux->info16[x]);
+ if (c->aux->info16[x] == 0) break;
+ }
+}
+
+
+
+