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
              Advanced Methods for the Detection of Infectious Disease
              by seqadmin




              The recent pandemic caused worldwide health, economic, and social disruptions with its reverberations still felt today. A key takeaway from this event is the need for accurate and accessible tools for detecting and tracking infectious diseases. Timely identification is essential for early intervention, managing outbreaks, and preventing their spread. This article reviews several valuable tools employed in the detection and surveillance of infectious diseases.
              ...
              11-27-2023, 01:15 PM
            • seqadmin
              Strategies for Investigating the Microbiome
              by seqadmin




              Microbiome research has led to the discovery of important connections to human and environmental health. Sequencing has become a core investigational tool in microbiome research, a subject that we covered during a recent webinar. Our expert speakers shared a number of advancements including improved experimental workflows, research involving transmission dynamics, and invaluable analysis resources. This article recaps their informative presentations, offering insights...
              11-09-2023, 07:02 AM

            ad_right_rmr

            Collapse

            News

            Collapse

            Topics Statistics Last Post
            Started by seqadmin, 12-05-2023, 02:24 PM
            0 responses
            14 views
            0 likes
            Last Post seqadmin  
            Started by seqadmin, 12-05-2023, 07:37 AM
            0 responses
            25 views
            0 likes
            Last Post seqadmin  
            Started by seqadmin, 12-04-2023, 08:23 AM
            0 responses
            9 views
            0 likes
            Last Post seqadmin  
            Started by seqadmin, 12-01-2023, 09:55 AM
            0 responses
            26 views
            0 likes
            Last Post seqadmin  
            Working...
            X