Find SNPs common to mutants, absent in WT (vcf files)
0
0
Entering edit mode
9 months ago
charlesh • 0

Hi,

I have 3 vcf files from WT and mut1, mut2 I would like to get a file with an intersection of common SNPs to the 2 Mutants, but absent from the WT sample.

I've run bcftools isec on all pairs of samples, but am hung up on using the resulting 3 sets of vcf files (0000, 00011,0002,0003) to identify SNPs common to the mutants, but absent from the WT.

Suggestions? Charles

vcf bcftools • 203 views
ADD COMMENT
0
Entering edit mode

merge the 3 vcf files with bcftools and select the variants with gatk SelectVariants --select <expr> https://gatk.broadinstitute.org/hc/en-us/articles/360036726131-SelectVariants

ADD REPLY

Login before adding your answer.

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