This is a test version of Biostars. For the public version, visit https://www.biostars.org.
Finding the maximum intron length from aligned exomic data

Hello,

I'm currently working on some plant exomic data and I need to calculate the maximum intron length, I've done some searching and found answers for RNA-Seq but not exome seq and I was wondering if anyone could help me with some suggestions regarding any methods or point me towards any literature calculating maximum intron length in exomic data? The annotations for the reference genome for the species I am working with does not contain any information on introns at all.

Thank you.

alignment exome analysis intron size
• 2,895 views • link
written by GLFrey • GLFrey

maximumintronlength makes no sense as a tag. Maybe add better, more relevant tags?

• 0 views • link
written by Ram Ram

1 answer

edit: I pointed to a fork of the old GFFutils. The repository for this old version is https://github.com/daler/GFFutils_old, and the repository for the current version is https://github.com/daler/gffutils.

end of edit.

What is the genome you are working with? What is the format of the annotation? If the annotation contains information on exons, it does contain on introns as well: introns are delimited by same-gene exons boundaries.

Some pointers for GTF / GFF files:

https://github.com/seandavi/GFFutils (see the Imputing introns section)

Intron size from gtf/gff file ?

http://bioops.info/2012/11/intron-size-gff3-perl/

And this gist with a perl script.

• 0 views • link
written by h.mon h.mon

Hello h.mon,

Once again, thank you for your help. The genome annotation is for wheat and the annotation files are in GTF and GFF3 format. I will have a look through the links you have posted.

Thank you!

• 0 views • link
written by GLFrey • GLFrey

In particular, have a look at the GFFutils repository, there is a section titled Imputing introns

• 0 views • link
written by h.mon h.mon

Thank you so much for this, h.mon it looks perfect.

• 0 views • link
written by GLFrey • GLFrey

Log in to answer this question.

More posts like this