X-Git-Url: http://woldlab.caltech.edu/gitweb/?p=samtools.git;a=blobdiff_plain;f=bam_maqcns.h;h=291ae533c3c38cd152b3c9c8fa8e820386a3718d;hp=2a82aeee8057f5992570d0896a1621b9116057cd;hb=016b50ab60e879a0b8f81cb76ce11ea360a03d4a;hpb=d363084f0412f3bcdeb0304aeb0974c9a10c7649 diff --git a/bam_maqcns.h b/bam_maqcns.h index 2a82aee..291ae53 100644 --- a/bam_maqcns.h +++ b/bam_maqcns.h @@ -3,11 +3,15 @@ #include "glf.h" +#define BAM_ERRMOD_MAQ2 0 +#define BAM_ERRMOD_MAQ 1 +#define BAM_ERRMOD_SOAP 2 + struct __bmc_aux_t; typedef struct { float het_rate, theta; - int n_hap, cap_mapQ; + int n_hap, cap_mapQ, errmod, min_baseQ; float eta, q_r; double *fk, *coef; @@ -16,8 +20,9 @@ typedef struct { } bam_maqcns_t; typedef struct { - int q_indel; - float r_indel; + int q_indel; // indel sequencing error, phred scaled + float r_indel; // indel prior + float r_snp; // snp prior // hidden parameters, unchangeable from command line int mm_penalty, indel_err, ambi_thres; } bam_maqindel_opt_t;