Get Introns from a gft annotation file
2
1
Entering edit mode
7.7 years ago
xd_d ▴ 110

Hey all,

i have a question how i get the intron coordinates from Hg38 in bed file or bam . Have anybody a Workflow ?

thanks all :)

RNA-Seq Intron • 8.6k views
ADD COMMENT
1
Entering edit mode

Hi xd_d,

It's always best to show what you tried to solve your issue, people will be more eager to solve your bugs and help you if you show some effort from your side.

Cheers,
Wouter

ADD REPLY
6
Entering edit mode
6.9 years ago
Saber HQ ▴ 60

If you have the gff3 file, you may use the GenomeTools with the following command to include intron coordinates in the gff3 file format. Later you may also convert it to gtf if you want.

gt -gff3 -retainids -addintrons input_gff3 > output_gff3

ADD COMMENT
2
Entering edit mode

Seems the command has changed slightly. It now is as follows (the "-gff3" flag is now its own tool "gff3":

gt gff3 -retainids -addintrons input_gff3 > output_gff3

ADD REPLY
0
Entering edit mode

This works really well - exactly what I was looking for!

ADD REPLY
0
Entering edit mode
7.7 years ago
xd_d ▴ 110

I found this script on https://github.com/riverlee/IntronGTF and used that to get the introns.

thanks all

ADD COMMENT

Login before adding your answer.

Traffic: 1088 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6