vcf files comparison to find unique variants
I have 2 vcf files for a single sample that has been sequenced twice, once before and again after developing resistance to a certain antibiotic. Is there a command I can use to extract only newly developed variants in the (post resistance) sample? (I mean to exclude mutual variants in both files)
1 answer
Use bcftools isec as explained here: How bcftools isec works ?
Log in to answer this question.
More posts like this
-
how to extract DNA seq from a file which contain both protein and DNA
written by Neel •
-
BCFTOOLS filter per sample per variant
written by Filago •
-
VCF using Python
written by V Vera •
-
Where can I find VCF files of the same source
written by ManuelDB
-
Performing so-called "SNP confirmation" for reads assigned to their corresponding Antibiotic Resist…
written by Matteo Schiavinato
-
population normalized read depth
written by sumudu_rangika •
-
Remove positions that are non-variant in a subset of samples from a vcf file
written by Aliaksei Holik
-
how to extract the variants using Bed file
written by k.manojkumar2991 •
-
combining sequencing variants
written by hmk123 •
-
Merging Two Fastq Files From Exome Sequencing
written by User 1933 •
This has been asked before : Best way to compare two samples in a VCF file