How to filter vcf file by MAF using bcftools?
I saw this thread (https://github.com/samtools/bcftools/issues/357) and it seems sometimes bcftools recomputes MAF and sometimes it takes it from INFO. Which one is recommended to use and how is this calculation different than what stated on INFO?
Also, what is the difference between:
bcftools view -q 0.01:minor data.vcf.gz -Ov -o out.vcf
bcftools view -i 'MAF > 0.01' data.vcf.gz -Ov -o out.vcf
bcftools isec -e'MAF<0.01' data.vcf.gz -Ov -o out.vcf
Thank you!
• 3,993 views
•
link
updated
by
Ram
•
written
by
_quantum_girl_ •
•
written
by
_quantum_girl_ •
0 answers
No answers yet.
Log in to answer this question.
More posts like this
-
Issues filtering VCF file with bcftools/vcftools
written by ekirsch •
-
bcftools=1.18 not filtering correcting MAF
written by AMARU •
-
Not returning any alleles when using bcftools view -q
written by kshaffman •
-
Error while using VCFtools (--maf) to filter by allele frequency
written by kshaffman •
-
Different estimates of nucleotide diversity (pi) from two pipelines: pixy vs vcftools
written by nitinra •
-
VCF parse error: Could not add dummy header for contig
written by User000
-
How to parallelize bcftools mpileup with GNU parallel?
written by beausoleilmo
-
Error with bcftools and vcftools and removing individuals
written by seta
-
Bcftools stats error: Incorrect allele ("1") at ...
written by msf •
-
CNVkit - Format of VCF file
written by fongchunchan •
Here is the
BCFtools cheat sheetthat will explain your questions and more.https://gist.github.com/elowy01/93922762e131d7abd3c7e8e166a74a0b
Thank you! Do you know when isec options should be used? What is the difference between just using view or isec?
isecoption is used to find overlapping variants in different vcf files.viewoption in used to subset, filter and convert VCF and BCF filesIMO that cheat sheet is quite confusing. The official manual is a lot clearer. OP needs to read the freaking manual, especially if they can't even tell
isecandviewapart.