library (ShortRead) genome <- readDNAStrings(Fasta file) sread(genome)[1] quality(genome)[1] ##PhredQuality instance pq <- PhredQuality (quality(fqsample)) # transform encoding into scores qs <- as(pq, "IntegerList") qs # Print score qaSummary <- qa(fqsample, lane = 1) #class : ShortReadQQA(1D) #Names accessible with the quality assessment summary names (qaSummary)
Preview:
downloadDownload PNG
downloadDownload JPEG
downloadDownload SVG
Tip: You can change the style, width & colours of the snippet with the inspect tool before clicking Download!
Click to optimize width for Twitter