How? I do not have any genomic coordinate?
• 1 views
•
link
written
by
abc.faisal •
How to find genomic coordinates of specific deletion in the HR gene?
We have found a novel deletion (429delC) in exon 2 of HR gene. I would like to find genomic coordinates of this location (may be based on GRCh37 assembly). Is there any easy solution for it?
• 2,964 views
•
link
updated
by
finswimmer
•
written
by
abc.faisal •
•
written
by
abc.faisal •
1 answer
- Ensembl's VEP
• 1 views
•
link
written
by
finswimmer
If it is the HR gene I posted a link for above you do have that info.
NCBI version or Ensembl (which you would use with VEP).
You just need to paste your gene and variant description into the data field. In your case the input looks like this: HR:c.429delC
That's all. The result includes the genomic region.
fin swimmer
• 1 views
•
link
written
by
finswimmer
Log in to answer this question.
More posts like this
-
RNAseq data for linked gene
written by gaom001 •
-
How to find read count for a newly discovered transcript
written by c_u
-
Is there a GRCh37 equivalent of EnsDb.Hsapiens.v86 or some other method for getting a complete list…
written by SpacemanSpiffo •
-
Extracting genomic coordinate, feature type and strand for a gene list
written by seta
-
How to find an inserted human gene in Drosophila RNA-seq data
written by williamsbrian5064
-
Localize the variant from the Crispr embryon editing paper
written by sacha
-
Find gene on genome from cds
written by guillaume.rbt
-
Difference between various GRCH 37 releases?
written by morteza.mahmoudisaber
-
Selecting and de-selecting exon-exon junctions based on depth or coverage from RNA-seq mapping
written by mjoyraj
-
Are Novel Junctions Only Give Novel Transcripts?
written by Ngsnewbie •
And what are your reference exons?
. The HR gene was amplified by specially designed oligonucleotide primers and the PCR products were purified by shrimp alkaline phosphatase and exonuclease enzymes and sequenced by using 3130 Genetic Analyzer (Applied Biosystems, Foster City, CA). I think it is RefSeqGene
Again, which RefSeq gene / sequence it refers to?
Do you have sequences of flanking regions? The do a BLAT https://genome.ucsc.edu/FAQ/FAQblat
OP is referring to HR gene.
A Gene can have multiple isoforms or transcripts. Moreover, the annotation of transcripts from different places can differ. So, one needs the reference sequence or transcript name to know exactly which transcript, as the coordinate is relative to that only.
You are absolutely right. To get an unambiguous result for sure, one needs the transcript, transcript version and reference genome.
But Ensembl's VEP accept also inputs, that might produce ambiguous results. It's the responsibility of the user to read these results correct.
fin swimmer
Good to know, thanks and +1 :) I did not expect this to be a known annotated variant!
And if you don't have the flanking regions, get the sequence of 2nd Exon from UCSC Table Browser (https://genome.ucsc.edu/cgi-bin/hgTables) and count the location manually
try: This annotation is for GRCh 37. At the end (hgvsg section), you will find all possible genomic coordinates. caution from ensembl:HR:c.429delC may resolve to multiple genomic locations"