I am using following scripts to read fasta files
>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?
Login before adding your answer.
Use of this site constitutes acceptance of our User Agreement and Privacy