BWA Aln Quality of 0 for perfect match read ?
0
0
Entering edit mode
2.0 years ago
ADDNOTHIING ▴ 10

Hi Everyone,

I have simulated short reads (<100bp) and trying to align to my reference using BWA_aln and I had a question when I saw my results. Most of my reads are mapping correctly and perfectly to the reference (90 matches out of 90 bp for example) but the MapQ is still equal to 0. Do you have any idea why is that happening ?

As those reads are simulated I know that they are correctly map aswell though.

Thanks a lot,

ADDNOTHIING

alignment bwa bwa aln MapQ • 397 views
ADD COMMENT
1
Entering edit mode

Please add representative outputs (from the resulting SAM file). Textual descriptions are often not very useful. Please show data.

ADD REPLY
0
Entering edit mode

How did you simulate the reads?

ADD REPLY
1
Entering edit mode

Using a software called Mitty (https://github.com/sbg/Mitty) incorporating sequencing errors and variation through VCF from 1k genomes project

ADD REPLY

Login before adding your answer.

Traffic: 2565 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6