Barrick Lab :: breseq

<noautolink>
<!--
   * Set PAGETITLE = Barrick Lab :: breseq
-->
---+ <b><i>breseq</i></b>
<div style="float: right; margin:10px;"><img style="border-color: black;  border-style: solid; border-width: 1px;" src="%PUBURL%/%WEB%/ToolList/breseq_icon.gif"></div>

__breseq__ is a computational pipeline for finding mutations relative to a reference sequence in short-read DNA re-sequencing data for haploid microbial-sized genomes. __breseq__ is a command line tool implemented in C++ and R. It will compile and function on a variety of Unix platforms, including <nop>MacOSX and Cygwin. It reports single-nucleotide mutations, point insertions and deletions, large deletions, and new junctions supported by mosaic reads (such as those produced by new mobile element insertions) in an annotated HTML format.

%TABLE{columnwidths="100,250"}%
| *Learn more* | %ICON{viewtopic}% [[%ATTACHURL%/documentation][Online __breseq__ documentation]] |

%TABLE{columnwidths="100,250"}%
| *Download* | %ICON{download}% [[https://github.com/barricklab/breseq/releases/latest][Latest version from GitHub]] |

---++ Comprehensive Analysis of Mutations in Evolved Microbes Using _breseq_
*ASM Microbe 2021 On-Demand Workshop (June 20-24)* %BR%
Introduction: Identifying Mutations and Studying Microbial Genome Evolution with breseq %BR% 
<a href="https://utexas.hosted.panopto.com/Panopto/Pages/Viewer.aspx?id=78b03d53-a98f-4f4c-8343-ad6a01540017"> %ATTACHURL%/ASM_Recording_Poster.png </a>

%ICON{mov}% [[https://utexas.hosted.panopto.com/Panopto/Pages/Viewer.aspx?id=78b03d53-a98f-4f4c-8343-ad6a01540017][ Watch Video]] %ICON{pdf}% [[%ATTACHURL%/ASM_Microbe_2021_slides_Barrick.pdf][Download Slides]]

---+++ Zoom Workshop: Introductory Topics (July 20, 2021)

%ICON{pdf}% [[%ATTACHURL%/2021_Introductory_Workshop.pdf][Download Slides]]

Example 1a: Analyzing an evolved _E. coli_ clone with a high quality reference sequence for its ancestor (LTEE Ara+1 50,000 generations, Clone A)
%BR%<code>breseq -p -l 80 -r REL606.gbk SRR2584524.fastq.gz</code>
%BR%[[https://barricklab.org/twiki/pub/Lab/ToolsBacterialGenomeResequencing/IntroWorkshop/REL11392_Ara+1_50K_clone_A/][View Results]]

Example 1b: What the results look like if you run this same clonal sample in polymorphism mode (LTEE Ara+1 50,000 generations, Clone A)
%BR%<code>breseq -p -l 80 -r REL606.gbk SRR2584524.fastq.gz</code>
%BR%[[https://barricklab.org/twiki/pub/Lab/ToolsBacterialGenomeResequencing/IntroWorkshop/REL11392_Ara+1_50K_clone_A_polymorphism/][View Results]]

Example 2: Results for another evolved clone that was sequenced with longer reads (LTEE Ara+1 50,000 generations, Clone B)
%BR%<code>breseq -r REL606.gbk SRR2584534_1.fastq.gz SRR2584534_2.fastq.gz</code>
%BR%[[https://barricklab.org/twiki/pub/Lab/ToolsBacterialGenomeResequencing/IntroWorkshop/REL11393_Ara+1_50K_clone_B/][View Results]]

Example 3: Analyzing the mixed population that both of these clones were isolated from (LTEE Ara+1 50,000 generations, Population)
%BR%<code>breseq -j 8 -p -r REL606.gbk SRR6173952_1.fastq.gz SRR6173952_2.fastq.gz</code>
%BR%[[https://barricklab.org/twiki/pub/Lab/ToolsBacterialGenomeResequencing/IntroWorkshop/REL11383_Ara+1_50K_population/][View Results]]

Example 4: Results from mapping to reference genome of a closely related strain–many predictions (links removed to save disk space).
%BR%<code>breseq -r NC_000913.3.MG1655.gbk SRR2584534_1.fastq.gz SRR2584534_2.fastq.gz</code>
%BR%[[https://barricklab.org/twiki/pub/Lab/ToolsBacterialGenomeResequencing/IntroWorkshop/REL11393_Ara+1_50K_clone_B_K12/][View Results]]

Example 5: Analyzing an _E. coli_ cell that contains a plasmid
%BR%<code>breseq -r E._coli_W3110_NC_007779.1.gbk -r GFP_Plasmid_SKO4.gbk AR_E1_GTTTCG_L005_R2_001.fastq.gz AR_E1_GTTTCG_L005_R1_001_1.fastq.gz AR_E1_GTTTCG_L005_R1_001.fastq.gz AR_E1_GTTTCG_L005_R2_001_1.fastq.gz</code>
%BR%[[https://barricklab.org/twiki/pub/Lab/ToolsBacterialGenomeResequencing/IntroWorkshop/Ecoli_with_plasmid/][View Results]]

Example 6a: Locating the insertion site of an integration cassette in the _A. baylyi_ genome using a junction reference (best option)
%BR%<code>breseq --junction-only-reference pBTK622_tdk-kanR_cassette_for_Golden_Transformation.gbk -r Acinetobacter-baylyi-ADP1-WT.gff3 G2_CCGTCC_L007_R1_001.fastq.gz G2_CCGTCC_L007_R2_001.fastq.gz</code>
%BR%[[https://barricklab.org/twiki/pub/Lab/ToolsBacterialGenomeResequencing/IntroWorkshop/Abaylyi_with_cassette/][View Results]]

Example 6b: Same sample not using junction reference
%BR%<code>breseq -r pBTK622_tdk-kanR_cassette_for_Golden_Transformation.gbk -r Acinetobacter-baylyi-ADP1-WT.gff3 G2_CCGTCC_L007_R1_001.fastq.gz G2_CCGTCC_L007_R2_001.fastq.gz</code>
%BR%[[https://barricklab.org/twiki/pub/Lab/ToolsBacterialGenomeResequencing/IntroWorkshop/Abaylyi_with_cassette_normal/][View Results]]

Example 6c: Same sample using targeted mode
%BR%<code>breseq -t -r pBTK622_tdk-kanR_cassette_for_Golden_Transformation.gbk -r Acinetobacter-baylyi-ADP1-WT.gff3 G2_CCGTCC_L007_R1_001.fastq.gz G2_CCGTCC_L007_R2_001.fastq.gz</code>
%BR%[[https://barricklab.org/twiki/pub/Lab/ToolsBacterialGenomeResequencing/IntroWorkshop/Abaylyi_with_cassette_targeted/][View Results]]

Example 7a: Mapping to a de novo assembled reference
%BR%<code>breseq -r assembly.fasta SRR2584534_1.fastq.gz SRR2584534_2.fastq.gz</code>
%BR%[[https://barricklab.org/twiki/pub/Lab/ToolsBacterialGenomeResequencing/IntroWorkshop/REL11393_Ara+1_50K_clone_B_assembly/][View Results]]

Example 7b: Mapping to a de novo assembled reference using contig mode (best option)
%BR%<code>breseq -c assembly.fasta SRR2584534_1.fastq.gz SRR2584534_2.fastq.gz</code>
%BR%[[https://barricklab.org/twiki/pub/Lab/ToolsBacterialGenomeResequencing/IntroWorkshop/REL11393_Ara+1_50K_clone_B_assembly_contig/][View Results]]

---++ Feedback and Bug Reporting

Please report bugs on [[https://github.com/barricklab/breseq/issues][GitHub]].

---++ Tutorial input and results files for _breseq_ book chapter

*Deatherage, D.E.*, *Barrick, J.E.*. (2014) Identification of mutations in laboratory-evolved microbes from next-generation sequencing data using _breseq_. _Methods Mol. Biol._ *1151*:165-188. [[http://www.ncbi.nlm.nih.gov/sites/entrez?Db=pubmed&Cmd=ShowDetailView&TermToSearch=24838886][&laquo;PubMed&raquo;]]

| *Download Link* | *Size* | *Description* |
| [[%ATTACHURL%/Clonal_Sample.tgz][%ICON{download}% <literal>Clonal_Sample.tgz</literal>]] | 267M | Input files for clonal sample |
| [[%ATTACHURL%/Clonal_Output.tgz][%ICON{download}% <literal>Clonal_Output.tgz</literal>]] | 1.2M | Output files for clonal sample |
| [[%ATTACHURL%/Population_Sample.tgz][%ICON{download}% <literal>Population_Sample.tgz</literal>]] | 1.4G | Input files for time-course of population samples |
| [[%ATTACHURL%/Population_Output.tgz][%ICON{download}% <literal>Population_Output.tgz</literal>]] | 11M | Output files for time-course of population samples, including comparison |
| [[%ATTACHURL%/Mutation_Examples.tgz][%ICON{download}% <literal>Mutation_Examples.tgz</literal>]] | 199K | Examples of using unassigned evidence to predict complex mutations |
| [[%ATTACHURL%/Poor_Evidence_Examples.tgz][%ICON{download}% <literal>Poor_Evidence_Examples.tgz</literal>]] | 139K | Example output showing characteristics of poor evidence | 

---++ Other Resources

*Previous Versions:*
| %ICON{download}% [[http://barricklab.org/release/breseq][Previous versions from 0.24rc6 to 0.25d (barricklab.org)]] |
| %ICON{download}% [[http://code.google.com/p/breseq/downloads/list][Previous versions up to 0.24rc6 (Google Code)]] |

*Additional test data and an example of __breseq__ output:*
| [[%ATTACHURL%/REL8593A.fastq.gz][%ICON{download}% FASTQ reads for _E. coli_ strain REL8593A (200M)]] | 
| [[%ATTACHURL%/REL606.gbk.gz][%ICON{download}% Reference genome (REL606)]] |
| [[%ATTACHURL%/REL8593A_output][%ICON{external}% Example of __breseq__ output]] |
This topic: Lab > ToolList > ToolsBacterialGenomeResequencing
Topic revision: r54 - 2021-07-20 - JeffreyBarrick