[BioC] Error in readSeqFile of qrqc package?
Sang Chul Choi
schoi at cornell.edu
Wed Feb 15 00:00:46 CET 2012
Hi,
I am trying to use "qrqc" package. I had two fastQ files (Illumina), one of which created about a year ago, and another was created a few months ago. The earlier one was read in using readSeqFile, and the latter one was not. I attach the error message from it with R version. The problem is about quality score format. I tried quality="phred" because for earlier one quality scores are Illumina Q score offset (ascii 64), and for latter one is in Sanger FASTQ format, the offset is ascii 33. I am wondering if I could change readSeqFile options to read the second fastQ file that is in Sanger FASTQ format created by Illumina.
Thank you,
SangChul
Error messages.
=========================================================
> fq.file <- readSeqFile("1.fq.gz")
Error in readSeqFile("FASTQ019.fq.gz") :
base quality out of range (0 <= b <= 62) encountered: 35
> fq.file <- readSeqFile("1.fq.gz",quality="phred")
Error in readSeqFile("FASTQ019.fq.gz", :
base quality out of range (4 <= b <= 60) encountered: 61
=========================================================
R version and sessionInfo.
=========================================================
R Under development (unstable) (2012-02-14 r58341)
Copyright (C) 2012 The R Foundation for Statistical Computing
ISBN 3-900051-07-0
Platform: x86_64-apple-darwin10.8.0 (64-bit)
other attached packages:
[1] qrqc_1.3.0 ShortRead_1.13.12 latticeExtra_0.6-19
[4] RColorBrewer_1.0-5 Rsamtools_1.7.33 lattice_0.20-0
[7] Biostrings_2.23.6 GenomicRanges_1.7.24 IRanges_1.13.24
[10] BiocGenerics_0.1.4
=========================================================
More information about the Bioconductor
mailing list