Seqanswers Leaderboard Ad

Collapse

Announcement

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

  • New k-mers overpresentation by adapter trimming

    Hi
    I am working on an Illumina Hiseq paired end whole genome data.

    The data shows overrepresentation of adapter sequences such as TruSeq Adapter, Index 14 (97% over 40bp) according to the quality reported by the fastqc analysis.

    So I performed trimming of the reads by the TruSeq-PE adapter sequences in trimmomatic data and performed fastqc after trimming.

    The post-trimming data shows appearance of new k-mers overrepresented at various positions within 0-50 bp. The bbmap analysis for the percentage match still shows no remarkable mismatch problems, indicating that the k-mers could be genomic.

    I have attached the before and after trimming k-mer overpresentation plots herewith.

    This observation makes me doubt the need of trimming in paired end data.

    My question is :
    Do we indeed need to perform trimming in such cases ?
    Any suggestion on change in the trimming was needed such as using only the overrepresented adapter sequences and but not other adapters from the same library in the analysis.

    Thanks
    Attached Files
    Last edited by ty23991; 06-22-2015, 11:28 AM.

  • #2
    Considering the Y-axis is log-scale, it looks like trimming vastly reduced the rate of overrepresented kmers (though I've never been sure exactly how to interpret the graph), which would be a good thing. What kind of trim rate did you get?

    Also, note that if you have contamination by synthetic artifacts of set length, that will not affect BBMap's error rate histograms because the fully-artifact reads won't map to your reference.

    Comment


    • #3
      Hi Brian
      Thanks for the suggestion.

      Here is the trimming report:

      Trimming was performed by using 1 prefix pairs, 10 forward/reverse sequences

      Input Read Pairs: 532073495
      Both Surviving: 510614940 (~96 %)
      Forward Only Surviving: 13002991 (2.4%)
      Reverse Only Surviving: 3487288 (0.7%)
      Dropped: 4968275 (0.9%)

      Comment


      • #4
        It's hard to extrapolate from that to derive the number of reads/bases that were trimmed due to adapter sequence, but it looks like trimming was probably a good idea in that case.

        Comment


        • #5
          Thanks for the idea. I am thinking to perform the alignment on both with and without trimming and compare the results.

          Comment

          Latest Articles

          Collapse

          • seqadmin
            The Impact of AI in Genomic Medicine
            by seqadmin



            Artificial intelligence (AI) has evolved from a futuristic vision to a mainstream technology, highlighted by the introduction of tools like OpenAI's ChatGPT and Google's Gemini. In recent years, AI has become increasingly integrated into the field of genomics. This integration has enabled new scientific discoveries while simultaneously raising important ethical questions1. Interviews with two researchers at the center of this intersection provide insightful perspectives into...
            02-26-2024, 02:07 PM
          • seqadmin
            Multiomics Techniques Advancing Disease Research
            by seqadmin


            New and advanced multiomics tools and technologies have opened new avenues of research and markedly enhanced various disciplines such as disease research and precision medicine1. The practice of merging diverse data from various ‘omes increasingly provides a more holistic understanding of biological systems. As Maddison Masaeli, Co-Founder and CEO at Deepcell, aptly noted, “You can't explain biology in its complex form with one modality.”

            A major leap in the field has
            ...
            02-08-2024, 06:33 AM

          ad_right_rmr

          Collapse

          News

          Collapse

          Topics Statistics Last Post
          Started by seqadmin, 02-23-2024, 04:11 PM
          0 responses
          57 views
          0 likes
          Last Post seqadmin  
          Started by seqadmin, 02-21-2024, 08:52 AM
          0 responses
          67 views
          0 likes
          Last Post seqadmin  
          Started by seqadmin, 02-20-2024, 08:57 AM
          0 responses
          56 views
          0 likes
          Last Post seqadmin  
          Started by seqadmin, 02-14-2024, 09:19 AM
          0 responses
          65 views
          0 likes
          Last Post seqadmin  
          Working...
          X