Hi everyone,
I am trying to get the quality of the read from the bam file with the samtools c++ library.
I am using the following (partial) code:
but I am not getting correct values for the read quality.
Does anyone have any suggestion on what I am doing wrong?
Thanks in advance!
I am trying to get the quality of the read from the bam file with the samtools c++ library.
I am using the following (partial) code:
Code:
bam1_t *b = bam_init1();
...
while(samread(reads, b) >=0){
uint8_t *readQual = bam1_qual(b);
...
cout<<"quality:"<<endl;
for(int k=0;k<readLength;k++){
cout<<readQual[k]<< " ";
}
cout<<endl;
}
Does anyone have any suggestion on what I am doing wrong?
Thanks in advance!