Where to get Gene start position to end positon for all gene in an organism
2
0
Entering edit mode
7.9 years ago
Bioblazer ▴ 50

In a transcript file (.gtf) we can able to get exon wise start position and stop position, CDS and start codon region and stop codon region,but I want start position and end position of all genes (not exon wise) with its chromosome number in an organism.

gene genome • 5.1k views
ADD COMMENT
3
Entering edit mode
7.9 years ago

Normally there's a "gene" entry for each gene, so:

awk 'BEGIN{FS="\t"; OFS="\t"}{if($3 == "gene") print $1, $4, $5}' foo.gtf
ADD COMMENT
0
Entering edit mode

Thank you very much for your kind reply

ADD REPLY
2
Entering edit mode
7.9 years ago
EagleEye 7.5k

As Devon said if the GTF file is from Gencode, you will have "gene" entry. Extracting "gene" entry will give you the desired results. If you use Ensembl GTF annotation (hg19 upto GRCh37.74), it does not have "gene" entry. In this case use the following script.

The output will look similar to

enter image description here

ADD COMMENT
0
Entering edit mode

Thank you so much for your valuable reply.

ADD REPLY
0
Entering edit mode

Hello, I have this exact same question but I'm not seeing the code here. Anyway you could please resend the code?

ADD REPLY

Login before adding your answer.

Traffic: 2433 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