Seqanswers Leaderboard Ad

Collapse

Announcement

Collapse
No announcement yet.
X
 
  • Filter
  • Time
  • Show
Clear All
new posts

  • #16
    Originally posted by cliffbeall View Post
    I am fairly sure fr-firststrand means read one (aka forward read or fr) corresponds to the reference genome sequence (firststrand). fr-secondstrand means read one corresponds to the reverse complement (secondstrand) of the reference genome. FR and RF would be a similar concept - but they are just making a shorthand where if the read one or forward read corresponds to the reference strand it maps to the left of the reverse read on a linear representation of the genome (FR) and if it corresponds to the reverse complement of the reference it maps to the right of the reverse (RF).

    I hope that makes sense -- you might need to go back to DNA/genome structure and how Illumina works to really understand the above.
    I'm afraid that this is incorrect.

    First 'fr' does not mean 'forward read'. What it means is the the two reads are arranged in the FR (Foward/Reverse) as opposed to the RF (Reverse/Forward) orientation. FR means that the read pairs are oriented with the 3' end point toward each other when aligned to the reference and RF means they are pointing away from each other.

    Code:
    FR:
    -------------------------------------------------------------------------
       ----------------->                              <-------------------
    
    RF:
    -------------------------------------------------------------------------
       <-----------------                               ------------------->
    Firststrand means that read #1 matches the first strand of cDNA generated during conversion of the single stranded mRNA to ds-cDNA. In other words the R1 is the reverse complement of the mRNA. Note that all of these reference landmarks are relative to transcript, not the genome. For Illumina strand specific, paired end mRNA-Seq reads the orientation of the reads relative to their original RNA is:

    Code:
      mRNA 5'-----------------------------------------------------------AAAAAAAAAAAA
                       R2 ------------>                <------------ R1

    Comment


    • #17
      Looks like Jim Robinson had it backwards for the color legends. Actually

      red=positive,
      blue=negative

      He has clarified here:

      Comment


      • #18
        Originally posted by sowmyai View Post
        Looks like Jim Robinson had it backwards for the color legends. Actually

        red=positive,
        blue=negative

        He has clarified here:
        https://groups.google.com/forum/#!to...lp/YiVzwnUuOZM
        I always thought that blue=positive and red=negative. See screen shots of my data in IGV.

        I used Illumina dUTP-type stranded library prep and fr-firststrand with Tophat2. I find several places that say dUTP libraries should use fr-firststrand with Tophat, for example: http://onetipperday.sterding.com/201...pe-to-use.html


        Is anyone able to clear this up?

        Thanks
        Attached Files

        Comment


        • #19
          Tablet has a "read concordance" mode for viewing strand direction on reads. Both pairs of reads that are consistent with a mapped direction following the reference sequence are coloured green (customisable), and reads that are consistent with the reverse complement of the reference sequence are coloured red.

          Comment

          Latest Articles

          Collapse

          • seqadmin
            Exploring the Dynamics of the Tumor Microenvironment
            by seqadmin




            The complexity of cancer is clearly demonstrated in the diverse ecosystem of the tumor microenvironment (TME). The TME is made up of numerous cell types and its development begins with the changes that happen during oncogenesis. “Genomic mutations, copy number changes, epigenetic alterations, and alternative gene expression occur to varying degrees within the affected tumor cells,” explained Andrea O’Hara, Ph.D., Strategic Technical Specialist at Azenta. “As...
            07-08-2024, 03:19 PM
          • seqadmin
            Exploring Human Diversity Through Large-Scale Omics
            by seqadmin


            In 2003, researchers from the Human Genome Project (HGP) announced the most comprehensive genome to date1. Although the genome wasn’t fully completed until nearly 20 years later2, numerous large-scale projects, such as the International HapMap Project and 1000 Genomes Project, continued the HGP's work, capturing extensive variation and genomic diversity within humans. Recently, newer initiatives have significantly increased in scale and expanded beyond genomics, offering a more detailed...
            06-25-2024, 06:43 AM

          ad_right_rmr

          Collapse

          News

          Collapse

          Topics Statistics Last Post
          Started by seqadmin, 07-10-2024, 07:30 AM
          0 responses
          26 views
          0 likes
          Last Post seqadmin  
          Started by seqadmin, 07-03-2024, 09:45 AM
          0 responses
          201 views
          0 likes
          Last Post seqadmin  
          Started by seqadmin, 07-03-2024, 08:54 AM
          0 responses
          212 views
          0 likes
          Last Post seqadmin  
          Started by seqadmin, 07-02-2024, 03:00 PM
          0 responses
          193 views
          0 likes
          Last Post seqadmin  
          Working...
          X