X-Git-Url: https://git.donarmstrong.com/?a=blobdiff_plain;f=BamConverter.h;h=af984ac4968c24a6a0c297154374bdaa32dbcd9e;hb=4c51136f0b5ae7758258642e1183de5bfb4624af;hp=e7253ba9415b566442f3efbfe266ffe052bb2753;hpb=635ca2939cfb1f519f19e9dec072ddd05e9fb450;p=rsem.git diff --git a/BamConverter.h b/BamConverter.h index e7253ba..af984ac 100644 --- a/BamConverter.h +++ b/BamConverter.h @@ -50,6 +50,8 @@ BamConverter::BamConverter(const char* inpF, const char* outF, const char* chr_l in = samopen(inpF, "rb", NULL); assert(in != 0); + transcripts.buildMappings(in->header->n_targets, in->header->target_name); + bam_header_t *out_header = sam_header_read2(chr_list); refmap.clear(); for (int i = 0; i < out_header->n_targets; i++) { @@ -93,7 +95,7 @@ void BamConverter::process() { // at least one segment is not properly mapped if ((b->core.flag & 0x0004) || (isPaired && (b2->core.flag & 0x0004))) continue; - const Transcript& transcript = transcripts.getTranscriptAt(b->core.tid + 1); + const Transcript& transcript = transcripts.getTranscriptViaEid(b->core.tid + 1); convert(b, transcript); if (isPaired) {