Quality control workflow

Table of Contents Resources Quality control workflow Quality control workflow Setup 1. Path to first example Fastq file 2. Read in Fastq example 3. Show method 4. Class of rfq is an instance of ShortReadQ S4 class 5. Print all S4 class slotNames 6. sread accessor method 7. id accessor method 8. quality accessor method 9. Run Quality control script 10. Create hmtml report 11. Open html report 12. Snippet to Investigate qaSummary structure 13. Generate qa reports for all Fastq files

Resources

• Workshop code
• An Introduction to ShortRead (Quality control page 8)

Quality control workflow

The ShortRead package contains a function called qa() which creates a quality report from several input formats, (usually a Fastq file).

 type=c("SolexaExport", "SolexaRealign", "Bowtie", "MAQMap", "MAQMapShort", "fastq", "BAM")

This function produces similar outputs to FastqQC but executes slightly faster. Functionality has been written to perform quality control either;

1. on a ShortReadQ structure of NGS reads loading into R using readFastq()
2. Specify the directory to run qa() directly from the NGS files.

Quality control workflow

Setup

The first command clears out any existing R objects from your current workspace

 # Clear working directory objects
 rm(list=ls())

Any contributed libraries need to be installed once, but loaded each time they are used with either require(), or library();

require(ShortRead)

Loading help documentation vignette;

 # Getting vignette help
 if(interactive()) vignette("Overview") ## Quality assessment, page 8 

The next commands source some global variables in the file ~/VISG-course-2012/supplementary_QC/globals.R

# Globals
cat("[ Code to source() ]\n")
cat(readLines("globals.R"), sep="\n")

# source globals code
source("globals.R")

cat("[ Directory and Fastq filenames ]\n")
print(indir)
print(fileNames)

1. Path to first example Fastq file

Define the path fo the first Fastq file in fileNames object to the R object 'file1;

 file1 <- file.path(indir, fileNames[1])
 print(file)

2. Read in Fastq example

Read in the first Fastq file into R using ShortReads readFastq() which creates an instance of the ShorReadQ class. This structure contains the sequence data, quality information, and header information;

 rfq       <- readFastq(file1, qualityType="Auto")

3. Show method

Using print() or show() displays auxillary information about the ShortReadQ structure, in this Roche 454 example we can tell that there is a variable number of cycles from the ragged array length: 4936 reads; width: 71..1200 cycles information;

 print(rfq)

4. Class of rfq is an instance of ShortReadQ S4 class

Checking the class;

 class(rfq)

5. Print all S4 class slotNames

Checking the slot names for the S4 class, usually there are accessor methods of the same name as the slots, in this case there are function methods sread(), quality(), and id(), to access slot components of the ShortReadQ object;

 slotNames(rfq)

6. sread accessor method

 Accessor method for the underlying reads, note only a truncated summary of the reads are printed not the entire set of NGS reads;

 sread(rfq)
  A DNAStringSet instance of length 4936
       width seq
   [1]   263 GACTACACGTAGTATGGCTCGCTGTTTCCAAC...GTCTCTCAAGGCACACAGGGGAGTAGGNGNNN
   [2]   161 GACTACACGTAGTATGGCTCGCTGGTTGGCGA...TCAAGGCACACAGGGGATAGGNNNNNNNNNNN
   [3]   196 GACTACACGTAGTATGGGCTCGCTGGAAACCT...GTCGGCGTCTCTCAAGGCACACAGGGGATAGG
   [4]   318 GACTACACGTAGTATGGCTCGCTGCCCGTTCG...CAAGGCACACAGGAGGTAGGAGNNNNNNNNNN
   [5]   405 GACTACACGTAGTATGGGCTCGCTGTAAAGTT...NNNNNNNNNNNNNNTNNTNNNNNNNNNNNNNN
   [6]   231 GACTACACGTAGTATGGGCTCGCTGTCCTACT...GTCGGCGTCTCTCAAGGCACACAGGGGATAGG
   [7]   236 GACTACACGTAGTATGGGCTCGCTGGTAATTC...GTCGGCGTCTCTCAAGGCACACAGGGGATAGG
   [8]   509 GACTACACGTAGTATGGGCTCGCTGGGCCACA...CTCAAGGCACACAGGAGTGAGNNNNNNNNNNN
   [9]   516 GACTACACGTAGTATGGGCTCGCTGCATGGCA...GTACGACGAGNGNGNGNGNGNGNNNNNNNNNN
   ...   ... ...
 [4928]   637 GACTACACGTAGTATGGGCTCGCTGCATATAA...AAGTCGTAAAGTAGTAAAGGTTCGTNTAACGG
 [4929]   516 GACTACACGTAGTATGGGCTCGCTGCCAATTG...GATGGATGCCGAAAAGGGCGTAGGCTTTAAGG
 [4930]   525 GACTACACGTAGTATGGGCTCGCTGTAAATGT...TTTAGGTATTATNTTACATGTAATGCTAAAAC
 [4931]   533 GACTACACGTAGTATGGGCTCGCTGATGCACG...AGGTGGTTCGGTCGGTCGTCCTCGTTAACCGG
 [4932]   631 GACTACACGTAGTATGGGCTCGCTGGTTTTTC...CTACTACGTTATTACGTTTAGTCGTTTTAAGG 
 [4933]   777 GACTACACGTAGTATGGGCTCGCTGTCATTTT...GTACGTTCCGTTTCGTAACGGTCTTTACGGGG
 [4934]   449 GACTACACGTAGTATGGGCTCGCTGCTGAACT...GGTCTCTCAAGGCACACAAGGGGCTAGGNNNN
 [4935]   557 GACTACACGTAGTATGGGCTCGCTGCCCTGGG...AACGAAACTTATAGAAAATTACTNAACTTAAG
 [4936]   621 GACTACACGTAGTATGGGGCTCGCTGATTATT...TTAACGGTAACCGTAGAACCTTTTTTAAAAGG

7. id accessor method

id(rfq)

  A BStringSet instance of length 4936
       width seq
   [1]    14 GYSSWG201ANAAU
   [2]    14 GYSSWG201AKKAN
   [3]    14 GYSSWG201BC6S3
   [4]    14 GYSSWG201BK297
   [5]    14 GYSSWG201BF59S
   [6]    14 GYSSWG201APQ25
   [7]    14 GYSSWG201ARCVR
   [8]    14 GYSSWG201B2QSX
   [9]    14 GYSSWG201BUQXE
   ...   ... ...
 [4928]    14 GYSSWG201AO4TW
 [4929]    14 GYSSWG201BYFJ6
 [4930]    14 GYSSWG201AH034
 [4931]    14 GYSSWG201BKBW3
 [4932]    14 GYSSWG201A3DQB
 [4933]    14 GYSSWG201A6O0R
 [4934]    14 GYSSWG201BUNO4
 [4935]    14 GYSSWG201BTU4B
 [4936]    14 GYSSWG201ALL3I

8. quality accessor method

9. Run Quality control script

10. Create hmtml report

11. Open html report

12. Snippet to Investigate qaSummary structure

13. Generate qa reports for all Fastq files