Useful script for NGS analysis