script for NGS analysis