Question: how to change pred+33 to pred+44 in shortread
gravatar for wang peter
3.8 years ago by
wang peter2.0k
wang peter2.0k wrote:
library(ShortRead); reads <- readFastq(fastqfile); seqs <- sread(reads); qual <- quality(quality(reads)); score_sys = data.class(quality(reads));# cat("the quality score system (SFastqQuality=Phred+64,FastqQuality=Phred+33) is",score_sys,"\n") it is a Phred+64 system, how can i change the qual to Phred+33. then i can make new reads by reads <- ShortReadQ(sread=seqs, quality=qual, id= BStringSet(paste0(sample, "-", first_number:last_number) ) ); -- shan gao Room 231(Dr.Fei lab) Boyce Thompson Institute for Plant Research Cornell University Tower Road, Ithaca, NY 14853-1801 Office phone: 1-607-254-1267(day) Official Facebook: [[alternative HTML version deleted]]
ADD COMMENTlink written 3.8 years ago by wang peter2.0k
Please log in to add an answer.


Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.2.0
Traffic: 99 users visited in the last hour