More posts like this
-
Genome annotation: convert gtf to gff3 file
written by shinyjj •
-
Plot read density across a normalized gene model for its each exon/intron length
written by kashiff007
-
How to create custom gtf annotation file?
written by John
-
reliability for the usage of merged GTF from tophat-cufflinks as an novel annotation?
written by Gabby •
-
Genomic feature graph
written by jb.sequence •
-
Get all intronic regions from a generic GTF file
written by alessandro.darmiento1991 •
-
parsing through gtf file
written by a.rex
-
C++ parse GTF file
written by scchess
-
Script for getting summary statistic of any genome using GTF or GFF3 ?
written by vahapel
-
Easiest way to compute RNA-Seq mapping stats (exons, introns, intergenic)?
written by Nick
What kind of output do you exactly want ? A toy example could help too.
Have you searched Biostars, there are many intron related questions. However, calculating introns is not as straight forward as it might seem. You cannot just take the gaps between all annotated exons, because different introns can exist for different isoforms. If you oversimplify the task (1 transcript isoform per gene only, no overlapping genes) anyway you can get very close to a one liner in R/BioC using the IRanges, GRanges, and their interval arithmetics (e.g. gaps).