bool canR = (samread(sam_in, b) >= 0);
if (!canR) return -1;
- if (b->core.flag & 0x0001) { fprintf(stderr, "Find a paired end read in the file!\n"); exit(-1); }
- //(b->core.flag & 0x0100) && && !(b->core.flag & 0x0004)
+ general_assert(!(b->core.flag & 0x0001), "Find a paired end read in the file!");
int readType = getReadType(b);
std::string name = getName(b);
bool canR = (samread(sam_in, b) >= 0);
if (!canR) return -1;
- if (b->core.flag & 0x0001) { fprintf(stderr, "Find a paired end read in the file!\n"); exit(-1); }
- //assert(!(b->core.flag & 0x0001)); //(b->core.flag & 0x0100) && && !(b->core.flag & 0x0004)
+ general_assert(!(b->core.flag & 0x0001), "Find a paired end read in the file!");
int readType = getReadType(b);
std::string name = getName(b);
bool canR = ((samread(sam_in, b) >= 0) && (samread(sam_in, b2) >= 0));
if (!canR) return -1;
- if (!((b->core.flag & 0x0001) && (b2->core.flag & 0x0001))) {
- fprintf(stderr, "One of the mate is not paired-end! (RSEM assumes the two mates of a paired-end read should be adjacent)\n");
- exit(-1);
- }
- //assert((b->core.flag & 0x0001) && (b2->core.flag & 0x0001));
+ general_assert((b->core.flag & 0x0001) && (b2->core.flag & 0x0001), \
+ "One of the mate is not paired-end! (RSEM assumes the two mates of a paired-end read should be adjacent)");
bam1_t *mp1 = NULL, *mp2 = NULL;
bool canR = ((samread(sam_in, b) >= 0) && (samread(sam_in, b2) >= 0));
if (!canR) return -1;
- if (!((b->core.flag & 0x0001) && (b2->core.flag & 0x0001))) {
- fprintf(stderr, "One of the mate is not paired-end! (RSEM assumes the two mates of a paired-end read should be adjacent)\n");
- exit(-1);
- }
- //assert((b->core.flag & 0x0001) && (b2->core.flag & 0x0001));
+ general_assert((b->core.flag & 0x0001) && (b2->core.flag & 0x0001), \
+ "One of the mate is not paired-end! (RSEM assumes the two mates of a paired-end read should be adjacent)");
bam1_t *mp1 = NULL, *mp2 = NULL;
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.null(vec)) exit_with_error(paste("Unknown transcript detected,", tids[i], "is not included in RSEM's indices."))
if (is.na(vec[[2]])) wiggle <- rep(0, vec[[1]]) else wiggle <- as.numeric(unlist(strsplit(vec[[2]], split = " ")))
len <- length(wiggle)
if (!show_uniq) {
}
} else {
for (gene_id in ids) {
- if (is.null(t2gmap[[gene_id]])) exit_with_error(paste("Cannot find gene", gene_id, sep = ""))
+ if (is.null(t2gmap[[gene_id]])) exit_with_error(paste("Unknown gene detected,", gene_id, "is not included in RSEM's in indices."))
generate_a_page(t2gmap[[gene_id]], gene_id)
}
}
cat("Plots are generated!\n)
dev.off.output <- dev.off()
-
-
-
-
-
-