Skip to content

Commit dbf065b

Browse files
committed
fix the regression issue of clustering reversed UMI
1 parent 86af0cf commit dbf065b

File tree

2 files changed

+2
-46
lines changed

2 files changed

+2
-46
lines changed

src/cluster.cpp

Lines changed: 1 addition & 45 deletions
Original file line numberDiff line numberDiff line change
@@ -48,51 +48,7 @@ int Cluster::umiDiff(const string& umi1, const string& umi2) {
4848
diff++;
4949
}
5050

51-
if( diff == 0)
52-
return 0;
53-
54-
int underline1 = -1;
55-
int underline2 = -1;
56-
57-
for(int i=0; i<len1; i++) {
58-
if(umi1[i] == '_') {
59-
underline1 = i;
60-
break;
61-
}
62-
}
63-
64-
if(underline1 <= 0)
65-
return diff;
66-
67-
for(int i=0; i<len2; i++) {
68-
if(umi2[i] == '_') {
69-
underline2 = i;
70-
break;
71-
}
72-
}
73-
74-
if(underline2 <= 0)
75-
return diff;
76-
77-
int len11 = underline1;
78-
int len12 = len1 - underline1 - 1;
79-
int len21 = underline2;
80-
int len22 = len2 - underline2 - 1;
81-
82-
// reversed
83-
int d1 = abs(len11 - len22);
84-
for(int i=0; i<min(len11, len22); i++) {
85-
if(umi1[i] != umi2[underline2 + i + 1])
86-
d1++;
87-
}
88-
int d2 = abs(len12 - len21);
89-
for(int i=0; i<min(len12, len21); i++) {
90-
if(umi1[underline1 + i + 1] != umi2[i])
91-
d2++;
92-
}
93-
int revDiff = d1 + d2;
94-
95-
return min(diff, revDiff);
51+
return diff;
9652
}
9753

9854
vector<Pair*> Cluster::clusterByUMI(int umiDiffThreshold, Stats* preStats, Stats* postStats) {

src/common.h

Lines changed: 1 addition & 1 deletion
Original file line numberDiff line numberDiff line change
@@ -1,7 +1,7 @@
11
#ifndef COMMON_H
22
#define COMMON_H
33

4-
#define VERSION_NUMBER "0.10.0"
4+
#define VERSION_NUMBER "0.11.0"
55

66
#define _DEBUG false
77

0 commit comments

Comments
 (0)