Question: Generate a VCF File for Different Reference Sequences
gravatar for pld
5.6 years ago by
United States
pld4.8k wrote:

Is there a means for taking two different reference sequences (e.g. two different virus isolates) and generating a VCF mapping those SNPs? In the sense that I want to see the differences in B with respect to A. I'd like to be able to overlay (in IGV) the differences in a second isolate against the SNPs I found after assembling to the primary isolate of interest.


virus snp igv ngs assembly • 2.1k views
ADD COMMENTlink modified 5.6 years ago by Istvan Albert ♦♦ 82k • written 5.6 years ago by pld4.8k
gravatar for Istvan Albert
5.6 years ago by
Istvan Albert ♦♦ 82k
University Park, USA
Istvan Albert ♦♦ 82k wrote:

In general the VCF format is not suited for describing an alignment between genomes. Nor is IGV a good tool to visualize that. 

What you could do however if your really wanted to use both of these is to generate reads from one of your genomes with wgsim, align these reads to the other genome then call SNPs from the alignment.

ADD COMMENTlink modified 5.6 years ago • written 5.6 years ago by Istvan Albert ♦♦ 82k

I ended up writing a script that converts a clustal formatted alignment (clustal omega) to VCF. I realize VCF isn't the best tool for this, most of the features of VCF aren't used in this case (depth, quality, etc). The only reason I went with VCF is that it is a format supported by many visualization tools.

The goal was to have something where I can visualize the snps of my isolate of species A with respect to reference A, and visualize the SNPs between reference A and reference B at the same time. IGV is just for me, the end game is to make a figure to show that none of the SNPs in our isolate (species A), match up with SNPs found in species B. I'd like to include a sashimi plot in this to show that the positions where the A-B SNPs are sufficiently covered. There's some concern (unreasonably imo), that we have some chocolate in our peanut butter.


ADD REPLYlink modified 5.6 years ago • written 5.6 years ago by pld4.8k

It would be better to convert your alignment to SAM, which can be visualized, too.

ADD REPLYlink written 5.6 years ago by lh332k

Hi Joe. Do you have this script available somewhere? :')

ADD REPLYlink written 4.5 years ago by Coryza380
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: 1847 users visited in the last hour