@@ -421,9 +421,7 @@ def coverage_from_array(valuesArray, zones, binSize, avgType):
421421 index_start = int (pos - start )
422422 index_end = int (index_start + step_size )
423423 try :
424- counts_list .append (
425- heatmapper .my_average (valuesArray [index_start :index_end ],
426- avgType ))
424+ counts_list .append (heatmapper .my_average (valuesArray [index_start :index_end ], avgType ))
427425 except Exception as detail :
428426 sys .stderr .write ("Exception found. "
429427 "Message: {}\n " .format (detail ))
@@ -488,12 +486,12 @@ def coverage_from_bam(bamfile, chrom, zones, binSize, avgType, verbose=True):
488486 def coverage_from_big_wig (bigwig , chrom , zones , binSize , avgType , nansAsZeros = False , verbose = True ):
489487
490488 """
491- uses bigwig file reader from bx-python
492- to query a region define by chrom and zones.
489+ uses bigwig file reader from pyBigWib
490+ to query a region defined by chrom and zones.
493491 The output is an array that contains the bigwig
494492 value per base pair. The summary over bins is
495493 done in a later step when coverage_from_array is called.
496- This method is more reliable than quering the bins
494+ This method is more reliable than querying the bins
497495 directly from the bigwig, which should be more efficient.
498496
499497 By default, any region, even if no chromosome match is found
0 commit comments