[BioC] how to plot quality report
wang peter
wng.peter at gmail.com
Tue Aug 14 17:19:03 CEST 2012
i tried to use ShortRead to get quality report
but i am confused there is no parameter to set quality score type
phred+33 or phred+64
library(ShortRead)
fastqfile="P11_ACTTGA.fastq"
qa <- qa(dirPath=".", pattern=fastqfile, type=c("fastq"))
report(qa, dest="qcReport1", type="html")
--
shan gao
Room 231(Dr.Fei lab)
Boyce Thompson Institute
Cornell University
Tower Road, Ithaca, NY 14853-1801
Office phone: 1-607-254-1267(day)
Official email:sg839 at cornell.edu
Facebook:http://www.facebook.com/profile.php?id=100001986532253
More information about the Bioconductor
mailing list