How to convert transcript-relative coordinates to genomic coordinates?
0
0
Entering edit mode
2.8 years ago
adixon3 • 0

I have queried using Entrez Utilities (efetch: https://www.ncbi.nlm.nih.gov/books/NBK25499/) and obtained annotations for transcripts like the following:

>Feature ref|NM_152486.3|

1 2557 gene

        gene    SAMD11

        gene_syn    MRS

        gene_desc   sterile alpha motif domain containing 11

        db_xref GeneID:148398

        db_xref HGNC:HGNC:28706

        db_xref MIM:616765

How/what database should I query to get the proper genomic coordinates for this gene and convert the transcript coordinates to absolute genomic coordinates?

IMPORTANT!: many of these transcripts are out of date, so just getting a table from UCSC is out of the question (unless transcript ID doesn't matter for getting genomic coordinates...)

Thanks!

refseq transcript • 1.1k views
ADD COMMENT
0
Entering edit mode

You will need to find the annotation release when the RefSeq was last annotated and extract that information from the GFF3/GTF files for that annotation. Since this is a human gene, which version of the genome are you interested: GRCh38/hg38 or GRCh37/hg19?

ADD REPLY

Login before adding your answer.

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