@follows(mkdir("taxonomy.dir"))
@merge(GENOMES, r"taxonomy.dir/gi_accessions.tsv")
def buildGiAccessionNumbers(infiles, outfile):
    # first line in each file is the contig name
    outf = open(outfile, "w")
    for inf in infiles:
        outf.write(open(inf).readline().split("|")[1] + "\n")
    outf.close()


######################################################
######################################################
######################################################
# AH: added default for taxonomy_taxdir
build a list of gi accession numbers from the genomes
dataset/ETHPy150Open CGATOxford/cgat/obsolete/pipeline_metagenomebenchmark.py/buildGiAccessionNumbers