@follows( liverTestesVenn )
@files( "liver_vs_testes/liver.testes.merge.bed", "liver_vs_testes/liver.testes.merge.transcript.tss.distance" )
def annotateLiverTestesMergedTranscriptTSSDistance( infile, outfile ):
    to_cluster = True
    annotation_file = os.path.join( PARAMS["geneset_dir"],PARAMS["geneset_transcript_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()

############################################################
Compute distance from CAPseq intervals to nearest transcript TSS
dataset/ETHPy150Open CGATOxford/cgat/obsolete/pipeline_proj007.py/annotateLiverTestesMergedTranscriptTSSDistance