This is a test version of Biostars. For the public version, visit https://www.biostars.org.
mapping quality for sanger reads

Hi

Anyone knows if BWA-mem automatically detect if the reads that I am aligning have their quality is fastq-sanger?. If not, this affect the mapping quality score?

Cheers

alignment
• 2,696 views • link
updated by Ram Ram • written by shinken123 shinken123

1 answer

BWA aln provides -I option that can be used to align fastq reads with Phred+64 quality scores. The reads in output bam file have base quality in Sanger or Phred+33 format. BWA-mem doesn't have that option (at least not stated in the manual). If you search on this forum there are multiple posts that explain in detail about how to test if reads belong to Phred+33 or Phred+64. Also, you can easily convert them from one to other. Base quality scores are not used in the calculation of mapping quality scores but BWA mem may throw an out of range error if it can't take Phred+64 fastq files.

• 0 views • link
updated by Ram Ram • written by Ashutosh Pandey Ashutosh Pandey

Thank you very much. So if my reads are from a sanger sequencing (Phred+33) and BWA-mem assumes that the sequences are in Phred+33, is going to work fine with my data.

Thanks again

Regards

• 0 views • link
written by shinken123 shinken123
Well if your reads are phred+33 then you don't need to do anything. Just use BWA mem.
• 0 views • link
written by Ashutosh Pandey Ashutosh Pandey

Log in to answer this question.

More posts like this