BWA Aln Quality of 0 for perfect match read ?
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
• 1,452 views
•
link
written
by
ADDNOTHIING •
0 answers
No answers yet.
Log in to answer this question.
More posts like this
-
Plant transcriptome analysis and annotation
written by martta95 •
-
Origin of unmapped reads
written by Kay •
-
How to not align reverse complementary reads?
written by dcheng1 •
-
BWA cigar not corresponding to reality
written by ADDNOTHIING •
-
How to blastn a very short query sequence? Identifying telomeres in assembly
written by mrj •
-
Fasta to FastQ with known qualities
written by ADDNOTHIING •
-
Bowtie2 - How can MAPQ be used to obtain uniquely mapped reads?
written by ej
-
Why alignment in "blat" and "nucmer" are different?
written by Yuyayuya •
-
bwa aln calculate SA coordinate issue
written by Crystal •
-
Are Some Sequences Uniquely Mapped By Bwa When Mapped As Paired-End That Would Not Be When Mapped A…
written by KCC
How did you simulate the reads?
Using a software called Mitty (https://github.com/sbg/Mitty) incorporating sequencing errors and variation through VCF from 1k genomes project
Please add representative outputs (from the resulting SAM file). Textual descriptions are often not very useful. Please show data.