--- a/dapars.py	Wed Oct 28 05:32:30 2015 -0400
+++ b/dapars.py	Wed Oct 28 05:56:52 2015 -0400
@@ -33,7 +33,7 @@
                         help="minimum coverage in each aligment to be considered for determining breakpoints")
     parser.add_argument("-b", "--breakpoint_bed", required=False, type=argparse.FileType('w'),
                         help="Write bedfile with coordinates of breakpoint positions to supplied path.")
-    parser.add_argument("-v", "--version", action='version', version='%(prog)s 0.1.4')
+    parser.add_argument("-v", "--version", action='version', version='%(prog)s 0.1.5')
     return parser.parse_args()
 
 
@@ -127,25 +127,6 @@
                     utr_dict[gene] = utr_d
         return utr_dict
 
-    def get_utr_coverage(self):
-        """
-        Returns a dict:
-        { UTR : [coverage_aligment1, ...]}
-        """
-        utr_coverages = {}
-        for utr, utr_d in self.utr_dict.iteritems():
-            if utr_d["chr"] in self.available_chromosomes:
-                if utr_d["strand"] == "+":
-                    is_reverse = False
-                else:
-                    is_reverse = True
-                utr_coverage = []
-                for bam in self.all_alignments:
-                    bp_coverage = get_bp_coverage(bam, utr_d["chr"], utr_d["new_start"], utr_d["new_end"], is_reverse)
-                    utr_coverage.append(bp_coverage)
-                utr_coverages[utr] = utr_coverage
-        return utr_coverages
-
     def get_coverage_weights(self):
         """
         Return weights for normalizing coverage.
@@ -230,8 +211,6 @@
     return res_control, res_treatment
 
 
-
-
 def breakpoint_to_result(res, utr, utr_d, breakpoint, breakpoint_type,
                          abundances, is_reverse, num_samples, num_control, num_treatment):
     """