Hello,
thank you for abra2, it is a nice realigner!
We just found an error in the realignment for a 3pb duplication, modifying cigar from
14M3I84M
to
11M3I87M
leading to an underestimation of the insertion and the creation of an additional mismatch during the variant calling.
If you want to have a look at the problem, I will be happy to send you small fastq/bam of the reads mapping in the region (or just carrying the duplication if you prefer).
Thank you
Camille