This past month I have made significant progress to my program. Aside from small modifications, it is near completion. First, it imports the global sequence from a given FASTA file. Next, it calculates the folding energy. Then, it prompts for input of the local sequence length and picks a random piece of that length from the global sequence. Following that, it prompts the user for a threshold and uses a generated Rfold file to generate various arrays that contain the positions in the local sequence that are above a certain threshold. Using these arrays, it makes synonymous changes in the positions that are above the given threshold and then inserts the modified local sequence back into the global sequence. Finally, it reruns mfold on the modified global sequence and compares the two output and keeps the modified sequence if it's ∆G is less the original global sequence's ∆G. Currently, I am running this program and creating graphs and other forms of visual data to document my results.
|
JAY VillariUndergraduate student studying Computer Science and Biology at The College of New Jersey Archives
May 2016
Categories |