X-Git-Url: https://git.donarmstrong.com/?a=blobdiff_plain;f=BamWriter.h;h=a4e7dc523a3cb8f0acbf3d903dd021142df64eab;hb=1d786dca8b4c9c1e2176a28a669d2d16cd93e395;hp=8c568f348f361fd009c964cd2b9ad021e50fc257;hpb=7bbb1f6e5c87e6a953789c1f3ef49d3ce891e31e;p=rsem.git

diff --git a/BamWriter.h b/BamWriter.h
index 8c568f3..a4e7dc5 100644
--- a/BamWriter.h
+++ b/BamWriter.h
@@ -55,9 +55,6 @@ BamWriter::BamWriter(char inpType, const char* inpF, const char* fn_list, const
 
 	//generate output's header
 	bam_header_t *out_header = bam_header_dwt(in->header);
-	for (int i = 0; i < out_header->n_targets; i++) {
-		out_header->target_len[i] = transcripts.getTranscriptViaEid(i + 1).getLength();  // transcript length without poly(A) tail
-	}
 
 	std::ostringstream strout;
 	strout<<"@HD\tVN:1.4\tSO:unknown\n@PG\tID:RSEM\n";
@@ -116,10 +113,10 @@ void BamWriter::work(HitWrapper<PairedEndHit> wrapper) {
 
 	while (samread(in, b) >= 0 && samread(in, b2) >= 0) {
 		cnt += 2;
-		if (verbose && cnt % 1000000 == 0) { std::cout<< cnt<< "alignment lines are loaded!"<< std::endl; }
+		if (verbose && cnt % 1000000 == 0) { std::cout<< cnt<< " alignment lines are loaded!"<< std::endl; }
 
 		assert((b->core.flag & 0x0001) && (b2->core.flag & 0x0001));
-		assert((b->core.flag & 0x0040) && (b2->core.flag & 0x0080) || (b->core.flag & 0x0080) && (b2->core.flag & 0x0040));
+		assert(((b->core.flag & 0x0040) && (b2->core.flag & 0x0080)) || ((b->core.flag & 0x0080) && (b2->core.flag & 0x0040)));
 
 		//unalignable reads, skip		
 		bool notgood = (b->core.flag & 0x0004) || (b2->core.flag & 0x0004);
@@ -141,24 +138,8 @@ void BamWriter::work(HitWrapper<PairedEndHit> wrapper) {
 			convert(b, hit->getConPrb());
 			convert(b2, hit->getConPrb());
 
-			b->core.mpos = b2->core.pos;
-			b2->core.mpos = b->core.pos;
-		}
-
-		/*
-		else {
-		  // if only one mate can be aligned, mask it as unaligned and put an additional tag Z0:A:!
-		  char exclamation = '!';
-		  if (!(b->core.flag & 0x0004)) { 
-		    b->core.flag |= 0x0004;
-		    bam_aux_append(b, "Z0", 'A', bam_aux_type2size('A'), (uint8_t*)&exclamation);
-		  }
-		  if (!(b2->core.flag & 0x0004)) {
-		    b2->core.flag |= 0x0004;
-		    bam_aux_append(b2, "Z0", 'A', bam_aux_type2size('A'), (uint8_t*)&exclamation);
-		  }
+			// omit "b->core.mpos = b2->core.pos; b2->core.mpos = b->core.pos" for the reason that it is possible that only one mate is alignable 
 		}
-		*/
 
 		samwrite(out, b);
 		samwrite(out, b2);