BQSR and MD fields
The GATK indel realigner appears to not write an MD field in realigned BAM file entries. I wonder how the base quality score recalibration
program deals with these entries. Does it consider them to be ambiguous and so ignore them, or does it use include them, using the
reference sequence to identify SNPs?