Unconfigured Ad

Collapse
X
 
  • Time
  • Show
Clear All
new posts
  • idedios
    Member
    • Mar 2014
    • 20

    #1

    bwa mis-alignment of TruSeq Amplicon reads



    I'm having the following issues using bwa mem on paired-end TruSeq Amplicon reads. I used the default settings with mem. The top sequence is from bwa mem. The bottom is from Illumina's BaseSpace analysis.
  • Bukowski
    Senior Member
    • Jan 2010
    • 388

    #2
    Cross posted to Biostars: https://www.biostars.org/p/153442/

    Comment

    • idedios
      Member
      • Mar 2014
      • 20

      #3
      Any comments?

      Comment

      • idedios
        Member
        • Mar 2014
        • 20

        #4
        So based on a reply from biostars it looks like the BaseSpace reads have been trimmed. Now I just need to understand how they trimmed it. This is desirable since my alignment has led to the variant caller calling a lot of false indels.

        It's possible that Illumina is trimming the reads based on the regions in their manifest file. If I could do the same then maybe I could get the reads I'm looking for.

        Comment

        • GenoMax
          Senior Member
          • Feb 2008
          • 7142

          #5
          You are likely doing an apples to oranges comparison (especially since you don't know the complete details of how BaseSpace analysis is being done). BaseSpace is also likely using isaac/isaac2 for the alignments (or do you know that they are using bwa).

          What is it that you are trying to do? Replicate BaseSpace analysis locally (with same software/parameters)?

          Comment

          • Brian Bushnell
            Super Moderator
            • Jan 2014
            • 2709

            #6
            Is it even the same data? The reads appear to be different lengths.

            Comment

            • idedios
              Member
              • Mar 2014
              • 20

              #7
              The reads are from the same fastqs. I'm just trying to figure out why my alignment is looking so messy when the TruSeq Amplicon stuff is supposed to have strictly sized blocks of reads. It is possible they are using the isaac aligner so I'll have to do more tinkering to figure out what is going on.
              I'm doing this to move my BaseSpace analyses in-house using bcbio-nextgen.

              Comment

              • GenoMax
                Senior Member
                • Feb 2008
                • 7142

                #8
                Originally posted by idedios View Post
                I'm doing this to move my BaseSpace analyses in-house using bcbio-nextgen.
                Is this what you are planning to use: https://github.com/chapmanb/bcbio-nextgen? Does not seem to include a trimming program at first glance (which you should plan to include).

                Comment

                • idedios
                  Member
                  • Mar 2014
                  • 20

                  #9
                  Originally posted by GenoMax View Post
                  Is this what you are planning to use: https://github.com/chapmanb/bcbio-nextgen? Does not seem to include a trimming program at first glance (which you should plan to include).
                  Yeah this is exactly what I am using. I could ask Brad Chapman about a trimming program.

                  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
                    ...
                    10-14-2026, 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, Yesterday, 11:13 AM
                  0 responses
                  17 views
                  0 reactions
                  Last Post SEQadmin2  
                  Started by SEQadmin2, Yesterday, 10:48 AM
                  0 responses
                  18 views
                  0 reactions
                  Last Post SEQadmin2  
                  Started by SEQadmin2, 10-15-2026, 10:52 AM
                  0 responses
                  36 views
                  0 reactions
                  Last Post SEQadmin2  
                  Started by SEQadmin2, 10-14-2026, 05:39 AM
                  0 responses
                  42 views
                  0 reactions
                  Last Post SEQadmin2  
                  Working...