![]() |
|
![]() |
||||
Thread | Thread Starter | Forum | Replies | Last Post |
Find unmapped read from sam/bam file | genelab | Bioinformatics | 9 | 03-18-2014 01:35 PM |
Base Quality values in BAM file | donniemarco | Bioinformatics | 2 | 06-22-2011 10:06 AM |
BWA: specifying SAM/BAM file header fields before read alignment? | nora | Bioinformatics | 3 | 12-04-2010 09:11 PM |
Threshold quality score to determine the quality read of ILLUMINA reads problem | edge | General | 1 | 09-13-2010 02:22 PM |
replacing read id in bam file | blu78 | Bioinformatics | 2 | 06-23-2010 11:14 AM |
![]() |
|
Thread Tools |
![]() |
#1 |
Member
Location: italy Join Date: Apr 2010
Posts: 20
|
![]()
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: 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! |
![]() |
![]() |
![]() |
Thread Tools | |
|
|