breseq  version 0.33.1  revision 8505477f25b3
mutation predictions | marginal predictions | summary statistics | genome diff | command line log

Predicted mutations
evidence position mutation freq annotation gene description
RA 32,807 C→A 6.0% T664N (ACC→AAC)  carB → carbamoyl‑phosphate synthase large subunit
RA 114,203 C→T 6.0% R254C (CGC→TGC)  guaC → GMP reductase
RA 129,076 G→A 5.0% A389T (GCT→ACT)  lpd → dihydrolipoyl dehydrogenase; E3 component of pyruvate and 2‑oxoglutarate dehydrogenases complexes; glycine cleavage system L protein; dihydrolipoamide dehydrogenase
JC JC 168,640 IS2 (+) +5 bp 100% coding (1157‑1161/2244 nt) fhuA → ferrichrome outer membrane transporter
RA 191,997 G→A 5.6% G48D (GGT→GAT)  pyrH → uridylate kinase
MC JC 257,908 Δ776 bp 100% [crl] [crl]
RA 342,393 C→T 8.2% D423D (GAC→GAT yahJ → putative metallo‑dependent hydrolase domain deaminase
RA 386,011 A→T 5.8% E260D (GAA→GAT tauA → taurine ABC transporter periplasmic binding protein
RA 405,463 G→A 5.3% A61V (GCG→GTG)  proC ← pyrroline‑5‑carboxylate reductase, NAD(P)‑binding
RA 414,144 T→G 6.2% N537T (AAC→ACC)  sbcC ← exonuclease, dsDNA, ATP‑dependent
RA 422,062 G→A 5.4% V359V (GTG→GTA proY → proline‑specific permease
RA 467,238 G→T 5.1% T25T (ACC→ACA ybaE ← putative ABC transporter periplasmic binding protein
RA 470,045 A→T 5.2% D392V (GAT→GTT)  mdlA → putative multidrug ABC transporter ATPase
RA 470,515 G→A 5.0% A549T (GCC→ACC)  mdlA → putative multidrug ABC transporter ATPase
RA 473,997 C→T 8.5% G344G (GGC→GGT amtB → ammonium transporter
RA 555,233 A→T 6.3% E208V (GAA→GTA)  cysS → cysteinyl‑tRNA synthetase
RA 596,222 G→T 7.1% S208I (AGC→ATC)  cusC → copper/silver efflux system, outer membrane component
RA 627,738 C→T 6.9% H15H (CAC→CAT entB → isochorismatase
RA 698,393 T→C 5.8% E262G (GAA→GGA)  asnB ← asparagine synthetase B
RA 716,326 C→T 8.0% pseudogene (272/651 nt) ybfG ← pseudogene
RA 787,359 C→T 6.1% L79L (CTG→CTA gpmA ← phosphoglyceromutase 1
RA 829,686 G→A 10.2% A96V (GCG→GTG)  ybhG ← putative membrane fusion protein (MFP) component of efflux pump, membrane anchor
RA 909,353 G→A 100% A566V (GCG→GTG)  poxB ← pyruvate dehydrogenase, thiamine triphosphate‑binding, FAD‑binding
JC JC 911,086 IS2 (–) +5 bp 100% intergenic (‑37/+92) poxB ← / ← hcr pyruvate dehydrogenase, thiamine triphosphate‑binding, FAD‑binding/HCP oxidoreductase, NADH‑dependent
RA 926,398 T→C 5.1% T16A (ACG→GCG)  infA ← translation initiation factor IF‑1
RA 951,524 T→C 6.3% K677K (AAA→AAG pflB ← formate C‑acetyltransferase 1, anaerobic; pyruvate formate‑lyase 1
RA 981,097 C→A 7.1% C17* (TGC→TGA ldtD → murein L,D‑transpeptidase
RA 999,964 T→C 5.2% V250A (GTA→GCA)  elfC → putative outer membrane fimbrial subunit export usher protein
RA 1,014,199 C→T 5.8% T314I (ACC→ATC)  pqiB → paraquat‑inducible, SoxRS‑regulated MCE domain protein
RA 1,046,089 G→A 7.3% S619S (AGC→AGT gfcD ← putative O‑antigen capsule production periplasmic protein
RA 1,126,716 A→T 5.4% D187V (GAT→GTT)  yceH → UPF0502 family protein
RA 1,130,038 C→T 5.1% A31T (GCA→ACA)  flgM ← anti‑sigma factor for FliA (sigma 28)
RA 1,158,773 T→C 5.9% F302S (TTC→TCC)  ptsG → fused glucose‑specific PTS enzymes: IIB component/IIC component
RA 1,158,881 C→T 100% S338L (TCG→TTG)  ptsG → fused glucose‑specific PTS enzymes: IIB component/IIC component
RA 1,161,067 C→T 6.3% R162H (CGT→CAT)  fhuE ← ferric‑rhodotorulic acid outer membrane transporter
RA 1,172,148 G→A 7.4% A606V (GCG→GTG)  mfd ← transcription‑repair coupling factor
RA 1,211,424 T→C 9.5% pseudogene (12/165 nt) icdC → pseudogene, isocitrate dehydrogenase C‑terminal gene fragment; idcC' is a 54 codon 3' gene fragment created during e14 prophage insertion
RA 1,258,185 C→T 6.0% G110D (GGC→GAC)  pth ← peptidyl‑tRNA hydrolase
RA 1,286,067 G→A 8.7% P73P (CCG→CCA narI → nitrate reductase 1, gamma (cytochrome b(NR)) subunit
JC JC 1,293,032 IS1 (–) +8 bp 100% intergenic (‑110/‑488) hns ← / → tdk global DNA‑binding transcriptional dual regulator H‑NS/thymidine kinase/deoxyuridine kinase
MC JC 1,299,499 Δ1,199 bp 100% intergenic (+254/‑485) ychE → / → oppA UPF0056 family inner membrane protein/oligopeptide ABC transporter periplasmic binding protein
RA 1,301,823 C→T 5.2% N214N (AAC→AAT oppA → oligopeptide ABC transporter periplasmic binding protein
RA 1,303,990 G→A 9.9% V53I (GTA→ATA)  oppC → oligopeptide ABC transporter permease
RA 1,426,842 A→G 5.1% pseudogene (389/933 nt) tmpR → Rac prophage; pseudogene, tail protein family;Phage or Prophage Related; putative alpha helix protein
RA 1,437,993 G→A 5.2% A931V (GCT→GTT)  pfo ← pyruvate‑flavodoxin oxidoreductase
RA 1,569,308 A→G 6.3% L394L (CTT→CTC gadC ← glutamate:gamma‑aminobutyric acid antiporter
RA 1,575,620 G→A 5.9% A667V (GCA→GTA)  yddB ← putative TonB‑dependent outer membrane receptor
RA 1,760,433 C→T 5.2% A31T (GCC→ACC)  sufS ← cysteine desulfurase, stimulated by SufE; selenocysteine lyase, PLP‑dependent
RA 1,830,829 T→C 8.3% E385G (GAG→GGG)  astC ← succinylornithine transaminase, PLP‑dependent
RA 1,841,764 C→T 8.6% T92M (ACG→ATG)  nudG → CTP pyrophosphohydrolase; also hydrolyzes 2‑hydroxy‑dATP, 8‑hydroxy‑dGTP, 5‑hydroxy‑CTP, dCTP and 5‑methyl‑dCTP
RA 1,872,748 G→A 5.8% V236V (GTG→GTA yeaJ → putative diguanylate cyclase
RA 1,890,091 G→C 6.9% H148D (CAT→GAT)  yeaY ← Slp family lipoprotein, RpoE‑regulated
RA 1,899,169 C→T 6.2% A248V (GCG→GTG)  yoaD → putative membrane‑anchored cyclic‑di‑GMP phosphodiesterase, regulator of cellulose production
RA 1,904,339 C→T 6.0% L149L (CTG→TTG)  manZ → mannose‑specific enzyme IID component of PTS
RA 1,916,484 T→A 7.4% M76K (ATG→AAG)  yebS → inner membrane PqiA domain protein
RA 1,927,030 G→T 5.1% Q604K (CAG→AAG)  ptrB ← protease II
RA 1,947,359 G→A 8.7% G6G (GGC→GGT ruvC ← component of RuvABC resolvasome, endonuclease
RA 1,973,263 T→C 8.7% D26G (GAT→GGT)  cheW ← purine‑binding chemotaxis protein
MC JC 1,978,503 Δ776 bp 100% insB1insA insB1, insA
JC JC 1,979,486 IS5 (+) +4 bp 100% intergenic (‑271/‑264) insA ← / → uspC IS1 repressor TnpA/universal stress protein
RA 2,021,984 Δ1 bp 6.0% coding (116/738 nt) fliP → flagellar biosynthesis protein
RA 2,075,061 T→C 5.9% S95P (TCT→CCT)  yeeR → CP4‑44 prophage; putative membrane protein
RA 2,085,179 T→C 6.8% Y116C (TAT→TGT)  yeeE ← UPF0394 family inner membrane protein
RA 2,086,469 C→T 6.0% A198A (GCG→GCA plaP ← putrescine importer, low affinity
RA 2,120,760 C→T 6.0% P265P (CCG→CCA wcaJ ← colanic biosynthesis UDP‑glucose lipid carrier transferase
RA 2,132,787 A→C 100% I204S (ATC→AGC)  wcaA ← putative glycosyl transferase
RA 2,173,361 Δ2 bp 100% pseudogene (1‑2/442 nt) gatC ← pseudogene, galactitol‑specific enzyme IIC component of PTS;transport; Transport of small molecules: Carbohydrates, organic acids, alcohols; PTS system galactitol‑specific enzyme IIC
JC JC 2,175,712 IS1 (–) +9 bp 100% coding (602‑610/1263 nt) gatZ ← D‑tagatose 1,6‑bisphosphate aldolase 2, subunit
RA 2,330,385 C→T 8.1% pseudogene (4018/4605 nt) yfaS ← pseudogene, bacterial alpha2‑macroglobulin YfaS variant family; putative membrane protein
RA 2,353,768 C→T 10.8% G374G (GGC→GGT glpA → anaerobic sn‑glycerol‑3‑phosphate dehydrogenase, large FAD/NAD(P)‑binding subunit
RA 2,391,376 C→T 5.2% G44R (GGG→AGG)  nuoN ← NADH:ubiquinone oxidoreductase, membrane subunit N
RA 2,413,461 G→A 100% intergenic (‑329/‑9) yfbV ← / → ackA UPF0208 family inner membrane protein/acetate kinase A and propionate kinase 2
RA 2,469,009 G→T 100% V266F (GTT→TTT)  gtrB → CPS‑53 (KpLE1) prophage; bactoprenol glucosyl transferase
RA 2,538,798 T→C 8.0% G262G (GGA→GGG cysM ← cysteine synthase B (O‑acetylserine sulfhydrolase B)
RA 2,551,174 G→C 7.1% L15V (CTG→GTG)  yfeY ← RpoE‑regulated lipoprotein
RA 2,583,110 G→A 7.3% G123G (GGC→GGT nudK ← GDP‑mannose pyrophosphatase
RA 2,584,549 A→G 5.6% V326A (GTG→GCG)  aegA ← putative oxidoreductase, FeS binding subunit/NAD/FAD‑binding subunit
RA 2,605,247 G→A 5.8% G146D (GGC→GAC)  hyfD → hydrogenase 4, membrane subunit
RA 2,614,970 C→T 5.9% A304A (GCG→GCA yfgO ← putative UPF0118 family inner membrane permease
RA 2,659,608 C→T 5.9% P93P (CCG→CCA iscA ← FeS cluster assembly protein
JC 2,722,123 Δ7 bp 100% coding (2171‑2177/2661 nt) pka → protein lysine acetyltransferase
RA 2,829,372 T→C 6.9% L109P (CTG→CCG)  srlR → sorbitol‑inducible srl operon transcriptional repressor
RA 2,834,784 C→T 5.6% S291S (AGC→AGT norW → NADH:flavorubredoxin oxidoreductase
RA 2,933,574 G→A 8.9% L39L (CTG→TTG)  fucA ← L‑fuculose‑1‑phosphate aldolase
RA 2,934,979 G→A 5.3% A249T (GCA→ACA)  fucP → L‑fucose transporter
RA 3,035,437 T→C 82.6% K282R (AAG→AGG)  prfB ← peptide chain release factor RF‑2
RA 3,081,526 G→A 6.3% A37V (GCC→GTC)  tktA ← transketolase 1, thiamine triphosphate‑binding
RA 3,134,891 A→T 8.2% I494K (ATA→AAA)  pitB ← phosphate transporter
RA 3,211,722 T→C 6.1% F206L (TTC→CTC)  dnaG → DNA primase
RA 3,259,683 G→A 5.7% intergenic (‑34/+38) tdcG ← / ← tdcF L‑serine dehydratase 3, anaerobic/putative reactive intermediate deaminase
RA 3,302,151 G→A 6.9% G223S (GGC→AGC)  yhbU → U32 peptidase family protein
RA 3,479,755 C→T 10.2% G285R (GGG→AGG)  kefB ← potassium:proton antiporter
RA 3,484,760 A→G 6.0% K91E (AAA→GAA)  prkB → putative phosphoribulokinase
RA 3,509,494 A→G 5.1% V107A (GTG→GCG)  yhfW ← phosphopentomutase‑related metalloenzyme superfamily protein
RA 3,560,455:1 +G 100% intergenic (‑2/+1) glpR ← / ← glpR pseudogene, DNA‑binding transcriptional repressor;regulator; Energy metabolism, carbon: Anaerobic respiration; repressor of the glp operon/pseudogene, DNA‑binding transcriptional repressor;regulator; Energy metabolism, carbon: Anaerobic respiration; repressor of the glp operon
RA 3,562,212 C→T 8.8% A67V (GCG→GTG)  glpD → sn‑glycerol‑3‑phosphate dehydrogenase, aerobic, FAD/NAD(P)‑binding
RA 3,594,232 T→C 5.9% N417D (AAC→GAC)  livM ← branched‑chain amino acid ABC transporter permease
RA 3,629,291 C→T 5.1% R83H (CGC→CAC)  rbbA ← ribosome‑associated ATPase: ATP‑binding protein/ATP‑binding membrane protein
RA 3,699,528 G→A 5.1% V439M (GTG→ATG)  bcsG → DUF3260 family cellulose production inner membrane protein
JC JC 3,771,709 IS3 (–) +3 bp 100% coding (34‑36/363 nt) yibI ← DUF3302 family inner membrane protein
RA 3,775,311 G→T 5.1% Q296H (CAG→CAT mtlD → mannitol‑1‑phosphate dehydrogenase, NAD‑dependent
RA 3,798,652 Δ1 bp 100% coding (661/1074 nt) waaU ← lipopolysaccharide core biosynthesis
RA 3,808,813 C→T 6.2% R92C (CGC→TGC)  waaA → 3‑deoxy‑D‑manno‑octulosonic‑acid transferase (KDO transferase)
MC JC 3,815,859 Δ82 bp 100% [rph][rph] [rph], [rph]
RA 3,816,244 T→C 7.8% pseudogene (306/669 nt) rph ← ribonuclease PH (defective);enzyme; Degradation of RNA; RNase PH
RA 3,834,765 C→A 6.1% G389* (GGA→TGA)  yicJ ← putative transporter
RA 3,850,593 C→T 6.2% G45E (GGG→GAG)  uhpA ← response regulator in two‑component regulatory system wtih UhpB
RA 3,926,829 C→T 5.3% V59I (GTC→ATC)  asnC ← transcriptional activator of asnA; autorepressor
RA 3,936,634 A→T 7.9% K119N (AAA→AAT rbsB → D‑ribose ABC transporter periplasmic binding protein; ribose chemotaxis receptor
RA 3,950,784 A→G 5.2% pseudogene (225/984 nt) ilvG → pseudogene, acetolactate synthase 2 large subunit, valine‑insensitive; acetolactate synthase II, large subunit, cryptic, interrupted
RA 3,956,462 T→A 5.1% C378S (TGC→AGC)  ilvA → l‑threonine dehydratase, biosynthetic; also known as threonine deaminase
JC 4,001,645 Δ5 bp 100% coding (220‑224/951 nt) corA → magnesium/nickel/cobalt transporter
RA 4,011,492 T→A 6.0% T139T (ACT→ACA yigM → putative inner membrane EamA‑like transporter
RA 4,021,557 T→C 6.5% P444P (CCT→CCC ubiB → regulator of octaprenylphenol hydroxylation, ubiquinone synthesis; regulator of 2'‑N‑acetyltransferase; putative ABC1 family protein kinase
RA 4,030,251 C→T 5.4% A241T (GCA→ACA)  fadB ← fused 3‑hydroxybutyryl‑CoA epimerase/delta(3)‑cis‑delta(2)‑trans‑enoyl‑CoA isomerase/enoyl‑CoA hydratase/3‑hydroxyacyl‑CoA dehydrogenase
RA 4,130,849 C→A 6.5% Q339K (CAA→AAA)  metL → Bifunctional aspartokinase/homoserine dehydrogenase 2
RA 4,184,543 C→A 100% P1100Q (CCG→CAG)  rpoB → RNA polymerase, beta subunit
RA 4,236,466 G→T 5.6% A69S (GCT→TCT)  yjbF → extracellular polysaccharide production lipoprotein
RA 4,296,060 C→T 20.4% intergenic (+266/+376) gltP → / ← yjcO glutamate/aspartate:proton symporter/Sel1 family TPR‑like repeat protein
RA 4,296,380:1 +CG 100% intergenic (+586/+56) gltP → / ← yjcO glutamate/aspartate:proton symporter/Sel1 family TPR‑like repeat protein
RA 4,320,735 G→A 5.1% G171G (GGC→GGT phnH ← ribophosphonate triphosphate synthase subunit
RA 4,325,804 C→T 7.3% intergenic (‑63/+595) yjdN ← / ← yjdM metalloprotein superfamily protein/zinc‑ribbon family protein
RA 4,422,372 G→A 5.8% A63T (GCC→ACC)  ulaD → 3‑keto‑L‑gulonate 6‑phosphate decarboxylase
RA 4,439,556 G→A 6.0% intergenic (‑294/‑31) ytfJ ← / → ytfK putative transcriptional regulator/DUF1107 family protein
RA 4,473,613 G→A 6.8% A147A (GCC→GCT bdcA ← c‑di‑GMP‑binding biofilm dispersal mediator protein
RA 4,509,643 A→G 6.2% pseudogene (165/315 nt) insO → pseudogene, IS911 transposase B;IS, phage, Tn; Transposon‑related functions; extrachromosomal; transposon related
RA 4,532,778 T→A 5.4% M136L (ATG→TTG)  yjhP ← putative methyltransferase
RA 4,548,579 A→G 7.0% K99E (AAA→GAA)  fimG → minor component of type 1 fimbriae
RA 4,626,708 A→G 6.3% E266G (GAG→GGG)  radA → DNA repair protein

Unassigned missing coverage evidence
   seq id start end size ←reads reads→ gene description
* * ÷ NC_000913 123061 127538 4478 30 [8] [18] 20 [aceE]–[aceF] [aceE],[aceF]
* * ÷ NC_000913 2066303–2066464 2066597–2066466 3–295 19 [17] [17] 20 [insH1] [insH1]
* * ÷ NC_000913 3423748–3424233 3424549–3424238 6–802 19 [18] [18] 20 [rrfD]–[rrlD] [rrfD],[rrlD]

Unassigned new junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 1207790 =26 (0.320)69 (0.970) 47/216 NT 73.2% coding (290/630 nt) stfP e14 prophage; uncharacterized protein
?NC_000913 1209619 = 28 (0.400)pseudogene (1/501 nt) stfE pseudogene, e14 prophage; side tail fiber protein fragment family;Phage or Prophage Related
* ? NC_000913 = 120780525 (0.310)57 (0.810) 42/216 NT 69.6% coding (305/630 nt) stfP e14 prophage; uncharacterized protein
?NC_000913 = 1209602 28 (0.400)pseudogene (18/501 nt) stfE pseudogene, e14 prophage; side tail fiber protein fragment family;Phage or Prophage Related
* ? NC_000913 1411899 =125 (1.540)19 (0.300) 18/196 NT 13.6% coding (50/936 nt) ttcA tRNA s(2)C32 thioltransferase, iron sulfur cluster protein
?NC_000913 = 1434958 142 (2.210)pseudogene (51/51 nt) ttcC pseudogene, prophage Rac integration site ttcA duplication;Phage or Prophage Related
* ? NC_000913 4542682 =14 (0.170)195 (2.590) 83/230 NT 94.2% intergenic (+49/‑433) fimE/fimA tyrosine recombinase/inversion of on/off regulator of fimA/major type 1 subunit fimbrin (pilin)
?NC_000913 4542996 = 11 (0.150)intergenic (+363/‑119) fimE/fimA tyrosine recombinase/inversion of on/off regulator of fimA/major type 1 subunit fimbrin (pilin)
* ? NC_000913 = 454269013 (0.160)202 (2.680) 97/230 NT 94.6% intergenic (+57/‑425) fimE/fimA tyrosine recombinase/inversion of on/off regulator of fimA/major type 1 subunit fimbrin (pilin)
?NC_000913 = 4542986 11 (0.150)intergenic (+353/‑129) fimE/fimA tyrosine recombinase/inversion of on/off regulator of fimA/major type 1 subunit fimbrin (pilin)