Seqanswers Leaderboard Ad

Collapse

Announcement

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

  • How do I find medium length somatic indels?

    I used IndelGenotyperV2 to find somatic indels. But the longest I find is 4bp. My reads are single 75bp exome reads. I am interested in indels around 20bp. Is that possible with my reads?

    I ran the mapping step with bwa-0.7.5 on default setting. Is that the cause? I heard that the -o and -e flags for bwa aln supposedly have something to do with indels. If so, how do I set it up to find indels around 20bp?

    Thanks a lot in advance!

  • #2
    try Pindel. there is no theoretical maximum deletion size. currently maximum insertion size for insertions is read length - 20.

    Comment


    • #3
      For single-end reads, try novoalign. You may also consider bwa-mem with (higher matching score):

      bwa mem -A2 -B3 -O5 -E1

      bwa-backtrack does not find a gap more than 4bp for single-end reads. Pindel would not work with single-end, as I remember (could be wrong).

      Comment


      • #4
        Strelka does a fantastic job of finding longer somatic indels: https://sites.google.com/site/strelk...variantcaller/

        For exome data, set the "isSkipDepthFilters" to 1 in order to skip depth filtration.

        Comment


        • #5
          I tried "bwa aln -e 50" and it seems to work. Is bwa_mem better than bwa for this? How so?

          Comment


          • #6
            Originally posted by ymc View Post
            I tried "bwa aln -e 50" and it seems to work. Is bwa_mem better than bwa for this? How so?
            Yes, that may work, but bwa-mem should be better, at least faster. Bwa and bwa-mem are completely different algorithms.

            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
            11 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
            51 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