I have a vcf compressed and indexed (.vcf.gz and .tbi)
I have a list of sites that I need to remove from this file (delete forever) while having the other information.
what would be the best way to go about this?
gunzip -c file.vcf.gz |\
awk -F '\t' '!($1=="chr8" && int($2)==125907085))' > out.vcf && \
bgzip out.vcf && \
tabix -p vcf out.vcf.gz
is there a more optimal way of achieving this, lets say i have like 100 positions to remove
See bedtools intersect
Login before adding your answer.
Use of this site constitutes acceptance of our User Agreement and Privacy