Seqanswers Leaderboard Ad

Collapse

Announcement

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

  • 0% duplicates in RNA-Seq/Drop-seq library

    Hi all,

    I recently sequenced a cDNA Drop-Seq library on a HiSeq2500 and later to get more reads per cell on a HiSeq4000. I am a bit confused when looking at the "Sequence dublication level" plots of the FastQC reports:
    Percent of seqs remaining if deduplicated: HiSeq2500 25% (about 25Mio reads in total) and HiSeq4000 100% (100Mio reads in total)

    Could some one explain me this? So for my understanding, it means that the HiSeq4000 run somehow doesn't have sequencing duplicates - which is super-unlikely given the fact that a Drop-Seq library is generated by quite some PCR cycles. Also, if anything I would have expected more duplicates as the same library was just sequenced deeper.
    I don't know whether it matters, but cycle length was 180 in HiSeq2500 and 70 in HiSeq4000. The Q30 mean seq quality was 85% in HiSeq2500 and 65% in HiSeq4000.
    Any ideas on this issue are more than welcome. Thanks
    Attached Files
    Last edited by Jaeb; 05-28-2017, 12:55 AM.

  • #2
    I agree with you that with deeper sequencing %duplicate should increase and read length differences are less likely to be the cause as FastQC uses initial 50 sequence of a subset of reads for duplicate calculation.

    It would be helpful if you could post the whole FastQC report for both runs as other plots might give some clues about the cause.

    Comment


    • #3
      Thanks for your fast reply. I did attach now the complete FastQC reports....

      Comment


      • #4
        I think HS4000 reads contain lots of errors due to positional lower quality so the sequences of duplicates do not match and they are reported as unique reads. Also lots of reads seems to have very low quality over the whole length of read. If you trim or filter low quality reads you should get similar duplication rate for both runs.

        Comment


        • #5
          I suggest that you run clumpify.sh from BBMap to get an exact idea of the duplication. You can allow for errors when doing the sequence match. FastQC does not look at the entire dataset for some of the modules (only a % of data is sampled).

          Even though there is a thread for clumpify here the one over at Biostars has the directions clearly defined on one page.

          Comment

          Latest Articles

          Collapse

          • seqadmin
            Genetic Variation in Immunogenetics and Antibody Diversity
            by seqadmin



            The field of immunogenetics explores how genetic variations influence immune responses and susceptibility to disease. In a recent SEQanswers webinar, Oscar Rodriguez, Ph.D., Postdoctoral Researcher at the University of Louisville, and Ruben Martínez Barricarte, Ph.D., Assistant Professor of Medicine at Vanderbilt University, shared recent advancements in immunogenetics. This article discusses their research on genetic variation in antibody loci, antibody production processes,...
            Yesterday, 07:24 PM
          • seqadmin
            Choosing Between NGS and qPCR
            by seqadmin



            Next-generation sequencing (NGS) and quantitative polymerase chain reaction (qPCR) are essential techniques for investigating the genome, transcriptome, and epigenome. In many cases, choosing the appropriate technique is straightforward, but in others, it can be more challenging to determine the most effective option. A simple distinction is that smaller, more focused projects are typically better suited for qPCR, while larger, more complex datasets benefit from NGS. However,...
            10-18-2024, 07:11 AM

          ad_right_rmr

          Collapse

          News

          Collapse

          Topics Statistics Last Post
          Started by seqadmin, 11-01-2024, 06:09 AM
          0 responses
          24 views
          0 likes
          Last Post seqadmin  
          Started by seqadmin, 10-30-2024, 05:31 AM
          0 responses
          21 views
          0 likes
          Last Post seqadmin  
          Started by seqadmin, 10-24-2024, 06:58 AM
          0 responses
          25 views
          0 likes
          Last Post seqadmin  
          Started by seqadmin, 10-23-2024, 08:43 AM
          0 responses
          56 views
          0 likes
          Last Post seqadmin  
          Working...
          X