Where to get Gene start position to end positon for all gene in an organism
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.
2 answers
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
• 0 views
•
link
written
by
Devon Ryan
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
Log in to answer this question.
More posts like this
-
How to add the coordinates of UTRs, CDS, start codon and stop codon to
written by Issac •
-
Downloading multiple gene sequence exon-wise using commandline
written by Saumya Gupta •
-
Copy Number Variations
written by Nyksubuz •
-
how to get distance of each SNPs to the transcripts' 5'UTR and 3'UTR?
written by mikysyc2016
-
SnpEff START/STOP codon not found error
written by thomasbersez
-
Getting transcriptomic coordinates for CDS exons from a GTF
written by rubic
-
Extracting entries from GTF file
written by EVR
-
should stop codon be counted as part of CDS
written by -_-
-
Does CDS, start_codon and stop_codon in gtf affect transcriptome assembly by StringTie?
written by syrttgump •
-
sorting merged GTFs for Tophat/Cufflinks pipeline
written by trakhtenberg
