# NOT RUN {
# }
# NOT RUN {
# read DNA from fasta file
fasta <- fasta.read("virus.fasta", "host.fasta")
fasta.v <- fasta[[1]]
enc.df.virus <- ENc.values(fasta.v)
# }
# NOT RUN {
gc.df <- GC.content(fasta.v)
ENc.GC3plot(enc.df.virus, gc.df)
# }
Run the code above in your browser using DataLab