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
                    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
                  • seqadmin
                    The Impact of AI in Genomic Medicine
                    by seqadmin



                    Artificial intelligence (AI) has evolved from a futuristic vision to a mainstream technology, highlighted by the introduction of tools like OpenAI's ChatGPT and Google's Gemini. In recent years, AI has become increasingly integrated into the field of genomics. This integration has enabled new scientific discoveries while simultaneously raising important ethical questions1. Interviews with two researchers at the center of this intersection provide insightful perspectives into...
                    02-26-2024, 02:07 PM

                  ad_right_rmr

                  Collapse

                  News

                  Collapse

                  Topics Statistics Last Post
                  Started by seqadmin, 03-14-2024, 06:13 AM
                  0 responses
                  34 views
                  0 likes
                  Last Post seqadmin  
                  Started by seqadmin, 03-08-2024, 08:03 AM
                  0 responses
                  72 views
                  0 likes
                  Last Post seqadmin  
                  Started by seqadmin, 03-07-2024, 08:13 AM
                  0 responses
                  81 views
                  0 likes
                  Last Post seqadmin  
                  Started by seqadmin, 03-06-2024, 09:51 AM
                  0 responses
                  68 views
                  0 likes
                  Last Post seqadmin  
                  Working...
                  X