Skip to content
Snippets Groups Projects
Commit 11aaf440 authored by Aaron's avatar Aaron
Browse files

[BUG] non-redundant cluster ids in stranded (-s) mode.

parent 2eeee4c3
No related branches found
No related tags found
No related merge requests found
......@@ -72,6 +72,8 @@ void BedCluster::ClusterBedStranded() {
// that we can easily compare "A" to it for overlaps
_bed->loadBedFileIntoMapNoBin();
uint32_t cluster_id = 0;
// loop through each chromosome and merge their BED entries
masterBedMapNoBin::const_iterator m = _bed->bedMapNoBin.begin();
masterBedMapNoBin::const_iterator mEnd = _bed->bedMapNoBin.end();
......@@ -85,7 +87,6 @@ void BedCluster::ClusterBedStranded() {
vector<string> strands(2);
strands[0] = "+";
strands[1] = "-";
uint32_t cluster_id = 0;
// do two passes, one for each strand.
for (unsigned int s = 0; s < strands.size(); s++) {
// cluster overlapping features for this chromosome.
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment