VCF - InDels and SNPs relation
1
0
Entering edit mode
2.5 years ago
misterie ▴ 90

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 • 692 views
ADD COMMENT
0
Entering edit mode
2.5 years ago

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
ADD COMMENT
0
Entering edit mode

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

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

ADD REPLY
0
Entering edit mode

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

yes

ADD REPLY

Login before adding your answer.

Traffic: 2718 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6