X-Git-Url: http://woldlab.caltech.edu/gitweb/?p=samtools.git;a=blobdiff_plain;f=bam_maqcns.h;h=291ae533c3c38cd152b3c9c8fa8e820386a3718d;hp=fa5489d13caad7c5c20489091302bfb5b2a5799d;hb=9f4bebab2e0917c676ae739b2d05cb22ad6c4ed5;hpb=317f5e8dd22cc9e1e5e05fbcaeb3b9aca7447351 diff --git a/bam_maqcns.h b/bam_maqcns.h index fa5489d..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, is_soap; + 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;