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 \
FORMAT/AD,FORMAT/DP,INFO/AD -f cmv3 scaffold_nofilt_PV002.bam | \
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
ADD COMMENT
0
Entering edit mode

Hello,

how low is "low"?

fin swimmer

ADD REPLY
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:

https://bioinformatics-core-shared-training.github.io/intro-to-IGV/InspectingVariantsInIGV.html

ADD REPLY

Login before adding your answer.

Traffic: 1766 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