Question: bimodal distribution of bam mapping quality
0
gravatar for irritable_phd_syndrom
3.9 years ago by
irritable_phd_syndrom70 wrote:

I have some RNA-Seq data that I analyzed with tophat2. The command that I used to generate it is

/path/to/tophat-2.1.0/tophat -p 20 -o ouputdir --library-type fr-firststrand /reference/homo_sapiens/GRCh38/ensembl/Sequence/Bowtie2Index/Homo_sapiens.GRC38 my_trimmed_data.fq.gz

This output a file, accepted_hits.bam with 50e6 aligned reads. When I plot histogram of the mapping quality scores, roughly 27e6 reads have a mapping quality value [0-3] and 23e6 reads have a mapping quality value of 50. There are _no_ values in between.

How should I interpret this bimodal distribution of mapping quality scores? This seems very strange to me.

rna-seq alignment tophat • 1.5k views
ADD COMMENTlink modified 3.9 years ago by dariober11k • written 3.9 years ago by irritable_phd_syndrom70
3
gravatar for John
3.9 years ago by
John12k
Germany
John12k wrote:

This is normal :)

Tophat2 Mapping Qualities

ADD COMMENTlink written 3.9 years ago by John12k
1
gravatar for dariober
3.9 years ago by
dariober11k
WCIP | Glasgow | UK
dariober11k wrote:

Besides Tophat, where intermediate mapq is avoided by design, the mapping quality from other aligners (e.g. bwa) tends to be markedly bimodal. I.e. reads tend to be either very unambiguously mapped somewhere or they can map equally well at multiple places.

ADD COMMENTlink written 3.9 years ago by dariober11k
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: 769 users visited in the last hour