magicDGS / thaplv

Tools for Haploid/Haplotype VCFs
BSD 3-Clause "New" or "Revised" License
0 stars 0 forks source link

Availables sites computed in PairwiseDifferencesWindow could be wrong #26

Open magicDGS opened 8 years ago

magicDGS commented 8 years ago

When counting the number of available sites in PairwiseDifferencesWindow, the current algorithm use interval.size() - counter.countNsRegion(interval) - 1 to maintain compatibility regarding previous (and not tested) implementations.

Thus, we should invertigate if this behaviour is what we expect. The available sites should be the same as the bases included in the window if the number of Ns within the window is 0. This is related with #23 because we will require to change this if it is not working.