How to generate variant DNA sequence from an HGVS-formatted variant?

I’m parsing clinvitae using the python hgvs module. For each variant, I want to generate the variant DNA sequence in a 500 bp radius around the variant position. I don’t see a way to generate the variant sequence with the python hgvs module. Am I missing something? Or is there another library I should be using instead? I can code it myself if necessary, but don’t want to reinvent any wheels. Thanks in advance!