Skip to content

Commit 8558e92

Browse files
committed
fixed bug of PCRdup removal of nonexistent unpaired alns
1 parent 308da50 commit 8558e92

File tree

2 files changed

+5
-5
lines changed

2 files changed

+5
-5
lines changed

Genrich.c

Lines changed: 3 additions & 3 deletions
Original file line numberDiff line numberDiff line change
@@ -4,7 +4,7 @@
44
55
Finding sites of enrichment from genome-wide assays.
66
7-
Version 0.6
7+
Version 0.6.1
88
*/
99

1010
#include <stdio.h>
@@ -3556,7 +3556,7 @@ void addHashPr(Read* r, HashAln** table,
35563556
0, 0, table, idx, dupsVerb ? r->name : NULL);
35573557

35583558
// also add both alignments as singletons to hashtable
3559-
if (tableSn != NULL) {
3559+
if (tableSn != NULL && hashSizeSn) {
35603560
checkAndAdd(tableSn, hashSizeSn, a->chrom, a->pos[0],
35613561
true, dupsVerb ? r->name : NULL);
35623562
checkAndAdd(tableSn, hashSizeSn, a->chrom, a->pos[1],
@@ -3680,7 +3680,7 @@ void addHashDc(Read* r, HashAln** table,
36803680
b->strand, table, idx, dupsVerb ? r->name : NULL);
36813681

36823682
// also add both alignments as singletons to hashtable
3683-
if (tableSn != NULL) {
3683+
if (tableSn != NULL && hashSizeSn) {
36843684
if (! j)
36853685
checkAndAdd(tableSn, hashSizeSn, a->chrom, pos,
36863686
a->strand, dupsVerb ? r->name : NULL);

Genrich.h

Lines changed: 2 additions & 2 deletions
Original file line numberDiff line numberDiff line change
@@ -4,9 +4,9 @@
44
55
Finding sites of enrichment from genome-wide assays.
66
7-
Version 0.6
7+
Version 0.6.1
88
*/
9-
#define VERSION "0.6"
9+
#define VERSION "0.6.1"
1010

1111
// macros
1212
#define MAX(a, b) ((a) > (b) ? (a) : (b))

0 commit comments

Comments
 (0)