I'm new to this community and looking for some help. I applied a variant calling pipeline with RNAseq data of four plant samples. Thus, after mapping each sample to a reference genome using Tophat, i used MPileup for variant calling, including all four samples. While I am content with the results, i wonder if there is any possibility to add the sequence information from the BAM files, to the VCF output of MPileup. I i get it right, the VCF output gives me information about the position of the SNP/Variant on the reference chromosome/contig. Since my goal is to create KASP assays from the SNP information, I'd like to add 100bp in before and after the called SNP. Is this possible, without having to find each position in the contig by hand?
Maybe someone had the same issue before?
Thanks in advance Cheers Helge