Seqanswers Leaderboard Ad

Collapse

Announcement

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

  • reads align but contigs do not??

    Here is something that puzzles me:
    I have 4000 reads that align to my reference genome. I can visualize them, they cover maybe 80 % or a little more of the reference sequence.
    I assemble these reads into contigs, get something like 160 contigs. ok.
    however, when I try to align these contigs, only 27 of them still match my reference.
    huh? how is that possible?
    Any thoughts?

  • #2
    Have you checked to see how much duplication you have in your data? 4000 data tends to have more duplicates if you don't have the right library insert size. Use clumpify.sh from BBMap to get an estimate of total duplicates and optical duplicates. You may not be getting a good assembly if there are too many duplicates.

    Comment


    • #3
      Additionally, lengths and coverage are important here. Possibly, 99% of the bases map to the 27 contigs that match the reference, and the other contigs are short, degenerate contigs from errors or chimerism. Or maybe you're assembling various adapter sequences (another form of chimerism, in a way). It depends on how you're doing the alignment and assembly, and what kind of data you have. Can you describe the situation in more detail? Also, have you BLASTed the unaligned sequences, or compared them to known adapter sequences?

      Comment


      • #4
        these are single end reads from an RNA sequencing project. adapter sequences have been removed.
        out of some 37 mil reads I get 4000 reads reads aligning to a just under 8000bp long virus that we know is present in the sample. But these are not siRNA in 20nt length range, these reads range from about 30 to 100 nt long.
        I am sure there is duplication, though I do not know how to find or eliminate that and would appreciate any hints as to where to read about this or find tools to work with this.
        out of about 160 contigs, 40 are very short, the same length as the shortest reads in fact.
        the reads matching the virus (bwa aligned, IGV visualized) seem to visually cover most of the virus sequence, the 27 contigs cover maybe 30 % of the virus sequence.
        Thanks for the help !

        Comment


        • #5
          Originally posted by sfh838t View Post
          I am sure there is duplication, though I do not know how to find or eliminate that and would appreciate any hints as to where to read about this or find tools to work with this.
          @Brian's post gives directions on how to use clumpify.sh from BBMap to identify, mark/eliminate duplicates. They can be optical or all (PCR+ optical).

          Comment


          • #6
            I did see Brians post, but did not check on that so far. Looks like I will have to take the plunge and install yet another toolbox. Let's see what will go wrong this time.....

            Comment


            • #7
              Nothing should go wrong. As long as you have Java 1.7 or greater available nothing else is needed.

              Comment


              • #8
                Ok, I did install and use this clumpify tool.
                and it did not make any difference. it found some 300 duplicates, I assemble slightly fewer reads, to similar n50 and total nt count and of the contigs produced I had 1 less contig match than before.....
                I have heard the term redundant vs non-redundant. just to make sure I understand: removing these duplicates would make my set of reads non-redundant? Or am I completely lost?

                Comment


                • #9
                  What setting did you use for dupedist? For a 4000 flowcell data you need to include
                  Code:
                  dedupe=t optical=t dupedist=2500 dupesubs=0 spantiles=f
                  to your clumpify command line. (dupesubs=0 will only look for perfectly matching reads without any errors).

                  If you did not use even one of those options then I suggest that you try again.

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