Unconfigured Ad

Collapse
X
 
  • Time
  • Show
Clear All
new posts
  • Saeideh
    Member
    • Aug 2015
    • 25

    #1

    Trimmomatic but no trimming

    Hii dear friends

    I have a question about Trimmomatic. I wrote this command:

    java -jar /usr/bin/trimmomatic-0.33.jar PE -phred64 72_L3_1.fq.gz 72_L3_2.fq.gz 71P_1.fq.gz 71UP_1.fq.gz 71P_2.fq.gz 71UP_2.fq.gz ILLUMINACLIP:Over.fasta:2:40:15


    72_L3_1.fq.gz
    and
    72_L3_2.fq.gz
    are my input files

    71P_1.fq.gz
    71UP_1.fq.gz
    71P_2.fq.gz
    and
    71UP_2.fq.gz
    are my output files. In my case P means paired and UP means unpaired.

    after I did fastqc I got overrepresented sequence (it is not adapter) and I stored it in Over.fasta and I want to trim it from my data.

    but when I run the command I see this:

    ILLUMINACLIP: Using 0 prefix pairs, 0 forward/reverse sequences, 0 forward only sequences, 0 reverse only sequences
    Input Read Pairs: 19759652 Both Surviving: 19759652 (100.00%) Forward Only Surviving: 0 (0.00%) Reverse Only Surviving: 0 (0.00%) Dropped: 0 (0.00%)

    It does not trim the sequence.

    Why it is like that?

    In addition, for Palindrome and Single mode U just selected 40 and 15 because the manual selected them. Actually I do not know which score I should select. How can I choose the correct score?

    Would you please help me in these two questions.

    -------------------------------------------------------------------
    I'm sorry for asking so many questions each day. I' working with RNA-seq data in these days. In each level, I should stop due to errors, warnings, blocks and confusions and this story continues.

    and thank you for helping me in solving all these problems.
  • yueluo
    Member
    • Aug 2013
    • 82

    #2
    Are you sure nothing was done the your sequences? 100% reads surviving does not necessarily mean no trimming was performed, especially when you don't specify MINLEN.

    Comment

    • Saeideh
      Member
      • Aug 2015
      • 25

      #3
      Yes, I'm almost sure because when I do fastqc, the number of reads and nucleotides in them have not changed (I mean "Total sequences" and "Sequence length").

      Comment

      • kmcarr
        Senior Member
        • May 2008
        • 1181

        #4
        Originally posted by Saeideh View Post
        ILLUMINACLIP: Using 0 prefix pairs, 0 forward/reverse sequences, 0 forward only sequences, 0 reverse only sequences
        The explanation is here (highlighted in red). Trimmomatic is not finding any valid sequences to use for trimming in the supplied clipping file 'Over.fasta' meaning you must have formatted it incorrectly.

        Comment

        • Saeideh
          Member
          • Aug 2015
          • 25

          #5
          I put following sequences in Over.fasta

          CTCCCATTTCGCTCGCCGCTACTACGGGAATCGCTTTTGCTTTCTTTTCC
          GGCTTGCGGTGGATACCTAGGTACCCAGAGACGAGGAAGGGCGTAGCAAG

          Isn't it correct?

          Comment

          • yueluo
            Member
            • Aug 2013
            • 82

            #6
            That is not a valid fasta file.

            Comment

            • Saeideh
              Member
              • Aug 2015
              • 25

              #7
              Alriiiiiight...

              I got it

              Thanks all

              Comment

              Latest Articles

              Collapse

              • SEQadmin2
                Multiomic Liquid Biopsy: What DNA/RNA Panels Reveal That DNA Alone Can’t
                by SEQadmin2



                Tumors contain complex mixtures of cell populations that evolve in response to treatment and host-specific factors. Tissue biopsies provide only static snapshots of a tumor’s state. But what if we could follow that evolution over time using a minimally invasive technique?

                Liquid biopsies can detect tumor-derived material circulating in the bloodstream, providing a way to study treatment response, detect emerging resistance, and discover new biomarkers. Circulating
                ...
                Today, 09:55 AM
              • SEQadmin2
                New Genomics Technologies Take Aim at Long-Standing Limits
                by SEQadmin2


                Researchers using sequencing and genomics tools often have to make trade-offs. They can choose between speed or scale, short reads or long-range information, or targeted panels or a view of the whole transcriptome. New technologies that have been released this year are built to address those tough choices.

                We asked six companies the same four questions to learn about their latest products. The new technologies bring a lot to the table, including rethinking sequencing
                ...
                09-28-2026, 10:25 AM

              ad_right_rmr

              Collapse

              News

              Collapse

              Topics Statistics Last Post
              Started by SEQadmin2, Today, 05:39 AM
              0 responses
              9 views
              0 reactions
              Last Post SEQadmin2  
              Started by SEQadmin2, Yesterday, 09:56 AM
              0 responses
              11 views
              0 reactions
              Last Post SEQadmin2  
              Started by SEQadmin2, 10-09-2026, 11:14 AM
              0 responses
              25 views
              0 reactions
              Last Post SEQadmin2  
              Started by SEQadmin2, 09-29-2026, 09:51 AM
              0 responses
              76 views
              0 reactions
              Last Post SEQadmin2  
              Working...