search

Altius / hotspot2

Implementation of hotspot2 by Eric Rynes
16 stars 3 forks source link

Total cut counts, and SPOT, should be restricted to mappable cut counts when possible #19

Closed erynes closed 8 years ago

erynes commented 8 years ago

cutcounts.bash writes the total number of cuts to an output file. Some of those cuts will lie within unmappable regions; cutcounts.bash currently doesn't receive a mappability file as additional input. When hotspot2.sh computes the SPOT score, the numerator is currently restricted to mappable cuts, but the denominator is the unrestricted tally from cutcounts.bash.

Correcting this should yield very little or no difference in the SPOT score (maybe in the 4th significant digit?), so this is not a high priority item.

erynes commented 8 years ago

Closed via commit 5f0db76724995dbc3501b80dab62fcf0aba6cef4.