BWA Aln Quality of 0 for perfect match read ?

0

Entering edit mode

5.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 • 969 views

ADD COMMENT • link 5.0 years ago by ADDNOTHIING ▴ 10

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 • link 5.0 years ago by ATpoint 85k

0

Entering edit mode

How did you simulate the reads?

ADD REPLY • link 5.0 years ago by Mark ★ 1.6k

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 • link 5.0 years ago by ADDNOTHIING ▴ 10

Login before adding your answer.