Entering edit mode
5.1 years ago
Mirza.Jawad • 0
I have one sequenced bacterial gene (three copies) taken from three different countries. I want to see the SNPs and variant calling (if possible for a single gene) to check out the mutation.
I am from NGS background but could not perform SNPs and variant calling analysis yet.
Can we do variant calling for a single gene? How can find the SNPs with having any reference? Later I will go for the structure prediction and the interactions/docking.
So then what will be the next procedure?
ADD COMMENT/ADD REPLYwhen responding to existing posts to keep threads logically organized.
What format is your data in? Fasta format files? You can use JalView to build multiple sequence alignments and see where the differences are. You can also use MEGA, ClustalW to do the alignments, if you have small number of sequences.
do global alignment with softwares like muscle, t-coffee as sequences are short in number.