SEQanswers

Go Back   SEQanswers > Bioinformatics > Bioinformatics



Similar Threads
Thread Thread Starter Forum Replies Last Post
Searching bases in reference genome. kasthuri Bioinformatics 2 01-12-2012 03:54 AM
VCF Reference in UnifiedGenotyper. kasthuri Bioinformatics 0 06-09-2011 04:31 AM
Retrieving single bases of the reference genome in R mixter Bioinformatics 1 10-18-2010 02:30 AM
Why do we use mapping programs instead of blast for mapping to a reference? thsuk1 Bioinformatics 6 08-27-2010 08:54 AM
reference mapping dina Bioinformatics 0 10-03-2009 06:18 AM

Reply
 
Thread Tools
Old 11-16-2011, 03:44 AM   #1
wswill
Junior Member
 
Location: china Beijing

Join Date: Nov 2011
Posts: 1
Default why the vcf has so little high_quality bases mapping on the reference?

hi,all
Now,i have met a problem when using samtools to get SNP call. In the vcf file generated from samtools, i find high_quality bases so little in col(DP4=0,0,21,15), although coverage is very high(DP=5773). Then i check the sam file, and find that reads' BQ are almost high.
So,i want to know why .
The command:
mpileup -DugS -C 50 -f scaffold.fa all.bam |bcftools view -bcvg -> all.var.bcf
bcftools view all.var.bcf | vcfutils.pl varFilter > all.var.flt.vcf
The vcf file(partical):
gi|158420570|ref|NC_009938.1| 115 . C T,A 168 . DP=5773;VDB=0.0472;AF1=1;AC1=2;DP4=0,0,21,15;MQ=36;FQ=-132 GT:PLP:SP:GQ 1/1:201,105,0,183,71,180:36:0:99
gi|158420570|ref|NC_009938.1| 784 . G A 177 . DP=6382;VDB=0.0444;AF1=1;AC1=2;DP4=2,1,145,11;MQ=36;FQ=-282;PV4=0.21,9.3e-10,0.29,1 GT:PLP:SP:GQ 1/1:210,255,0:159:7:99
wswill is offline   Reply With Quote
Old 12-12-2011, 09:16 PM   #2
tedtoal
Junior Member
 
Location: Davis, CA

Join Date: Dec 2011
Posts: 3
Default DP4 too small

I have the same problem that you described. DP4 values are small even though DP is high and quality is high. Did you ever find an answer to this?
tedtoal is offline   Reply With Quote
Old 12-13-2011, 09:22 AM   #3
swbarnes2
Senior Member
 
Location: San Diego

Join Date: May 2008
Posts: 912
Default

Quick thing to try, rerun mpileup with -B or -E. Sometimes, the BAQ calcautions will wrongly downgrade the heck out of quality scores at variations, make it hard for the caller to call them.

The other thing to do is to actually look at the .sam file of the reads that cross that locus. What is the mapping quality of those reads? What is the quality of the letter at that locus?
swbarnes2 is offline   Reply With Quote
Old 12-13-2011, 10:05 AM   #4
Heisman
Senior Member
 
Location: St. Louis

Join Date: Dec 2010
Posts: 535
Default

Also include -A. That has solved this problem for me in the past.
Heisman is offline   Reply With Quote
Reply

Tags
sam file, samtools mpileup, vcf format

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 01:12 PM.


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