fix snp beg/end flank for longer reads

huishenlab · May 2, 2024 · 12dce78 · 12dce78
1 parent 93313ea
commit 12dce78
Showing 1 changed file with 2 additions and 2 deletions.
diff --git a/src/epiread.c b/src/epiread.c
@@ -544,8 +544,8 @@ static void *process_func(void *data) {
         rec.tid = w.tid;
         char *chrm = header->target_name[w.tid];
 
-        uint32_t snp_beg = w.beg>1000 ? w.beg-1000 : 1; // start location of snps
-        uint32_t snp_end = w.end+1000;
+        uint32_t snp_beg = w.beg > flank ? w.beg-flank : 1; // start location of snps
+        uint32_t snp_end = w.end+flank;
 
         // Make snp lookup table (if supplied)
         uint8_t *snps = NULL;