簡介
基礎功能
（1）下載fastq.gz文件
wget ftp://ftp.sra.ebi.ac.uk/vol1/fastq/ERR127/ERR127302/ERR127302_1.fastq.gz
wget ftp://ftp.sra.ebi.ac.uk/vol1/fastq/ERR127/ERR127302/ERR127302_2.fastq.gz
下載Fastq文件
ps：下載中途卡死了落君，(⊙﹏⊙)b
補救： 1. we can use ctrl+z to stop it, Note: not ctrl+c.
2. resume it via wget -c url
又可以了鸭丛。吭产。坯钦。
（2）加載數據
library(ShortRead)
sampler <- FastqSampler('ERR127302_1.fastq.gz', 20000)
set.seed(123); ERR127302_1 <- yield(sampler)
sampler <- FastqSampler('ERR127302_2.fastq.gz', 20000)
出現在右上角Global Environment上
（3）基礎功能
fq <- yield(sampler)
fl <- system.file(package="ShortRead", "extdata", "E-MTAB-1147",
+ "ERR127302_1_subset.fastq.gz")
countFastq(fl)
fq <- readFastq(fl)
fq[1:5]
head(sread(fq), 3)
head(quality(fq), 3)
encoding(quality(fq))
基礎功能
