hw2 4 .docx - COMP.4290/5510 HW #2 DUE ON 2/7/22 G-C skew...

COMP.4290/5510 HW #2 DUE ON 2/7/22
G-C skew is examined in the textbook.
One can take this further into oligomers (multiple nucleotides).
Consider a genome sequence G, partitioned into 50 equal-length segments. In each segment, the
frequency of each k-mer m is obtained. Let m * be the reverse complement of m .
Now, define
S
i
= log ( f
m
/f
m*
),
as a measure of the skew. If frequency is equal to 0, it is set to 1.
Write a program (preferably in Python) that produce a table of S
i
for each k = 1, 2, and 3. Each k will have
4