Seqanswers Leaderboard Ad

Collapse

Announcement

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

  • False negative variant calling in haploids. Call variants using coverage (not stats)?

    Greetings,

    I am looking for the simplest possible way to call variants computationally, just using position coverage and the ratio of reference/alternative allele, for divergent haploid genomes. I.e. no need for prior probability or a prediction based on an allele frequency spectrum. There is good coverage (average depth of coverage is 316X), and the sequences are 100bp reads of high-quality.

    I am mapping with BWA and using bcftools to call variants using Bayesian inference in "mpileup" files. It seems bcftools is just too fancy for this to work among these divergent haploid genomes, and I am getting only half the variants I should be, even setting the quality scores as low as 5 for filtering.

    Looking at my BAMs, the variants are there but are getting removed somehow, probably because they don't fit with some expected allele-frequency spectrum due to the level of divergence (~10-20 variants/kb)

    Is there a simple way, using bcftools, a script that goes into the pileup, or another program to get better calls? FreeBayes is recommended here for haploid genomes, and I'll look into it, but I think it may end-up with the same problem.

    Thanks!
    Last edited by Genomics101; 06-05-2012, 10:23 AM.

  • #2
    I am assuming by variants you mean SNPs. FreeBayes is good. Also try JGIL.

    Comment


    • #3
      SNPs and small indels, yes. I'm looking into FreeBayes now. I also ran the analysis on each BAM individually, rather than using a population sample, with bcftools, and the results were much better.

      Comment


      • #4
        Here is another seqanswers thread that might give some ideas if you haven't already seen it:



        If your looking for a simple way to call SNPs from the base frequencies, you might check out



        where they call something a homozygous SNP if the non-reference allele is >80% (after filtering on quality Q20 and coverage 20X). They called hets if the alternate allele was between 20-80% (which wouldn't be pertinent for you since you are dealing with haploids). This could at least serve as a base line sanity check.

        The IGV browser employs a similar heuristic to show variants:

        If a nucleotide differs from the reference sequence in greater than 20% of quality weighted reads, IGV colors the bar in proportion to the read count of each base
        Justin

        Comment

        Latest Articles

        Collapse

        • seqadmin
          Exploring the Dynamics of the Tumor Microenvironment
          by seqadmin




          The complexity of cancer is clearly demonstrated in the diverse ecosystem of the tumor microenvironment (TME). The TME is made up of numerous cell types and its development begins with the changes that happen during oncogenesis. “Genomic mutations, copy number changes, epigenetic alterations, and alternative gene expression occur to varying degrees within the affected tumor cells,” explained Andrea O’Hara, Ph.D., Strategic Technical Specialist at Azenta. “As...
          07-08-2024, 03:19 PM
        • seqadmin
          Exploring Human Diversity Through Large-Scale Omics
          by seqadmin


          In 2003, researchers from the Human Genome Project (HGP) announced the most comprehensive genome to date1. Although the genome wasn’t fully completed until nearly 20 years later2, numerous large-scale projects, such as the International HapMap Project and 1000 Genomes Project, continued the HGP's work, capturing extensive variation and genomic diversity within humans. Recently, newer initiatives have significantly increased in scale and expanded beyond genomics, offering a more detailed...
          06-25-2024, 06:43 AM

        ad_right_rmr

        Collapse

        News

        Collapse

        Topics Statistics Last Post
        Started by seqadmin, Yesterday, 07:20 AM
        0 responses
        24 views
        0 likes
        Last Post seqadmin  
        Started by seqadmin, 07-16-2024, 05:49 AM
        0 responses
        38 views
        0 likes
        Last Post seqadmin  
        Started by seqadmin, 07-15-2024, 06:53 AM
        0 responses
        44 views
        0 likes
        Last Post seqadmin  
        Started by seqadmin, 07-10-2024, 07:30 AM
        0 responses
        41 views
        0 likes
        Last Post seqadmin  
        Working...
        X