breseq  version 0.32.0  revision 6ff6de7d1b87
mutation predictions | marginal predictions | summary statistics | genome diff | command line log

Predicted mutations
evidence position mutation freq annotation gene description
RA 33,407 T→G 15.8% V864G (GTG→GGG)  carB → carbamoyl‑phosphate synthase large subunit
RA 55,437 A→C 11.5% V558G (GTG→GGG)  lptD ← LPS assembly OM complex LptDE, beta‑barrel component
RA 217,141 T→A 2.6% E545D (GAA→GAT proS ← prolyl‑tRNA synthetase
RA 217,144 G→A 2.5% I544I (ATC→ATT proS ← prolyl‑tRNA synthetase
RA 246,879 G→A 3.1% L56L (CTG→CTA yafL → putative lipoprotein and C40 family peptidase
MC JC 257,908 Δ776 bp 100% [crl] [crl]
RA 281,089 A→C 22.2% Y299D (TAT→GAT)  yagA ← CP4‑6 prophage, putative DNA‑binding transcriptional regulator
RA 326,812 T→C 15.5% S146G (AGC→GGC)  betA ← choline dehydrogenase, a flavoprotein
RA 328,024 A→C 14.4% V237G (GTG→GGG)  betB ← betaine aldehyde dehydrogenase, NAD‑dependent
RA 337,463 A→G 20.4% E229G (GAG→GGG)  yahF → putative NAD(P)‑binding succinyl‑CoA synthase
RA 349,642 A→G 3.8% intergenic (+70/‑370) prpB → / → prpC 2‑methylisocitrate lyase/2‑methylcitrate synthase
RA 349,673 G→C 2.8% intergenic (+101/‑339) prpB → / → prpC 2‑methylisocitrate lyase/2‑methylcitrate synthase
RA 349,717 T→C 4.0% intergenic (+145/‑295) prpB → / → prpC 2‑methylisocitrate lyase/2‑methylcitrate synthase
RA 366,509 C→G 13.0% T334T (ACG→ACC lacI ← lactose‑inducible lac operon transcriptional repressor
RA 373,637 T→G 19.9% G239G (GGT→GGG mhpF → acetaldehyde‑CoA dehydrogenase II, NAD‑binding
RA 433,649 T→G 19.9% G65G (GGT→GGG ribD → fused diaminohydroxyphosphoribosylaminopyrimidine deaminase and 5‑amino‑6‑(5‑phosphoribosylamino) uracil reductase
RA 557,156 G→T 3.7% L196M (CTG→ATG)  folD ← bifunctional 5,10‑methylene‑tetrahydrofolate dehydrogenase/ 5,10‑methylene‑tetrahydrofolate cyclohydrolase
RA 599,581 A→G 19.6% I290V (ATC→GTC)  cusA → copper/silver efflux system, membrane component
RA 766,694 T→A 16.2% N237K (AAT→AAA mngA → fused 2‑O‑a‑mannosyl‑D‑glycerate specific PTS enzymes: IIA component/IIB component/IIC component
RA 811,721 T→G 30.1% V67G (GTG→GGG)  bioC → malonyl‑ACP O‑methyltransferase, SAM‑dependent
RA 825,233 C→A 12.9% M292I (ATG→ATT ybhR ← putative ABC transporter permease
RA 921,243 G→T 6.3% M299I (ATG→ATT macB → macrolide ABC transporter peremase/ATPase
RA 954,360 T→A 11.8% Y36F (TAT→TTT)  focA ← formate channel
RA 954,374 C→G 9.8% P31P (CCG→CCC focA ← formate channel
RA 989,945 A→G 17.2% S138P (TCA→CCA)  pncB ← nicotinate phosphoribosyltransferase
RA 1,126,895:1 +G 5.5% coding (90/924 nt) yceM → putative oxidoreductase
JC 1,195,194 Δ3 bp 2.6% coding (72‑74/1251 nt) icd → isocitrate dehydrogenase, e14 prophage attachment site, tellurite reductase
RA 1,196,131 G→A 31.0% A337T (GCC→ACC)  icd → isocitrate dehydrogenase, e14 prophage attachment site, tellurite reductase
RA 1,196,305 C→T 49.6% R395C (CGT→TGT)  icd → isocitrate dehydrogenase, e14 prophage attachment site, tellurite reductase
RA 1,238,523:1 +G 3.3% coding (953/1299 nt) dadA → D‑amino acid dehydrogenase
RA 1,287,149 G→A 5.9% noncoding (88/171 nt)
R10* (CGA→TGA) 
rttR ←
tpr ←
rtT sRNA, processed from tyrT transcript
protamine‑like protein
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,466,201 T→C 5.3% pseudogene (810/2513 nt) ydbA → pseudogene, autotransporter homolog, interrupted by IS2 and IS30
RA 1,466,432 A→G 9.1% pseudogene (1041/2513 nt) ydbA → pseudogene, autotransporter homolog, interrupted by IS2 and IS30
RA 1,466,438 T→G 7.8% pseudogene (1047/2513 nt) ydbA → pseudogene, autotransporter homolog, interrupted by IS2 and IS30
RA 1,466,579 G→C 11.1% pseudogene (1188/2513 nt) ydbA → pseudogene, autotransporter homolog, interrupted by IS2 and IS30
RA 1,466,729 G→A 4.5% pseudogene (1338/2513 nt) ydbA → pseudogene, autotransporter homolog, interrupted by IS2 and IS30
RA 1,708,758 G→C 8.8% E614Q (GAA→CAA)  rsxC → SoxR iron‑sulfur cluster reduction factor component, putative membrane‑associated NADH oxidoreductase of electron transport complex
RA 1,708,814 A→G 11.0% A632A (GCA→GCG rsxC → SoxR iron‑sulfur cluster reduction factor component, putative membrane‑associated NADH oxidoreductase of electron transport complex
RA 1,708,844 C→G 12.4% A642A (GCC→GCG rsxC → SoxR iron‑sulfur cluster reduction factor component, putative membrane‑associated NADH oxidoreductase of electron transport complex
RA 2,040,143 A→G 21.8% V222V (GTA→GTG yedY → membrane‑anchored, periplasmic TMAO, DMSO reductase
RA 2,132,787 A→C 100% I204S (ATC→AGC)  wcaA ← putative glycosyl transferase
RA 2,151,689 G→T 28.8% intergenic (+43/+22) yegJ → / ← yegK uncharacterized protein/ser/thr phosphatase‑related protein
RA 2,165,068 G→C 2.9% intergenic (+70/‑121) baeR → / → yegP response regulator in two‑component regulatory system with BaeS/UPF0339 family protein
RA 2,165,215 G→T 4.1% K9N (AAG→AAT yegP → UPF0339 family protein
RA 2,167,707 Δ1 bp 3.8% pseudogene (108/213 nt) yegZ ← pseudogene, gpD phage P2‑like protein D
RA 2,228,227 G→C 17.8% pseudogene (381/1287 nt) mdtQ ← pseudogene, putative channel/filament proteins
RA 2,367,397 T→C 22.6% I110T (ATC→ACC)  arnC → undecaprenyl phosphate‑L‑Ara4FN transferase
RA 2,480,378 A→G 16.1% G393G (GGA→GGG dsdA → D‑serine dehydratase
RA 2,509,878 A→C 19.0% D84A (GAC→GCC)  yfeO → putative ion channel protein
RA 2,643,098 Δ1 bp 2.7% intergenic (‑256/+29) rodZ ← / ← rlmN MreB assembly cytoskeletal protein/dual specificity 23S rRNA m(2)A2503, tRNA m(2)A37 methyltransferase, SAM‑dependent
RA 2,643,098:1 +T 2.9% intergenic (‑256/+29) rodZ ← / ← rlmN MreB assembly cytoskeletal protein/dual specificity 23S rRNA m(2)A2503, tRNA m(2)A37 methyltransferase, SAM‑dependent
RA 2,643,101 A→T 3.2% intergenic (‑259/+26) rodZ ← / ← rlmN MreB assembly cytoskeletal protein/dual specificity 23S rRNA m(2)A2503, tRNA m(2)A37 methyltransferase, SAM‑dependent
RA 2,643,103 A→T 3.2% intergenic (‑261/+24) rodZ ← / ← rlmN MreB assembly cytoskeletal protein/dual specificity 23S rRNA m(2)A2503, tRNA m(2)A37 methyltransferase, SAM‑dependent
RA 2,643,105 A→T 3.2% intergenic (‑263/+22) rodZ ← / ← rlmN MreB assembly cytoskeletal protein/dual specificity 23S rRNA m(2)A2503, tRNA m(2)A37 methyltransferase, SAM‑dependent
RA 3,053,583 T→G 8.6% N419T (AAC→ACC)  pepP ← proline aminopeptidase P II
RA 3,206,434 C→A 13.3% intergenic (‑180/‑27) ttdR ← / → ttdA transcriptional activator of ttdABT/L‑tartrate dehydratase, alpha subunit
RA 3,206,436 G→A 18.2% intergenic (‑182/‑25) ttdR ← / → ttdA transcriptional activator of ttdABT/L‑tartrate dehydratase, alpha subunit
RA 3,206,454 A→G 8.8% intergenic (‑200/‑7) ttdR ← / → ttdA transcriptional activator of ttdABT/L‑tartrate dehydratase, alpha subunit
RA 3,206,458 G→T 13.8% intergenic (‑204/‑3) ttdR ← / → ttdA transcriptional activator of ttdABT/L‑tartrate dehydratase, alpha subunit
RA 3,206,465 T→G 8.3% M2R (ATG→AGG)  ttdA → L‑tartrate dehydratase, alpha subunit
RA 3,206,469 C→G 16.6% S3R (AGC→AGG ttdA → L‑tartrate dehydratase, alpha subunit
RA 3,206,474 G→A 9.3% S5N (AGT→AAT) ‡ ttdA → L‑tartrate dehydratase, alpha subunit
RA 3,206,475 T→A 18.6% S5R (AGT→AGA) ‡ ttdA → L‑tartrate dehydratase, alpha subunit
RA 3,206,481 G→A 23.7% K7K (AAG→AAA ttdA → L‑tartrate dehydratase, alpha subunit
RA 3,206,482 C→A 13.2% Q8K (CAA→AAA)  ttdA → L‑tartrate dehydratase, alpha subunit
RA 3,206,492 T→G 12.9% V11G (GTG→GGG) ‡ ttdA → L‑tartrate dehydratase, alpha subunit
RA 3,206,493 G→A 13.5% V11V (GTG→GTA) ‡ ttdA → L‑tartrate dehydratase, alpha subunit
RA 3,206,504 C→A 13.7% T15K (ACA→AAA)  ttdA → L‑tartrate dehydratase, alpha subunit
RA 3,234,702 G→A 2.8% intergenic (‑249/+35) higB ← / ← rlmG mRNA interferase toxin of the HigB‑HigA toxin‑antitoxin system/23S rRNA m(2)G1835 methyltransferase, SAM‑dependent
RA 3,364,342 C→A 8.0% S513* (TCA→TAA)  yhcD → putative outer membrane fimbrial subunit usher protein
RA 3,520,794 G→T 11.0% T75K (ACA→AAA)  hofP ← DNA catabolic protein
RA 3,528,943 G→A 2.7% V93I (GTA→ATA)  yrfG → GMP/IMP nucleotidase
RA 3,560,493 G→A 8.5% pseudogene (111/149 nt) glpR ← pseudogene, DNA‑binding transcriptional repressor,regulator, Energy metabolism, carbon: Anaerobic respiration, repressor of the glp operon
RA 3,714,370 G→T 8.1% A675D (GCC→GAC)  bisC ← biotin sulfoxide reductase
RA 3,760,368 T→C 16.6% D294G (GAC→GGC)  selA ← selenocysteine synthase
RA 3,782,901 T→G 17.3% T254P (ACC→CCC)  gpsA ← glycerol‑3‑phosphate dehydrogenase (NAD+)
RA 3,782,912 A→C 14.9% V250G (GTG→GGG)  gpsA ← glycerol‑3‑phosphate dehydrogenase (NAD+)
RA 3,813,071 A→C 21.1% T114T (ACA→ACC dfp → fused 4'‑phosphopantothenoylcysteine decarboxylase/phosphopantothenoylcysteine synthetase, FMN‑binding
RA 3,815,499 T→G 12.0% D90A (GAC→GCC)  pyrE ← orotate phosphoribosyltransferase
MC JC 3,815,858 Δ82 bp 100% [rph][rph] [rph], [rph]
RA 3,821,098 G→A 5.1% A24A (GCC→GCT ligB ← DNA ligase, NAD(+)‑dependent
RA 3,862,660 T→G 29.6% pseudogene (943/1617 nt) glvC ← pseudogene, arbutin specific enzyme IIC component of PTS,enzyme, Transport of small molecules: Carbohydrates, organic acids, alcohols, PTS system, arbutin‑like IIB component, PTS system, arbutin‑like IIC component
RA 3,862,661 T→G 20.7% pseudogene (942/1617 nt) glvC ← pseudogene, arbutin specific enzyme IIC component of PTS,enzyme, Transport of small molecules: Carbohydrates, organic acids, alcohols, PTS system, arbutin‑like IIB component, PTS system, arbutin‑like IIC component
RA 3,862,669 T→A 17.9% pseudogene (934/1617 nt) glvC ← pseudogene, arbutin specific enzyme IIC component of PTS,enzyme, Transport of small molecules: Carbohydrates, organic acids, alcohols, PTS system, arbutin‑like IIB component, PTS system, arbutin‑like IIC component
RA 3,872,494 T→G 6.4% A168A (GCA→GCC dgoD ← D‑galactonate dehydratase
RA 3,971,091 T→G 14.5% G320G (GGT→GGG wecB → UDP‑N‑acetyl glucosamine‑2‑epimerase
RA 3,971,096 T→G 16.0% V322G (GTG→GGG)  wecB → UDP‑N‑acetyl glucosamine‑2‑epimerase
JC 4,001,644 Δ5 bp 100% coding (220‑224/951 nt) corA → magnesium/nickel/cobalt transporter
RA 4,027,131 A→G 14.2% G202G (GGA→GGG fre → NAD(P)H‑flavin reductase
RA 4,145,075 A→C 14.6% N361T (AAC→ACC)  pflD → putative glycine radical domain‑containing pyruvate formate‑lyase
RA 4,151,328 T→G 19.2% A590A (GCA→GCC ppc ← phosphoenolpyruvate carboxylase
RA 4,184,542 C→A 100% P1100Q (CCG→CAG)  rpoB → RNA polymerase, beta subunit
RA 4,219,538 C→T 11.0% A315V (GCG→GTG)  aceK → isocitrate dehydrogenase kinase/phosphatase
RA 4,256,196 A→C 16.9% R90R (CGT→CGG plsB ← glycerol‑3‑phosphate O‑acyltransferase
RA 4,274,842 T→G 9.2% T67P (ACC→CCC)  yjcB ← putative inner membrane protein
RA 4,296,059 C→T 5.4% intergenic (+266/+378) gltP → / ← yjcO glutamate/aspartate:proton symporter/Sel1 family TPR‑like repeat protein
RA 4,296,267 T→C 10.9% intergenic (+474/+170) gltP → / ← yjcO glutamate/aspartate:proton symporter/Sel1 family TPR‑like repeat protein
RA 4,307,560 T→G 12.8% T81P (ACC→CCC)  alsK ← D‑allose kinase
RA 4,342,595 G→A 3.3% A228A (GCG→GCA melA → alpha‑galactosidase, NAD(P)‑binding
RA 4,421,184 A→C 17.2% T402P (ACC→CCC)  ulaA → L‑ascorbate‑specific enzyme IIC permease component of PTS
RA 4,421,215 A→C 13.1% D412A (GAC→GCC)  ulaA → L‑ascorbate‑specific enzyme IIC permease component of PTS
RA 4,457,130 C→G 16.1% L449V (CTG→GTG)  mpl → UDP‑N‑acetylmuramate:L‑alanyl‑gamma‑D‑glutamyl‑ meso‑diaminopimelate ligase
RA 4,465,779 T→C 11.1% I135V (ATC→GTC)  treB ← trehalose‑specific PTS enzyme: IIB and IIC component
RA 4,590,880 T→C 11.8% D134G (GAC→GGC)  yjiY ← putative transporter

Unassigned missing coverage evidence
   seq id start end size ←reads reads→ gene description
* * ÷ NC_000913 2033585 2034051 467 42 [39] [36] 48 rseX–[yedS] rseX, [yedS]
* * ÷ NC_000913 4233762 4235410 1649 43 [41] [41] 43 [pgi] [pgi]

Unassigned new junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 = 217135148 (1.150)5 (0.040) 5/410 NT 3.5% coding (1641/1719 nt) proS prolyl‑tRNA synthetase
?NC_000913 = 217144 138 (1.110)coding (1632/1719 nt) proS prolyl‑tRNA synthetase
* ? NC_000913 = 251816124 (0.960)5 (0.040) 4/412 NT 4.2% coding (919/1056 nt) dinB DNA polymerase IV
?NC_000913 = 251864 109 (0.880)coding (967/1056 nt) dinB DNA polymerase IV
* ? NC_000913 257908 =NA (NA)109 (0.840) 67/428 NT 99.1% intergenic (+9/‑768) crl/crl pseudogene, sigma factor‑binding protein, RNA polymerase holoenzyme formation stimulator,regulator, Surface structures, transcriptional regulator of cryptic csgA gene for curli surface fibers/pseudogene, sigma factor‑binding protein, RNA polymerase holoenzyme formation stimulator,regulator, Surface structures, transcriptional regulator of cryptic csgA gene for curli surface fibers
?NC_000913 = 1293039 1 (0.010)intergenic (‑117/‑488) hns/tdk global DNA‑binding transcriptional dual regulator H‑NS/thymidine kinase/deoxyuridine kinase
* ? NC_000913 458787 =115 (0.890)51 (0.400) 28/422 NT 31.0% intergenic (+87/‑101) clpX/lon ATPase and specificity subunit of ClpX‑ClpP ATP‑dependent serine protease/DNA‑binding ATP‑dependent protease La
?NC_000913 608011 = NA (NA)intergenic (+23/‑54) hokE/insL1 toxic polypeptide, small/IS186 transposase
* ? NC_000913 = 458799108 (0.840)54 (0.430) 41/420 NT 33.8% intergenic (+99/‑89) clpX/lon ATPase and specificity subunit of ClpX‑ClpP ATP‑dependent serine protease/DNA‑binding ATP‑dependent protease La
?NC_000913 = 609347 NA (NA)intergenic (+170/+112) insL1/entD IS186 transposase/phosphopantetheinyltransferase component of enterobactin synthase multienzyme complex
* ? NC_000913 1207790 =33 (0.260)120 (1.000) 75/396 NT 82.3% coding (290/630 nt) stfP e14 prophage, uncharacterized protein
?NC_000913 1209619 = 21 (0.180)pseudogene (1/501 nt) stfE pseudogene, e14 prophage, side tail fiber protein fragment family,Phage or Prophage Related
* ? NC_000913 = 120780533 (0.260)92 (0.770) 64/396 NT 78.1% coding (305/630 nt) stfP e14 prophage, uncharacterized protein
?NC_000913 = 1209602 21 (0.180)pseudogene (18/501 nt) stfE pseudogene, e14 prophage, side tail fiber protein fragment family,Phage or Prophage Related
* ? NC_000913 1293032 =2 (0.020)138 (1.070) 78/428 NT 98.6% intergenic (‑110/‑495) hns/tdk global DNA‑binding transcriptional dual regulator H‑NS/thymidine kinase/deoxyuridine kinase
?NC_000913 = 1979270 NA (NA)intergenic (‑55/‑483) insA/uspC IS1 repressor TnpA/universal stress protein
* ? NC_000913 1729156 =115 (0.890)3 (0.020) 3/410 NT 2.6% coding (70/4617 nt) lhr putative ATP‑dependent helicase
?NC_000913 1729218 = 118 (0.950)coding (132/4617 nt) lhr putative ATP‑dependent helicase
* ? NC_000913 = 19785020 (0.000)38 (0.310) 30/412 NT 100% intergenic (‑305/+16) flhD/insB1 flagellar class II regulon transcriptional activator, with FlhC/IS1 transposase B
?NC_000913 1979279 = 0 (0.000)intergenic (‑64/‑474) insA/uspC IS1 repressor TnpA/universal stress protein
* ? NC_000913 1979486 =0 (0.000)48 (0.370) 25/428 NT 100% intergenic (‑271/‑267) insA/uspC IS1 repressor TnpA/universal stress protein
?NC_000913 2101749 = NA (NA)intergenic (‑5/+146) wbbL/insH1 pseudogene, lipopolysaccharide biosynthesis protein/IS5 transposase and trans‑activator
* ? NC_000913 = 19794890 (0.000)104 (0.800) 64/428 NT 100% intergenic (‑274/‑264) insA/uspC IS1 repressor TnpA/universal stress protein
?NC_000913 = 2290111 NA (NA)intergenic (‑32/+1) insH1/yejO IS5 transposase and trans‑activator/pseudogene, autotransporter outer membrane homology,putative transport, Not classified, putative ATP‑binding component of a transport system
* ? NC_000913 2514275 =NA (NA)7 (0.050) 7/422 NT 6.0% intergenic (+33/‑54) nupC/insL1 nucleoside (except guanosine) transporter/IS186 transposase
?NC_000913 = 2536318 111 (0.860)coding (487/510 nt) crr glucose‑specific enzyme IIA component of PTS
* ? NC_000913 = 2515612NA (NA)12 (0.090) 12/424 NT 9.9% intergenic (+171/+29) insL1/yfeA IS186 transposase/putative diguanylate cyclase
?NC_000913 2536308 = 110 (0.850)coding (477/510 nt) crr glucose‑specific enzyme IIA component of PTS
* ? NC_000913 = 2643098128 (0.990)5 (0.040) 5/400 NT 3.8% intergenic (‑256/+29) rodZ/rlmN MreB assembly cytoskeletal protein/dual specificity 23S rRNA m(2)A2503, tRNA m(2)A37 methyltransferase, SAM‑dependent
?NC_000913 = 2643108 136 (1.130)intergenic (‑266/+19) rodZ/rlmN MreB assembly cytoskeletal protein/dual specificity 23S rRNA m(2)A2503, tRNA m(2)A37 methyltransferase, SAM‑dependent
* ? NC_000913 = 2713827140 (1.080)4 (0.030) 4/406 NT 2.8% coding (934/1335 nt) srmB ATP‑dependent RNA helicase
?NC_000913 = 2713906 143 (1.170)coding (1013/1335 nt) srmB ATP‑dependent RNA helicase
* ? NC_000913 2791460 =117 (0.910)3 (0.020) 3/414 NT 2.8% coding (190/1449 nt) gabD succinate‑semialdehyde dehydrogenase I, NADP‑dependent
?NC_000913 2791491 = 99 (0.790)coding (221/1449 nt) gabD succinate‑semialdehyde dehydrogenase I, NADP‑dependent
* ? NC_000913 3220235 =102 (0.790)5 (0.040) 5/408 NT 4.1% coding (744/1380 nt) patA putrescine:2‑oxoglutaric acid aminotransferase, PLP‑dependent
?NC_000913 3220407 = 139 (1.130)coding (916/1380 nt) patA putrescine:2‑oxoglutaric acid aminotransferase, PLP‑dependent
* ? NC_000913 = 3424748148 (1.170)5 (0.040)
+TCTTCTTC
5/392 NT 4.1% noncoding (2034/2904 nt) rrlD 23S ribosomal RNA of rrnD operon
?NC_000913 = 4039546 94 (0.740)noncoding (2029/2905 nt) rrlA 23S ribosomal RNA of rrnA operon
* ? NC_000913 4160851 =178 (1.380)13 (0.100) 13/426 NT 6.8% intergenic (‑62/‑272) sthA/fabR pyridine nucleotide transhydrogenase, soluble/transcriptional repressor of fabA and fabB
?NC_000913 = 4499511 NA (NA)intergenic (+10/+83) insD1/yjgX IS2 transposase TnpB/pseudogene fragment
* ? NC_000913 = 4160857180 (1.390)18 (0.140) 17/426 NT 9.1% intergenic (‑68/‑266) sthA/fabR pyridine nucleotide transhydrogenase, soluble/transcriptional repressor of fabA and fabB
?NC_000913 4498183 = NA (NA)intergenic (+242/‑90) intB/insC1 pseudogene, integrase homology,IS, phage, Tn, Phage‑related functions and prophages, KpLE2 phage‑like element, P4‑like integrase/IS2 repressor TnpA