better error with missing fasta on bed.gz files

brentp · Feb 2, 2021 · 91d92de · 91d92de
1 parent 5681ed6
commit 91d92de
Showing 1 changed file with 2 additions and 0 deletions.
diff --git a/src/seqcoverpkg/typeenum.nim b/src/seqcoverpkg/typeenum.nim
@@ -63,6 +63,8 @@ proc chromosomes*(c:var Cover, fai:Fai): OrderedTableRef[string, uint32] =
     result = c.d4.chromosomes
   of CoverType.BED:
     result = newOrderedTable[string, uint32]()
+    if fai == nil:
+      raise newException(OSError, "--fasta index required when .bed.gz files are given")
     for i in 0..<fai.len:
       var cname = fai[i]
       result[cname] = fai.chrom_len(cname).uint32