Include the regionblocks even if the total genome is < 1 million

src/covcalc.rs

@@ -1323,6 +1323,9 @@ pub fn region_divider(regs: &Vec<Region>) -> Vec<Vec<Region>> {
             bplen = 0;
         }
     }
+    if tempregionvec.len() > 0 {
+        blocks.push(tempregionvec.clone());
+    }
     blocks
 }
 

src/multibamsummary.rs

@@ -136,6 +136,8 @@ pub fn r_mbams(
 
     // Divide up the regions into regionBlocks
     let regionblocks = region_divider(&regions);
+
+    assert!(regionblocks.len() > 0, "No regions to process. Exiting.");
     if verbose {
         println!("Regions divided into {} parallel blocks", regionblocks.len());
         println!("Start coverage calculation");