Skip to content
Snippets Groups Projects
Commit 9628eca3 authored by Aaron's avatar Aaron
Browse files

Tweaked Makefile and added notboth option to pairToPair.

parent 82617370
No related branches found
No related tags found
No related merge requests found
......@@ -21,10 +21,8 @@ SUBDIRS = $(SRC_DIR)/annotateBed \
$(SRC_DIR)/closestBed \
$(SRC_DIR)/complementBed \
$(SRC_DIR)/coverageBed \
$(SRC_DIR)/cuffToTrans \
$(SRC_DIR)/fastaFromBed \
$(SRC_DIR)/flankBed \
$(SRC_DIR)/fjoin \
$(SRC_DIR)/genomeCoverageBed \
$(SRC_DIR)/intersectBed \
$(SRC_DIR)/linksBed \
......
......@@ -124,7 +124,7 @@ void PairToPair::FindOverlaps(const BEDPE &a) {
if ((hitsA2B1.size() > 0) || (hitsA1B2.size() > 0))
found2 = FindHitsOnBothEnds(a, hitsA2B1, hitsA1B2);
if (!found1 && !found2)
if (!found1 || !found2)
_bedA->reportBedPENewLine(a);
}
else if (_searchType == "either") {
......
......@@ -65,7 +65,7 @@ private:
void FindQualityHitsBetweenEnds(CHRPOS start, CHRPOS end,
const vector<MATE> &hits, vector<MATE> &qualityHits, int &numOverlaps);
void FindHitsOnBothEnds(const BEDPE &a, const vector<MATE> &qualityHitsEnd1,
bool FindHitsOnBothEnds(const BEDPE &a, const vector<MATE> &qualityHitsEnd1,
const vector<MATE> &qualityHitsEnd2);
void FindHitsOnEitherEnd(const BEDPE &a, const vector<MATE> &qualityHitsEnd1,
......
......@@ -122,8 +122,8 @@ int main(int argc, char* argv[]) {
showHelp = true;
}
if (haveSearchType && (searchType != "neither") && (searchType != "both") && (searchType != "either")) {
cerr << endl << "*****" << endl << "*****ERROR: Request \"both\" or \"neither\"" << endl << "*****" << endl;
if (haveSearchType && (searchType != "neither") && (searchType != "both") && (searchType != "either") && (searchType != "notboth")) {
cerr << endl << "*****" << endl << "*****ERROR: Request \"both\",\"neither\",\"either\",or \"notboth\"" << endl << "*****" << endl;
showHelp = true;
}
......@@ -160,11 +160,11 @@ void ShowHelp(void) {
cerr << "\t-type \t" << "Approach to reporting overlaps between A and B." << endl << endl;
cerr << "\t\tneither\tReport overlaps if neither end of A overlaps B." << endl;
cerr << "\t\teither\tReport overlaps if either ends of A overlap B." << endl;
cerr << "\t\tboth\tReport overlaps if both ends of A overlap B." << endl;
cerr << "\t\t\t- Default = both." << endl << endl;
cerr << "\t\tnotboth\tReport overlaps if one or neither of ends of A overlap B." << endl;
cerr << "\t\t- Default = both." << endl << endl;
cerr << "\t-slop \t" << "The amount of slop (in b.p.). to be added to each footprint." << endl;
cerr << "\t\t*Note*: Slop is subtracted from start1 and start2 and added to end1 and end2." << endl << endl;
......
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