Seqanswers Leaderboard Ad

Collapse

Announcement

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

  • How to trim dual index adapters with Trim Galore!

    Hi everybody,
    I am trying to trim some adapters with Trim Galore from my bisulfite sequencing data, but I don't know what command to use.
    It has been really difficult to take the information of library construction and I only know that it has been done using IDT unique dual indexing adapters in paired end data, detailed as follows:

    1 - P5 ATATGCGC, P7 CTGATCGT
    30 - P5 AACGTCTG, P7 GCGTCATT
    31 - P5 TATCGGTC, P7 ATCCGGTA
    32 - P5 CGCTCTAT, P7 CGTTGCAA

    When doing the FASTQC analysis I have found some overrepresented sequences such as:
    - In one of the sequence of the pair
    AAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAA
    - In the other sequence of the pair
    GGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGG
    TTTTTTTTTTTTTTTTTTTTTTTTTTTTTTTTTTTTTTTTTTTTTTTTTT

    An also in the adapter content section the Illumina Universal adapters are identified

    Does all this make sense?
    So I need to add something more to those P5 and P7 adapter sequences when running Trim Galore?

    Will this command be correct for the first example of unique dual index adapters?

    > trim_galore --paired -a ATATGCGC -a2 CTGATCGT -q 15 --stringency 5 -e 0.05 --length 50 --fastqc --gzip file_1.fastq.gz file_2.fastq.gz

    Thanks in advance!

  • #2
    Index sequences are never part of an actual read in Illumina sequencing so you are never going to find index sequences you noted above in R1/R2 reads.

    As for the poly-A and ploy-G (no signal, two color chemistry) you just need to remove them during trimming.

    Comment


    • #3
      Thanks for your answer!
      As the only information of the sequencing service was that about the dual index adapters I though that those were the sequences that I need to trim.
      So for the command how do I need to add this poly A, T and Gs?
      One by one in each sequence file? or combining the paired end sequences somehow?

      Thanks again!

      Comment


      • #4
        I suggest you try "bbduk.sh" from BBMap suite instead. A guide is available. You can trim out all these things in one pass.

        I am not a trim galore user but you may be able to additional specify poly-A and poly-G there as well.

        Comment

        Latest Articles

        Collapse

        • seqadmin
          Understanding Genetic Influence on Infectious Disease
          by seqadmin




          During the COVID-19 pandemic, scientists observed that while some individuals experienced severe illness when infected with SARS-CoV-2, others were barely affected. These disparities left researchers and clinicians wondering what causes the wide variations in response to viral infections and what role genetics plays.

          Jean-Laurent Casanova, M.D., Ph.D., Professor at Rockefeller University, is a leading expert in this crossover between genetics and infectious...
          09-09-2024, 10:59 AM
        • seqadmin
          Addressing Off-Target Effects in CRISPR Technologies
          by seqadmin






          The first FDA-approved CRISPR-based therapy marked the transition of therapeutic gene editing from a dream to reality1. CRISPR technologies have streamlined gene editing, and CRISPR screens have become an important approach for identifying genes involved in disease processes2. This technique introduces targeted mutations across numerous genes, enabling large-scale identification of gene functions, interactions, and pathways3. Identifying the full range...
          08-27-2024, 04:44 AM

        ad_right_rmr

        Collapse

        News

        Collapse

        Topics Statistics Last Post
        Started by seqadmin, Today, 06:25 AM
        0 responses
        13 views
        0 likes
        Last Post seqadmin  
        Started by seqadmin, Yesterday, 01:02 PM
        0 responses
        12 views
        0 likes
        Last Post seqadmin  
        Started by seqadmin, 09-18-2024, 06:39 AM
        0 responses
        14 views
        0 likes
        Last Post seqadmin  
        Started by seqadmin, 09-11-2024, 02:44 PM
        0 responses
        14 views
        0 likes
        Last Post seqadmin  
        Working...
        X