Convert classify_snps function to R script
Update USAGE doc