How can I see the low frequency SNP

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

• 1,556 views • link

written 7.3 years ago by apl00028

Hello,

how low is "low"?

fin swimmer

• 0 views • link

written 7.3 years ago by finswimmer

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

• 0 views • link

written 7.3 years ago by apl00028

0 answers

No answers yet.

More posts like this

AD sum about half of DP in vcf file
written 3.1 years ago by MH85 •

In my vcf file, the sum of total allelic depth (AD) is about half of the raw read depth (DP) for all SNPs. Is it …
(Bcftools) SNP allele frequency in DNA pools
written 6.0 years ago by gwyneth0011 •

I tried to use Bcftools to calculate allele frequency by SNP variants after bowtie 2 alignment. bcftools mpilup --redo-BAQ --min-BQ 30 --per-sample-mF --annotate FORMAT/AD, FORMAT/ADF, …
samtools having to use -A input option
written 7.0 years ago by williamsbrian5064

Hi, I am using samtools to make a bcf file and then a vcf file. I am using the following code samtools mpileup --redo-BAQ --min-BQ …
How to set up my vcf with a minimum depth?
written 7.1 years ago by apl00028

I am trying to set up the minimum depth which samtools mpileup considers a polymorphism for RNAseq samples. For that I did this: I did …
How can I detect low frequent polymorphisms ?
written 7.1 years ago by apl00028

Hi everybody, I am trying to detect polymorphisms (and their proportions) to sanger sequences of the genes that codify for Watermelon Mosaic Virus Coat protein: …
getting the HOM_REF genotypes using bcftools mpileup+call ?
written 7.5 years ago by Pierre Lindenbaum

I've been asked to find the status of some variants (= overlaping a bed file) in a '1000 genomes' data set, **including** the **HOM_REF** genotypes …
Samtools mpileup taking a long time to finish
written 8.2 years ago by williamsbrian5064

Hi, I am using samtools for variant calling. The issue I am having is that samtools mpileup seems to be taking a very long time. …
How can I add PV4 or BQ or NS to my VCF file
written 8.7 years ago by sahar.voghoei •

I am trying to extract maximum features for my machine learning input, and I need PV4, BQ, G3, NS to be added to my VCF …
STAR aligner samtools mpileup C50
written 9.7 years ago by zeyuchen2016 •

I use STAR align different RNA-SEQ sample to ref genome，then use samtools mpileup call snp![enter image description here][1] samtools mpileup -u -q 20 -Q 25 …
direct call of indels/snps with bcftools?
written 10.0 years ago by cristina_sabiers

I am wondering... With vcftools I can call directly my snps and indels from my vcf file, but with bcftools I only found this way …