Getting the mRNA sequence from genome annotation
From my understanding, if an annotated gene from Ensembl, UCSC, etc. is marked as being on the positive strand, that means its sense strand is the exact sequence from the assembly fasta file, and so the mRNA sequence is a copy of the genome assembly sequence for the exon regions (except with uracil replacement). For a gene on the negative strand, then, the reverse complement of the fasta sequence needs to be used to get the mRNA sequence. Is this correct, or do I have this backwards?
• 1,257 views
•
link
updated
by
swbarnes2
•
written
by
colin.kern
•
written
by
colin.kern
2 answers
Yes that is correct assuming you want the exonic sequence in 5' to 3' order relative to the mRNA. Also note RNA-editing or alternative splicing may change the sequence from the genome reference.
• 0 views
•
link
written
by
benformatics
Log in to answer this question.
More posts like this
-
Strand in gff file
written by blz •
-
Extracting promoter sequence with GenomicFeatures ?
written by lokraj2003
-
In RSAT retrieve_Ensembl_seq do FASTA sequences of genes on reverse strands take this into account?
written by Monika515 •
-
Distinguishing coding DNA using NCBI and UCSC databases
written by pairwiseseq •
-
Confusion in distinguishing between (+) and (-), cDNA and mRNA using RefSeq and UCSC genome
written by pairwiseseq •
-
Genes on negative strand?
written by christylynn002 •
-
Bowtie2 local alignment SAM file SEQ output - reverse complement or not ?
written by Roger •
-
Reference sequence and gtf annotation.
written by juncheng
-
Is the mRNA sequence given in the nucleotide database of NCBI is strand specific?
written by bedantamohanty •
-
Negative Sense Exon Strand Coordinates In Ucsc
written by Max
