## Load a FASTQ file
s.fastq <- readSeqFile(system.file('extdata', 'test.fastq', package='qrqc'))
## Plot sequence lengths
seqlenPlot(s.fastq)
## Plot sequence legnths before and after trimming
s.trimmed.fastq <- readSeqFile(system.file('extdata',
'test-trimmed.fastq', package='qrqc'))
seqlenPlot(list("not trimmed"=s.fastq, "trimmed"=s.trimmed.fastq))
Run the code above in your browser using DataLab