Header Leaderboard Ad

Collapse

How do I find medium length somatic indels?

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
              Improved Targeted Sequencing: A Comprehensive Guide to Amplicon Sequencing
              by seqadmin



              Amplicon sequencing is a targeted approach that allows researchers to investigate specific regions of the genome. This technique is routinely used in applications such as variant identification, clinical research, and infectious disease surveillance. The amplicon sequencing process begins by designing primers that flank the regions of interest. The DNA sequences are then amplified through PCR (typically multiplex PCR) to produce amplicons complementary to the targets. RNA targets...
              Today, 01:49 PM
            • seqadmin
              Targeted Sequencing: Choosing Between Hybridization Capture and Amplicon Sequencing
              by seqadmin




              Targeted sequencing is an effective way to sequence and analyze specific genomic regions of interest. This method enables researchers to focus their efforts on their desired targets, as opposed to other methods like whole genome sequencing that involve the sequencing of total DNA. Utilizing targeted sequencing is an attractive option for many researchers because it is often faster, more cost-effective, and only generates applicable data. While there are many approaches...
              03-10-2023, 05:31 AM
            • seqadmin
              Expert Advice on Automating Your Library Preparations
              by seqadmin



              Using automation to prepare sequencing libraries isn’t a new concept, and most researchers are aware that there are numerous benefits to automating this process. However, many labs are still hesitant to switch to automation and often believe that it’s not suitable for their lab. To combat these concerns, we’ll cover some of the key advantages, review the most important considerations, and get real-world advice from automation experts to remove any lingering anxieties....
              02-21-2023, 02:14 PM

            ad_right_rmr

            Collapse

            News

            Collapse

            Topics Statistics Last Post
            Started by seqadmin, 03-17-2023, 12:32 PM
            0 responses
            12 views
            0 likes
            Last Post seqadmin  
            Started by seqadmin, 03-15-2023, 12:42 PM
            0 responses
            18 views
            0 likes
            Last Post seqadmin  
            Started by seqadmin, 03-09-2023, 10:17 AM
            0 responses
            67 views
            1 like
            Last Post seqadmin  
            Started by seqadmin, 03-03-2023, 12:03 PM
            0 responses
            64 views
            0 likes
            Last Post seqadmin  
            Working...
            X