Question: Best way to look at nonsynonymous changes in a collection of sequences?
gravatar for nchuang
4.2 years ago by
United States
nchuang210 wrote:

Hi guys,

I have a collection of 5kb human sequences that I want to compare their differences, particularly if base changes result in nonsynonymous changes. I know there are tools such as VEP, ANNOVAR, etc that annotates changes, but they all seem to require a VCF first. 

What is the simplest way to generate a VCF of their differences. I figure since there is no reference sequence for this region, I did a multiple sequence alignment. Perhaps I can make a consensus sequence and use that as reference for variant calling?

Since what I have are just Fasta sequences I am not sure how to further proceed. I did try jvarkit tool, biostar94573 from a previous thread, but it does not seem to work for me.



snp • 813 views
ADD COMMENTlink written 4.2 years ago by nchuang210

What about MEGA?

ADD REPLYlink written 4.2 years ago by Biomonika (Noolean)3.1k

for reference:

ADD REPLYlink modified 10 weeks ago by RamRS25k • written 4.2 years ago by Pierre Lindenbaum126k
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: 735 users visited in the last hour