@follows(indexGenome,loadExpressionData)
@transform( buildPileups, suffix(".pileup.gz"), ".validated.gz" )
def createSNPValidationData( infile, outfile ):

    track = infile[:-len(".pileup.gz")]
    strain = track[len("mouse"):]
    
    pattern = "%s/*%s*.bam" % (PARAMS["rnaseq_data"], strain)
    
    to_cluster = True

    statement = '''gunzip
        < %(infile)s
        | python %(scriptsdir)s/snp2snp.py 
              --method=validate
              --filename-reference="%(pattern)s"
              --min-coverage=%(rnaseq_min_coverage)i
              --filename-genome=%(genome)s.fa
              --log=%(outfile)s.log
        | gzip
        > %(outfile)s'''

    P.run()
Remove node n.
dataset/ETHPy150Open CGATOxford/cgat/obsolete/pipeline_snps.py/createSNPValidationData--crossover with--dataset/ETHPy150Open gkno/gkno_launcher/src/networkx/classes/graph.py/Graph.remove_node