command line script for querying a pan genome