This is a test version of Biostars. For the public version, visit https://www.biostars.org.
Retrieving RNA sequence from corresponding genomic coordinates

I have some genomic coordinates and I want to retrieve their RNA sequence especially the ones for reverse complement ones (since RNA sequence would not be strand specific). Where could I find such tool or package?

rna-seq
• 4,180 views • link
updated by Ram Ram • written by Aishwarya Kulkarni Aishwarya Kulkarni

What organism?

• 0 views • link
updated by Ram Ram • written by Neilfws Neilfws

2 answers

Hello,

you can also go to the UCSC Table Browser, upload your input file as a custom track ("add custom track" button), select your genome, choose "sequence" as output format, press "get output" and choose "mRNA" from the following menu. This would output your mRNA sequences in FASTA format.

• 0 views • link
updated by Ram Ram • written by Anima Mundi Anima Mundi

Thanks this is most useful! It's just that I cannot see the mRNA option once I click get sequence, all I get is:

  • All upper case.
  • All lower case.
  • Mask repeats: to lower case to N
  • Reverse complement (get '-' strand sequence)
• 0 views • link
updated by Ram Ram • written by Aishwarya Kulkarni Aishwarya Kulkarni

I see, then you can create an intersection with "mRNA and EST" (track "Human RNAs"), either to retrieve just the sequences overlapping between your data and human mRNAs, or also to get the complete mRNA sequences with at least partial overlap with your custom track.

• 0 views • link
updated by Ram Ram • written by Anima Mundi Anima Mundi

See this thread for details: Extract Cds Fastas From A Gff Annotation + Reference Sequence

• 0 views • link
updated by Ram Ram • written by arnstrm arnstrm

Log in to answer this question.

More posts like this