Advertisement
Not a member of Pastebin yet?
Sign Up,
it unlocks many cool features!
- 7 303 0.207756232686981
- 16 23 0.208562019758507
- 6 57 0.208727272727273
- 7 80 0.209065354884048
- 11 124 0.209500609013398
- chr7 303 292526
- chr16 23 169805
- chr6 57 62822
- chr11 124 320564
- chr7 80 300291
- 6 57 0.208727272727273
- for line in inputfile:
- line = line.rstrip()
- fields = line.split("t")
- #!/usr/bin/env python
- f2 = open( '/mnt/genotyping/CT/GreatApes/HKA/callability/callable_sites_per_region_500Kb.txt', 'r')
- d2 = {}
- print f2
- for line in f2:
- line = line.rstrip()
- fields = line.split("t")
- key = (fields[0].replace('chr', ''), fields[1])
- d2[key] = int(fields[2])
- f1 = open( '/mnt/genotyping/CT/GreatApes/HKA/Barcelona_approach/500kb/cov_5/Homo-Gorilla/R_plots/Gorilla_genome_dist_cov5.txt', 'r')
- for line in f1:
- line = line.rstrip()
- fields = line.split("t")
- if 'region' not in line:
- key = (fields[0], fields[1])
- if d2[key] >= 100000:
- print line
- >>> f2 = """chr7t303t292526
- chr16t23t169805
- chr6t57t62822
- chr11t124t320564
- chr7t80t300291"""
- >>> d2 = {}
- >>> for line in f2.split('n'):
- line = line.rstrip()
- fields = line.split("t")
- key = (fields[0].replace('chr', ''), fields[1])
- d2[key] = int(fields[2])
- >>> d2
- {('7', '303'): 292526, ('7', '80'): 300291, ('16', '23'): 169805, ('6', '57'): 62822, ('11', '124'): 320564}
- >>> f1 = """7t303t0.207756232686981
- 16t23t0.208562019758507
- 6t57t0.208727272727273
- 7t80t0.209065354884048
- 11t124t0.209500609013398"""
- >>> for line in f1.split('n'):
- line = line.rstrip()
- fields = line.split("t")
- key = (fields[0], fields[1])
- if d2[key] >= 100000:
- print line
- 7 303 0.207756232686981
- 16 23 0.208562019758507
- 7 80 0.209065354884048
- 11 124 0.209500609013398
- >>>
Advertisement
Add Comment
Please, Sign In to add comment
Advertisement