Seqanswers Leaderboard Ad

Collapse

Announcement

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

  • Calling SNP frequencies in BWA alignments

    Namaste all.

    I am using BWA to align my reads. I need to not only call the variants but also their frequencies, for example, if I have 30X coverage at a chromosomal position, the reference has A and of the 30 reads I have 20 As and 10 Ts, how do I do that? any tools available for this?

    Any help would be highly appreciated.

  • #2
    Have a look at samtools mpileup and GATK UnifiedGenotyper, I think they can do what you want.

    Comment


    • #3
      Thanks YiLong. Both the tools you mentioned can only call allele frequencies if multiple aligned .bam files are present. They do not call allele frequencies from different reads of the same sample. What I need is: if base (say at chr2 position 135487) is heterozygous A/T, say reference is A, and I have 30 reads at that position. I want to know how many reads are As and how many are Ts. I want to this for every position in the genome for every sample. I am not sure whether mpileup or GATK do such a thing.

      Comment


      • #4
        In mpileup in the INFO field, there's the AF1 subfield that'll tell you the frequency.

        In GATK UnifiedGenotyper the same info can be found in the last field in the AD subfield.

        At least this is how I have understood it though I haven't read this from anywhere so use with caution... Also you might have to play around with the parameters to get these fields into the output - I don't remember what parameters I used for generating these .vcf files.

        Comment

        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
        44 views
        0 likes
        Last Post seqadmin  
        Started by seqadmin, 05-09-2024, 02:46 PM
        0 responses
        54 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