Seqanswers Leaderboard Ad

Collapse

Announcement

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

  • Number of Reads Increasing

    Hi there,
    I trimmed my reads and then merged them using FLASh. However, I found that the number of reads before processing (trim and merge) is lower than the number of reads after processing.
    Is there an explanation for this?
    Thanks

  • #2
    When you merge a read pair, you turn 2 reads into 1 read. Therefore, you should expect the number of reads to decrease.

    Comment


    • #3
      i know i should expect a decrease but i'm getting an increase. do you know why?

      Comment


      • #4
        Perhaps you can post your exact commands used at each step, as well as how you are calculating the number of reads?

        Comment


        • #5
          First step I'm doing is using Trimmomatic to trim and remove my adapters using the PE setting:

          java -jar trimmomatic-0.36.jar PE File1_L001_R1_001.fastq File1_L001_R2_001.fastq /home/rachelle/File1_R1paired.fastq File1_R1unpaired.fastq File1_R2paired.fastq File1_R2unpaired.fastq ILLUMINACLIP:/usr/local/trimmomatic/adapters/NexteraPE-PE.fa:2:30:10:6 LEADING:3 TRAILING:3 SLIDINGWINDOW:4:25 MINLEN:36

          Afterwards, I take R1paired and R2 paired and merge them using FLASh with a minum of overlap of 15.
          I then concatenate the out.extended frags with the two notcombined files for a final file that is then concatenated with R1UNpaired and R2UNpaired.

          That final file has more reads than the beginning.
          I looked into it some more and saw that my reads are increasing after FLASh. I am assuming it has to do with concatenation of the two notCombined files. Does anyone know what those two files are? And if I should be concatenating them with my out.extendedfrags file?

          Thanks

          Comment


          • #6
            Best option would be to merge first (use bbmerge.sh from BBMap) and then trim (bbduk.sh from BBMap).

            The only way number of reads is going to increase is if you double dip into the read pool. You must be somehow doing that. Since your merged reads become "single-end" you should not merge them with remaining PE reads.
            Last edited by GenoMax; 01-19-2017, 10:43 AM.

            Comment


            • #7
              Paired-end reads are typically counted once (unless it's the Illumina marketing team). Merged reads will also be counted once after FLASH, but unmerged reads now become two single-end reads and are counted twice.

              Comment


              • #8
                that makes perfect sense!!! Thank you so much

                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
                67 views
                0 likes
                Last Post seqadmin  
                Working...
                X