Entering edit mode
I am using following scripts to read fasta files
>library(Biostrings)
>dna <- readDNAStringSet("<<PATH TO FASTA FILE>>")
But, further I would like to extract SNPs from these alignment file, but I don't know how extract the SNPs.
Does anyone know?