[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