# Report version p.print_version() # Internal flags fix_chromosome_name = True correct_stop_position = True # Track name, description and output file name track_name = os.path.splitext(os.path.basename(infile))[0] track_description = track_name outfile = os.path.splitext(os.path.basename(infile))[0] + ".bed" print "Output file: %s" % outfile # Read in data data = BedMaker(infile,column_names=('chr','start','stop','strand','transcript', 'fold_change','p_value')) # Fix chromosome name if fix_chromosome_name: print "Prepending 'chr' to chromosome names where it's missing" prependChromosomeName(data,'chr') # Subtract one from end position if correct_stop_position: print "Correcting 'stop' position by subtracting one base" adjustStopPosition(data) # Set name and RBG values data.computeColumn('name',lambda line: "%s_fc%s" % (line['transcript'],line['fold_change'])) data.computeColumn('RGB',lambda line: computeRGB(line['p_value']))
# Report version p.print_version() # Internal flags fix_chromosome_name = True correct_stop_position = True # Track name, description and output file name track_name = os.path.splitext(os.path.basename(infile))[0] track_description = track_name outfile = os.path.splitext(os.path.basename(infile))[0] + ".bed" print "Output file: %s" % outfile # Read in data data = BedMaker(infile,column_names=('chr','start','stop','sample_id','length', 'average_coverage')) # Fix chromosome name if fix_chromosome_name: print "Prepending 'chr' to chromosome names where it's missing" prependChromosomeName(data,'chr') # Subtract one from end position if correct_stop_position: print "Correcting 'stop' position by subtracting one base" adjustStopPosition(data) # Set name, strand, score and RBG values data.computeColumn('name',lambda line: "%s_%sbp" % (line['sample_id'],line['length'])) data.computeColumn('strand',lambda line: "+") data.computeColumn('score',lambda line: min(int(line['average_coverage']),1000))