This is a test version of Biostars. For the public version, visit https://www.biostars.org.
Remove indels from VCF file- not working

Hi,

I am trying to remove indels sites from vcf file. I used following command line;

vcftools --vcf inputfile.vcf --remove-indels --recode --recode-INFO-all --out outputfile

It still keeps all sites without removing indels. Please anybody has a solution for this? Thanks in advance.

snp software error
• 10,759 views • link
written by gaurab.bhattarai46 • gaurab.bhattarai46

Hi Kevin, Thank you so much for this information. I am kind of new to the bioinformatics. So, if you can suggest me the bcf command line to remove indels, it will help me a lot. In my SNPs data, there are both '-' indels as well as '0' i,e -/+ are present (Based on IUPAC). Thanks in advance.

• 0 views • link
written by gaurab.bhattarai46 • gaurab.bhattarai46

, there are both '-' indels as well as '0' i,e -/+ are present (Based on IUPAC)

?

can you please show us such line please

• 0 views • link
written by Pierre Lindenbaum Pierre Lindenbaum

1 answer

VCFtools is more or less deprecated, and even one of the co-authors of VCFtools (with whom I had the courtesy of having lunch one day) suggest to use BCFtools.

So, tis will do it for you:

bcftools view --types snps

Kevin

• 0 views • link
written by Kevin Blighe Kevin Blighe

Log in to answer this question.

More posts like this