Bug Bulletin: we have identified a bug that affects indexing when producing gzipped VCFs. This will be fixed in the upcoming 3.2 release; in the meantime you need to reindex gzipped VCFs using Tabix.

Errors while running PrintReads walker of BQSR step

kssrkssr Posts: 6Member
edited October 2013 in Ask the team

I ran BaseRecalibrator on some of my realigned bam files, there were no errors reported while running it and I was to able to generate the "recal.grp" files. However, while running the PrintReads walker to generate the recalibrated BAM files, I get the following error message for some of the BAMs.

##### ERROR MESSAGE: Exception when processing alignment for BAM index WTCHG_35305_101:4:1105:10332:129428#TGTTAACT 2/2 100b aligned read.

So, I tried to validate my realigned BAM as well as the original BAM before realignment using Picard's ValidateSAMFile and I get the following:

original BAM: Mate unmapped flag does not match read unmapped flag of mate, Mate alignment does not match alignment start of mate, Mate negative strand flag does not match read negative strand flag of mate

Aditionally these errors on realigned BAMs: Mate reference index (MRNM) does not match reference index of mate, Mate not found for paired read

Do I need to worry about initial alignments? I read on the forums that using -rf MateSameStrand Filter should help me work around this, what does this filter exactly do? Any other approaches to solve this problem would be appreciated.

Post edited by kssr on
Sign In or Register to comment.