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 96,240 G→A 10.7% intergenic (+19/‑152) trnJ‑Ala → / → rrnW‑16S tRNA‑Ala/ribosomal RNA‑16S
RA BS168_reference_genome 143,659 C→T 100% S100F (TCT→TTT)  rpsE → ribosomal protein S5
RA BS168_reference_genome 146,055 C→T 7.0% P60L (CCT→CTT)  adk → adenylate kinase
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 186,061 T→G 19.5% K126N (AAA→AAC ybbC ← conserved protein of unknown function
RA BS168_reference_genome 193,897 G→A 5.3% G43G (GGC→GGT ybbK ← conserved protein of unknown function
RA BS168_reference_genome 209,957 G→T 5.8% G109C (GGT→TGT)  ybcF → putative enzyme; prophage 1 region
RA BS168_reference_genome 278,734 G→A 6.8% A120T (GCA→ACA)  ycbL → two‑component response regulator [YcbM]
RA BS168_reference_genome 303,830 G→A 5.9% P160S (CCC→TCC)  cwlK ← murein L,D:‑endopeptidase
RA BS168_reference_genome 330,016 C→A 5.6% C81* (TGC→TGA lctE → L‑lactate dehydrogenase
RA BS168_reference_genome 333,932 G→A 81.0% L16L (CTC→CTT mdr ← multidrug‑efflux transporter
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 390,006 A→C 19.2% T755P (ACA→CCA)  srfAB → surfactin synthetase
RA BS168_reference_genome 434,334 G→A 5.5% G27S (GGC→AGC)  pbtP → petrobactin iron‑siderophore ABC transporter (ATP‑binding protein)
RA BS168_reference_genome 485,993 G→T 13.3% intergenic (+37/‑439) ydaN → / → kimA putative regulator of exopolysaccharide synthesis/cyclic di‑AMP controlled transporter of potassium
RA BS168_reference_genome 501,094 C→T 5.6% T310I (ACT→ATT)  dctP → C4‑dicarboxylate transport protein
RA BS168_reference_genome 507,361 T→G 10.1% F166V (TTT→GTT)  ydbO → putative cation efflux system
RA BS168_reference_genome 515,053 G→C 38.7% S188W (TCG→TGG)  ydcA ← putative rhomboid protease
RA BS168_reference_genome 557,865:1 +T 100% intergenic (+416/‑8) yddT → / → ydzN putative exported phage protein/hypothetical protein
RA BS168_reference_genome 577,462 A→T 6.5% H224Q (CAT→CAA ydeR ← putative efflux transporter
RA BS168_reference_genome 792,925 G→A 6.7% A82T (GCA→ACA)  cypD → bifunctional P‑450/NADPH‑P450 reductase 1
RA BS168_reference_genome 859,991 G→A 8.1% D83N (GAC→AAC)  sufLB → general stress protein 18; deglycase
RA BS168_reference_genome 892,022 A→T 100% N196I (AAT→ATT)  malR → transcriptional activator of the Mal operon
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 1,045,961 T→C 6.1% V215A (GTT→GCT)  bmrC → efflux ABC transporter (ATP‑binding protein)
RA BS168_reference_genome 1,073,374 C→T 92.3% W115* (TGG→TAG)  scoC ← transcriptional regulator of extracellular protease production, sporulation and bacilysin production (MarR family)
RA BS168_reference_genome 1,073,426 C→T 9.4% D98N (GAT→AAT)  scoC ← transcriptional regulator of extracellular protease production, sporulation and bacilysin production (MarR family)
RA BS168_reference_genome 1,271,830 G→A 6.3% R400R (CGG→CGA manR → transcriptional DNA‑binding activator
RA BS168_reference_genome 1,291,923 T→G 40.4% F194V (TTT→GTT)  cypE → cytochrome P450 CYP109B1, monooxygenase
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,359,996 T→C 7.4% F181F (TTT→TTC proCC → redundant pyrroline‑5‑carboxylate reductase
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,556,359 G→T 5.6% E725D (GAG→GAT pycA → pyruvate carboxylase
RA BS168_reference_genome 1,565,798 C→A 16.8% intergenic (‑89/‑51) ylbA ← / → ylbB conserved hypothetical protein/putative enzyme
RA BS168_reference_genome 1,648,989 G→C 24.2% D351H (GAC→CAC)  rsmB → RNA‑binding Sun protein; 16S rRNA m5C967 methyltransferase, S‑adenosyl‑L‑methionine‑dependent
RA BS168_reference_genome 1,690,758 C→G 74.0% R214G (CGT→GGT)  codY → transcriptional regulator, GTP and BCAA‑dependent
RA BS168_reference_genome 1,708,535 A→G 5.8% A274A (GCA→GCG flhA → component of the flagellar export machinery
RA BS168_reference_genome 1,944,411 A→G 5.1% I415T (ATT→ACT)  xynD ← arabinoxylan arabinofuranohydrolase
RA BS168_reference_genome 1,953,412 T→C 6.8% E368E (GAA→GAG ldeHA ← biotin carboxylase for subunit LdeHB of methylcrotonyl‑CoA carboxylase
RA BS168_reference_genome 1,988,649 A→G 7.8% P533P (CCT→CCC ppsB ← nonribosomal plipastatin synthetase B involved in synthesis of plipastatin
RA BS168_reference_genome 2,003,223 G→T 30.6% intergenic (+41/+53) iseA → / ← trnSL‑Arg1 inhibitor of DL‑endopeptidases involved in cell‑separation/tRNA‑Arg
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,098,463 T→A 6.4% K110N (AAA→AAT yocM ← putative chaperone
RA BS168_reference_genome 2,128,521 C→T 8.3% Q182Q (CAG→CAA yodD ← putative hydrolase
RA BS168_reference_genome 2,141,680 Δ1 bp 6.1% coding (46/828 nt) kamB ← epsilon‑amino‑beta‑lysine acetyl transferase
RA BS168_reference_genome 2,143,167 C→T 10.4% C170Y (TGT→TAT)  yodR ← putative acyloate‑acetoacetate CoA‑transferase
RA BS168_reference_genome 2,148,040 G→A 9.7% A21V (GCG→GTG)  cgeD ← putative glycosyltransferase involved in maturation of the outermost layer of the spore
RA BS168_reference_genome 2,178,088 A→T 5.4% S476T (TCA→ACA)  dpoL ← DNA polymerase with 3'‑5' exonuclease activity; phage SPbeta
RA BS168_reference_genome 2,253,230 G→C 6.0% Q1359H (CAG→CAC cwlP → lytic transglycosylase; SPbeta phage 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,345,473 T→C 5.3% Y220C (TAT→TGT)  dnaD ← DNA‑remodelling primosomal protein
RA BS168_reference_genome 2,369,488 T→G 6.4% G293G (GGA→GGC tyrA ← prephenate dehydrogenase
RA BS168_reference_genome 2,374,282 A→G 7.8% pseudogene (607/750 nt) trpC ← indole‑3‑glycerol phosphate synthase
RA BS168_reference_genome 2,374,283 G→T 7.8% pseudogene (606/750 nt) trpC ← indole‑3‑glycerol phosphate synthase
RA BS168_reference_genome 2,374,296 A→G 6.1% pseudogene (593/750 nt) trpC ← indole‑3‑glycerol phosphate synthase
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,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,518,776 G→A 5.7% S17S (AGC→AGT spo0A ← response regulator, phosphorylated in response to complex YlbF/YmcA/YaaT
RA BS168_reference_genome 2,566,977 G→T 14.8% Q108K (CAA→AAA)  yqgU ← putative lipoprotein
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,739,967 T→G 14.5% L161R (CTC→CGC)  yrpD → putative lipoprotein
RA BS168_reference_genome 2,829,602 T→C 7.1% L13L (CTT→CTC spoVB → involved in spore cortex synthesis (stage V sporulation); translocase with flippase function for peptidoglycan synthesis
RA BS168_reference_genome 2,908,248 C→A 5.4% A164S (GCT→TCT)  sdhC ← succinate dehydrogenase (cytochrome b558 subunit)
RA BS168_reference_genome 3,031,556 C→A 7.8% intergenic (‑139/‑58) ezrA ← / → hisK negative regulator of FtsZ ring formation/histidinol phosphate phosphatase
RA BS168_reference_genome 3,041,446 C→T 5.5% Q19* (CAG→TAG)  acuC → protein deacetylase (regulates AcsA activity)
RA BS168_reference_genome 3,069,029 T→G 28.4% V41G (GTT→GGT)  pbuO → hypoxanthine/guanine permease
RA BS168_reference_genome 3,089,254 A→C 24.2% L387R (CTT→CGT)  bioI ← cytochrome P450 for pimelic acid formation for biotin biosynthesis
RA BS168_reference_genome 3,125,061 A→G 7.8% F194S (TTT→TCT)  ytnA ← putative amino acid permease
RA BS168_reference_genome 3,212,872 C→T 6.4% Y94Y (TAC→TAT tgl → protein‑glutamine gamma‑glutamyltransferase (transglutaminase)
RA BS168_reference_genome 3,305,521 C→A 16.1% intergenic (‑260/+78) paiA ← / ← sufA polyamine N‑acetyltransferase/sulfur carrier chaperone involved in Fe‑S cluster assembly
RA BS168_reference_genome 3,315,315 T→A 7.6% H272L (CAT→CTT)  hom ← homoserine dehydrogenase
RA BS168_reference_genome 3,373,796 T→A 100% D18E (GAT→GAA yuzM → conserved protein of unknown function
RA BS168_reference_genome 3,442,493 T→A 6.8% E346V (GAA→GTA)  copA ← copper transporter ATPase
RA BS168_reference_genome 3,482,162 T→G 5.7% T182P (ACA→CCA)  gapA ← glyceraldehyde‑3‑phosphate dehydrogenase (NAD‑dependent, glycolytic)
RA BS168_reference_genome 3,504,684 T→G 19.9% T134P (ACA→CCA)  ganA ← short chain beta‑1,4‑galacto‑oligosaccharides beta‑galactosidase (beta‑galacto‑pyranoside hydrolase)
RA BS168_reference_genome 3,521,702 C→T 6.8% K148K (AAG→AAA epsH ← putative glycosyltransferase involved in biofilm formation
RA BS168_reference_genome 3,531,109 C→G 26.1% L159V (CTT→GTT)  pnbA → para‑nitrobenzyl esterase (intracellular esterase B)
RA BS168_reference_genome 3,542,242 C→A 6.4% L150L (CTG→CTT yvdQ ← conserved protein of unknown function
RA BS168_reference_genome 3,626,895 A→C 22.4% S40A (TCC→GCC)  yvjA ← putative integral inner membrane protein
RA BS168_reference_genome 3,658,838 G→T 24.6% pseudogene (55/486 nt) tuaAn ← putative undecaprenyl‑phosphate N‑acetylgalactosaminyl‑1‑phosphate transferase; N‑terminal part of TuaA
RA BS168_reference_genome 3,672,798 G→T 7.1% intergenic (‑1176/‑131) yvzI ← / → yvzE putative CDP‑glycerol:poly(glycerophosphate) glycerophosphotransferase (fragment)/putative UTP‑glucose‑1‑phosphate uridylyltransferase (UDP‑glucose pyrophosphorylase) (UDPGP) (fragment)
RA BS168_reference_genome 3,702,947 C→T 5.5% H185H (CAC→CAT rbsK → ribokinase
RA BS168_reference_genome 3,731,793 A→G 6.4% intergenic (‑24/+29) ptpZ ← / ← ptkA protein tyrosine‑phosphatase/maintenance protein tyrosine kinase involved in biofilm formation
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,834,805 C→T 6.1% A172A (GCG→GCA argS ← arginyl‑tRNA synthetase
RA BS168_reference_genome 3,837,955 A→C 24.3% N39T (AAC→ACC)  albC → putative transporter involved in subtilosin production
RA BS168_reference_genome 3,865,646 T→G 20.2% K227N (AAA→AAC pta ← phosphotransacetylase
RA BS168_reference_genome 3,919,423 C→A 9.8% L404L (CTG→CTT 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 3,966,634 G→A 9.6% P546P (CCG→CCA katX → major catalase in spores
RA BS168_reference_genome 3,973,838 T→G 12.0% K418N (AAA→AAC cydD ← ABC membrane transporter (ATP‑binding protein) required for cytochrome bb' function (reductant efflux pump)
RA BS168_reference_genome 4,021,185:1 +T 5.8% coding (122/330 nt) yxiH ← conserved protein of unknown function
RA BS168_reference_genome 4,072,534 C→T 5.4% R147H (CGT→CAT)  tpeJ ← two‑component response regulator [TpeK(YxdK)]
RA BS168_reference_genome 4,090,006 T→G 12.9% Y435S (TAC→TCC)  htpG ← class III heat‑shock protein (ATP‑dependent molecular chaperone HSP90)
RA BS168_reference_genome 4,093,357 A→C 16.7% S161A (TCT→GCT)  yxbF ← putative transcriptional regulator
RA BS168_reference_genome 4,099,811 T→G 5.2% F296V (TTT→GTT)  asnH → asparagine synthetase (glutamine‑hydrolyzing)
RA BS168_reference_genome 4,108,412 T→G 24.4% S119A (TCG→GCG)  yxnA → putative oxidoreductase
RA BS168_reference_genome 4,121,077 G→T 12.1% intergenic (+21/+89) ahpF → / ← bglA alkyl hydroperoxide reductase (large subunit)/aryl‑6‑phospho‑beta‑glucosidase
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 pSIJ663_pHT315_araE_xy 121 G→C 5.3% intergenic (–/‑2857)  / → eryR –/–
RA pSIJ663_pHT315_araE_xy 131 G→C 7.0% intergenic (–/‑2847)  / → eryR –/–
RA pSIJ663_pHT315_araE_xy 137 A→C 5.0% intergenic (–/‑2841)  / → eryR –/–
RA pSIJ663_pHT315_araE_xy 2,179 Δ1 bp 100% intergenic (–/‑799)  / → eryR –/–
JC pSIJ663_pHT315_araE_xy 2,655 Δ3 bp 100% intergenic (–/‑321)  / → eryR –/–
RA pSIJ663_pHT315_araE_xy 3,729 T→G 7.2% intergenic (+14/‑142) eryR → / → ampR –/–
RA pSIJ663_pHT315_araE_xy 3,735 A→T 9.8% intergenic (+20/‑136) eryR → / → ampR –/–
RA pSIJ663_pHT315_araE_xy 3,741 G→T 9.9% intergenic (+26/‑130) eryR → / → ampR –/–
RA pSIJ663_pHT315_araE_xy 3,743 A→T 5.1% intergenic (+28/‑128) eryR → / → ampR –/–
RA pSIJ663_pHT315_araE_xy 3,749 A→T 5.2% intergenic (+34/‑122) eryR → / → ampR –/–
RA pSIJ663_pHT315_araE_xy 3,783 G→T 5.3% intergenic (+68/‑88) eryR → / → ampR –/–
JC pSIJ663_pHT315_araE_xy 3,843 Δ5 bp 100% intergenic (+128/‑24) eryR → / → ampR –/–
RA pSIJ663_pHT315_araE_xy 4,625 A→C 6.4% K252T (AAG→ACG)  ampR →
RA pSIJ663_pHT315_araE_xy 4,648 T→A 10.6% Y260N (TAC→AAC) ‡ ampR →
RA pSIJ663_pHT315_araE_xy 4,649 A→C 8.3% Y260S (TAC→TCC) ‡ ampR →
RA pSIJ663_pHT315_araE_xy 5,049 C→T 100% intergenic (+318/–) ampR → /  –/–
RA pSIJ663_pHT315_araE_xy 5,173 T→G 6.9% intergenic (+442/–) ampR → /  –/–
RA pSIJ663_pHT315_araE_xy 5,302 A→C 6.1% intergenic (+571/–) ampR → /  –/–
RA pSIJ663_pHT315_araE_xy 5,464 T→G 6.3% intergenic (+733/–) ampR → /  –/–
RA pSIJ663_pHT315_araE_xy 5,467 T→G 5.9% intergenic (+736/–) ampR → /  –/–
RA pSIJ663_pHT315_araE_xy 5,612 T→A 71.5% intergenic (+881/–) ampR → /  –/–

Unassigned missing coverage evidence
   seq id start end size ←reads reads→ gene description
* * ÷ BS168_reference_genome 160952–161238 165598–161669 432–4647 16 [14] [6] 18 [rrnI‑16S]–[rrnI‑5S] [rrnI‑16S],rrnI‑23S,[rrnI‑5S]
* * ÷ BS168_reference_genome 166445–166813 171219–170393 3581–4775 15 [12] [13] 15 rrnH‑16S–[rrnH‑5S] rrnH‑16S,rrnH‑23S,[rrnH‑5S]
* * ÷ BS168_reference_genome 171470–171677 176208–172977 1301–4739 18 [11] [14] 15 rrnG‑16S–[rrnG‑5S] rrnG‑16S,rrnG‑23S,[rrnG‑5S]
* * ÷ BS168_reference_genome 3485361–3485463 3486758 1296–1398 16 [14] [0] 42 [araE]–[araR] [araE],[araR]

Unassigned new junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? BS168_reference_genome 9457 =96 (1.550)3626 (61.670) 8/222 NT 97.5% coding (2464/2466 nt) gyrA DNA gyrase (subunit A)
?BS168_reference_genome 3484075 = 94 (1.600)coding (1392/1395 nt) araE arabinose‑related compounds permease
* ? BS168_reference_genome 3486759 =0 (0.000)19 (0.010)
+50 bp
16/134 NT 91.8% intergenic (+1/+48) araR/yvbT transcriptional repressor of the ara regulon (LacI family)/putative alkanal monooxygenase
?pSIJ663_pHT315_araE_xy 9781 = 6 (0.000)intergenic (+5050/–) ampR/– –/–
* ? pSIJ663_pHT315_araE_xy = 555118 (0.000)10379 (1.110) 204/224 NT 99.8% intergenic (+820/–) ampR/– –/–
?pSIJ663_pHT315_araE_xy 9889 = NA (NA)intergenic (+5158/–) ampR/– –/–
* ? pSIJ663_pHT315_araE_xy = 56144 (0.000)13 (0.000) 11/228 NT 79.0% intergenic (+883/–) ampR/– –/–
?pSIJ663_pHT315_araE_xy 9532 = 3 (0.000)intergenic (+4801/–) ampR/– –/–