Post date: Aug 26, 2016 2:59:49 PM
The previous matlab script was having problems coping with two transposons of opposite orientations inserted very close to each other. This caused few transposons to be counted twice. This did not affect any of the conclusion but we anyway produce a corrected script here for the future.
Now the script first sorts every read according to its orientations, processes all reads in direct orientation, then all reads in reverse orientation, so that no more confusion happens.