bam1_t sequence setting

set_base.c

/*
 * :abstract: set the base in pSeq at index i in the read to the base in bSeq at index i in the read.
 * These 8-bit integers are the points in memory returned by bam_get_seq(p) and bam_get_seq(b), respectively.
 * :param: pSeq [uint8_t *]
 * :param: bSeq [uint8_t *]
 * :param: i [int/integral type]
*/
#define set_base(pSeq, bSeq, i) (pSeq)[(i)>>1] = ((bam_seqi(bSeq, i) << (((~i) & 1) << 2)) | (((pSeq)[(i)>>1]) & (0xf0U >> (((~i) & 1) << 2))))
/*
 * :abstract: set the base at index i in the read to "N".
 * :param: pSeq [uint8_t *]
 * :param: i [int/integral type]
*/
#define n_base(pSeq, i) pSeq[(i)>>1] |= (0xf << ((~(i) & 1) << 2));