1. mirnylab
  2. hiclib
  3. Issues
Issue #3 resolved

enzyme-specific fragmentLarge

Anonymous created an issue

this might be better as enzyme-specific, current parameters are good for HindIII (ie. 6-cutters)

Comments (2)

  1. mimakaev

    Perhaps, the low cutoff will need to be adjusted a bit, and a high cutoff may be adjusted too.

    Reads with fragments which are much less than 100 bp would likely spawn multiple restriction sites, and 50 is already the read length.

    Do you have any suggestion on what should these numbers be for a 4-cutters?

    I can do something like high = 10 * MedianFragentDistance, low = max(50, 0.05 * medianFragmentDistance)

  2. Log in to comment