[BioC] about the quality score

wang peter wng.peter at gmail.com
Wed Jan 11 18:11:01 CET 2012


i have data, like this
@HWI-ST397_0000:5:1101:1110:2156#NTNTCG/1
CAAGGAGCAGAGTGAGGTCTCTCCGCAGAGCGCCAAGGCAAAGCAGATGGTGGAGCTGCTCTACGAGACGGCGCACATTGCGTCTGGATTCCCGCCTG
+HWI-ST397_0000:5:1101:1110:2156#NTNTCG/1
gdgggegc_gffdeeggdgfggggdgbfadaeeeeefaefgcebebdffebbbbcg`gdgfgeeebgbddcbccbfdffcdcadbW_]]`]dd_bb\^
reads <- readFastq(fastqfile);
seqs <- sread(reads)
it seems phred+64 system

can i use the function:
qual <- PhredQuality(quality(quality(reads)))
i think can only read phred+32 system

-- 
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