@iskandr is changing a few things around. Also, his comment from before:
My only suggestion would be to maybe increase the minimum MAPQ to 1 and if you're ultimately using the counts of neoepitopes with >= 3 spanning reads then this will be very sensitive to the degree of exonic coverage. Maybe normalize by number of reads mapping to exons?
@iskandr is changing a few things around. Also, his comment from before: