Seqanswers Leaderboard Ad

Collapse

Announcement

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

  • Allele frequency calculation in SNP calling

    Hi all,

    I was wondering how the allele frequency is calculated by SAMTOOLS (I'm using mpileup and then bcftools).
    It seems that my SNPs always have a frequency of either 1 or 0.5, and do not always match the DP4 values.

    For example:
    Code:
    chr01   11961   .       G       C       52      .       DP=38;AF1=0.5;CI95=0.5,0.5;DP4=12,16,3,7;MQ=42;FQ=55;PV4=0.71,0.26,0.052,0.41   GT:PL:GQ        0/1:82,0,255:85
    Here it seems that the allele frequency should be more 0.3 than 0.5, doesn't it?

    Another example:
    Code:
    chr02   667170  .       C       T       14.2    .       DP=135;AF1=1;CI95=1,1;DP4=13,1,79,38;MQ=60;FQ=-71;PV4=0.064,6.2e-46,1,1 GT:PL:GQ        1/1:47,44,0:55
    In this case the frequency is closer to 0.9 than to 1...

    How is the frequency calculated?

    Thanks!
    Rachelly.

  • #2
    I have the same doubt!!!!!!!!

    Comment


    • #3
      I am running into the same question with GATK unified genotyper output. Any help? I want the exact frequency using the counts. It seems to always just give a .5 or 1

      Comment


      • #4
        GATK and samtools assume samples are all diploid. You should look at DP4 or pileup instead.

        Also note that samtools does not work with pooled samples.

        Comment


        • #5
          Originally posted by lh3 View Post
          GATK and samtools assume samples are all diploid. You should look at DP4 or pileup instead.

          Also note that samtools does not work with pooled samples.

          Could you clarify what DP4 or pileup is? I wasn't able to find a software package called either of those things. Could you point me in the right direction?

          Comment


          • #6
            Originally posted by shawpa View Post
            Could you clarify what DP4 or pileup is? I wasn't able to find a software package called either of those things. Could you point me in the right direction?
            Ignore what I just wrote. I realized now I was being very stupid because you were referring to a format. I guess I am just trying to find an easy way to extract those counts so that I can calculate allele frequency from those. I can't do that with GATK.

            Comment

            Latest Articles

            Collapse

            • seqadmin
              Current Approaches to Protein Sequencing
              by seqadmin


              Proteins are often described as the workhorses of the cell, and identifying their sequences is key to understanding their role in biological processes and disease. Currently, the most common technique used to determine protein sequences is mass spectrometry. While still a valuable tool, mass spectrometry faces several limitations and requires a highly experienced scientist familiar with the equipment to operate it. Additionally, other proteomic methods, like affinity assays, are constrained...
              04-04-2024, 04:25 PM
            • seqadmin
              Strategies for Sequencing Challenging Samples
              by seqadmin


              Despite advancements in sequencing platforms and related sample preparation technologies, certain sample types continue to present significant challenges that can compromise sequencing results. Pedro Echave, Senior Manager of the Global Business Segment at Revvity, explained that the success of a sequencing experiment ultimately depends on the amount and integrity of the nucleic acid template (RNA or DNA) obtained from a sample. “The better the quality of the nucleic acid isolated...
              03-22-2024, 06:39 AM

            ad_right_rmr

            Collapse

            News

            Collapse

            Topics Statistics Last Post
            Started by seqadmin, 04-11-2024, 12:08 PM
            0 responses
            15 views
            0 likes
            Last Post seqadmin  
            Started by seqadmin, 04-10-2024, 10:19 PM
            0 responses
            21 views
            0 likes
            Last Post seqadmin  
            Started by seqadmin, 04-10-2024, 09:21 AM
            0 responses
            16 views
            0 likes
            Last Post seqadmin  
            Started by seqadmin, 04-04-2024, 09:00 AM
            0 responses
            45 views
            0 likes
            Last Post seqadmin  
            Working...
            X