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

1 Expert Approved Answer
Step: 1 Unlock blur-text-image
Question Has Been Solved by an Expert!

Get step-by-step solutions from verified subject matter experts

Step: 2 Unlock
Step: 3 Unlock

Students Have Also Explored These Related Databases Questions!