SEQanswers

Go Back   SEQanswers > Bioinformatics > Bioinformatics



Similar Threads
Thread Thread Starter Forum Replies Last Post
what is DP4 in VCF (Variant Call Format) and AF1? gcrdb Bioinformatics 8 01-20-2015 09:37 AM
Samtools mpileup - Call somatic mutations from a pair of samples - detail steps ? swapnil2188 General 2 09-23-2013 04:49 PM
Doubts regarding samtools mpileup - Call somatic mutations from a pair of samples ? swapnil2188 RNA Sequencing 2 05-18-2012 06:54 AM
strand bias variant call yinshe Bioinformatics 2 10-10-2011 06:01 AM
Variant call using BAM from Tophat emilyjia2000 Bioinformatics 5 08-01-2011 04:37 PM

Reply
 
Thread Tools
Old 09-23-2012, 07:18 PM   #1
lchang1101
Junior Member
 
Location: St. Louis

Join Date: Apr 2012
Posts: 2
Default samtools mpileup did not call variant

I used samtools mpileup and bcftools to call variants in my RNA-Seq data. At one locus, I have 10 reads that are different from the ref sequence (DP4=32,21,5,5), but bcftools did not call this variant.

3 13812321 . T . 47 . DP=65;VDB=0.0190;AF1=0;AC1=0;DP4=32,21,5,5;MQ=20;FQ=-44;PV4=0.73,1,1,1 PL: DP:SP 0:63:1

I tried to change the -p value from the default 0.5 to 0.6, 0.7, 0.8 and 0.9 but still none of them called this variant. I understand that 10 reads in 65 reads is a low fraction, but is there anyway to make bcftools less stringent? I know I can just print out info for every locus and do my own filtering but I was just wondering if there are other parameters of bcftool tat I can modify. Thanks!
lchang1101 is offline   Reply With Quote
Old 09-24-2012, 06:58 AM   #2
vivek_
PhD Student
 
Location: Denmark

Join Date: Jul 2012
Posts: 164
Default

Isn't Mapq=20 rather too low to confidently take the alignment at that position?
vivek_ is offline   Reply With Quote
Old 09-24-2012, 07:22 AM   #3
lchang1101
Junior Member
 
Location: St. Louis

Join Date: Apr 2012
Posts: 2
Default

Thank you for the reply. No, I dont think MQ=20 was the problem. Some of the called variants in my run had MQ as low as 16. Any other thoughts?
lchang1101 is offline   Reply With Quote
Reply

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 09:07 PM.


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