Seqanswers Leaderboard Ad

Collapse

Announcement

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

  • alignment reads in homologous region

    Hi NGS users,
    I'm analyzing SureSelect sequencing
    I noticed that BWA/SAMTOOLS in mapping step, align reads in target genes and also in regions that are homologous to my target gene but the are not in baits.
    does anyone knows how correct this bias? This unappropriate alignment could affect the SNP calling.
    Thanx a lot,
    bye!

  • #2
    Hi. I have never done capture, but we will very soon, so I though about these problems. I hope my contribution is useful.

    You could use, as "reference genome", a file that only contains your targets.

    however, the pull down is probably not 100% on target, so this might introduce artifacts

    The baits, on the other hand, might actually pull down a region very similar to their desired target, so it might be that your sequence comes from a homologous region, and the aligner get it right.

    Comment


    • #3
      Thanx a lot Stefano.
      I have already tried to align my reads against my target sequence (using sure select design as reference genome) but for experience, I can say that it forces too much the alignment with high risk of false positive.

      I don't know if my problem is at an experimental level (in the enrichment protocol) or in the mapping, because it's not a problem of repetitive regions but a problem in enrichment for homology.

      Comment


      • #4
        just to know, what sort of similarity is there between a homologous region and your target?

        Comment


        • #5
          Limiting your reference is a bad idea as similar regions differing by a single base will look like a SNP if both regions are not included in the reference.

          A common approach to dealing with this is to scale the quality of your alignment based on the uniqueness of the mapping (following the GATK best practices will do this).

          While some the the off-target effect you are seeing is likely due to incorrect mapping, because you are enriching for these regions by hybridization, there is likely off-target capture as well.

          Comment


          • #6
            about 94% of homology

            Comment


            • #7
              Adamdeluca, you are right!
              in my alignment pipeline, I used always GATk realign and recalibration score. But I can observe reads mapped in homologuos region in the final bam.
              So...the question is: it is an experimental problem of enrichment for homology (and in this case I cannot do anything) or there is a way (parameter-script...) to understand when I map reads off target in homologuos genes?

              And, if the problem is an incorrect mapping (I use BWA and Samtools for this step), how can I limit the mapping to other region? In some case, mapping quality of the reads is very low, so I can discriminate off target region, in other cases I have a high mapping quality and a high coverage, so I can't differentiate target and off target (with the risk of false positive in SNP calling)

              Comment


              • #8
                Originally posted by m_elena_bioinfo View Post
                Adamdeluca, you are right!
                in my alignment pipeline, I used always GATk realign and recalibration score. But I can observe reads mapped in homologuos region in the final bam.
                So...the question is: it is an experimental problem of enrichment for homology (and in this case I cannot do anything) or there is a way (parameter-script...) to understand when I map reads off target in homologuos genes?

                And, if the problem is an incorrect mapping (I use BWA and Samtools for this step), how can I limit the mapping to other region? In some case, mapping quality of the reads is very low, so I can discriminate off target region, in other cases I have a high mapping quality and a high coverage, so I can't differentiate target and off target (with the risk of false positive in SNP calling)
                great questions. Let me know if you figure them out.

                By manually inspecting reads in these types of regions I can observe obvious miss-mapping, and obvious off-target capture.

                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
                12 views
                0 likes
                Last Post seqadmin  
                Started by seqadmin, Yesterday, 06:07 PM
                0 responses
                10 views
                0 likes
                Last Post seqadmin  
                Started by seqadmin, 03-22-2024, 10:03 AM
                0 responses
                52 views
                0 likes
                Last Post seqadmin  
                Started by seqadmin, 03-21-2024, 07:32 AM
                0 responses
                68 views
                0 likes
                Last Post seqadmin  
                Working...
                X