}
std::string getName(const bam1_t* b) {
- return std::string((char*)bam1_qname(b));
+ return std::string(bam1_qname(b));
}
std::string getReadSeq(const bam1_t*);
mp1 = b2; mp2 = b;
}
+ general_assert(!strcmp(bam1_qname(mp1), bam1_qname(mp2)), "Detected a read pair whose two mates have different names: " + getName(mp1) + " , " + getName(mp2) + " !");
+
int readType = getReadType(mp1, mp2);
std::string name = getName(mp1);
mp1 = b2; mp2 = b;
}
+ general_assert(!strcmp(bam1_qname(mp1), bam1_qname(mp2)), "Detected a read pair whose two mates have different names: " + getName(mp1) + " , " + getName(mp2) + " !");
+
int readType = getReadType(mp1, mp2);
std::string name = getName(mp1);
record_read.write(n_os, cat[record_val]);
++N[record_val];
}
+
+ general_assert(record_val == 1 || hits.getNHits() == 0, "Read " + record_read.getName() + " is both unalignable and alignable according to the input SAM/BAM file!");
+
// flush out previous read's hits if the read is alignable reads
if (record_val == 1) {
hits.updateRI();