Seqanswers Leaderboard Ad

Collapse

Announcement

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

  • GATK indel raligner after bfast alignement

    Hello,

    i'm attempting to perform a pipeline to find out somatic mutations of tumor-normal samples sequenced by Illumina paired end reads (100 + 100). After the alignement, I want to do some additional steps (according for instance to http://www.broadinstitute.org/gsa/wi...th_the_GATK_v2), as RemoveDuplicates, Recalibration and Indel_Realignement.

    The point is that I want to try the alignement by using bfast, which includes a localign step, in which as far as I know each read is aligned by Smith_Waterman algorithm (allowing gaps) for each CAL. So my question is: when using bfast alignement, is it necessary the GATK_indel_realigner step? I do not know if it is a non wise question, but I want to be sure of performing such GATK step since it's quite demanding and my cluster is not awesome at all.

    Thank you for your feedback!

    cheers
    david

  • #2
    When aligning you will only look at one read at a time, while with IndelRealignment all reads around a predicted indel will be looked at simultaneously. I would go ahead and try both and see what differences you get in your variant calls. If there's say a 1-2% difference, then maybe that wouldn't be enough to care.

    Comment


    • #3
      thank you for your explanation, Heisman.

      Then, I guess I have to try it. If the cluster and I are alive to get some results about this issue, I will post them in case anyone more is interested.

      Thank you again!
      cheers
      David

      Comment


      • #4
        Using the most recent GATK you do not have to do the complete SW realignment if your mapper already produced gapped alignments. See the documentation, there is a shortcut option that saves a little time.

        Comment


        • #5
          I've just seen it, I will try the indel_realignement command of the latest GATK by using the "--consensusDeterminationModel USE_READS" argument, it seems exactly what I was looking for.

          Many thanks!

          Comment

          Latest Articles

          Collapse

          • seqadmin
            Best Practices for Single-Cell Sequencing Analysis
            by seqadmin



            While isolating and preparing single cells for sequencing was historically the bottleneck, recent technological advancements have shifted the challenge to data analysis. This highlights the rapidly evolving nature of single-cell sequencing. The inherent complexity of single-cell analysis has intensified with the surge in data volume and the incorporation of diverse and more complex datasets. This article explores the challenges in analysis, examines common pitfalls, offers...
            06-06-2024, 07:15 AM
          • seqadmin
            Latest Developments in Precision Medicine
            by seqadmin



            Technological advances have led to drastic improvements in the field of precision medicine, enabling more personalized approaches to treatment. This article explores four leading groups that are overcoming many of the challenges of genomic profiling and precision medicine through their innovative platforms and technologies.

            Somatic Genomics
            “We have such a tremendous amount of genetic diversity that exists within each of us, and not just between us as individuals,”...
            05-24-2024, 01:16 PM

          ad_right_rmr

          Collapse

          News

          Collapse

          Topics Statistics Last Post
          Started by seqadmin, Today, 07:24 AM
          0 responses
          9 views
          0 likes
          Last Post seqadmin  
          Started by seqadmin, Yesterday, 08:58 AM
          0 responses
          11 views
          0 likes
          Last Post seqadmin  
          Started by seqadmin, 06-12-2024, 02:20 PM
          0 responses
          15 views
          0 likes
          Last Post seqadmin  
          Started by seqadmin, 06-07-2024, 06:58 AM
          0 responses
          182 views
          0 likes
          Last Post seqadmin  
          Working...
          X