bcftools mpileup giving more wildtypes <*> than expected
1
0
Entering edit mode
3.4 years ago
gt ▴ 30

Hi there, I am having trouble with interpreting bcftools mpileup vcf output. My file has almost all <*> listed for the alternate field.

enter image description here

I am aligning my paired-end reads to an E.coli reference genome. Anyone have any intuition as to why this is happening?

bcftools DNA-seq mpileup • 921 views
ADD COMMENT
1
Entering edit mode
3.4 years ago

bcftools mpileup must be followed by bcftools call to call the variants. see "Variant calling" in http://www.htslib.org/workflow/

ADD COMMENT
0
Entering edit mode

Yep, I performed variant calling on the mpileup. I am just wondering why I have no ALT equal to "." in my vcf file at all. So when I call variants, my vcf file only has a very small number of mutations. I understand that "<*>" is a wild type but I would like to know if any of those could be possible indels.

ADD REPLY

Login before adding your answer.

Traffic: 1586 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6