Question: mpileup mapping quality
0
gravatar for hpapoli
2.0 years ago by
hpapoli90
Sweden
hpapoli90 wrote:

Hello,

When using samtools mpileup with -s option, I get a 7th column with mapping quality. I realise the mapping quality in SAM format is calculated for a given read. I wonder what does a mapping quality mean for a single base position as it is the case in mpileup? Below is an example:

NW_008793873.1  12      A       1       .       @       H
NW_008793873.1  13      G       2       .^S.    AB      HS
NW_008793873.1  14      A       2       ..      @B      HS
NW_008793873.1  15      C       2       ..      @B      HS
NW_008793873.1  16      g       2       ..      BC      HS
NW_008793873.1  17      a       2       ..      A@      HS

How can I interpret and use the mapping quality field for filtering bases?

Thank you!

mapq samtools • 1.1k views
ADD COMMENTlink modified 2.0 years ago by Kevin Blighe63k • written 2.0 years ago by hpapoli90
0
gravatar for Kevin Blighe
2.0 years ago by
Kevin Blighe63k
Kevin Blighe63k wrote:

These are Sanger-encoded quality scores. I cannot confirm if the given score is the max MAPQ over the position, the average, the median, or the min (or something else).

Sanger encoding goes from characters 33 to 126 on the ASCII scale, covering the MAPQ range 0 to 93, as per Table 1 from: The Sanger FASTQ file format for sequences with quality scores, and the Solexa/Illumina FASTQ variants.

You can look up your ASCII codes here: https://pt.wikipedia.org/wiki/ASCII

Kevin

ADD COMMENTlink modified 2.0 years ago • written 2.0 years ago by Kevin Blighe63k
Please log in to add an answer.

Help
Access

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 1239 users visited in the last hour