breseq  version 0.35.4  revision f352f80f4bc9
mutation predictions | marginal predictions | summary statistics | genome diff | command line log

Predicted mutations
evidence seq id position mutation freq annotation gene description
RA BS168_reference_genome 20,055 G→C 8.7% V30L (GTA→CTA)  pdxT → glutamine amidotransferase for pyridoxal phosphate synthesis; pyridoxal 5'‑phosphate synthase complex, glutamine amidotransferase subunit PdxT
RA BS168_reference_genome 113,389 A→G 7.7% E197G (GAA→GGA)  cysE → serine O‑acetyltransferase
RA BS168_reference_genome 143,659 C→T 100% S100F (TCT→TTT)  rpsE → ribosomal protein S5
RA BS168_reference_genome 165,749 Δ1 bp 100% intergenic (+42/‑5) rrnI‑5S → / → trnI‑Asn ribosomal RNA‑5S/tRNA‑Asn
RA BS168_reference_genome 165,751 Δ1 bp 100% intergenic (+44/‑3) rrnI‑5S → / → trnI‑Asn ribosomal RNA‑5S/tRNA‑Asn
RA BS168_reference_genome 165,825:1 +C 100% noncoding (72/72 nt) trnI‑Asn → tRNA‑Asn
RA BS168_reference_genome 166,037:1 +T 100% intergenic (+4/‑27) trnI‑Gly → / → trnI‑Arg tRNA‑Gly/tRNA‑Arg
RA BS168_reference_genome 166,344 Δ1 bp 100% intergenic (+16/‑156) trnI‑Ala → / → rrnH‑16S tRNA‑Ala/ribosomal RNA‑16S
RA BS168_reference_genome 376,018 Δ1 bp 100% intergenic (‑220/‑14) hxlA ← / → hxlR 3‑hexulose‑6‑phosphate synthase (HPS)/positive regulator of hxlAB expression (formaldehyde sensing)
RA BS168_reference_genome 380,714 G→A 7.3% P1249P (CCG→CCA srfAA → surfactin synthetase
RA BS168_reference_genome 442,499 C→T 10.8% T310M (ACG→ATG)  gabT → 4‑aminobutyrate aminotransferase
RA BS168_reference_genome 557,865:1 +T 100% intergenic (+416/‑8) yddT → / → ydzN putative exported phage protein/hypothetical protein
RA BS168_reference_genome 592,831 G→T 5.5% H125Q (CAC→CAA cesA ← carboxylesterase NP
RA BS168_reference_genome 727,505 C→G 10.3% G272R (GGC→CGC)  opuE ← proline transporter
RA BS168_reference_genome 892,022 A→T 100% N196I (AAT→ATT)  malR → transcriptional activator of the Mal operon
RA BS168_reference_genome 934,016 C→T 13.8% T713I (ACC→ATC)  yfhO → conserved membrane protein
RA BS168_reference_genome 961,032:1 +A 100% intergenic (‑46/‑362) katA ← / → ssuB vegetative catalase 1/aliphatic sulfonate ABC transporter (ATP‑binding protein)
RA BS168_reference_genome 996,888 G→T 13.8% K82N (AAG→AAT yhcU → conserved protein of unknown function
RA BS168_reference_genome 1,006,868 A→C 12.5% Q32P (CAG→CCG)  pgcA → alpha‑phosphoglucomutase
RA BS168_reference_genome 1,044,003 T→C 7.8% I82V (ATC→GTC)  nhaC ← Na+/H+ antiporter
RA BS168_reference_genome 1,073,040 C→T 100% A91A (GCC→GCT yhaI → conserved protein of unknown function
RA BS168_reference_genome 1,073,489 C→T 100% A77T (GCA→ACA)  scoC ← transcriptional regulator of extracellular protease production, sporulation and bacilysin production (MarR family)
RA BS168_reference_genome 1,074,835 C→T 100% G110G (GGG→GGA trpP ← tryptophan transporter or sensor
RA BS168_reference_genome 1,159,412 A→G 5.3% H93H (CAT→CAC yizA ← conserved hypothetical protein (sporulation related)
RA BS168_reference_genome 1,213,393 C→T 8.7% P312S (CCA→TCA)  appF → oligopeptide ABC transporter (ATP‑binding protein)
RA BS168_reference_genome 1,241,555 G→T 16.7% E400D (GAG→GAT kefQ → gated K+/H+ antiporter
RA BS168_reference_genome 1,253,605 T→C 8.1% F74S (TTT→TCT)  spoVIF → sporulation‑specific protein needed for heat resistance
RA BS168_reference_genome 1,255,382 T→C 7.0% Y204C (TAC→TGC)  yjcD ← putative ATP‑dependent DNA helicase (sporulation associated)
RA BS168_reference_genome 1,317,152:1 +GT 100% intergenic (+23/+5) phrA → / ← yjpA secreted inhibitor of the activity of phosphatase RapA (quorum sensing)/putative enzyme
RA BS168_reference_genome 1,317,153:1 +T 100% intergenic (+24/+4) phrA → / ← yjpA secreted inhibitor of the activity of phosphatase RapA (quorum sensing)/putative enzyme
RA BS168_reference_genome 1,364,966 T→C 10.0% F272F (TTT→TTC dppE → dipeptide ABC transporter (dipeptide‑binding lipoprotein)
RA BS168_reference_genome 1,473,747 T→C 11.1% V48A (GTG→GCG)  cheV → coupling protein and response regulator for CheA activity in response to attractants (chemotaxis)
RA BS168_reference_genome 1,474,890 C→T 100% L45L (TTG→TTA kre ← regulator of transcription factor ComK function via modulation of mRNA stability
RA BS168_reference_genome 1,494,589 T→C 9.3% E260E (GAA→GAG sppO ← spore protein cse15
RA BS168_reference_genome 1,690,758 C→A 100% R214S (CGT→AGT)  codY → transcriptional regulator, GTP and BCAA‑dependent
RA BS168_reference_genome 1,742,383 G→A 8.9% D256N (GAT→AAT)  ylxY → putative sugar deacetylase
RA BS168_reference_genome 1,904,636 T→A 8.5% V38V (GTA→GTT ynzK ← putative membrane protein of unknown function (phage origin)
RA BS168_reference_genome 1,910,707 C→T 20.0% P181S (CCG→TCG)  yndF → putative spore germination lipoprotein
RA BS168_reference_genome 1,910,754 C→T 19.9% S196S (TCC→TCT yndF → putative spore germination lipoprotein
RA BS168_reference_genome 2,097,080:1 +A 100% intergenic (‑104/+46) azoJ ← / ← yocK FMN‑dependent NADH‑azoreductase/putative general stress protein
RA BS168_reference_genome 2,271,424 T→C 100% R78R (AGA→AGG uvrX ← lesion bypass phage DNA polymerase; phage SPbeta
RA BS168_reference_genome 2,271,505 C→T 100% V51V (GTG→GTA uvrX ← lesion bypass phage DNA polymerase; phage SPbeta
RA BS168_reference_genome 2,271,523 A→C 100% D45E (GAT→GAG uvrX ← lesion bypass phage DNA polymerase; phage SPbeta
RA BS168_reference_genome 2,387,878:1 +A 100% coding (711/735 nt) yphF ← putative lipoprotein of unknown function
RA BS168_reference_genome 2,403,857 C→G 33.6% K78N (AAG→AAC mecB ← adaptor to ClpC; regulator of competence and sporulation
RA BS168_reference_genome 2,412,532 G→T 7.4% V483F (GTT→TTT)  serA → 3‑phosphoglycerate dehydrogenase
RA BS168_reference_genome 2,480,646 2 bp→AT 100% intergenic (+21/+103) zwf → / ← gndA glucose‑6‑phosphate 1‑dehydrogenase (NADP‑dependent)/NADP+‑dependent 6‑P‑gluconate dehydrogenase
RA BS168_reference_genome 2,480,654 Δ1 bp 100% intergenic (+29/+96) zwf → / ← gndA glucose‑6‑phosphate 1‑dehydrogenase (NADP‑dependent)/NADP+‑dependent 6‑P‑gluconate dehydrogenase
RA BS168_reference_genome 2,480,667 Δ1 bp 100% intergenic (+42/+83) zwf → / ← gndA glucose‑6‑phosphate 1‑dehydrogenase (NADP‑dependent)/NADP+‑dependent 6‑P‑gluconate dehydrogenase
RA BS168_reference_genome 2,581,726:1 +T 100% intergenic (‑109/+45) pstS ← / ← pbpA phosphate ABC transporter (phosphate binding lipoprotein)/transpeptidase (penicillin‑binding protein 2A)
RA BS168_reference_genome 2,633,093 C→T 100% G43E (GGG→GAG)  yqxA ← conserved protein of unknown function
RA BS168_reference_genome 2,770,869 C→T 100% R174C (CGT→TGT)  rsiV → anti‑sigma(V) factor
RA BS168_reference_genome 2,831,221 A→G 12.9% L187S (TTG→TCG)  yrbG ← conserved hypothetical protein
RA BS168_reference_genome 2,842,959 A→G 6.8% Q450R (CAA→CGA)  alaP → sodium/proton‑dependent alanine transporter
RA BS168_reference_genome 2,861,703 C→T 15.3% A16T (GCG→ACG)  mreB ← cell‑shape determining protein
RA BS168_reference_genome 2,907,559 C→A 10.5% L179L (CTG→CTT sdhA ← succinate dehydrogenase (flavoprotein subunit)
RA BS168_reference_genome 3,104,817 A→G 12.1% M76T (ATG→ACG)  leuS ← leucyl‑tRNA synthetase
RA BS168_reference_genome 3,136,949 G→A 6.8% H143Y (CAT→TAT)  ytkA ← putative lipoprotein
RA BS168_reference_genome 3,287,650 A→G 9.7% P2P (CCT→CCC dhbF ← siderophore 2,3‑dihydroxybenzoate‑glycine‑threonine trimeric ester bacillibactin synthetase
RA BS168_reference_genome 3,376,306 A→G 11.1% M451V (ATG→GTG)  mdtP → multidrug‑efflux transporter
RA BS168_reference_genome 3,471,869 T→C 12.5% K203R (AAA→AGA)  yvbG ← putative integral membrane protein
RA BS168_reference_genome 3,503,708 A→G 8.3% V459A (GTC→GCC)  ganA ← short chain beta‑1,4‑galacto‑oligosaccharides beta‑galactosidase (beta‑galacto‑pyranoside hydrolase)
RA BS168_reference_genome 3,770,058:1 +A 100% intergenic (‑250/+46) ywzF ← / ← csbD conserved protein of unknown function/stress response protein
RA BS168_reference_genome 3,935,823 Δ1 bp 100% intergenic (+35/+1) ywbE → / ← rlmI conserved protein of unknown function/large subunit ribosomal RNA m5C1962 methyltransferase
RA BS168_reference_genome 4,153,115 G→T 100% L192I (CTT→ATT)  walK ← two‑component sensor histidine kinase [YycG]
RA BS168_reference_genome 4,155,390:1 +A 100% intergenic (‑181/+43) trnY‑Lys ← / ← purA tRNA‑Lys/adenylosuccinate synthetase
RA BS168_reference_genome 4,174,602 A→T 5.7% S65C (AGT→TGT)  yybL → putative transporter permease subunit
RA BS168_reference_genome 4,191,127 A→G 10.8% intergenic (+1/‑71) yyaM → / → yyaL putative efflux transporter/conserved hypothetical protein
RA Exported 121 G→C 5.4% intergenic (–/‑2857)  / → eryR –//vntifkey=4
RA Exported 131 G→C 6.9% intergenic (–/‑2847)  / → eryR –//vntifkey=4
RA Exported 2,039 C→T 18.6% intergenic (–/‑939)  / → eryR –//vntifkey=4
RA Exported 2,205 C→T 8.0% intergenic (–/‑773)  / → eryR –//vntifkey=4
RA Exported 2,206 C→T 8.0% intergenic (–/‑772)  / → eryR –//vntifkey=4
JC Exported 2,210 +TCTCG 12.5% intergenic (–/‑768)  / → eryR –//vntifkey=4
RA Exported 2,214 C→T 11.0% intergenic (–/‑764)  / → eryR –//vntifkey=4
JC Exported 2,655 Δ3 bp 100% intergenic (–/‑321)  / → eryR –//vntifkey=4
RA Exported 3,729 T→G 6.8% intergenic (+14/‑142) eryR → / → ampR /vntifkey=4//vntifkey=4
RA Exported 3,735 A→T 9.3% intergenic (+20/‑136) eryR → / → ampR /vntifkey=4//vntifkey=4
RA Exported 3,741 G→T 9.8% intergenic (+26/‑130) eryR → / → ampR /vntifkey=4//vntifkey=4
RA Exported 3,783 G→T 5.1% intergenic (+68/‑88) eryR → / → ampR /vntifkey=4//vntifkey=4
JC Exported 3,843 Δ5 bp 100% intergenic (+128/‑24) eryR → / → ampR /vntifkey=4//vntifkey=4
RA Exported 4,625 A→C 6.6% K252T (AAG→ACG)  ampR → /vntifkey=4
RA Exported 4,648 T→A 11.3% Y260N (TAC→AAC) ‡ ampR → /vntifkey=4
RA Exported 4,649 A→C 9.3% Y260S (TAC→TCC) ‡ ampR → /vntifkey=4
RA Exported 5,049 C→T 100% intergenic (+318/‑839) ampR → / → lacZalpha /vntifkey=4//vntifkey=4
RA Exported 5,173 T→G 6.5% intergenic (+442/‑715) ampR → / → lacZalpha /vntifkey=4//vntifkey=4
RA Exported 5,302 A→C 5.4% intergenic (+571/‑586) ampR → / → lacZalpha /vntifkey=4//vntifkey=4
RA Exported 5,464 T→G 5.4% intergenic (+733/‑424) ampR → / → lacZalpha /vntifkey=4//vntifkey=4
RA Exported 6,006 T→G 6.1% V40G (GTT→GGT)  lacZalpha → /vntifkey=4
RA Exported 6,012 A→C 7.3% Q42P (CAA→CCA)  lacZalpha → /vntifkey=4
RA Exported 6,031 A→C 6.4% A48A (GCA→GCC lacZalpha → /vntifkey=4
RA Exported 6,033 A→C 5.0% H49P (CAT→CCT)  lacZalpha → /vntifkey=4
RA Exported 6,256 A→C 5.4% intergenic (+45/–) lacZalpha → /  /vntifkey=4/–
RA Exported 6,277 A→C 5.3% intergenic (+66/–) lacZalpha → /  /vntifkey=4/–

Unassigned missing coverage evidence
   seq id start end size ←reads reads→ gene description
* * ÷ BS168_reference_genome 160976–161238 165590–164903 3666–4615 5 [4] [2] 5 [rrnI‑16S]–rrnI‑23S [rrnI‑16S],rrnI‑23S
* * ÷ BS168_reference_genome 166450–166702 171216–170393 3692–4767 5 [4] [3] 5 rrnH‑16S–[rrnH‑5S] rrnH‑16S,rrnH‑23S,[rrnH‑5S]
* * ÷ BS168_reference_genome 171474–171678 176227–172977 1300–4754 5 [4] [3] 6 rrnG‑16S–[rrnG‑5S] rrnG‑16S,rrnG‑23S,[rrnG‑5S]
* * ÷ BS168_reference_genome 616339 616412 74 5 [3] [4] 5 ydhC putative transcriptional regulator (GntR family)

Unassigned new junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? BS168_reference_genome 2653329 =22 (0.700)5 (0.170) 5/224 NT 26.3% pseudogene (337/471 nt) sigKn RNA polymerase sporulation‑specific sigma‑K factor precursor (Sigma‑27) (N‑terminal half)
?BS168_reference_genome = 2701361 7 (0.230)pseudogene (24/417 nt) sigKc RNA polymerase sporulation‑specific sigma‑K factor precursor (Sigma‑27) (C‑terminal fragment)
* ? BS168_reference_genome = 265333322 (0.700)3 (0.100) 3/224 NT 17.6% pseudogene (341/471 nt) sigKn RNA polymerase sporulation‑specific sigma‑K factor precursor (Sigma‑27) (N‑terminal half)
?BS168_reference_genome 2701367 = 7 (0.230)pseudogene (30/417 nt) sigKc RNA polymerase sporulation‑specific sigma‑K factor precursor (Sigma‑27) (C‑terminal fragment)
* ? Exported = 183878 (0.000)28746 (1.490) 225/230 NT 99.8% intergenic (–/‑1140) –/eryR –//vntifkey=4
?Exported 2572 = 25 (0.000)intergenic (–/‑406) –/eryR –//vntifkey=4