minor fix

Author: mfiers

hagfish_coverage_combine

@@ -59,6 +59,9 @@ if __name__ == '__main__':
 
     no_seqs_without_coverage = 0
     for seqId in seqInfo.keys():
+        l.debug("processig seqid %s" % seqId)
+
+        seqLen = seqInfo[seqId]['length']
 
         r_ok = None
         r_low = None
@@ -109,10 +112,8 @@ if __name__ == '__main__':
                 r_low_ends += d_low_ends
 
         if r_ok != None:
-            l.info("Processed %s ln %d min/max for ok:%d/%d low:%d/%d high: %d/%d" % (
-                    bamBase, len(d_ok), max(d_ok), min(d_ok),
-                    max(d_low), min(d_low), max(d_high),
-                    min(d_high)))
+            l.info("Processed %s ln %d max for ok: %d low: %d high: %d" % (
+                    bamBase, len(d_ok), max(d_ok), max(d_low), max(d_high)))
             np_savez(os.path.join('combined', seqId),
                      r_ok = r_ok,
                      r_high = r_high,
@@ -122,7 +123,7 @@ if __name__ == '__main__':
                      r_high_ends = r_high_ends,
                      r_low_ends = r_low_ends)
         else:
-            l.info("PROBLEM - FOUND NOTHING FOR %s" % bamBase)
+            l.info("PROBLEM - FOUND NOTHING FOR %s/%s" % (bamBase,seqId))
             no_seqs_without_coverage += 1
 
     l.info("ignored %d sequences (no coverage??)" % no_seqs_without_coverage)