cDNA to genomic position

[cmccabe]

Is there a tool or script that will convert cDNA cordinates to genomic positions?


c.2T>C to chr15:38545390
c.7G>T to chr15: 38545393
c.7_20del to chr15:38545393
c.26A>T to chr15:38545412
c.30C>A to chr:38545416
c.42T>C to chr15:38591583

Thank you.

[TiborNagy]

I have never heared a tool like that. You need to align cDNA to the genome.

[cmccabe]

cDNA to genomic position

The hg19 reference was used for alingnment. Thanks.

[dpryan]

It turns out that those sorts of coordinates are a bit ambiguous (I've reported those sorts of coordinates myself in the past, so I too am guilty here). They only really make sense if a gene has only one splice form. Otherwise, you end up needing to first know which form the coordinates are for. If you know that, then it wouldn't be very tough to put together a little algorithm in R using GenomicRanges and an annotation file.

[tobias.mann]

Mutalyzer works pretty well

But it's not scriptable.

https://mutalyzer.nl/positionConverter

[cmccabe]

cDNA to genomic position

So the NM_ transcript that was used and coordinate is needed? Thanks.

[dpryan]

Yeah, otherwise you wouldn't know what transcript to use (the example is of SDHD, which has 4 splice forms).

[cmccabe]

cDNA to genomic position

So if NM_002880.3 if the transcript for RAF1, what coordinate information is needed? Thank you very much.

[dpryan]

For SNPs, you would just need the position (e.g., the 2 in c.2T>C). For deletions or insertion, you would need the start and width/length. That combined with a transcript designator is enough to reconstruct the genomic coordinates.