I'm currently analyzing my paired-end RNA-seq data using the sorted-position.bam files and a reference.genome.gtf file.
I used the default settings of the featureCounts tool to generate a count.matrix file, which generated count.matrix contains only gene-IDs without any corresponding actual gene names or feature.
While this isn't necessarily a problem, at this stage, since the gene-IDs are unique, I would like to extract the actual names from the gtf file and merge them into the featureCounts.CSV file before proceeding with edgeR/deseq2 analysis.
When using edgeR and Deseq2 to generate a CSV file for differential gene expression analysis, only the gene-IDs are included in the output. As a result, when visualizing the results using ggplot2, the gene-IDs are displayed instead of the actual gene names.
If anyone has any suggestions or solutions, I would greatly appreciate your help. Thank you.