Seqanswers Leaderboard Ad

Collapse

Announcement

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

  • Quality trimming for RNA Seq

    When I trim sequences for RNASeq, can I trim by sliding window on quality score or will the fact that the Illumina reads are now of different lengths affect the FPKM values (using cufflinks)? Should I just do a 3' and 5' trim of a specific length and leave all lengths equal?

  • #2
    Not sure if there are any papers on this yet, but I don't think so.

    I tried the window based trimmer ea-utils on my RNAseq datasets and realigned, getting extremely similar results to my untrimmed datasets.
    I should mention this is on very short read data (35-40bp).

    For longer reads with low quality ends this is of course more of an issue.

    Comment


    • #3
      This is unless the aligner takes base-call quality scores into account, in which case trimming should be unnecessary anyway.

      Comment


      • #4
        Can you please elaborate on that? Why would the trimming be unnecessary?

        Comment


        • #5
          The point of trimming low-quality base-calls is to remove bases which might be erroneous. A naive aligner would be deterred from a valid mapping by mismatches that are only due to sequencing errors. A quality-aware aligner, however, knows that a mismatch of a low-quality base is no reason to reject an otherwise good mapping and will report it. Many of the currently popular aligners work this way, i.e., they pay less attention to the low-quality ends of reads. By trimming, you basically take away from the aligner the chance to make use of this feature. However, trimming makes a hard-cut decision (everything below a quality threshold is removed) while a well-designed quality-aware aligner may make a more sophisticated, gradual, decision.
          Last edited by Simon Anders; 03-01-2012, 09:57 AM. Reason: reworded

          Comment


          • #6
            Thanks for the answer- do you know which aligners work in this way? (I am working with bowtie and cufflinks)
            The reason we were thinking of trimming is because we are working with mixed communities (co-cultures) so we may have two bacteria in the same Illumina. We then align all the data to each of two genomes, allowing zero mismatches.
            Given this, would your answer still be not to trim?
            And just for my general knowledge- do you know if different length reads if we did trim would affect the final FPKM calculations?
            Thanks!

            Comment

            Latest Articles

            Collapse

            • seqadmin
              Strategies for Sequencing Challenging Samples
              by seqadmin


              Despite advancements in sequencing platforms and related sample preparation technologies, certain sample types continue to present significant challenges that can compromise sequencing results. Pedro Echave, Senior Manager of the Global Business Segment at Revvity, explained that the success of a sequencing experiment ultimately depends on the amount and integrity of the nucleic acid template (RNA or DNA) obtained from a sample. “The better the quality of the nucleic acid isolated...
              03-22-2024, 06:39 AM
            • seqadmin
              Techniques and Challenges in Conservation Genomics
              by seqadmin



              The field of conservation genomics centers on applying genomics technologies in support of conservation efforts and the preservation of biodiversity. This article features interviews with two researchers who showcase their innovative work and highlight the current state and future of conservation genomics.

              Avian Conservation
              Matthew DeSaix, a recent doctoral graduate from Kristen Ruegg’s lab at The University of Colorado, shared that most of his research...
              03-08-2024, 10:41 AM

            ad_right_rmr

            Collapse

            News

            Collapse

            Topics Statistics Last Post
            Started by seqadmin, Yesterday, 06:37 PM
            0 responses
            8 views
            0 likes
            Last Post seqadmin  
            Started by seqadmin, Yesterday, 06:07 PM
            0 responses
            8 views
            0 likes
            Last Post seqadmin  
            Started by seqadmin, 03-22-2024, 10:03 AM
            0 responses
            49 views
            0 likes
            Last Post seqadmin  
            Started by seqadmin, 03-21-2024, 07:32 AM
            0 responses
            66 views
            0 likes
            Last Post seqadmin  
            Working...
            X