@files( [( [(x,y) % "%s%s.norm.bam" for y in REPLICATES], "%s.readcorrelations.gz" % x) for x in TRACKS_EXPERIMENTS ] )
def makeReadCorrelation( infiles, outfile ):
    '''compute correlation between reads.
    '''

    to_cluster = True
    job_options = "-l mem_free=8000M"

    infiles = " ".join(infiles)

    statement = '''
    python %(scriptsdir)s/bam_correlation.py \
    --log=%(outfile)s.log \
    --genome=%(genome)s \
    %(infiles)s | gzip > %(outfile)s
    ''' 

    P.run()

############################################################
############################################################
############################################################
dataset/ETHPy150Open CGATOxford/cgat/obsolete/pipeline_vitaminD.py/makeReadCorrelation/swapped operands of `%`