I am trying to filter a VCF file to only include variants that are within another file, which is a txt file with VCF formatted columns (CHR POS REF ALT). I have been having a hard time finding a way to filter by REF and ALT, in addition to CHR POS (which can be accomplished in plink). Anyone have tips?
•
written
by
another_bioinformatician •
2 answers
The usual way to do this with plink2 is to first use --set-all-var-ids to make all your variant IDs based on position+alleles, and then use --extract for the filtering step. Note that you can use --recover-var-ids at a later time to restore your original variant IDs when you need them.
convert your tab delimited file to VCF by adding a dummy VCF header and some empty columns (.). You can then use tools like bcftools isec --collapse none .
Log in to answer this question.
-
looking for 1000g hg38 reference data to anchor with my genomic data to plot the principal compone…
written by nicole •
-
Filtering VCF based on codon-transcript
written by avelarbio46 •
-
Reomve the first occurance of overalapping variants from VCF file
written by Frieda •
-
Change column order in vcf file
written by gms03223 •
-
Merge two files with multiple columns by 4 common columns
written by User000
-
Calculating coverage from ref and alt columns?
written by PsFeathers •
-
How to seperate vcf annoation info field
written by raya.girish •
-
Comparing samtools,gatk and varscan against gold standard 1000G snps and need to decide the best me…
written by bioinforesearchquestions
-
Best way to query VCF for specific variants
written by Sean
-
Replace fields CHROM and POS in a vcf file
written by hosseinv •
