Answers.Body: <p>Yes, the mRNA results should include CDSs (not only, don't forget you also have UTR sequences which are non-coding so are not part of the CDS), but there is no reason why they should all be from genes on the + strand of the genome. Remember that the distinction between &quot;sense/+&quot; and &quot;antisense/-&quot; is essentially a human invention; it doesn't imply that the cell considers one strand to be in any way superior than the other, and genes are pretty randomly distributed across both strands.</p>
  <p>This means that when you align to a human reference genome, you will always have some sequences aligning to the + strand, if the RNAs were transcribed from genes that happen to be on the sense strand, and others will align to the - strand, if they were transcribed from genes found in the antisense strand.</p>
  <p>This is completely normal and expected, there is nothing wrong with your data. For example, human <a href="https://genome-euro.ucsc.edu/cgi-bin/hgTracks?db=hg38&amp;lastVirtModeType=default&amp;lastVirtModeExtraState=&amp;virtModeType=default&amp;virtMode=0&amp;nonVirtPosition=&amp;position=chr17%3A7668421%2D7687490&amp;hgsid=296456572_3scecBBFuWBIs3UAbZqJhMdmIY8A" rel="nofollow noreferrer">tP53</a> is on the antisense strand (note how the arrows are all pointing to the left):</p>
  <p><a href="https://i.sstatic.net/RaA7e.png" rel="nofollow noreferrer"><img src="https://i.sstatic.net/RaA7e.png" alt="screnshot of the UCSC genome browser showing tp53" /></a></p>
  
Answers.ContentLicense: CC BY-SA 4.0
Answers.CreationDate: 2023-04-25T15:37:36.977Z
Answers.Text: Can you help me understand this: does the bam file sequence always correspond to the +strand of the reference genome, regardless of which mate it came from?
Answers.UserDisplayName: geom_na
Answers.ContentLicense: CC BY-SA 4.0
Answers.CreationDate: 2023-04-25T15:49:22.143Z
Answers.Text: @Nomad420 yes, exactly. Basically, we arbitrarily decided that "this is the 'sense' strand" and then all coordinates are given with respect to that strand.
Answers.UserDisplayName: terdon
Answers.ContentLicense: CC BY-SA 4.0
Answers.CreationDate: 2023-04-25T16:06:33.797Z
Answers.Text: Thank you. This clears up a lot of my confusion!
Answers.UserDisplayName: geom_na
Answers.ContentLicense: CC BY-SA 4.0
Answers.CreationDate: 2023-04-25T20:30:37.877Z
Answers.Text: Sequences in a bam ar always in the forward direction, with respect to the reference, even if the original read was in the opposite orientation.  The flag indicates if the read as you see it in the bam is just like the original read, or if the bam is displaying the rev-comp.
Answers.UserDisplayName: swbarnes2
Answers.ContentLicense: CC BY-SA 4.0
Answers.CreationDate: 2023-04-25T14:59:09.827Z
Answers.LastActivityDate: 2023-04-25T14:59:09.827Z
Answers.UserDisplayName: terdon
Body: <p>I'm looking at RNAseq data from CCLE. The data is paired-end.</p>
  <p>Take the cell line Hs578T and the gene <em>HRAS</em> as an example.</p>
  <p>The cell line carries a G12D mutation (c.35G&gt;A), so the change in cds is:</p>
  <pre><code>ggc ggtgtgggca agagtgcgct g - Wildtype CDS
  gAc ggtgtgggca agagtgcgct g - Mutant CDS
   ^
  </code></pre>
  <p>My question is, when I <code>grep</code> the mutant CDS <code>gAcggtgtgggcaagagtgcgctg</code>, I do not get a match in my .bam file. But when I <code>grep</code> the reverse complement (i.e., anti-sense sequence), I get matches coming from both mates.</p>
  <p>My question is, why are there reads in the .bam file that correspond to the anti-sense strand of the DNA? Shouldn't all the reads of mRNA be the CDS?</p>
  <hr />
  <p>I got the answer -
  The punch line is that the sequence in .bam files is ALWAYS '+' strand of the reference, no matter what mate the read came from, or which strand the gene sits on.</p>
  
Tags: <rna-seq><strandedness>
Title: Why is there antisense sequence in RNAseq data
ThreadText: # Why is there antisense sequence in RNAseq data
  
  (Asked by: username_0 on 2023-04-25)
  
  I'm looking at RNAseq data from CCLE. The data is paired-end.
  
  Take the cell line Hs578T and the gene *HRAS* as an example.
  
  The cell line carries a G12D mutation (c.35G>A), so the change in cds is:
  
  ```
  ggc ggtgtgggca agagtgcgct g - Wildtype CDS
  gAc ggtgtgggca agagtgcgct g - Mutant CDS
   ^
  
  ```
  
  My question is, when I `grep` the mutant CDS `gAcggtgtgggcaagagtgcgctg`, I do not get a match in my .bam file. But when I `grep` the reverse complement (i.e., anti-sense sequence), I get matches coming from both mates.
  
  My question is, why are there reads in the .bam file that correspond to the anti-sense strand of the DNA? Shouldn't all the reads of mRNA be the CDS?
  
  ---
  
  I got the answer -
  The punch line is that the sequence in .bam files is ALWAYS '+' strand of the reference, no matter what mate the read came from, or which strand the gene sits on.
  
  **Tags:** rna-seq, strandedness
  
  
  ---
  
  ## 1 answers
  
  ### Answer by username_1 on 2023-04-25. Score: 3 | Accepted answer
  
  Yes, the mRNA results should include CDSs (not only, don't forget you also have UTR sequences which are non-coding so are not part of the CDS), but there is no reason why they should all be from genes on the + strand of the genome. Remember that the distinction between "sense/+" and "antisense/-" is essentially a human invention; it doesn't imply that the cell considers one strand to be in any way superior than the other, and genes are pretty randomly distributed across both strands.
  
  This means that when you align to a human reference genome, you will always have some sequences aligning to the + strand, if the RNAs were transcribed from genes that happen to be on the sense strand, and others will align to the - strand, if they were transcribed from genes found in the antisense strand.
  
  This is completely normal and expected, there is nothing wrong with your data. For example, human [tP53](https://genome-euro.ucsc.edu/cgi-bin/hgTracks?db=hg38&lastVirtModeType=default&lastVirtModeExtraState=&virtModeType=default&virtMode=0&nonVirtPosition=&position=chr17%3A7668421%2D7687490&hgsid=296456572_3scecBBFuWBIs3UAbZqJhMdmIY8A) is on the antisense strand (note how the arrows are all pointing to the left):
  
  [![screnshot of the UCSC genome browser showing tp53](https://i.sstatic.net/RaA7e.png)](https://i.sstatic.net/RaA7e.png)
  
  > **Comments:**
  > **username_0**: Can you help me understand this: does the bam file sequence always correspond to the +strand of the reference genome, regardless of which mate it came from?  
  > **username_1**: @Nomad420 yes, exactly. Basically, we arbitrarily decided that "this is the 'sense' strand" and then all coordinates are given with respect to that strand.  
  > **username_0**: Thank you. This clears up a lot of my confusion!  
  > **username_2**: Sequences in a bam ar always in the forward direction, with respect to the reference, even if the original read was in the opposite orientation.  The flag indicates if the read as you see it in the bam is just like the original read, or if the bam is displaying the rev-comp.  
  
  ---
  
  