This is a test version of Biostars. For the public version, visit https://www.biostars.org.
Tophat2 Mapping Qualities

When I run Tophat2 on human genome, I am using default parameters and with a GTF file from ensemble.

The accepted hits bam file has only three Mapping quality scores in 5th column in my file i.e 0, 3, 50. I would like to know is this correct or I have done something wrong.

As per my knowledge, mapping quality of 0 means the read mapped at multiple places but for other reads, I could see only 3 and 50. The data is from Illumina HiSeq 1000, paired end library.

I would be happy if anybody could help me in figuring out why I am getting only three figures .

tophat2 bowtie2
• 11,390 views • link
updated by Devon Ryan Devon Ryan • written by GouthamAtla GouthamAtla

you may also be interested in checking different option for these in tophat2 like -g 1

• 0 views • link
written by GR GR

1 answer

Your results are normal, the MAPQ scores reported by tophat2 are not related to -10*log10(probability the mapping is wrong). It's 50 for uniquely mapped, and then 0-3 for various degrees of multiple mapping.

• 284 views • link
written by Devon Ryan Devon Ryan

Is this documented anywhere?

• 0 views • link
written by Fedor Gusev Fedor Gusev

Not that I know of. To make life slightly more complicated, the scores used to be 0-3 and 255, so don't be surprised if you see that if you have older datasets.

• 0 views • link
written by Devon Ryan Devon Ryan

I see quality scores of mostly 3 and 50 in recent tophat (2.0.13) but also of 41, 42, 44 and 24, 28

• 0 views • link
written by brentp brentp

Gotta love undocumented changes.

• 0 views • link
written by Devon Ryan Devon Ryan

Some addition to the above answer:

255 = unique mapping

3 = maps to 2 locations in the target

2 = maps to 3 locations

1 = maps to 4-9 locations

0 = maps to 10 or more locations.

• 139 views • link
written by Ashutosh Pandey Ashutosh Pandey

50 means unique mapping for tophat2 ??

• 0 views • link
written by GouthamAtla GouthamAtla

dpryan79 may be right with that. I haven't used the latest version of tophat2. May be now they score uniquely aligned reads with 50 MAPQ. One reason could be that some downstream tools like GATK complain when they see a MAPQ of 255.

• 0 views • link
written by Ashutosh Pandey Ashutosh Pandey

Yeah, they changed from 255 to 50 at some point. I have no clue which release had the change.

• 314 views • link
written by Devon Ryan Devon Ryan

Log in to answer this question.

More posts like this
  • Htseq Count
    written by sethugunja • sethugunja
    <p>Hi guys,</p> <p>I want to use htseq-edgeR to validate my differentially expression cuffdiff results.</p> <p>For this, first htseq-count have to be used to count the …
  • Extract Regions From A Bam To A Bed File Where The Average Mapping Quality Is Below 30.
    written by William William
    <p>How can I extract the regions with a mapping quality below 30 from a <a href='http://samtools.sourceforge.net/SAM1.pdf'>bam</a> file to a bed file ?</p> <p>I am not …
  • Error In Converting Sam To Bam By Samtools
    written by Raghav Raghav
    <p>hello every one,</p> <p>after mapping,i encountered an error during converting my <a href='http://samtools.sourceforge.net/SAM1.pdf'>sam</a> out file into <a href='http://samtools.sourceforge.net/SAM1.pdf'>bam</a> out put file.</p> <p>[cdac@nbri <a href='http://samtools.sourceforge.net/'>samtools</a>-0.1.19]$ <a …
  • How To Extract Consensus Sequence Of Mapped Regions From A Bam File?
    written by lin.barnum • lin.barnum
    <p>I would like to extract consensus sequence of only the mapped regions from a <a href='http://samtools.sourceforge.net/SAM1.pdf'>BAM</a> file. I have a <a href='http://samtools.sourceforge.net/SAM1.pdf'>bam</a> file which contains …
  • How To Interpret The Flag Field In The Sam Specification
    written by deepthithomaskannan deepthithomaskannan
    <p>Hi all ,</p> <p>I created a <a href='http://samtools.sourceforge.net/SAM1.pdf'>sam</a> file and converted it to <a href='http://samtools.sourceforge.net/SAM1.pdf'>bam</a> file using <a href='http://samtools.sourceforge.net/'>samtools</a>.</p> <p>One of the record in <a …
  • Problem With Sorting Bam File: Samtools
    written by lyfsa • lyfsa
    <p>Hi!</p> <p>I have three aligned pair end reads in <a href='http://samtools.sourceforge.net/SAM1.pdf'>BAM</a> format. Two of them are of same sample run in different lanes. I first …
  • Htseq-Count Not Working With Bam From Tophat
    written by siddharth.sethi5 siddharth.sethi5
    <p>Hi everyone <br /> I have RNA seq data. I used Tophat2 for the alignment, got the junctions successfully.</p> <p>When I use HTSeq-count (latest version) …
  • Filtering A Sam File For Quality Scores
    written by Varun Gupta Varun Gupta
    <p>Hi I have a <a href='http://samtools.sourceforge.net/SAM1.pdf'>sam</a> file. The 5th column of the <a href='http://samtools.sourceforge.net/SAM1.pdf'>sam</a> file tells us the quality score. In my <a href='http://samtools.sourceforge.net/SAM1.pdf'>sam</a> alignment …
  • Sam To Bam Conversion Reference Fasta File
    written by Atom Smasher • Atom Smasher
    <p>Hello,</p> <p>I had some data in "eland_export.txt" format and I want to convert it to the <a href='http://samtools.sourceforge.net/SAM1.pdf'>SAM</a> format first and then later to the …
  • Soap Result Converted To Bam Not Working In Samtool
    written by Sarwar • Sarwar
    <p>I converted the output of soap2 which is map file to <a href='http://samtools.sourceforge.net/SAM1.pdf'>sam</a> file using map2sam.pl and then converted this <a href='http://samtools.sourceforge.net/SAM1.pdf'>sam</a> to <a href='http://samtools.sourceforge.net/SAM1.pdf'>bam</a> …