Orient allele counts toward the predicted stronger-binding allele

Given reference / alternative counts and a per-site PWM binding delta (delta = score(ALT) - score(REF)), returns the read counts on the predicted "high-binding" (H) and "low-binding" (L) allele. Sites with delta > 0 bind better on the ALT allele, so H = alt; sites with delta < 0 bind better on REF, so H = ref.

ascOrientReads(refV, altV, delta)

Arguments

refV, altV

integer vectors of reference / alternative counts (aligned).

delta

numeric vector, PWM score difference alt - ref.

Value

list(H, L) oriented count vectors.

Author

Irem B. GUNDUZ