Question: How to add sequence information to vcf from MPileup
gravatar for helge.fluess
21 months ago by
helge.fluess0 wrote:

Hello everyone,

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

rna-seq mpileup snp bam vcf • 533 views
ADD COMMENTlink modified 21 months ago • written 21 months ago by helge.fluess0

Thanks a lot, this should provide exactly what i need. I should've mentioned though, that i work within Galaxy and that I don't have very much experience with python or coding in general. I'll give the java tool a try and hope i can get it to work. Wish me luck :)

ADD REPLYlink written 21 months ago by helge.fluess0

May the force be with you..Check if your plant genome is present in the pull down list and upload the vcf file here to get flanking regions:

ADD REPLYlink modified 21 months ago • written 21 months ago by cpad011213k

You may want to post this over at Galaxy biostars to see if they have any other suggestions related to doing this is galaxy.

ADD REPLYlink written 21 months ago by genomax83k
Please log in to add an answer.


Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 1783 users visited in the last hour