deev

hitbc · Mar 17, 2016 · 53082ce · 53082ce
1 parent 56a8cac
commit 53082ce
Show file tree

Hide file tree

Showing 2 changed files with 4 additions and 3 deletions.
diff --git a/src/lamsa_aln.c b/src/lamsa_aln.c
@@ -246,7 +246,8 @@ int split_seed(const char *prefix, lamsa_aln_para AP, seed_msg *s_msg)
     m_read = s_msg->read_m;
     while (kseq_read(seq) >= 0)
     {
-        seed_all = (1+ (seq->seq.l - seed_len) / seed_step);
+        if (seq->seq.l < seed_len) seed_all = 0;
+        else seed_all = (1+ (seq->seq.l - seed_len) / seed_step);
         seed_seq[seed_len] = '\n';
         if (seed_all > s_msg->seed_max) s_msg->seed_max = seed_all;
         if ((int)seq->seq.l > s_msg->read_max_len) s_msg->read_max_len = seq->seq.l;

diff --git a/src/lamsa_aln.h b/src/lamsa_aln.h
@@ -63,10 +63,10 @@
 #define ON_SEED_LEN 50
 #define ON_SEED_STEP 25
 
-#define ON_ED_RATE 0.4 // 0.45
+#define ON_ED_RATE 0.25 // 0.45
 #define ON_ID_RATE 0.1 
 #define ON_MAT_RATE 0.6
-#define ON_MIS_RATE 0.04
+#define ON_MIS_RATE 0.06
 #define ON_MISMATCH_THD 10 //XXX
 
 #define ON_MAT_SCORE 1