diff --git a/small_RNA-seq/small_RNA-seq.snakefile b/small_RNA-seq/small_RNA-seq.snakefile index c4616e3ef9fae16f7eab187864cbc4320ba8160e..107e1eda7b8c7bbe4ee3d17c7b80652a6f1e29c3 100644 --- a/small_RNA-seq/small_RNA-seq.snakefile +++ b/small_RNA-seq/small_RNA-seq.snakefile @@ -2788,6 +2788,12 @@ rule merge_bigwig_reps: warn(f"{msg}:\n") #raise warn(f"The bigwig files without {chrom} will be skipped.\n") + except AssertionError as e: + chrom_lens = [] + for bw in bws: + chrom_lens.append(bw.chroms()[chrom]) + warn(f"chrom_len for {chrom}: {chrom_len}\nchrom_lens: {chrom_lens}") + raise to_use = [bw for bw in bws if chrom in bw.chroms()] if to_use: means = np.nanmean(np.vstack([bw.values(chrom, 0, chrom_len, numpy=True) for bw in to_use]), axis=0)