Smith - Programs

GC Skew

I modified an R script that calculates GC skew using a sliding window method from the seqinR manual.  

To run this script you first need to load the seqinR package.  When working with R programs I use Rstudio.

The fasta file must then be read into the environment and stored in a variable.  The variable I use for the reminder of my program is "myseq".

 myseq <- read.fasta(file = "Desktop/pMP.fasta", as.string = FALSE, forceDNAtolower = TRUE, set.attributes = FALSE, seqonly = TRUE, strip.desc = TRUE)

After that I simply load my R script into Rstudio, cmd+A, and click the Run function. 

The result: