Mercurial > repos > gregory-minevich > ems_variant_density_mapping
comparison EMS_VariantDensityMapping.py @ 4:d2158e476320 draft
Uploaded
| author | gregory-minevich |
|---|---|
| date | Thu, 14 Jun 2012 21:22:09 -0400 |
| parents | a43cb9a57a9a |
| children |
comparison
equal
deleted
inserted
replaced
| 3:2e236e914815 | 4:d2158e476320 |
|---|---|
| 1 #!/usr/bin/python | 1 #!/usr/bin/python |
| 2 | 2 |
| 3 import re | |
| 3 import sys | 4 import sys |
| 4 import optparse | 5 import optparse |
| 5 import csv | 6 import csv |
| 6 from rpy import * | 7 from rpy import * |
| 7 | 8 |
| 53 iv_position_list = [] | 54 iv_position_list = [] |
| 54 v_position_list = [] | 55 v_position_list = [] |
| 55 x_position_list = [] | 56 x_position_list = [] |
| 56 | 57 |
| 57 for row in reader: | 58 for row in reader: |
| 58 chromosome = row[0] | 59 chromosome = row[0].upper() |
| 60 chromosome = re.sub("chr", "", chromosome, flags = re.IGNORECASE) | |
| 61 chromosome = re.sub("CHROMOSOME_", "", chromosome, flags = re.IGNORECASE) | |
| 62 | |
| 59 position = row[1] | 63 position = row[1] |
| 60 ref_allele = row[3] | 64 ref_allele = row[3] |
| 61 alt_allele = row[4] | 65 alt_allele = row[4] |
| 62 | 66 |
| 63 if (ems=='true'): | 67 if (ems=='true'): |
