# #——ShortRead eval = FALSE,整洁= FALSE - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - # # # #使用“ShortRead”包# #库(ShortRead) # # # #从一个文件中创建一个对象来表示一个示例# #取样器< - FastqSampler (ERR127302_1.fastq.gz) # # # #应用方法产生随机样本# # fq < -收益率(取样器)# # # #访问的采样序列读取使用的sread()“# # # #总结核苷酸循环使用# # # #“abc”是一个核苷酸x循环矩阵的计数# # abc < - alphabetByCycle (sread (fq)有趣的核苷酸# # # # # #子集abc < - abc (c (“a”、“c”、“G”、“T”,“N”),) # #——ShortRead-vis, eval = FALSE,整洁= FALSE - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - # # # #创建一个情节从# # # # # # matplot矩阵(T (abc)类型=“l”, # # lty = 1, lwd = 3, # # xlab =“循环”,# # ylab =“计数”,# # cex.lab = 2) # # # #添加一个传说# #传说(“topright”,传说# # = rownames (abc), # # lty = 1, lwd = 3,坳= 1:5,# # cex = 1.8)