This is a test version of Biostars. For the public version, visit https://www.biostars.org.
VCF - InDels and SNPs relation

Hi, I have files with my SNPs and InDels. Do you have any idea how to check whether indels are called mainly in the vicinity of SNPs or not? How can I check it and calculate this relation? Thank you for your all proposition.

snp indel vcf
• 1,591 views • link
updated by Pierre Lindenbaum Pierre Lindenbaum • written by misterie misterie

1 answer

you could used bedtools closest to get a distribution of the distances indel/snp

for T in snps indels 
do
bcftools view -G -v $T  inputvcf.gz |  bcftools query --format "%CHROM\t%POS\t%END\n" | sort -T. -t $'\t' -k1,1 -k2,2n > $T.bed
done

bedtools closest -a indels.bed -b snps.bed -d | cut -f 7 | sort -n
• 0 views • link
written by Pierre Lindenbaum Pierre Lindenbaum

wow! Thank you, I did not know this site of bedtools.

How can I understand zeros? That''s mean that SNP is inside indel?

• 0 views • link
written by misterie misterie

How can I understand zeros? That''s mean that SNP is inside indel?

yes

• 0 views • link
written by Pierre Lindenbaum Pierre Lindenbaum

Log in to answer this question.

More posts like this