@transform(copyCapseqReplicatedBedFiles, suffix(".bed"), ".lncrna.tss.distance" )
def getCapseqlncRNATSSDistance( infile, outfile ):
    to_cluster = True
    annotation_file = os.path.join( PARAMS["geneset_dir"], PARAMS["geneset_lncrna_tss"] )
    statement = """cat < %(infile)s 
                   | python %(scriptsdir)s/bed2gff.py --as-gtf 
                   | python %(scriptsdir)s/gtf2table.py 
		               --counter=distance-tss 
		               --log=%(outfile)s.log 
                       --filename-gff=%(annotation_file)s 
                       --filename-format="bed" 
                   > %(outfile)s"""
    P.run()

############################################################
Verify that the C{write} method of the transport of a protocol
dataset/ETHPy150Open CGATOxford/cgat/obsolete/pipeline_proj007.py/getCapseqlncRNATSSDistance--crossover with--dataset/ETHPy150Open kuri65536/python-for-android/python-modules/twisted/twisted/test/test_stdio.py/StandardInputOutputTestCase.test_write