X-Git-Url: http://woldlab.caltech.edu/gitweb/?p=pysam.git;a=blobdiff_plain;f=samtools%2Fbam.c.pysam.c;h=28cdbd74b0d751246c6726f5f88a1499eb87a30f;hp=ad263395de51eb3b38d56ffdcd9041a8dac1864e;hb=e1756c41e7a1d7cc01fb95e42df9dd04da2d2991;hpb=ca46ef4ba4a883c57cea62d5bf1bc021f1185109 diff --git a/samtools/bam.c.pysam.c b/samtools/bam.c.pysam.c index ad26339..28cdbd7 100644 --- a/samtools/bam.c.pysam.c +++ b/samtools/bam.c.pysam.c @@ -34,7 +34,7 @@ int32_t bam_cigar2qlen(const bam1_core_t *c, const uint32_t *cigar) int32_t l = 0; for (k = 0; k < c->n_cigar; ++k) { int op = cigar[k] & BAM_CIGAR_MASK; - if (op == BAM_CMATCH || op == BAM_CINS || op == BAM_CSOFT_CLIP || op == BAM_CEQUAL || op == BAM_CDIFF) + if (op == BAM_CMATCH || op == BAM_CINS || op == BAM_CSOFT_CLIP) l += cigar[k] >> BAM_CIGAR_SHIFT; } return l; @@ -270,7 +270,7 @@ char *bam_format1_core(const bam_header_t *header, const bam1_t *b, int of) else { for (i = 0; i < c->n_cigar; ++i) { kputw(bam1_cigar(b)[i]>>BAM_CIGAR_SHIFT, &str); - kputc("MIDNSHP=X"[bam1_cigar(b)[i]&BAM_CIGAR_MASK], &str); + kputc("MIDNSHP"[bam1_cigar(b)[i]&BAM_CIGAR_MASK], &str); } } kputc('\t', &str);