Masking sites in a vcf file
I am new to Bioinformatics therefore, I am sorry if my question is too basic but, I am following this article to filter my vcf file from 1000Genomes:
see method section 'Application to human chromosome 18 from the 1000 Genomes CEU sample'
They say in the section that they masked all sites flagged by the 1000 Genomes Project as being unfit for population genetic analyses. 1000Genomes has those sites in bed format and Fasta as well. The link is as follows:
My question is how can I use the bed file or Fasta files provided to mask those sites in a vcf file. Is there a tool that does it? Insights will be appreciated.
0 answers
No answers yet.
Log in to answer this question.
More posts like this
-
Reference mouse phased vcf
written by machoo •
-
Why is NA18498 sample missing from 1000 genome GRCh38 vcf file?
written by octpus616 •
-
Concensus from 1000 genome project
written by Peerzada •
-
Allele Specific Expression (ASE) from RNA seq data
written by basuanubhav •
-
Using 1000 genomes phase 3 variants for BQSR
written by prasundutta87
-
1000 genomes missing SNPs on chr X in GRCh38
written by b_ha •
-
Mate mismatch errors in 1000genomes CRAM files
written by Mehulsharma.253 •
-
Report haplotype frequencies from a region within a phased vcf file
written by Krisr
-
How to get individual chromosome sequence in fasta format from vcf.gz and its vcf.gz.tbi file of 10…
written by David Simmons •
-
News: 1000 Genomes data on GRCh38
written by fairley •
I think that you do the opposite i.e. mask a fasta file using a VCF file. You mask a fasta file (usually a reference genome) based on the variant positions in a VCF file so that any downstream analysis software avoid those sites. This should help: https://bedtools.readthedocs.io/en/latest/content/tools/maskfasta.html
Thank you @prasundutta87 for you reply! I came across vcftools that I think I can use. I'll try to use that first.