Merge pull request #52 from lpantano/devel

fix bug in report making profile with long tail
lpantano · Mar 4, 2021 · e4ef358 · e4ef358
2 parents 38297c9 + 334197d
commit e4ef358
Show file tree

Hide file tree

Showing 3 changed files with 6 additions and 2 deletions.
diff --git a/HISTORY.md b/HISTORY.md
@@ -1,3 +1,7 @@
+- 1.2.8
+
+  * Fix bug when detecting the end of the cluster due to big gaps in biopython function
+
 - 1.2.7
 
   * Fix bug when writing files for debug of big meta-clusters: 

diff --git a/seqcluster/libs/read.py b/seqcluster/libs/read.py
@@ -103,7 +103,7 @@ def _align(x, y, local = False):
         sorted_alignments = sorted(aligned_x, key=operator.itemgetter(2))
         e = enumerate(sorted_alignments[0][0])
         nts = [i for i,c in e if c != "-"]
-        return [min(nts), max(nts)]
+        return [min(nts), min(nts)+len(x)]
 
 def map_to_precursor_biopython(seqs, names, loci, args):
     """map the sequences using biopython package"""

diff --git a/setup.py b/setup.py
@@ -13,7 +13,7 @@ def readme():
 
 
 setup(name='seqcluster',
-      version='1.2.7',
+      version='1.2.8',
       description='Small RNA-seq pipeline',
       long_description=readme(),
       long_description_content_type="text/markdown",