The CircularMapper method just
modifies the header of the resulting "rewrapped" BAM file to match the
original, unmodified reference genome. Apparently, this is not done
entirely correct in some cases which is why mapDamage and GATK complain
about this.
The CircularMapper method just modifies the header of the resulting "rewrapped" BAM file to match the original, unmodified reference genome. Apparently, this is not done entirely correct in some cases which is why mapDamage and GATK complain about this.