diff dante.py @ 15:3151a72a6671 draft

Uploaded
author petr-novak
date Tue, 03 Sep 2019 05:20:02 -0400
parents a6c55d1bdb6c
children 039c45c01b47
line wrap: on
line diff
--- a/dante.py	Wed Aug 28 08:08:47 2019 -0400
+++ b/dante.py	Tue Sep 03 05:20:02 2019 -0400
@@ -586,10 +586,10 @@
             if count_region == len(indices_plus):
                 strand_gff = "-"
             if strand_gff == "+":
-                feature_start = min(start_hit[regions_above_threshold])-1
+                feature_start = min(start_hit[regions_above_threshold]) + 1
                 feature_end = max(end_hit[regions_above_threshold])
             else:
-                feature_end = seq_len[region][0] - min(start_hit[regions_above_threshold]) - 1
+                feature_end = seq_len[region][0] - min(start_hit[regions_above_threshold])
                 feature_start = seq_len[region][0] - max(end_hit[regions_above_threshold]) + 1
             create_gff3(domain_type, ann_substring, unique_annotations,
                         ann_pos_counts, feature_start,feature_end,