SEQanswers

Go Back   SEQanswers > Bioinformatics > Bioinformatics



Similar Threads
Thread Thread Starter Forum Replies Last Post
How to represent reference and observed alleles in BED files for variants on negative LauraSmith Bioinformatics 0 02-17-2012 12:47 PM
Using UnifiedGenotyper to call variants from haploid data? oiiio Bioinformatics 0 01-24-2012 08:29 PM
samtools false positive/negative dnusol Bioinformatics 0 11-08-2011 07:13 AM
Variant calling for high-coverage Illumina data dgmacarthur Bioinformatics 19 04-08-2011 06:34 AM
TopHat results-- negative coverage values seqseqseq Bioinformatics 0 01-21-2010 12:54 PM

Reply
 
Thread Tools
Old 06-05-2012, 07:09 AM   #1
Genomics101
Member
 
Location: Maryland, USA

Join Date: May 2012
Posts: 60
Default 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 at 10:23 AM.
Genomics101 is offline   Reply With Quote
Old 06-05-2012, 07:29 AM   #2
adaptivegenome
Super Moderator
 
Location: US

Join Date: Nov 2009
Posts: 437
Default

I am assuming by variants you mean SNPs. FreeBayes is good. Also try JGIL.
adaptivegenome is offline   Reply With Quote
Old 06-07-2012, 07:39 AM   #3
Genomics101
Member
 
Location: Maryland, USA

Join Date: May 2012
Posts: 60
Default

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.
Genomics101 is offline   Reply With Quote
Old 06-07-2012, 11:29 AM   #4
BAMseek
Senior Member
 
Location: St. Louis, MO, USA

Join Date: Apr 2011
Posts: 124
Default

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

http://seqanswers.com/forums/showthread.php?t=10182

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

http://protocols.netlab.uky.edu/~liu...rs/Vijay_1.pdf

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:

Quote:
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
BAMseek is offline   Reply With Quote
Reply

Tags
bcftools, false negatives, haploid, mapping, variant analysis

Thread Tools

Posting Rules
You may not post new threads
You may not post replies
You may not post attachments
You may not edit your posts

BB code is On
Smilies are On
[IMG] code is On
HTML code is Off




All times are GMT -8. The time now is 06:07 AM.


Powered by vBulletin® Version 3.8.9
Copyright ©2000 - 2019, vBulletin Solutions, Inc.
Single Sign On provided by vBSSO