Seqanswers Leaderboard Ad

Collapse

Announcement

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

  • david.tamborero
    replied
    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!

    Leave a comment:


  • adaptivegenome
    replied
    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.

    Leave a comment:


  • david.tamborero
    replied
    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

    Leave a comment:


  • Heisman
    replied
    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.

    Leave a comment:


  • david.tamborero
    started a topic GATK indel raligner after bfast alignement

    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

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