SEQanswers

Go Back   SEQanswers > Bioinformatics > Bioinformatics



Similar Threads
Thread Thread Starter Forum Replies Last Post
Understanding bcftools call output evoevo Bioinformatics 0 08-31-2017 11:05 PM
Empty VCF file with bcftools call AP38 Bioinformatics 8 05-10-2016 09:24 AM
bcftools issue willgrand Bioinformatics 11 12-03-2015 01:20 PM
bcftools doesn't call suspected indel rwmills Bioinformatics 3 10-05-2015 09:43 AM
Getting FQ from bcftools 1.2 Greg_owens Bioinformatics 2 06-18-2015 08:18 AM

Reply
 
Thread Tools
Old 03-06-2020, 02:42 AM   #1
samhokin
Member
 
Location: Santa Fe, NM

Join Date: Nov 2013
Posts: 20
Question bcftools call: why is this a 1/1?

I can't for the life of me figure out why bcftools call genotypes this position (and many, many others like it) as 1/1:

Here's the input VCF line:

Code:
10	227179	.	C	T,<*>	0	.	DP=191;I16=91,87,7,2,6502,238474,333,12321,235,3775,420,25200,3675,85899,97,1569;QS=0.733,0.267,0;VDB=0.00520152;SGB=-0.662043;RPB=0.000875247;MQB=0.0198712;MQSB=0.892843;BQB=0.980427;MQ0F=0.0209424	PL:DP:ADF:ADR:AD	64,255,0,255,27,95:187:91,7,0:87,2,0:178,9,0
and here's the output from bcftools call -m :

Code:
10	227179	.	C	T	34.7149	.	DP=191;VDB=0.00520152;SGB=-0.662043;RPB=0.000875247;MQB=0.0198712;MQSB=0.892843;BQB=0.980427;MQ0F=0.0209424;AC=2;AN=2;DP4=91,87,7,2;MQ=3	GT:PL:DP:ADF:ADR:AD	1/1:64,255,0:187:91,7:87,2:178,9
Note that DP4=91,87,7,2 -- there are 178 REF reads and 9 ALT reads used to make the call. How on earth is that ALT homozygous "1/1"??

I hugely appreciate any explanation of this call!
__________________
Sam Hokin
Computational Scientist, Carnegie and NCGR
samhokin is offline   Reply With Quote
Reply

Tags
bcftools, diploid calling, variants calling

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 04:02 PM.


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