calculate_fst function

Calculate FST