Mercurial > repos > alenail > chipsequtil

comparison chipsequtil/map_to_known_genes.py @ 16:4c645d02a428 draft

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Uploaded
author alenail
date Mon, 04 Apr 2016 16:20:07 -0400
parents 425c48c48dde
children 435adb6000f1
comparison
equal deleted inserted replaced
15:d7e08e5fa1a7 16:4c645d02a428
118 peak_loc = (peak[start_field]+peak[end_field])/2 118 peak_loc = (peak[start_field]+peak[end_field])/2
119 119
120 chrom_genes = gene_ref[peak[chr_field]] 120 chrom_genes = gene_ref[peak[chr_field]]
121 121
122 if len(chrom_genes) == 0 : 122 if len(chrom_genes) == 0 :
123 sys.stderr.write('WARNING: peak chromosome %s not found in gene reference, skipping: %s\n'%(peak[chr_field],peak)) 123 sys.stdout.write('WARNING: peak chromosome %s not found in gene reference, skipping: %s\n'%(peak[chr_field],peak))
124 continue 124 continue
125 125
126 mapped = False 126 mapped = False
127 127
128 # walk through the genes for this chromosome 128 # walk through the genes for this chromosome