I am using Rsubread to analyze human FASTQ files from the ENA. I downloaded a GRCh38 FASTA from ENCODE to create and index. I used the featureCounts(bam.files, annot.inbuilt = "hg38") function to count mapped reads for genomic features. Everything seems to work fine, but the count matrices are annotated to Entrez IDs, and I want the output using Ensembl IDs. Is there a way to do that?