X-Git-Url: http://woldlab.caltech.edu/gitweb/?p=erange.git;a=blobdiff_plain;f=siteintersects.py;fp=siteintersects.py;h=58d9e3059ad8a587c10badd3fb8c7cc820462047;hp=e012bb49436c92f50a1d8722182d0cd72a754897;hb=cfc5602b26323ad2365295145e3f6c622d912eb4;hpb=c4561c55cfa9726530c6777b6515c4ef66306b2f diff --git a/siteintersects.py b/siteintersects.py index e012bb4..58d9e30 100755 --- a/siteintersects.py +++ b/siteintersects.py @@ -126,6 +126,8 @@ def compareSites(siteFilename, rejectFilename, outfilename, siteDict, rejectSite continue siteNotCommon = True + #TODO: as handed to me this is O(n^m) - probably could be better... + #TODO: QForAli - if there are multiple matches will count them but only write the first match to the file for (rstart, rstop, rline) in siteDict[chrom]: if regionsOverlap(start, stop, rstart, rstop): commonSites += 1 @@ -162,4 +164,4 @@ def writeRejectSiteFile(siteRejectFilename, rejectSiteDict): if __name__ == "__main__": - main(sys.argv) + main(sys.argv) \ No newline at end of file