This is a test version of Biostars. For the public version, visit https://www.biostars.org.
why gtf file has 0 length of CDS region

I retrieve CDS region from gtf file (ensembl v75). When I check the region length of CDS, I am very curious about the following result

library(data.table)

annotation <- fread('Homo_sapiens.GRCh37.75.gtf')

> summary(annotation[annotation$V3=="CDS",]$V5 - annotation[annotation$V3=="CDS",]$V4)    

   Min. 1st Qu.  Median    Mean 3rd Qu.    Max. 
    0.0    77.0   115.0   149.2   162.0 21692.0

More than 100 CDS regions have 0 length. Why? Please help me if you know, thanks a lot.

r genome annotation
• 1,937 views • link
updated by Devon Ryan Devon Ryan • written by wangshx • wangshx

1 answer

GTF is 1-based, so those are 1 base long microexons.

• 0 views • link
written by Devon Ryan Devon Ryan

It is new for me that exists 1 base long microexons.

• 0 views • link
written by wangshx • wangshx

Log in to answer this question.

More posts like this