[Genome] Converting coordinates on a mRNA to coordinates on the genome
Rachel Harte
hartera at soe.ucsc.edu
Wed Jul 4 17:31:56 PDT 2007
Hello Archie,
We have a tool called pslMap which is in the directory,
src/hg/pslMap/, in the Genome Browser source tree that will help you
do this calculation. If you type pslMap at the command line, you will get
help for using the program.
The input files should be in PSL format - see
http://genome.ucsc.edu/FAQ/FAQformat#format2
You should create a PSL file of coordinates for your mRNAs as the inPsl
and then the mapPsl should be a PSL file of mRNA alignments to a
genome downloaded from our downloads server:
http://hgdownload.cse.ucsc.edu
Click on the organism of interest, find the correct assembly, and then
click on the "Annotation database" link.
all_mrna.txt.gz is the PSL file of mRNA alignments and refSeqAli.txt.gz is
the PSL file for RefSeq mRNA alignments.
I hope that this helps you. Please let us know if you have further
questions.
Rachel
Rachel Harte
UCSC Genome Bioinformatics Group
http://genome.ucsc.edu
On Tue, 3 Jul 2007, Russell, Archie wrote:
>
> Hi,
>
> I have some coordinates relative to mRNAs (transcript features) that I'd
> like to turn into coordinates on the genome. I can blat the mRNA to
> the genome, but doing the math for the translation is a little tricky,
> in particular on the reverse strand. Do you have any code or
> recommendations for how to do this?
>
> Thanks,
> Archie
>
>
> ------------------------------------------------------------------------------
> Notice: This e-mail message, together with any attachments, contains
> information of Merck & Co., Inc. (One Merck Drive, Whitehouse Station,
> New Jersey, USA 08889), and/or its affiliates (which may be known
> outside the United States as Merck Frosst, Merck Sharp & Dohme or MSD
> and in Japan, as Banyu - direct contact information for affiliates is
> available at http://www.merck.com/contact/contacts.html) that may be
> confidential, proprietary copyrighted and/or legally privileged. It is
> intended solely for the use of the individual or entity named on this
> message. If you are not the intended recipient, and have received this
> message in error, please notify us immediately by reply e-mail and then
> delete it from your system.
>
> ------------------------------------------------------------------------------
>
> _______________________________________________
> Genome maillist - Genome at soe.ucsc.edu
> http://www.soe.ucsc.edu/mailman/listinfo/genome
>
More information about the Genome
mailing list