I am analysing bulk RNAseq data of 86 samples by DESeq2. I really liked one feature of Seurat in scRNAseq analysis "mitochondrial content percentage" . Is there a simple way to extract "mt percentage" information from deseq2 object. Like
plotCounts(dds, gene="MT-", intgroup="Group")
FYI: My deseq2 object have gene symbols, so all mitochondrial genes have MT- as suffix.