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 87,944 T→A 6.0% V104E (GTG→GAG)  dusB → tRNA‑dihydrouridine synthase B
RA BS168_reference_genome 104,346 G→T 5.8% D259Y (GAT→TAT)  clpC → class III stress response‑related ATPase, AAA+ superfamily
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 196,543 A→G 9.3% K111E (AAA→GAA)  cdaA → diadenylate cyclase
RA BS168_reference_genome 238,265 G→T 21.1% pseudogene (184/285 nt) ybfEc ← conserved protein of unknown function (C‑terminal end)
RA BS168_reference_genome 291,913 A→G 11.5% Q333Q (CAA→CAG ansZ → L‑asparaginase 2 (putative lipoprotein)
RA BS168_reference_genome 353,749 A→G 13.0% I298V (ATT→GTT)  ffoR → Fur‑regulated NADPH:ferredoxin oxidoreductase
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 404,950 T→C 5.4% E22E (GAA→GAG ycxB ← hypothetical protein
RA BS168_reference_genome 419,079 G→A 7.6% D363N (GAC→AAC)  yclG → putative uronase
RA BS168_reference_genome 434,025 C→T 9.2% V237V (GTC→GTT pbtO → petrobactin iron‑siderophore ABC transporter (permease)
RA BS168_reference_genome 452,995 G→A 6.3% G56S (GGC→AGC)  ycsA → putative tartrate dehydrogenase
RA BS168_reference_genome 501,751 G→A 5.8% S58N (AGC→AAC)  ydbI → putative integral inner membrane protein
RA BS168_reference_genome 557,865:1 +T 100% intergenic (+416/‑8) yddT → / → ydzN putative exported phage protein/hypothetical protein
RA BS168_reference_genome 851,659 C→T 5.5% G431G (GGC→GGT treP → phosphotransferase system (PTS) trehalose‑specific enzyme IIBC component
RA BS168_reference_genome 874,963 T→C 7.7% intergenic (+161/‑465) yfjO → / → yfzA putative RNA methyltransferase/conserved hypothetical protein
RA BS168_reference_genome 892,022 A→T 100% N196I (AAT→ATT)  malR → transcriptional activator of the Mal operon
RA BS168_reference_genome 900,276 C→T 7.3% A66V (GCG→GTG)  yfiG → putative sugar/proton transporter
RA BS168_reference_genome 922,964 Δ1 bp 7.5% coding (347/882 nt) yfhB → putative isomerase
RA BS168_reference_genome 944,526 A→G 92.1% K14E (AAG→GAG)  perR → transcriptional regulator (Fur family)
RA BS168_reference_genome 960,972 T→G 10.8% K5N (AAA→AAC katA ← vegetative catalase 1
RA BS168_reference_genome 966,494 G→C 26.1% A59G (GCG→GGG)  ygaO ← putative integral inner membrane protein
RA BS168_reference_genome 981,799 G→A 6.7% V66I (GTT→ATT)  yhcG → putative ABC transporter ATP‑binding protein
RA BS168_reference_genome 1,023,150 C→A 29.3% intergenic (+50/‑200) yhdF → / → bcaP putative NAD(P)‑dependent dehydrogenase/branched‑chain amino acid transporter or methylthioribose‑degradation metabolite exporter
RA BS168_reference_genome 1,067,071 G→A 6.1% W332* (TGG→TAG)  sbcE → ATPase involved in DNA double strand break repair and recombination
RA BS168_reference_genome 1,073,367:1 +A 6.9% coding (351/612 nt) scoC ← transcriptional regulator of extracellular protease production, sporulation and bacilysin production (MarR family)
RA BS168_reference_genome 1,073,367:2 +A 6.9% coding (351/612 nt) scoC ← transcriptional regulator of extracellular protease production, sporulation and bacilysin production (MarR family)
RA BS168_reference_genome 1,204,963 T→A 22.8% intergenic (‑47/‑202) yjzD ← / → yjaU forespore targeted protein/conserved hypothetical protein
RA BS168_reference_genome 1,220,375 A→C 16.7% N176T (AAC→ACC)  oppA → oligopeptide ABC transporter (binding lipoprotein)
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,320,397 T→G 17.4% H236Q (CAT→CAG xpdC → phage PBSX; manganese‑containing peroxidase
RA BS168_reference_genome 1,338,342 A→C 15.3% N1126T (AAC→ACC)  xkdO → phage PBSX; putative tape measure protein
RA BS168_reference_genome 1,349,495 G→T 100% A325A (GCC→GCA pitB ← putative low‑affinity inorganic phosphate transporter
RA BS168_reference_genome 1,450,050 A→G 5.8% A267A (GCA→GCG ykvU → spore membrane protein involved in germination
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,511,878 C→A 24.8% Q191K (CAA→AAA)  sipT → type I signal peptidase
RA BS168_reference_genome 1,511,930 G→T 23.2% S88* (TCA→TAA)  ykoA ← putative membrane protein of unknown function
RA BS168_reference_genome 1,522,146 T→C 10.8% F266L (TTT→CTT)  adeC → adenine deaminase
RA BS168_reference_genome 1,616,534 A→C 18.4% T90P (ACG→CCG)  ylyA → sporulation‑related RNA polymerase‑binding factor
RA BS168_reference_genome 1,620,410 T→C 7.4% intergenic (+80/‑66) pyrP → / → pyrB uracil permease/aspartate carbamoyltransferase
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,719,477 T→C 10.2% F235L (TTC→CTC)  tsf → elongation factor Ts
RA BS168_reference_genome 1,747,662 A→G 6.8% V180V (GTA→GTG dapG → aspartokinase I (alpha and beta subunits)
RA BS168_reference_genome 1,749,907 G→T 100% D164Y (GAC→TAC)  rnjB → dual activity 5' exo‑and endoribonuclease J2
RA BS168_reference_genome 1,915,114 T→C 15.4% intergenic (+119/‑107) yndK → / → pghL putative phage/plasmid replication protein/phage gamma‑polyglutamate hydrolase
RA BS168_reference_genome 1,923,465 Δ1 bp 5.6% coding (232/708 nt) ccdA → cytochrome c‑type biogenesis protein CcdA; thiol‑disulfide oxido‑reductase
RA BS168_reference_genome 1,999,126 A→T 24.6% A230A (GCT→GCA dacC ← D‑alanyl‑D‑alanine carboxypeptidase
RA BS168_reference_genome 2,070,110 T→C 7.5% intergenic (+135/+134) bsrE → / ← yobJ type I toxin (BsrE/AsrE)/conserved protein of unknown function; putative defective prophage 6
RA BS168_reference_genome 2,096,912 A→G 8.1% L22P (CTT→CCT)  azoJ ← FMN‑dependent NADH‑azoreductase
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,128,562 C→A 7.1% G169W (GGG→TGG)  yodD ← putative hydrolase
RA BS168_reference_genome 2,232,019 T→A 5.5% Y375* (TAT→TAA yonF → putative P‑loop ATPase; phage SPbeta
RA BS168_reference_genome 2,247,591 A→G 9.9% intergenic (‑48/‑298) youB ← / → yomK conserved phage protein of unknown function; phage SPbeta/conserved protein of unknown function; phage SPbeta
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,296,578 C→A 8.7% V9F (GTT→TTT)  ypkP ← putative 1‑acyl‑sn‑glycerol‑3‑phosphate acyltransferase
RA BS168_reference_genome 2,306,179 G→A 6.9% G268S (GGC→AGC)  metAA → homoserine O‑acetyltransferase
RA BS168_reference_genome 2,318,201 C→A 27.7% V415L (GTA→TTA)  pbuX ← xanthine permease
MC JC BS168_reference_genome 2,344,013 Δ39 bp 100% coding (2570‑2608/2745 nt) ponA → peptidoglycan glycosyltransferase (penicillin‑binding proteins 1A and 1B)
RA BS168_reference_genome 2,375,679 A→G 15.9% I73I (ATT→ATC) ‡ trpD ← anthranilate phosphoribosyltransferase
RA BS168_reference_genome 2,375,681 T→C 15.3% I73V (ATT→GTT) ‡ trpD ← anthranilate phosphoribosyltransferase
RA BS168_reference_genome 2,375,710 T→C 15.9% H63R (CAC→CGC)  trpD ← anthranilate phosphoribosyltransferase
RA BS168_reference_genome 2,375,721 T→C 14.7% A59A (GCA→GCG trpD ← anthranilate phosphoribosyltransferase
RA BS168_reference_genome 2,375,727 C→T 16.0% V57V (GTG→GTA trpD ← anthranilate phosphoribosyltransferase
RA BS168_reference_genome 2,375,733 A→G 18.2% G55G (GGT→GGC trpD ← anthranilate phosphoribosyltransferase
RA BS168_reference_genome 2,375,738 C→T 18.1% A54T (GCG→ACG)  trpD ← anthranilate phosphoribosyltransferase
RA BS168_reference_genome 2,375,739 A→T 20.0% L53L (CTT→CTA trpD ← anthranilate phosphoribosyltransferase
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,394,731 T→A 9.5% K361M (AAG→ATG)  rpfA ← RNA degradation presenting factor (ribosomal protein S1 homolog)
RA BS168_reference_genome 2,421,950 C→T 6.5% A75T (GCG→ACG)  rluB ← 23S rRNA pseudouridine 2633 (=2605 standard) pseudouridine synthase
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,555,006 A→G 10.0% I81T (ATA→ACA)  tapA ← lipoprotein for biofilm formation
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,595,353 T→C 9.9% T50A (ACG→GCG)  cshB ← ATP‑dependent RNA helicase; cold shock
RA BS168_reference_genome 2,749,256 A→T 6.5% intergenic (‑439/‑4) csn ← / → yraL chitosanase/conserved protein of unknown function
RA BS168_reference_genome 2,758,999 A→G 9.7% S360P (TCT→CCT)  sacC ← levanase
RA BS168_reference_genome 2,839,474 C→A 29.2% K133N (AAG→AAT yrbE ← putative inositol‑related oxidoreductase
RA BS168_reference_genome 2,912,810 T→G 15.4% T27P (ACT→CCT)  uvrC ← excinuclease ABC (subunit C)
RA BS168_reference_genome 2,987,550 T→C 7.1% intergenic (‑3/+181) pfkA ← / ← accA 6‑phosphofructokinase/acetyl‑CoA carboxylase (carboxyltransferase alpha subunit)
RA BS168_reference_genome 3,021,860 C→A 8.2% L210M (CTG→ATG)  ppnKB → inorganic polyphosphate/ATP‑NAD kinase
RA BS168_reference_genome 3,036,709 C→T 6.3% R388H (CGC→CAC)  tyrS ← tyrosyl‑tRNA synthetase
RA BS168_reference_genome 3,096,417 C→T 8.9% A50T (GCG→ACG)  ytaP ← putative hydrolase
RA BS168_reference_genome 3,141,615 A→G 6.4% A270A (GCA→GCG ythB → putative cytochrome bd menaquinol oxidase subunit II
RA BS168_reference_genome 3,197,776 A→T 6.5% K291* (AAA→TAA)  iolU → scyllo‑inositol dehydrogenase (NADP+‑dependent); biofilm formation
RA BS168_reference_genome 3,290,481 C→T 5.5% M335I (ATG→ATA dhbC ← isochorismate synthase (siderophore‑specific)
RA BS168_reference_genome 3,392,796 G→T 32.7% G200* (GGA→TGA)  gerAB → component of the germination receptor GerA; putative transporter
RA BS168_reference_genome 3,496,398 C→T 9.9% M27I (ATG→ATA yvfU ← two‑component response regulator [YvfT]
RA BS168_reference_genome 3,695,698 C→T 6.7% D112D (GAC→GAT ribZB → 5‑amino‑6‑ribitylamino‑2,4(1H, 3H)‑pyrimidinedione 5'‑phosphate phosphatase (promiscuous activity)
RA BS168_reference_genome 3,701,226 T→A 11.4% intergenic (‑599/‑185) capB ← / → rbsR capsular polyglutamate synthetase (ATP‑dependent amide ligase)/ribose operon repressor (LacI family, D‑ribose)
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,920,358 G→A 7.7% P93S (CCT→TCT)  galT ← galactose‑1‑phosphate uridyltransferase
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,010,739 A→T 6.6% G112G (GGA→GGT citH → secondary transporter of divalent metal ions/citrate complexes
RA BS168_reference_genome 4,048,670 T→G 20.0% V378G (GTC→GGC)  hutM → histidine degradation metabolites permease
RA BS168_reference_genome 4,073,940 G→T 6.5% S5* (TCG→TAG)  iolJ ← 2‑deoxy‑5‑keto‑D‑gluconic acid 6‑phosphate aldolase
RA BS168_reference_genome 4,076,868 T→G 6.7% K431Q (AAA→CAA)  iolF ← inositol transport protein
RA BS168_reference_genome 4,077,654 T→C 6.9% N169D (AAT→GAT)  iolF ← inositol transport protein
RA BS168_reference_genome 4,087,669 C→T 6.3% H72Y (CAC→TAC)  yxcD → conserved protein of unknown function
RA BS168_reference_genome 4,087,723 G→C 14.2% G90R (GGC→CGC)  yxcD → conserved protein of unknown function
RA BS168_reference_genome 4,102,595 G→A 5.3% T356M (ACG→ATG)  yxaL ← membrane associated protein kinase with beta‑propeller domain
RA BS168_reference_genome 4,137,538 G→C 20.6% H25D (CAT→GAT)  yycQ ← conserved membrane protein
RA BS168_reference_genome 4,151,742 G→A 7.2% S41F (TCT→TTT)  walH ← regulator of WalRK(YycFG)
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,165,884 G→A 7.0% G235G (GGC→GGT yybS ← putative integral inner membrane protein
RA BS168_reference_genome 4,191,258 C→T 8.5% H21Y (CAC→TAC)  yyaL → conserved hypothetical protein
RA Exported 78 G→T 5.3% intergenic (–/‑2900)  / → eryR –//vntifkey=4
RA Exported 121 G→C 5.9% intergenic (–/‑2857)  / → eryR –//vntifkey=4
RA Exported 131 G→C 7.0% intergenic (–/‑2847)  / → eryR –//vntifkey=4
RA Exported 137 A→C 6.0% intergenic (–/‑2841)  / → eryR –//vntifkey=4
RA Exported 2,206 C→T 5.1% intergenic (–/‑772)  / → eryR –//vntifkey=4
RA Exported 2,224 T→C 15.1% intergenic (–/‑754)  / → eryR –//vntifkey=4
RA Exported 2,236 A→G 7.5% intergenic (–/‑742)  / → eryR –//vntifkey=4
RA Exported 2,239 A→G 7.8% intergenic (–/‑739)  / → eryR –//vntifkey=4
JC Exported 2,655 Δ3 bp 100% intergenic (–/‑321)  / → eryR –//vntifkey=4
RA Exported 3,729 T→G 6.2% intergenic (+14/‑142) eryR → / → ampR /vntifkey=4//vntifkey=4
RA Exported 3,735 A→T 9.0% intergenic (+20/‑136) eryR → / → ampR /vntifkey=4//vntifkey=4
RA Exported 3,741 G→T 9.0% intergenic (+26/‑130) 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.9% K252T (AAG→ACG)  ampR → /vntifkey=4
RA Exported 4,649 A→C 8.4% 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.7% 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.5% V40G (GTT→GGT)  lacZalpha → /vntifkey=4
RA Exported 6,012 A→C 7.4% Q42P (CAA→CCA)  lacZalpha → /vntifkey=4
RA Exported 6,031 A→C 6.7% A48A (GCA→GCC lacZalpha → /vntifkey=4
RA Exported 6,033 A→C 5.2% H49P (CAT→CCT)  lacZalpha → /vntifkey=4
RA Exported 6,047 A→C 5.2% S54R (AGC→CGC)  lacZalpha → /vntifkey=4

Unassigned missing coverage evidence
   seq id start end size ←reads reads→ gene description
* * ÷ BS168_reference_genome 160964–161238 165598–164903 3666–4635 12 [10] [9] 11 [rrnI‑16S]–[rrnI‑5S] [rrnI‑16S],rrnI‑23S,[rrnI‑5S]
* * ÷ BS168_reference_genome 166441–166703 171220–170393 3691–4780 11 [9] [10] 13 rrnH‑16S–[rrnH‑5S] rrnH‑16S,rrnH‑23S,[rrnH‑5S]
* * ÷ BS168_reference_genome 171470–171679 176232–172978 1300–4763 13 [5] [10] 12 rrnG‑16S–[rrnG‑5S] rrnG‑16S,rrnG‑23S,[rrnG‑5S]
* * ÷ BS168_reference_genome 1073639 1074757 1119 11 [10] [10] 11 [scoC]–[trpP] [scoC],yhaH,yhzF,[trpP]

Unassigned new junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? BS168_reference_genome = 10734584 (0.080)34 (0.680) 22/230 NT 87.2% coding (260/612 nt) scoC transcriptional regulator of extracellular protease production, sporulation and bacilysin production (MarR family)
?BS168_reference_genome 1074991 = 6 (0.120)coding (174/519 nt) trpP tryptophan transporter or sensor