errmod_t *errmod_init(float depcorr);
void errmod_destroy(errmod_t *em);
+
+/*
+ n: number of bases
+ m: maximum base
+ bases[i]: qual:6, strand:1, base:4
+ q[i*m+j]: phred-scaled likelihood of (i,j)
+ */
int errmod_cal(const errmod_t *em, int n, int m, uint16_t *bases, float *q);
#endif