scripts for pulling out consensus SNPs from reference and dealing with strand alignment