Question: Write a program that takes as input an alignment, and outputs the PWM and the information content of the PWM. To prevent prob = 0
Write a program that takes as input an alignment, and outputs the PWM and the information content of the PWM. To prevent prob = 0 from happening, add a pseducount 1 for each nucleotide. i.e., p_A = (c_A + 1) / (c_A + c_C + c_G + c_T + 4), where c_A is the number of A's observed in the column. (It is best to implement the main functionality as a function so that you can reuse it in Problem 4 or possibly also Problem 3.) Apply your program on the following alignment and report the results. Also use the web tool enoLogos to compute the sequence logo. Make sure you set %GC to "equiprobable". Click on textout to check the probability matrix and information content there. Do they match your results?
GTGTAGC GTGTTGT GGGTTGC GTGTAGC GTGCTGC GTGTAGC TTGTTGC GTGTAGC ATGTTGC GTGGTGC GTGTTCC GCGTTGC GTGTTAC GCGTTGC GTGCTGC GTTTTGC TTGTTGC GTCTTGC GTGATGC GTGTTGA
Step by Step Solution
There are 3 Steps involved in it
Get step-by-step solutions from verified subject matter experts
