Rare variant gene burden test annotation

Hi,

I am currently conducting rare variant gene burden test, but before that, I would like to generate an annotation file that showing the variant, gene and annotation (or the most deleterious).
variant | gene_symbol | csq
| “chr1:17379:G:A” | “MIR6859-1”|“mature_miRNA_variant”
| “chr1:95068:G:A” | “AL627309.1”|“intron_variant”
| “chr1:111735:C:A” | “AL627309.1”|“intron_variant”

I would like to appreciate any suggestion by using vep() or other tools. Thank you.

Best,
Yuyang