X-Git-Url: http://woldlab.caltech.edu/gitweb/?p=samtools.git;a=blobdiff_plain;f=bam_color.c;h=ce637f708db56de777e86f2f0ad50478fb1a6657;hp=75aedd603a1798e85a40c982dd9392ee2a3e5b3f;hb=4a17fa7e1f91b2fe04ad334a63fc2b0d5e859d8a;hpb=b27e00385f41769d03a8cca4dbd71275fc9fa906 diff --git a/bam_color.c b/bam_color.c index 75aedd6..ce637f7 100644 --- a/bam_color.c +++ b/bam_color.c @@ -100,8 +100,8 @@ char bam_aux_getCEi(bam1_t *b, int i) cs_i = strlen(cs) - 1 - i; // get current color cur_color = cs[cs_i]; - // get previous base - prev_b = (0 == cs_i) ? cs[0] : bam_nt16_rev_table[bam1_seqi(bam1_seq(b), i+1)]; + // get previous base. Note: must rc adaptor + prev_b = (cs_i == 1) ? "TGCAN"[(int)bam_aux_nt2int(cs[0])] : bam_nt16_rev_table[bam1_seqi(bam1_seq(b), i+1)]; // get current base cur_b = bam_nt16_rev_table[bam1_seqi(bam1_seq(b), i)]; }