A question about the MAPQ (mapping quality) in SAM file

[kaixinsjtu]

hey everyone,

I am looking into a SAM file now and found the MAPQ field all filled with 255. 255 means mapping quality not available, right? So I want to use bowtie to remap the tags in Fastq format, and hopefully get the mapping quality scores instead of 255.

My question is, with bowtie, how to calculate the MAPQ for each tag? I have read the bowtie manual and got no clue at all... If I use "--MAPQ <int>", that means I set the quality by myself. But I don't know it. Shouldn't bowtie be able to calculate quality scores for each tag based on the fastq file??

Thanks for your help in advance!

[pole_sana]

Peak calling for TSS detection

...........

[Bukowski]

hey everyone,

I am looking into a SAM file now and found the MAPQ field all filled with 255. 255 means mapping quality not available, right? So I want to use bowtie to remap the tags in Fastq format, and hopefully get the mapping quality scores instead of 255.

My question is, with bowtie, how to calculate the MAPQ for each tag? I have read the bowtie manual and got no clue at all... If I use "--MAPQ <int>", that means I set the quality by myself. But I don't know it. Shouldn't bowtie be able to calculate quality scores for each tag based on the fastq file??

Thanks for your help in advance!

Despite what it says in the SAM spec, I think for Bowtie this just means "alignment is non-repetitive" if you're using the default parameters.

[kaixinsjtu]

Despite what it says in the SAM spec, I think for Bowtie this just means "alignment is non-repetitive" if you're using the default parameters.

Thanks! But where is the mapping quality in SAM file? I want to calculate the probability that the tag is mapped incorrectly.

Thanks again.

[Bukowski]

If you're looking for MAQ stype mapping qualities I don't think it does them, but I'm happy to be proven otherwise.

[kaixinsjtu]

If you're looking for MAQ stype mapping qualities I don't think it does them, but I'm happy to be proven otherwise.

Thanks! I will look into MAQ

[Bukowski]

Thanks! I will look into MAQ

If you're going to take that step back, look at BWA instead which largely supercedes MAQ and still has sensible mapping qualities generated for variant calling.

[kaixinsjtu]

If you're going to take that step back, look at BWA instead which largely supercedes MAQ and still has sensible mapping qualities generated for variant calling.

Thanks!!