X-Git-Url: https://git.donarmstrong.com/?a=blobdiff_plain;f=extractRef.cpp;h=0a46de0950534ea4c22f3acb8f239de66ef8415f;hb=1b2999c4407ef8419fb89b66b843b7141caff4da;hp=3173d6517c1fccd6edcf03fa15502b64faadb28d;hpb=f7d64c1502c10b0de87415343e9d6bdac726f155;p=rsem.git diff --git a/extractRef.cpp b/extractRef.cpp index 3173d65..0a46de0 100644 --- a/extractRef.cpp +++ b/extractRef.cpp @@ -134,14 +134,14 @@ void parse_gtf_file(char* gtfF) { else { if (hasMappingFile) { tid = item.getTranscriptID(); - mi_iter = mi_table.find(tid); - if (mi_iter == mi_table.end()) { - fprintf(stderr, "Mapping Info is not correct, cannot find %s's gene_id!\n", tid.c_str()); - exit(-1); - } - //assert(iter != table.end()); - gid = mi_iter->second; - item.setGeneID(gid); + mi_iter = mi_table.find(tid); + if (mi_iter == mi_table.end()) { + fprintf(stderr, "Mapping Info is not correct, cannot find %s's gene_id!\n", tid.c_str()); + exit(-1); + } + //assert(iter != table.end()); + gid = mi_iter->second; + item.setGeneID(gid); } items.push_back(item); } @@ -304,8 +304,7 @@ int main(int argc, char* argv[]) { for (int i = 1; i <= M; i++) { if (seqs[i] == "") { const Transcript& transcript = transcripts.getTranscriptAt(i); - fprintf(stderr, "Cannot extract transcript %s's sequence from chromosome %s, whose information might not be provided! \" - "Please check if the chromosome directory is set correctly or the list of chromosome files is complete.\n", \ + fprintf(stderr, "Cannot extract transcript %s's sequence from chromosome %s, whose information might not be provided! Please check if the chromosome directory is set correctly or the list of chromosome files is complete.\n", \ transcript.getTranscriptID().c_str(), transcript.getGeneID().c_str()); exit(-1); }