Seqanswers Leaderboard Ad

Collapse

Announcement

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

  • Jon_Keats
    replied
    You need to be careful of calling so aggressively that the "not-so-likely" variants are called. What will most likely happen is your will end up with a bunch of recurrent variants that are all noise as some regions will come up again and again with such relaxed settings and you will waste your time and money validating false positives. If you goal is what you say try screening samples with a stringent filter, then generate a pileup of the individual positions with samtools to see if the same variant exists in the other samples. Assuming the variant if frequent is should show up with a stringent filter in at least one sample if not a couple/most of the samples. This should limit the number of false positives detected and eliminate a horrible realization that your bioinformatics filtering strategy resulted in a validation rate of 2%. Tends to make PIs and pipet jockeys chase you around the institute with anger in their eyes.

    Leave a comment:


  • Yilong Li
    replied
    Yes, I fully agree that this particular position has so many dubious factors that it is unlikely to be a true variant .

    However, currently we would like to find recurrent variants present in multiple samples. For example, if a high confidence variant was called at this position in another sample, we would like to know whether there is even the weakest evidence that a variant exists also in the current sample. That's why we want to get the not-so-likely variants called as well.

    Leave a comment:


  • nilshomer
    replied
    Also observe that all the Cs are from reads on the same strand. Not a good thing!

    Leave a comment:


  • Yilong Li
    replied
    I have done realignment but not base quality recalibration.

    If I understood correctly though, GATK UG calls variants from a bam file and does not know whether the realignment or recal has been done, so even if the base qualities I use were recalibrated, GATK UG wouldn't make any differing calls.

    Leave a comment:


  • jstjohn
    replied
    Hi Yilong,
    Did you follow the full GATK protocol with realignment and base quality recalibration prior to genotyping? I agree that a "C" call does not seem very likely in that position though.
    -John

    Leave a comment:


  • Yilong Li
    replied
    Sorry, I didn't provide enough information regarding the situation. I ran GATK UG using --min_base_quality 1 (so all the grey C bases should have been included in the calling). Data is from Illumina.

    There are several C bases in the position with decent quality, although I must say those reads do seem a bit dubious. I agree that this could be a false positive (haven't Sanger sequenced), but I am still interested what happens internally in GATK UG at this position and how I could call this or similar position as variant loci. I am trying to find recurrent mutation within a set of samples and therefore want to achieve maximal a relatively high sensitivity for single sample variant calls.

    Leave a comment:


  • nilshomer
    replied
    Two things to observe:

    1. The "C" bases are shaded or greyed out, meaning they are of low quality.
    2. They are towards the ends of reads (SOLID?), and are therefore more ambiguous. Also, he top C read has another improper C to the left (same for some of the others), so how much should we trust this/these reads?

    I think the genotyper did the absolute correct thing. Have you tried validating this with Sanger Sequencing?

    I would not trust a C here.

    Leave a comment:


  • "Obvious variant" missed by GATK UnifiedGenotyper

    Hi all,

    We have noticed in our lab that GATK UnifiedGenotyper misses some variants (as one shown in the IGV figure below) that seem reasonably obvious with good alt. allele counts.

    I ran GATK UG -all_bases, and this site was called as
    Code:
    1       43308444        .       A       .       35.64   .       AC=0;AF=0.00;AN=2;DP=38;Dels=0.00;MQ=57.26;MQ0=0        GT:DP:GQ:PL     0/0:24:5.65:0,6,677
    Does anybody have any ideas why GATK UG called this site as ref (with such high confidence!) and how to loosen up the GATK UG parameters to have this site called as a variant?

    Yilong

    Last edited by Yilong Li; 04-05-2011, 06:23 AM. Reason: updated image url

Latest Articles

Collapse

  • seqadmin
    Recent Advances in Sequencing Analysis Tools
    by seqadmin


    The sequencing world is rapidly changing due to declining costs, enhanced accuracies, and the advent of newer, cutting-edge instruments. Equally important to these developments are improvements in sequencing analysis, a process that converts vast amounts of raw data into a comprehensible and meaningful form. This complex task requires expertise and the right analysis tools. In this article, we highlight the progress and innovation in sequencing analysis by reviewing several of the...
    05-06-2024, 07:48 AM
  • seqadmin
    Essential Discoveries and Tools in Epitranscriptomics
    by seqadmin




    The field of epigenetics has traditionally concentrated more on DNA and how changes like methylation and phosphorylation of histones impact gene expression and regulation. However, our increased understanding of RNA modifications and their importance in cellular processes has led to a rise in epitranscriptomics research. “Epitranscriptomics brings together the concepts of epigenetics and gene expression,” explained Adrien Leger, PhD, Principal Research Scientist...
    04-22-2024, 07:01 AM

ad_right_rmr

Collapse

News

Collapse

Topics Statistics Last Post
Started by seqadmin, 05-14-2024, 07:03 AM
0 responses
19 views
0 likes
Last Post seqadmin  
Started by seqadmin, 05-10-2024, 06:35 AM
0 responses
43 views
0 likes
Last Post seqadmin  
Started by seqadmin, 05-09-2024, 02:46 PM
0 responses
53 views
0 likes
Last Post seqadmin  
Started by seqadmin, 05-07-2024, 06:57 AM
0 responses
42 views
0 likes
Last Post seqadmin  
Working...
X