/scripts/microbes/create_bacteria_table.py

https://bitbucket.org/cistrome/cistrome-harvard/ · Python · 79 lines · 60 code · 14 blank · 5 comment · 22 complexity · b263c164cd15bffd375bdffbed7f014b MD5 · raw file

  1. #!/usr/bin/env python
  2. #Dan Blankenberg
  3. import sys, os
  4. assert sys.version_info[:2] >= ( 2, 4 )
  5. def __main__():
  6. base_dir = os.path.join( os.getcwd(), "bacteria" )
  7. try:
  8. base_dir = sys.argv[1]
  9. except:
  10. pass
  11. #print "using default base_dir:", base_dir
  12. organisms = {}
  13. for result in os.walk(base_dir):
  14. this_base_dir,sub_dirs,files = result
  15. for file in files:
  16. if file[-5:] == ".info":
  17. dict = {}
  18. info_file = open(os.path.join(this_base_dir,file),'r')
  19. info = info_file.readlines()
  20. info_file.close()
  21. for line in info:
  22. fields = line.replace("\n","").split("=")
  23. dict[fields[0]]="=".join(fields[1:])
  24. if 'genome project id' in dict.keys():
  25. name = dict['genome project id']
  26. if 'build' in dict.keys():
  27. name = dict['build']
  28. if name not in organisms.keys():
  29. organisms[name] = {'chrs':{},'base_dir':this_base_dir}
  30. for key in dict.keys():
  31. organisms[name][key]=dict[key]
  32. else:
  33. if dict['organism'] not in organisms.keys():
  34. organisms[dict['organism']] = {'chrs':{},'base_dir':this_base_dir}
  35. organisms[dict['organism']]['chrs'][dict['chromosome']]=dict
  36. orgs = organisms.keys()
  37. for org in orgs:
  38. if 'name' not in organisms[org]:del organisms[org]
  39. orgs = organisms.keys()
  40. #need to sort by name
  41. swap_test = False
  42. for i in range(0, len(orgs) - 1):
  43. for j in range(0, len(orgs) - i - 1):
  44. if organisms[orgs[j]]['name'] > organisms[orgs[j + 1]]['name']:
  45. orgs[j], orgs[j + 1] = orgs[j + 1], orgs[j]
  46. swap_test = True
  47. if swap_test == False:
  48. break
  49. print "||'''Organism'''||'''Kingdom'''||'''Group'''||'''Links to UCSC Archaea Browser'''||"
  50. for org in orgs:
  51. org = organisms[org]
  52. at_ucsc = False
  53. #if no gpi, then must be a ncbi chr which corresponds to a UCSC org, w/o matching UCSC designation
  54. try:
  55. build = org['genome project id']
  56. except: continue
  57. if 'build' in org:
  58. build = org['build']
  59. at_ucsc = True
  60. out_str = "||"+org['name']+"||"+org['kingdom']+"||"+org['group']+"||"
  61. if at_ucsc:
  62. out_str = out_str + "Yes"
  63. out_str = out_str + "||"
  64. print out_str
  65. if __name__ == "__main__": __main__()