Mercurial > repos > mvdbeek > dapars

changeset 4:73b932244237 draft

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
planemo upload for repository https://github.com/mvdbeek/dapars commit c2344ee1b8c6371007b7484b02d17fb056eb2427-dirty
author mvdbeek
date Wed, 28 Oct 2015 06:22:18 -0400
parents aca85eb3eb5b
children a5d8b08af089
files dapars.py
diffstat 1 files changed, 5 insertions(+), 1 deletions(-) [+]
line wrap: on
line diff
--- a/dapars.py	Wed Oct 28 06:05:21 2015 -0400
+++ b/dapars.py	Wed Oct 28 06:22:18 2015 -0400
@@ -185,7 +185,7 @@
 
     def write_bed(self):
         w = csv.writer(self.bed_output, delimiter='\t')
-        bed = [(result.chr, result.breakpoint, result.breakpoint+1, result.gene+"_"+result.breakpoint_type, 0, result.strand) for result in self.result_d.itervalues()]
+        bed = [(result.chr, result.breakpoint, int(result.breakpoint)+1, result.gene+"_"+result.breakpoint_type, 0, result.strand) for result in self.result_d.itervalues()]
         w.writerows(bed)
 
 
@@ -208,6 +208,10 @@
     if treatment_breakpoint:
         breakpoint_to_result(res_treatment, utr, utr_d, treatment_breakpoint, "treatment_breakpoint", treatment_abundance, is_reverse,
                              num_samples, num_control, num_treatment)
+    if res_control == dict(zip(result_tuple_fields, result_tuple_fields)):
+        res_control = False
+    if res_treatment == dict(zip(result_tuple_fields, result_tuple_fields)):
+        res_treatment == False
     return res_control, res_treatment