Seqanswers Leaderboard Ad

Collapse

Announcement

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

  • adapter trimming miRNA sequencing

    Below is the output from FastQC. Any idea what to trim for? Tried the Illumina Single End Adapter 1 but got no mapping. What are the A, and how do I get rid of them?


    ATCTCGTATGCCGTCTTCTGCTTGAAAAAAAAAAAA 4805262 44.91006094070466 Illumina Single End Adapter 1 (95% over 24bp)
    GTATGCCGTCTTCTGCTTGAAAAAAAAAAAAAAAAA 1107781 10.353340196425243 No Hit
    GTATGCCGTCTTCTGCTTGAAAAAAAAAAAAAAACA 200679 1.8755493705691118 No Hit
    CGATCATCTCGTATGCCGTCTTCTGCTTGAAAAAAA 144914 1.3543687256098158 Illumina PCR Primer Index 7 (96% over 29bp)
    TAGCTTATCAGACTGATGTTGACATCTCGTATGCCG 140931 1.3171435394021074 No Hit
    GTATGCCGTCTTCTGCTTGAAAAAAAAAAAAAACAA 124510 1.1636725921972906 No Hit
    TGTAAACATCCTTGACTGGAAGCATCTCGTATGCCG 114871 1.0735863331322382 No Hit
    TCGTATGCCGTCTTCTGCTTGAAAAAAAAAAAAAAA 112198 1.048604429357896 Illumina Single End Adapter 1 (100% over 21bp)
    TCTCGTATGCCGTCTTCTGCTTGAAAAAAAAAAAAA 110008 1.0281366518547872 Illumina Single End Adapter 1 (95% over 23bp)
    ATCTCGTATGCCGTCTTCTGCTTGAAAAAAAAAAAC 94442 0.8826565492915952 Illumina Single End Adapter 1 (95% over 24bp)
    GTATGCCGTCTTCTGCTTGAAAAAAAAAAAAAAAAC 82271 0.7689061748667841 No Hit
    TGAGGTAGTAGGTTGTATAGTTATCTCGTATGCCGT 77257 0.7220452450035024 No Hit
    GTAGTGTTTCCTACTTTATGGAATCTCGTATGCCGT 62069 0.5800979369134498 No Hit
    GTATGCCGTCTTCTGCTTGAAAAAAAAAAAAACAAA 58659 0.5482280185181984 No Hit

  • #2
    Hi Palgrave,
    Have you visited here -
    Oligonucleotide (oligo) sequences of Illumina adapters used in AmpliSeq, Nextera, TruSeq, and TruSight library prep kits.

    The pdf file enlists all adapter sequences. If you copy the first entry in your list (without the polyA tail) and search in the pdf, you will see that its the TruSeq adapter.

    I wouldn't bother about the polyA and just give the remaining sequence as a contaminant to any adapter trimmer tool. When the tool chops, the polyA's get removed too. I use Trimmomatic (can use multiple threads).

    Comment


    • #3
      I've used trimmomatic with the adapters file that comes along with fastqc converted to a fasta file. Worked good. The As might be after it reads all the way through the adapter and isn't reading anything, maybe?

      Comment


      • #4
        try skewer

        Originally posted by Palgrave View Post
        Below is the output from FastQC. Any idea what to trim for? Tried the Illumina Single End Adapter 1 but got no mapping. What are the A, and how do I get rid of them?


        ATCTCGTATGCCGTCTTCTGCTTGAAAAAAAAAAAA 4805262 44.91006094070466 Illumina Single End Adapter 1 (95% over 24bp)
        GTATGCCGTCTTCTGCTTGAAAAAAAAAAAAAAAAA 1107781 10.353340196425243 No Hit
        GTATGCCGTCTTCTGCTTGAAAAAAAAAAAAAAACA 200679 1.8755493705691118 No Hit
        CGATCATCTCGTATGCCGTCTTCTGCTTGAAAAAAA 144914 1.3543687256098158 Illumina PCR Primer Index 7 (96% over 29bp)
        TAGCTTATCAGACTGATGTTGACATCTCGTATGCCG 140931 1.3171435394021074 No Hit
        GTATGCCGTCTTCTGCTTGAAAAAAAAAAAAAACAA 124510 1.1636725921972906 No Hit
        TGTAAACATCCTTGACTGGAAGCATCTCGTATGCCG 114871 1.0735863331322382 No Hit
        TCGTATGCCGTCTTCTGCTTGAAAAAAAAAAAAAAA 112198 1.048604429357896 Illumina Single End Adapter 1 (100% over 21bp)
        TCTCGTATGCCGTCTTCTGCTTGAAAAAAAAAAAAA 110008 1.0281366518547872 Illumina Single End Adapter 1 (95% over 23bp)
        ATCTCGTATGCCGTCTTCTGCTTGAAAAAAAAAAAC 94442 0.8826565492915952 Illumina Single End Adapter 1 (95% over 24bp)
        GTATGCCGTCTTCTGCTTGAAAAAAAAAAAAAAAAC 82271 0.7689061748667841 No Hit
        TGAGGTAGTAGGTTGTATAGTTATCTCGTATGCCGT 77257 0.7220452450035024 No Hit
        GTAGTGTTTCCTACTTTATGGAATCTCGTATGCCGT 62069 0.5800979369134498 No Hit
        GTATGCCGTCTTCTGCTTGAAAAAAAAAAAAACAAA 58659 0.5482280185181984 No Hit
        Hi Palgrave,

        From the pdf file provided by illimina, you will find ATCTCGTATGCCGTCTTCTGCTTG is actually the v1.5 Small RNA 3' Adapter. Never mind those tailing poly-A as they are fake base-calls by the illumina sequencer if there is no base detected. A mature adapter trimmer will trim those tailing poly-A next to adapter sequence.

        Though trimmomatic seems to be specifically optimized for small RNA adapter trimming, I also recommend you using the skewer (http://sourceforge.net/projects/skewer ) tool which is not optimized for any specific applications but performs equally well or even better for small RNA adapter trimming.

        For your case, the command is very simple like this:
        Code:
        skewer -x ATCTCGTATGCCGTCTTCTGCTTG -l 16 -L 30 yourData.fastq

        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
        10 views
        0 likes
        Last Post seqadmin  
        Started by seqadmin, Yesterday, 06:07 PM
        0 responses
        10 views
        0 likes
        Last Post seqadmin  
        Started by seqadmin, 03-22-2024, 10:03 AM
        0 responses
        51 views
        0 likes
        Last Post seqadmin  
        Started by seqadmin, 03-21-2024, 07:32 AM
        0 responses
        67 views
        0 likes
        Last Post seqadmin  
        Working...
        X