cerebis / qc3C

Reference-free quality assessment for Hi-C sequencing data
GNU Affero General Public License v3.0
13 stars 1 forks source link

Expected fraction at 50% GC #13

Closed cerebis closed 5 years ago

cerebis commented 5 years ago

@koadman

For this calculation, the reported value seems quite small. Is this correct?

What are we trying to communicate by mentioning this? That we are seeing far more than what would be random chance?

Eg. For a Hi-C library made using a 4-mer with 4 bp overhang.

The equation is: 1 / np.power(4, site_size / 2)

Which produces an expected fraction at 50% GC of 0.00391.