src/obitools/utils/bioseq.py

@@ -7,11 +7,15 @@ def uniqSequence(seqIterator):
         s = str(seq)
         if s in uniques:
             s = uniques[s]
-            s['count']+=1
+            if 'count' in seq:
+                s['count']+=seq['count']
+            else:
+                s['count']+=1
             s['merged'].append(seq.id)
         else:
             uniques[s]=seq
-            seq['count']=1
+            if 'count' not in seq:
+                seq['count']=1
             seq['merged']=[seq.id]
             uniqSeq.append(seq)