...
The optional tutorial from day 2 (Advanced variant calling tutorial (GVA14)) detailed the use of the breseq pipeline to call variants on clonal samples from an evolving E. Colicoli population. While this tutorial does not require completion of the previous tutorial, many of the finer points breseq are better covered there. This tutorial will focus on the use of the polymorphic mode of breseq to identify variants from a mixed population.
...
All fastq files necessary for this tutorial can be found inside the $BI/gva_course/mixed_population
folder folder. Copy all fastq files to a new folder named "fastq
", and REL606.6.gbk
to to a new folder named "reference
".
Expand | ||
---|---|---|
| ||
|
By default breseq preforms several statistical tests to rule out false positives. To make use of these tests, simply add a -p
flag flag to any breseq command. To highlight what breseq is normally doing by default we will run the same fastq files with and without several of the statistical tests. Specifically, base quality scores, polymorphism scores, polymorphism bias, and minimum strand coverage will be ignored. All 4 of these arguments can be found in the breseq -h output and their values should be set to 0. While running breseq in polymorphism mode is a fairly simple, due to the complexity of the command with turning off all the additional options, it is recommended that you copy paste these commands into a commands file or an idev session.
Expand | ||
---|---|---|
| ||
Be sure to make a new folder named
|
...