Search
Question: how to change pred+33 to pred+44 in shortread
0
gravatar for wang peter
3.5 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 email:sg839@cornell.edu Facebook:http://www.facebook.com/profile.php?id=100001986532253 [[alternative HTML version deleted]]
ADD COMMENTlink written 3.5 years ago by wang peter2.0k
Please log in to add an answer.

Help
Access

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