@transform( TRACKS + [exportIntrons,], 
            suffix(".gtf.gz"), 
            ".fasta")
def exportSequences( infile, outfile ):
    
    prefix = outfile[:-len(".fasta")]

    to_cluster = True
    statement = '''gunzip 
        < %(infile)s
	| %(scriptsdir)s/gff_sort gene  
	| python %(scriptsdir)s/gff2fasta.py 
		--is-gtf 
		--genome-file=genome 
		--log=%(outfile)s.log 
	| python %(toolsdir)s/index_fasta.py --force %(prefix)s - 
        > %(outfile)s.log'''

    P.run()

############################################################
Marks notifications as either read or unread depending of POST parameters.
dataset/ETHPy150Open CGATOxford/cgat/obsolete/pipeline_kamilah.py/exportSequences--crossover with--dataset/ETHPy150Open v1k45/django-notify-x/notify/views.py/mark_all