Hello,
I came across a bam file with one read having flag 103, which means:
read paired
read mapped in proper pair
read unmapped
mate reverse strand
first in pair
I think there is something weird as a read cannot be "mapped in proper pair" and "read unmapped" at the same time. It wasn't me doing the alignment, but as far as I know it was done with bwa 0.6.1 against hg19. Everything looks fine in the bam file (sort & indexed ok) apart from this read.
This is the incriminated read pair:
I think the problem comes from chr9_gl000201_random being 36148bp long and the read being aligned to the end of it (36144)?
Am I missing something or is this a bug in bwa? Any thought appreciated!
Dario
I came across a bam file with one read having flag 103, which means:
read paired
read mapped in proper pair
read unmapped
mate reverse strand
first in pair
I think there is something weird as a read cannot be "mapped in proper pair" and "read unmapped" at the same time. It wasn't me doing the alignment, but as far as I know it was done with bwa 0.6.1 against hg19. Everything looks fine in the bam file (sort & indexed ok) apart from this read.
This is the incriminated read pair:
Code:
1:488660:43 103 chr9_gl000201_random 36144 60 36M chrM 267 0 CGATGGATCACAGGTCTATCACCCTATTAACCACTC 1:BDEFDHHHHABGH@>DHEEHJJJFHIFBFI>DD9 XT:A:U NM:i:2 SM:i:25 AM:i:25 X0:i:1 X1:i:0 XM:i:2 XO:i:0 XG:i:0 MD:Z:0G3C31 1:488660:43 147 chrM 267 60 36M chr9_gl000201_random 36144 0 TCCACACAGACATCATAACAAAAAATTTCCACCAAA HCD:<BDFGGE:?AAGFFIGHF>BFCA<A=2?B?@? XT:A:U NM:i:0 SM:i:37 AM:i:25 X0:i:1 X1:i:0 XM:i:0 XO:i:0 XG:i:0 MD:Z:36
Am I missing something or is this a bug in bwa? Any thought appreciated!
Dario
Comment