Advertisement
Guest User

Untitled

a guest
Feb 20th, 2019
71
0
Never
Not a member of Pastebin yet? Sign Up, it unlocks many cool features!
text 1.74 KB | None | 0 0
  1. 7 303 0.207756232686981
  2. 16 23 0.208562019758507
  3. 6 57 0.208727272727273
  4. 7 80 0.209065354884048
  5. 11 124 0.209500609013398
  6.  
  7. chr7 303 292526
  8. chr16 23 169805
  9. chr6 57 62822
  10. chr11 124 320564
  11. chr7 80 300291
  12.  
  13. 6 57 0.208727272727273
  14.  
  15. for line in inputfile:
  16. line = line.rstrip()
  17. fields = line.split("t")
  18.  
  19. #!/usr/bin/env python
  20.  
  21.  
  22.  
  23. f2 = open( '/mnt/genotyping/CT/GreatApes/HKA/callability/callable_sites_per_region_500Kb.txt', 'r')
  24. d2 = {}
  25. print f2
  26. for line in f2:
  27. line = line.rstrip()
  28. fields = line.split("t")
  29. key = (fields[0].replace('chr', ''), fields[1])
  30. d2[key] = int(fields[2])
  31.  
  32.  
  33.  
  34.  
  35.  
  36. f1 = open( '/mnt/genotyping/CT/GreatApes/HKA/Barcelona_approach/500kb/cov_5/Homo-Gorilla/R_plots/Gorilla_genome_dist_cov5.txt', 'r')
  37. for line in f1:
  38. line = line.rstrip()
  39. fields = line.split("t")
  40. if 'region' not in line:
  41. key = (fields[0], fields[1])
  42. if d2[key] >= 100000:
  43. print line
  44.  
  45. >>> f2 = """chr7t303t292526
  46. chr16t23t169805
  47. chr6t57t62822
  48. chr11t124t320564
  49. chr7t80t300291"""
  50. >>> d2 = {}
  51. >>> for line in f2.split('n'):
  52. line = line.rstrip()
  53. fields = line.split("t")
  54. key = (fields[0].replace('chr', ''), fields[1])
  55. d2[key] = int(fields[2])
  56.  
  57.  
  58. >>> d2
  59. {('7', '303'): 292526, ('7', '80'): 300291, ('16', '23'): 169805, ('6', '57'): 62822, ('11', '124'): 320564}
  60.  
  61. >>> f1 = """7t303t0.207756232686981
  62. 16t23t0.208562019758507
  63. 6t57t0.208727272727273
  64. 7t80t0.209065354884048
  65. 11t124t0.209500609013398"""
  66. >>> for line in f1.split('n'):
  67. line = line.rstrip()
  68. fields = line.split("t")
  69. key = (fields[0], fields[1])
  70. if d2[key] >= 100000:
  71. print line
  72.  
  73.  
  74. 7 303 0.207756232686981
  75. 16 23 0.208562019758507
  76. 7 80 0.209065354884048
  77. 11 124 0.209500609013398
  78. >>>
Advertisement
Add Comment
Please, Sign In to add comment
Advertisement