Seqanswers Leaderboard Ad

Collapse
X
 
  • Filter
  • Time
  • Show
Clear All
new posts

  • BAMseek
    replied
    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

    Leave a comment:


  • Genomics101
    replied
    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.

    Leave a comment:


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

    Leave a comment:


  • 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.

Latest Articles

Collapse

  • seqadmin
    Pathogen Surveillance with Advanced Genomic Tools
    by seqadmin




    The COVID-19 pandemic highlighted the need for proactive pathogen surveillance systems. As ongoing threats like avian influenza and newly emerging infections continue to pose risks, researchers are working to improve how quickly and accurately pathogens can be identified and tracked. In a recent SEQanswers webinar, two experts discussed how next-generation sequencing (NGS) and machine learning are shaping efforts to monitor viral variation and trace the origins of infectious...
    03-24-2025, 11:48 AM
  • seqadmin
    New Genomics Tools and Methods Shared at AGBT 2025
    by seqadmin


    This year’s Advances in Genome Biology and Technology (AGBT) General Meeting commemorated the 25th anniversary of the event at its original venue on Marco Island, Florida. While this year’s event didn’t include high-profile musical performances, the industry announcements and cutting-edge research still drew the attention of leading scientists.

    The Headliner
    The biggest announcement was Roche stepping back into the sequencing platform market. In the years since...
    03-03-2025, 01:39 PM

ad_right_rmr

Collapse

News

Collapse

Topics Statistics Last Post
Started by seqadmin, 03-20-2025, 05:03 AM
0 responses
49 views
0 reactions
Last Post seqadmin  
Started by seqadmin, 03-19-2025, 07:27 AM
0 responses
57 views
0 reactions
Last Post seqadmin  
Started by seqadmin, 03-18-2025, 12:50 PM
0 responses
50 views
0 reactions
Last Post seqadmin  
Started by seqadmin, 03-03-2025, 01:15 PM
0 responses
201 views
0 reactions
Last Post seqadmin  
Working...