Using gtf file to get gene positions
Hi, I am wanting to extract gene positions from a gtf file using R. I have uploaded the gtf file as a df in R, with the headings 'seqnames', 'start', 'end', 'strand' and 'gene_id'. I am wanting to use this genome annotation to get the gene positions. Does anyone know a way I can do this on R.
My desired output would produce the different gene positions from the gtf file.
1 answer
Hi
You should get subset of the data frame in another data frame, try like following for 3 column of it.
df2 <- unique( df[,c("start","end","gene_id")])
• 0 views
•
link
written
by
ahmad mousavi
Log in to answer this question.
More posts like this
-
Extracting Sequences from a FASTA File Using Exon Coordinates from a GTF File
written by Milica •
-
How to use GFF3 annotation to split genome fasta into gene sequence fasta in R
written by Chilly •
-
Extract transcript ID and gene ID from ITAG4.1_gene_models.gff
written by Tian •
-
Intersecting transcriptome bam file with GTF file
written by Jalil Sharif •
-
Converting bedgraph file to wig file
written by margo •
-
Extracting fasta sequence in R
written by margo •
-
Writing a function in R which will identify peaks within a data frame
written by margo •
-
Finding the read ratio of transcription termination data using R
written by margo •
-
How to locate and extract data from a df in R
written by margo •
-
How to edit a GTF file
written by margo •
