NGS-Var2018 Exercise.1

From BITS wiki
Jump to: navigation, search


[ Main_Page | Hands-on_introduction_to_NGS_variant_analysis-2018 | NGS-Var2018 Exercise.2 ]


QC paired end reads using fastQC


ex01_wf.png

You should have learned how to use FastQC during the Introductory training session. We now apply it to the sample data for the sake of time and to build this as a reflex in your Trainee's mind.

Apply FastQC to the FastQ data

You usually receive FastQ data from your sequence provider. It is always good to inspect it prior to starting analyses as it may be biased by artefacts and 'may' in that case need cleaning.

Handicon.png We could also use the fastQC GUI version to do this and will return to GenePattern after that


gpjob_fastqc.png

Q: What if you got a unique fastq file containing interleaved paired reads?

You can ignore this fact and run FastQC on the interleaved data but generally reverse reads of the pair are less pretty than forward reads and you will get the mixture of both in the QC plots. Otherwise, you can easily convert the unique fastQ file to two files, one for each read of a pair by using a number of tools found on the web then run fastQC twice. We recommend a simple bash tool developed by Nathan S. Watson-Haigh that will do the job very fast and can be found HERE. Just copy paste the code to a script under unix terminal then give it your fastq as input.

review fastQC results

The FastQC tool stores all results in a html file and linked pictures. You can review these results after decompressing the zip files created above. The QC results for the first read in pair data (10% sample) ares used to exemplify this.

You should get something like this (read1, read2).


HG001.GRCh38_chr22_0.1_1_fastqc.png


HG001.GRCh38_chr22_0.1_2_fastqc.png

 

Handicon.png Have a look at https://sequencing.qcfail.com/software/fastqc/ to find examples of QC failures

download exercise files

Download exercise files here

Use the right application to open the files present in ex1-files

References:

[ Main_Page | Hands-on_introduction_to_NGS_variant_analysis-2018 | NGS-Var2018 Exercise.2 ]