![](https://blogger.googleusercontent.com/img/b/R29vZ2xl/AVvXsEhNj1DJuSXAGumyJ8O_rw0fTt5OMa2OLDqTD7n94hEqkgNIUYKH2-mXMeold4kn44DDlVzwMHIvUs6g3eSvGaz9vqUnbXPTL3I8FfEeg0SHKh4BesP_i9ot5yQF8zW2zoNs6TnY1FPPV9pT/s320/Screen+shot+2011-03-05+at+4.02.03+PM.png)
There is one more thing to say about the HITS project. Using bowtie worked well enough that I feel comfortable posting a plot similar to Fig 2 of the paper. Each point represents one of 1148 genes plotted. The x-axis is saturation (fraction of TA sites in that gene which were hit by the transposon in the library) and the y-axis is the selection index (ratio of total HITS in the lung compared to the library).
It doesn't look exactly like the paper, but it's pretty good.
I couldn't get matplotlib to do the semi-log plot, so I used R:
And that really is it for this. Here is hemR again:
![](https://blogger.googleusercontent.com/img/b/R29vZ2xl/AVvXsEiR4u9t3od9lhO1DCpHnGfRppFb8wM67dmPoAJpUrVbpapyqOux8a20Oc-jHGmbv9iVH8mzG1uOytIP9D4ix1Oa-per1KDgdyIZmyT1uJcMPCbfztmMaKh0RCn29lNnOkzHhb2_iPHYX9po/s320/hemR.png)