How can I see the low frequency SNP
0
0
Entering edit mode
2.3 years ago
apl00028 ▴ 70

I have created a vcf file from bam file using bcftools:

 bcftools mpileup \
--redo-BAQ \
--min-BQ 30 \
--per-sample-mF \
--annotate \
bcftools call --multiallelic-caller --variants-only -Ov > \
PV002_nofilt.vcf


After that I visualized this file using R, but It had lower number of records in comparation to IGV program results using the same bam file. I think that it is due to bcftools take only SNP with high frequency.

What modification should do add to the bcftools program to see all records although they have a low frequency?

SNP ngs • 481 views
0
Entering edit mode

Hello,

how low is "low"?

fin swimmer

0
Entering edit mode

With a Coveage allele-fraction threshold of 0.05. I see them when I run IGV but I can not with the vcf file: