if (is_gene) par(oma = c(0, 0, 3, 0))
for (i in 1:n) {
vec <- readdepth[[tids[i]]]
if (is.null(vec)) exit_with_error(paste("Cannot find transcript", tids[i], sep = ""))
if (is_gene) par(oma = c(0, 0, 3, 0))
for (i in 1:n) {
vec <- readdepth[[tids[i]]]
if (is.null(vec)) exit_with_error(paste("Cannot find transcript", tids[i], sep = ""))