Remapping genomic coordinates to account for indels
I think that you need a LiftOver Chain file to transform your coordinates. You can obtain such a file using bcftools consensus with the -c parameter:

-c, --chain &lt;file&gt;         write a chain file for liftover


Then you can use it to transform coordinates in various genomic formats using CrossMap.