docs: def readBimFile(basefilename):
      """
      Helper fuinction that reads bim files
      """
      # read bim file
      bim_fn = basefilename+'.bim'
      rv = SP.loadtxt(bim_fn,delimiter='\t',usecols = (0,3),dtype=int)
      return rv