Dear All,
Now I would like to change the cigar value in the alignment in the bam
file
Below is the source code.
samfile_t *fp;
bam1_t *b=bam_init1();
bam_header_t *header;
bam1_core_t *c;
// go through each alignment
while(samread(fp,b) >= 0)
{
// Extract relevant information concerning current read
header = fp->header;
c=&b->core;
uint32_t *cigar = bam1_cigar(b)
string new_cigar="70M6S"
// how can I change the cigare in b using this new_cigar value
}
The question is that how can I change the cigare in b using this new_cigar
value in the above source code?
Thanks
------------------------------------------------------------------------------
_______________________________________________
Samtools-help mailing list
[email protected]
https://lists.sourceforge.net/lists/listinfo/samtools-help
