get flanking sequences of a SNP in a vcf
1
1
Entering edit mode
6 weeks ago
eric.londin ▴ 60

Hi All,

I think something like this has been asked before, but I having trouble finding a simple solution. Essentially, I'm trying to take a VCF, and for each variant in it, I would like to get the flanking ~25nts of the site on either side of the SNP. I would like there to be 2 sequences, 1 with the reference allele in the middle and a second sequence with the the alternate allele in the middle. I would think that there was some application to do this, without me having to re-invente the wheel!

Thanks for any help

vcf sequences • 266 views
ADD COMMENT
0
Entering edit mode

It will be a pain but you could iterate variant by variant, use each variant as a VCF format entry somehow with bcftools consensus that picks the appropriate region based on POS of that variant. May I ask why you wish to do this?

ADD REPLY

Login before adding your answer.

Traffic: 2448 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