Also added number of matches in bed score column and switched from pyGATB to mappy (no more segmentation faults observed).