Seqanswers Leaderboard Ad

Collapse

Announcement

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

  • Discrepancy in DP value in vcf file and reads visualization in IGV

    A sample variant from my vcf file :

    I 27000 . A G 37.88 . AC=2;AF=1.00;AN=2;DP=51;Dels=0.00;FS=0.000;HRun=1;HaplotypeScore=0.0000;MQ=6.53;MQ0=49;QD=0.74 GT:ADP:GQ:PL 1/1:37,14:51:6.01:69,6,0

    which gives DP=51 but while visualizing the variant on IGV i'm getting
    Total Count as 56. Seems like some reads have been disqualified for calling variant and the number is lesser in vcf than the bam file.

    I need to understand this discrepancy. I would like to view same counts while in vcf and in viewer.
    Vaibhav Kulkarni

  • #2
    What software produced the output?

    Does the meta-information line match this?
    ##FORMAT=<ID=DP,Number=1,Type=Integer,Description="Read Depth">

    If the vcf metainfo matches the line below refer to documentation at https://www.broadinstitute.org/gatk

    ##FORMAT=<ID=DP,Number=1,Type=Integer,Description="Approximate read depth (reads with MQ=255 or with bad mates are filtered)">

    If you are using varscan see:

    Comment


    • #3
      It's also quite likely that whatever made the VCF file was accounting for overlapping paired-end reads, whereas IGV doesn't.

      Comment


      • #4
        Originally posted by m_two View Post
        What software produced the output?

        Does the meta-information line match this?
        ##FORMAT=<ID=DP,Number=1,Type=Integer,Description="Read Depth">

        If the vcf metainfo matches the line below refer to documentation at https://www.broadinstitute.org/gatk

        ##FORMAT=<ID=DP,Number=1,Type=Integer,Description="Approximate read depth (reads with MQ=255 or with bad mates are filtered)">

        If you are using varscan see:
        http://varscan.sourceforge.net/suppo...unts-different
        As per your quote above my vcf file header shows:

        My vcf file has been generated through UnifiedGenotyper.

        ##INFO=<ID=DP,Number=1,Type=Integer,Description="Filtered Depth">
        Vaibhav Kulkarni

        Comment


        • #5
          I see something similar in my outputs from Illumina's somatic variant caller. My relevant header line is this:

          ##INFO=<ID=DP,Number=1,Type=Integer,Description="Total Depth">

          and my allele depth values don't add up to the total, although it's usually only off by fewer than 10 reads.

          Comment

          Latest Articles

          Collapse

          • seqadmin
            The Impact of AI in Genomic Medicine
            by seqadmin



            Artificial intelligence (AI) has evolved from a futuristic vision to a mainstream technology, highlighted by the introduction of tools like OpenAI's ChatGPT and Google's Gemini. In recent years, AI has become increasingly integrated into the field of genomics. This integration has enabled new scientific discoveries while simultaneously raising important ethical questions1. Interviews with two researchers at the center of this intersection provide insightful perspectives into...
            02-26-2024, 02:07 PM
          • seqadmin
            Multiomics Techniques Advancing Disease Research
            by seqadmin


            New and advanced multiomics tools and technologies have opened new avenues of research and markedly enhanced various disciplines such as disease research and precision medicine1. The practice of merging diverse data from various ‘omes increasingly provides a more holistic understanding of biological systems. As Maddison Masaeli, Co-Founder and CEO at Deepcell, aptly noted, “You can't explain biology in its complex form with one modality.”

            A major leap in the field has
            ...
            02-08-2024, 06:33 AM

          ad_right_rmr

          Collapse

          News

          Collapse

          Topics Statistics Last Post
          Started by seqadmin, Yesterday, 06:12 AM
          0 responses
          19 views
          0 likes
          Last Post seqadmin  
          Started by seqadmin, 02-23-2024, 04:11 PM
          0 responses
          67 views
          0 likes
          Last Post seqadmin  
          Started by seqadmin, 02-21-2024, 08:52 AM
          0 responses
          74 views
          0 likes
          Last Post seqadmin  
          Started by seqadmin, 02-20-2024, 08:57 AM
          0 responses
          66 views
          0 likes
          Last Post seqadmin  
          Working...
          X