@transform( getChipseqOverlap, suffix(".chipseq"), ".chipseq.load")
def loadChipseqIntervals(infile, outfile):
    track = P.snip( os.path.basename( infile ), ".chipseq" ).replace(".","_").replace("-","_")
    header = "track,overlap"
    statement = '''cat %(infile)s | python %(scriptsdir)s/csv2db.py
                      --database=%(database)s
                      --table=%(track)s_chipseq
                      --header=%(header)s
                      --allow-empty
                   > %(outfile)s '''
    P.run()

############################################################
## Compare CAPseq intervals with external CAPseq intervals
Fetch newsgroups posts, using the given nntp object s.
dataset/ETHPy150Open CGATOxford/cgat/obsolete/pipeline_proj007.py/loadChipseqIntervals--crossover with--dataset/ETHPy150Open GoogleCloudPlatform/appengine-sockets-python-java-go/python_socket_demo/main.py/fetch_posts